Jung-Il Choi

An immersed boundary method for complex incompressible flows

An immersed boundary method for time-dependent, three-dimensional, incompressible flows is presented in this paper. The incompressible Navier-Stokes equations are discretized using a low-diffusion flux splitting method for the inviscid fluxes and second-order central-differences for the viscous components. Higher-order accuracy achieved by using weighted essentially non-oscillatory (WENO) or total variation diminishing (TVD) schemes. An implicit method based on artificial compressibility and dual-time stepping is used for time advancement. The immersed boundary surfaces are defined as clouds of points, which may be structured or unstructured. Immersed-boundary objects are rendered as level sets in the computational domain, and concepts from computational geometry are used to classify points as being outside, near, or inside the immersed boundary. The velocity field near an immersed surface is determined from separate interpolations of the components tangent and normal to the surface. The tangential velocity near the surface is constructed as a power-law function of the local wall normal distance. Appropriate choices of the power law enable the method to approximate the energizing effects of a turbulent boundary layer for higher Reynolds number flows. Five different flow problems (flow over a circular cylinder, an in-line oscillating cylinder, a NACA0012 airfoil, a sphere, and a stationary mannequin) are simulated using the present immersed boundary method, and the predictions show good agreement with previous computational and experimental results. Finally, the flow induced by realistic human walking motion is simulated as an example of a problem involving multiple moving immersed objects.

Drag Reduction by Spanwise Wall Oscillation in Wall-Bounded Turbulent Flows

Drag reduction in turbulent channel and pipe e ows by spanwise (circumferential) wall oscillations is studied numerically. The ine uence of the wall oscillation on near-wall streamwise vortices is examined. By the use of the Stokes second problem, a wall-normal distanceparameter and an acceleration parameterare obtained toestimate the drag reduction rate. A simple equation is derived for expressing the drag reduction rate by spanwise wall oscillations. The relation between near-wall streamwise vortices and low- and high-speed e uids is scrutinized to extract the key parameters. The drag reduction mechanism is analyzed in terms of the attenuation of Reynolds shear stress.

Large Eddy/Reynolds-Averaged Navier-Stokes Simulation of a Mach 5 Compression-Corner Interaction

Simulations of Mach 5 turbulent flow over a 28-deg compression corner are performed using a hybrid large-eddy/ Reynolds-averaged Navier-Stokes method. The model captures the mean-flow structure of the interaction reasonably well, with observed deficiencies relating to an underprediction of the displacement effects of the shock-induced separation region. The computational results provide some support for a recent theory concerning the underlying causes of low-frequency shock-wave oscillation. In the simulations, the sustained presence of a collection of streaks of fluid with lower/higher momentum than the average induces a low-frequency undulation of the separation front. Power spectra obtained at different streamwise stations are in good agreement with experimental results. Downstream of reattachment, the simulations capture a three-dimensional mean-flow structure, dominated by counter-rotating vortices that produce wide variations in the surface skin friction. Predictions of the structure of the reattaching boundary layer agree well with experimental pitot pressure measurements. In comparison with Reynolds-averaged model predictions, the hybrid large-eddy/Reynolds-averaged Navier-Stokes model predicts more amplification of the Reynolds stresses and a broadening of the Reynolds stress distribution within the boundary layer that is probably due to reattachment-shock motion.

Numerical Simulations of Effects of Micro Vortex Generators Using Immersed-Boundary Methods

This work presents an immersed-boundary technique for compressible, turbulent flows and applies the technique to simulate the effects of micro vortex generators in controlling oblique-shock/turbulent boundary-layer interactions. The Reynolds-averaged Navier-Stokes equations, closed using the Menter k-ω turbulence model, are solved in conjunction with the immersed-boundary technique. The approach is validated by comparing solutions obtained using the immersed-boundary technique with solutions obtained on a body-fitted mesh and with experimental laser Doppler anemometry data collected at Cambridge University for Mach 2.5 flow over single micro vortex generators. Simulations of an impinging oblique-shock boundary-layer interaction at Mach 2.5 with and without micro vortex-generator flow control are also performed, considering the development of the flow in the entire wind tunnel. Comparisons are made with experimental laser Doppler anemometry data and surface-pressure measurements from Cambridge University and an analysis of the flow structure is performed. The results show that three dimensional effects initiated by the interaction of the oblique shock with the sidewall boundary layers significantly influence the flow patterns in the actual experiment. The general features of the interactions with and without the micro vortex-generator array are predicted to good accord by the Reynolds-averaged Navier-Stokes/ immersed-boundary model.

