Human thermophysiological models: Quantification of uncertainty in the output quantities of the passive system due to uncertainties in the control equations of the active system via the Monte Carlo method.

Uncertainty propagation analysis in the Fiala thermophysiological model is performed by the Monte Carlo Method. The uncertainties of the output quantities of the passive system, due to imported uncertainties in the coefficients of the control equations of the active system, caused by the variation of the experimental data, are computed. The developed and implemented in-house code is accordingly validated. The effect of the input uncertainties, in each of the four main responses (shivering, vasodilatation, vasoconstriction, sweating) of the active system, is separately examined by simulating the human exposure from neutral conditions to cold and hot environments. It is predicted that the maximum output uncertainties of the response mechanisms may be of the same order of magnitude as the imported ones, while the corresponding maximum uncertainties in core and skin temperatures always remain less than 2%. The maximum absolute deviations of the rectal (core) temperatures from their estimated mean values may be up to 0.72 °C and 0.22 °C, due to input uncertainties in shivering and sweating respectively, while the corresponding deviations due to uncertainties in vasomotion processes are negligible. The deviations, particularly the ones due to shivering, are significant, since differences of a few tenths of a degree may have large impact in human health. The maximum absolute deviations of the skin temperatures are 0.42 °C in the hands due to uncertainties in shivering and 0.69 °C in the feet due to uncertainties in vasodilatation. These deviations are less significant than the core ones, but they may still affect human thermal sensation and comfort. The present analysis provides a better insight in the dynamic response of the model and indicates which response mechanism needs to be further investigated by more accurate estimates in order to improve model reliability. It can be also applied in other human thermophysiological models.

Oscillatory pressure-driven rarefied binary gas mixture flow between parallel plates.

The rarefied, oscillatory, pressure-driven binary gas mixture flow between parallel plates is computationally investigated in terms of the mixture molar fraction and molecular mass ratio of the species, in a wide range of gas rarefaction and oscillation frequency. Modeling is based on the McCormack kinetic model. The output quantities are in dimensionless form and include the flow rate, wall shear stress and pumping power of the mixture, as well as the velocity and shear stress distributions and flow rates of the species. The presented results are for He-Xe and Ne-Ar. The heavier species are affected more drastically than the lighter ones from the inertial forces, resulting to large differences between the flow rate amplitudes of the species, which are increased as the flow becomes less rarefied, provided that the oscillation frequency is adequately high. At very high frequencies the ratio of the flow rate amplitudes of the light over the heavy species tends to the inverse of their molecular mass ratio in the whole range of gas rarefaction. The velocity overshooting effect becomes more pronounced as the molecular mass is increased. The mixture flow rate amplitude is larger, while its phase angle is smaller, than the corresponding ones of single gas, and they both vary nonmonotonically with the molar fraction. The effect of the mixture composition on the wall shear stress and pumping power is small. The present work may be useful in the design of gas separation devices, operating at moderate and high frequencies in rarefied and dense atmospheres.

Design Guidelines for Thermally Driven Micropumps of Different Architectures Based on Target Applications via Kinetic Modeling and Simulations.

The manufacturing process and architecture of three Knudsen type micropumps are discussed and the associated flow performance characteristics are investigated. The proposed fabrication process, based on the deposition of successive dry film photoresist layers with low thermal conductivity, is easy to implement, adaptive to specific applications, cost-effective, and significantly improves thermal management. Three target application designs, requiring high mass flow rates (pump A), high pressure differences (pump B), and relatively high mass flow rates and pressure differences (pump C), are proposed. Computations are performed based on kinetic modeling via the infinite capillary theory, taking into account all foreseen manufacturing and operation constraints. The performance characteristics of the three pump designs in terms of geometry (number of parallel microchannels per stage and number of stages) and inlet pressure are obtained. It is found that pumps A and B operate more efficiently at pressures higher than 5 kPa and lower than 20 kPa, respectively, while the optimum operation range of pump C is at inlet pressures between 1 kPa and 20 kPa. In all cases, it is advisable to have the maximum number of stages as well as of parallel microchannels per stage that can be technologically realized.

Gas Mixing and Final Mixture Composition Control in Simple Geometry Micro-mixers via DSMC Analysis.

The mixing process of two pressure driven steady-state rarefied gas streams flowing between two parallel plates was investigated via DSMC (Direct Simulation Monte Carlo) for different combinations of gases. The distance from the inlet, where the associated relative density difference of each species is minimized and the associated mixture homogeneity is optimized, is the so-called mixing length. In general, gas mixing progressed very rapidly. The type of gas surface interaction was clearly the most important parameter affecting gas mixing. As the reflection became more specular, the mixing length significantly increased. The mixing lengths of the HS (hard sphere) and VHS (variable hard sphere) collision models were higher than those of the VSS (variable soft sphere) model, while the corresponding relative density differences were negligible. In addition, the molecular mass ratio of the two components had a minor effect on the mixing length and a more important effect on the relative density difference. The mixture became less homogenous as the molecular mass ratio reduced. Finally, varying the channel length and/or the wall temperature had a minor effect. Furthermore, it was proposed to control the output mixture composition by adding in the mixing zone, the so-called splitter, separating the downstream flow into two outlet mainstreams. Based on intensive simulation data with the splitter, simple approximate expressions were derived, capable of providing, once the desired outlet mixture composition was specified, the correct position of the splitter, without performing time consuming simulations. The mixing analysis performed and the proposed approach for controlling gas mixing may support corresponding experimental work, as well as the design of gas micro-mixers.

Predicting the Knudsen paradox in long capillaries by decomposing the flow into ballistic and collision parts.

The well-known Knudsen paradox observed in pressure driven rarefied gas flows through long capillaries is quantitatively explored by decomposing the particle distribution function into its ballistic and collision parts. The classical channel, tube, and duct Poiseuille flows are considered. The solution is obtained by a typical direct simulation Monte Carlo algorithm supplemented by a suitable particle decomposition indexation process. It is computationally confirmed that in the free-molecular and early transition regimes the reduction rate of the ballistic flow is larger than the increase rate of the collision flow deducing the Knudsen minimum of the overall flow. This description interprets in a precise, quantitative manner the appearance of the Knudsen minimum and verifies previously reported qualitative physical arguments.