Ozone reacts with skin lipids such as squalene, generating an array of organic compounds, some of which can act as respiratory or skin irritants. Thus, it is important to quantify and predict the formation of these products under different conditions in indoor environments. We developed the kinetic multilayer model that explicitly resolves mass transport and chemical reactions at the skin and in the gas phase (KM-SUB-Skin). It can reproduce the concentrations of ozone and organic compounds in previous measurements and new experiments. This enabled the spatial and temporal concentration profiles in the skin oil and underlying skin layers to be resolved. Upon exposure to ~30 ppb ozone, the concentrations of squalene ozonolysis products in the gas phase and in the skin reach up to several ppb and on the order of ~10 mmol m-3 . Depending on various factors including the number of people, room size, and air exchange rates, concentrations of ozone can decrease substantially due to reactions with skin lipids. Ozone and dicarbonyls quickly react away in the upper layers of the skin, preventing them from penetrating deeply into the skin and hence reaching the blood.