Lagrangian statistics in turbulent channel flow

The Lagrangian dispersion of fluid particles in inhomogeneous turbulence is investigated by a direct numerical simulation of turbulent channel flow. Lagrangian velocity and acceleration along a particle trajectory are computed by employing several interpolation schemes. Among the schemes tested, the four-point Hermite interpolation in the homogeneous directions combined with Chebyshev polynomials in the wall-normal direction seems to produce most reliable Lagrangian statistics. Inhomogeneity of Lagrangian statistics in turbulent boundary layer is investigated by releasing many particles at several different wall-normal locations and tracking those particles. The fluid particle dispersion and Lagrangian structure function of velocity are investigated for the Kolmogorov similarity. The behavior of the Lagrangian integral time scales, Kolmogorov constants a0 and C0 of the velocity structure function near the wall are discussed. The intermittent behavior of the fluid particle acceleration is also examined by ...

Large‐eddy simulation of human‐induced contaminant transport in room compartments

UNLABELLED A large-eddy simulation is used to investigate contaminant transport owing to complex human and door motions and vent-system activity in room compartments where a contaminated and clean room are connected by a vestibule. Human and door motions are simulated with an immersed boundary procedure. We demonstrate the details of contaminant transport owing to human- and door-motion-induced wake development during a short-duration event involving the movement of a person (or persons) from a contaminated room, through a vestibule, into a clean room. Parametric studies that capture the effects of human walking pattern, door operation, over-pressure level, and vestibule size are systematically conducted. A faster walking speed results in less mass transport from the contaminated room into the clean room. The net effect of increasing the volume of the vestibule is to reduce the contaminant transport. The results show that swinging-door motion is the dominant transport mechanism and that human-induced wake motion enhances compartment-to-compartment transport. PRACTICAL IMPLICATIONS The effect of human activity on contaminant transport may be important in design and operation of clean or isolation rooms in chemical or pharmaceutical industries and intensive care units for airborne infectious disease control in a hospital. The present simulations demonstrate details of contaminant transport in such indoor environments during human motion events and show that simulation-based sensitivity analysis can be utilized for the diagnosis of contaminant infiltration and for better environmental protection.

Large eddy simulation and zonal modeling of human‐induced contaminant transport

UNLABELLED An immersed boundary method for particulate flow in an Eulerian framework is utilized to examine the effects of complex human motion on the transport of trace contaminants. The moving human object is rendered as a level set in the computational domain, and realistic human walking motion is implemented using a human kinematics model. A large eddy simulation (LES) technique is used to simulate the fluid and particle dynamics induced by human activity. Parametric studies are conducted within a Room-Room and a Room-Hall configuration, each separated by an open doorway. The effects of the average walking speed, initial proximity from the doorway, and the initial mass loading on room-to-room contaminant transport are examined. The rate of mass transport increases as the walking speed increases, but the total amount of material transported is more influenced by the initial proximity of the human from the doorway. The Room-Hall simulations show that the human wake transports material over a distance of about 8 m. Time-dependent data extracted from the simulations is used to develop a room-averaged zonal model for contaminant transport due to human walking motion. The model shows good agreement with the LES results. PRACTICAL IMPLICATIONS The effect of human activity on contaminant transport may be important in applications such as clean or isolation room design for biochemical production lines, in airborne infection control, and in entry/exit into collective protection or decontamination systems. The large eddy simulations (LES) performed in this work allow precise capturing of the local wakes generated by time-dependent human motion and thus provide a means of quantifying contaminant transport due to wake effects. The LES database can be used to develop zonal models for the bulk effects of human-induced contaminant transport. These may be incorporated into multi-zone infiltration models for use in threat-response and exposure mitigation studies.

