The epidermis is the outmost skin tissue. It operates as a first defense system to process inflammatory signals and responds by producing inflammatory mediators that promote the recruitment of immune cells. Various skin diseases such as atopic dermatitis occur as a result of the defect of proper skin barrier function and successive impaired inflammatory responses. The onset of such a skin disease links to the disturbed epidermal homeostasis regulated by appropriate self-renewal and differentiation of epidermal stem cells. The theory of physiologically structured population models provides a versatile framework to formulate mathematical models which describe the growth dynamics of a cell population such as the epidermis. In this paper, we develop an algorithm to implement stochastic simulation for a class of physiologically structured population models. We demonstrate that the developed algorithm is applicable to several cell population models and typical age-structured population models. On the basis of the developed algorithm, we investigate stochastic dynamics of skin cell populations and spread of inflammation. It is revealed that demographic stochasticity can bring considerable impact on the outcome of inflammation spread at the tissue level.