Mathematical Analysis of Particle Deposition in Human Lungs: An Improved Single Path Transport Model

A dynamic single-path mathematical model was developed that is capable of analyzing detailed deposition patterns of inhaled particles in human lungs. Weibels symmetric lung morphology was adopted as the basic lung structure, and detailed transport processes were evaluated numerically using the fully implicit procedure. Deposition efficiencies by specific mechanisms were individually examined for accuracy and new empirical formulas were incorporated whenever appropriate. Deposition in the alveolar region was divided into deposition fractions in the alveolar duct and alveoli, considering active transport processes between the two regions. The deposition fractions were obtained for each airway generation, serial lung volumetric compartments, and conventional three-compartment anatomic lung regions. In addition, the surface dose and cumulative deposition with time were analyzed. The results showed excellent agreement with available experimental data. The present model provides an improvement from the previously reported models and can be used as a tool in assessing internal dose of inhaled particles under various inhalation conditions.

Compressible-Flow Simulations Using a New Large-Eddy Simulation/Reynolds-Averaged Navier-Stokes Model

A new hybrid large-eddy simulation/Reynolds-averaged Navier–Stokes simulation (LES/RANS) method is presented in this work. In this approach, the resolved turbulence kinetic energy, ensemble-averaged modeled turbulence kinetic energy and turbulence frequency, and time-resolved turbulence frequency are used to form an estimate of an outer-layer turbulence length scale that is nearly Reynolds-number-independent. The ratio of this outer-layer scale with an inner-layer length scale (proportional to the wall distance) is used to construct a blending function that facilitates the shift between an unsteady RANS formulation near solid surfaces and a LES formulation away from the wall. The new model is tested through simulations of compressible flat-plate boundary layers over a widerangeofReynoldsnumbersandMach2.86 flowoverasmoothcompressionramp.Theresultsshowthatthenew modelpredicts mean andsecond-moment statistics that arein goodagreement withexperiment andare comparable with those obtained using an earlier model (Edwards, J. R., Choi, J-I., and Boles, J. A., “Hybrid Large-Eddy/ Reynolds-Averaged Navier–Stokes Simulation of a Mach-5 Compression Corner Interaction,” AIAA Journal, Vol. 464, 2008, pp. 977–991.) which required a case-by-case calibration of a model constant.

Evolution of Ultrafine Particle Size Distributions Following Indoor Episodic Releases: Relative Importance of Coagulation, Deposition and Ventilation

Indoor ultrafine particles (UFP, <100 nm) undergo aerosol processes such as coagulation and deposition, which alter UFP size distribution and accordingly the level of exposure to UFP of different sizes. This study investigates the decay of indoor UFP originated from five different sources: a gas stove and an electric stove, a candle, a hair dryer, and power tools in a residential test building. An indoor aerosol model was developed to investigate differential effects of coagulation, deposition, and ventilation. The coagulation model includes Brownian, van der Waals, and viscosity forces, and also fractal geometry for particles of >24 nm. The model was parameterized using different values of the Hamaker constant for predicting the coagulation rate. Deposition was determined for two different conditions: central fan on versus central fan off. For the case of a central fan running, deposition rates were measured by using a nonlinear solution to the mass balance equation for the whole building. For the central fan off case, an empirical model was used to estimate deposition rates. Ventilation was measured continuously using an automated tracer gas injection and sampling system. The study results show that coagulation is a significant aerosol process for UFP dynamics and the primary cause for the shift of particle size distribution following an episodic high-concentration UFP release with no fans operating. However, with the central mechanical fan on, UFP deposition loss is substantial and comparable to the coagulation loss. These results suggest that coagulation should be considered during high concentration periods (>20,000 cm−3), while particle deposition should be treated as a major loss mechanism when air recirculates through ductwork or mechanical systems. Copyright 2012 American Association for Aerosol Research