Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

An atomistic-continuum hybrid simulation of fluid flows over superhydrophobic surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths which cannot be obtained by molecular dynamics simulation alone.

Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

Bubble Behavior and Heat Transfer in Quasi-Steady Pool Boiling in Microgravity

Pool boiling of degassed FC-72 on a plane plate heater has been studied experimentally in microgravity. A quasi-steady heating method is adopted, in which the heating voltage is controlled to increase exponentially with time. Compared with terrestrial experiments, bubble behaviors are very different, and have direct effect on heat transfer. Small, primary bubbles attached on the surface seem to be able to suppress the activation of the cavities in the neighborhoods, resulting in a slow increase of the wall temperature with the heat flux. For the high subcooling, the coalesced bubble has a smooth surface and a small size. It is difficult to cover the whole heater surface, resulting in a special region of gradual transitional boiling in which nucleate boiling and local dry area can co-exist. No turning point corresponding to the transition from nucleate boiling to film boiling can be observed. On the contrary, the surface oscillation of the coalesced bubble at low subcooling may cause more activated nucleate sites, and then the surface temperature may keep constant or even fall down with the increasing heat flux. Furthermore, an abrupt transition to film boiling can also be observed. It is shown that heat transfer coefficient and CHF increase with the subcooling or pressure in microgravity, as observed in normal gravity.

Boiling heat transfer enhancement by using micro-pin-finned surface for electronics cooling

For efficiently cooling electronic components with high heat flux, experiments were conducted to study the flow boiling heat transfer performance of FC-72 over square silicon chips with the dimensions of 10 × 10 × 0.5 mm3. Four kinds of micro-pin-fins with the dimensions of 30 × 60, 30 × 120, 50 × 60, 50 × 120 μm2 (thickness, t × height, h) were fabricated on the chip surfaces by the dry etching technique for enhancing boiling heat transfer. A smooth surface was also tested for comparison. The experiments were made at three different fluid velocities (0.5, 1 and 2 m/s) and three different liquid subcoolings (15, 25 and 35 K). The results were compared with the previous published data of pool boiling. All micro-pin-fined surfaces show a considerable heat transfer enhancement compared with a smooth surface. Flow boiling can remarkably decrease wall superheat compared with pool boiling. At the velocities lower than 1 m/s, the micro-pin-finned surfaces show a sharp increase in heat flux with increasing wall superheat. For all surfaces, the maximum allowable heat flux, qmax, for the normal operation of LSI chips increases with fluid velocity and subcooling. For all micro-pin-finned surfaces, the wall temperature at the critical heat flux (CHF) is less than the upper limit for the reliable operation of LSI chips, 85◦C. The largest value of qmax can reach nearly 148 W/cm2 for micro-pin-finned chips with the fin height of 120 μm at the fluid velocity of 2 m/s and the liquid subcooling of 35 K. The perspectives for the boiling heat transfer experiment of the prospective micro-pin-finned sur- faces, which has been planned to be made in the Drop Tower Beijing/NMLC in the future, are also presented.

Fluid flow induced calcium response in bone cell network

In our previous work, bone cell networks with controlled spacing and functional intercellular gap junctions had been successfully established by using microcontact printing and self assembled monolayers technologies [Guo, X. E., E. Takai, X. Jiang, Q. Xu, G. M. Whitesides, J. T. Yardley, C. T. Hung, E. M. Chow, T. Hantschel, and K. D. Costa. Mol. Cell. Biomech. 3:95-107, 2006]. The present study investigated the calcium response and the underlying signaling pathways in patterned bone cell networks exposed to a steady fluid flow. The glass slides with cell networks were separated into eight groups for treatment with specific pharmacological agents that inhibit pathways significant in bone cell calcium signaling. The calcium transients of the network were recorded and quantitatively evaluated with a set of network parameters. The results showed that 18 alpha-GA (gap junction blocker), suramin (ATP inhibitor), and thapsigargin (depleting intracellular calcium stores) significantly reduced the occurrence of multiple calcium peaks, which were visually obvious in the untreated group. The number of responsive peaks also decreased slightly yet significantly when either the COX-2/PGE(2) or the NOS/nitric oxide pathway was disrupted. Different from all other groups, cells treated with 18 alpha-GA maintained a high concentration of intracellular calcium following the first peak. In the absence of calcium in the culture medium, the intracellular calcium concentration decreased slowly with fluid flow without any calcium transients observed. These findings have identified important factors in the flow mediated calcium signaling of bone cells within a patterned network.

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included groundbased tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.

Experimentally Determined Flame Properties near Flammability Limits Under Gravity and Microgravity Conditions

Flammability limits for flames propagating in a rich propane/air mixture under gravity conditions appeared to be 6.3% C3H8 for downward propagation and 9.2% C3H8 for upward propagation. Different limits might be explained by the action of preferential diffusion of the deficient reactant (Le < 1) on the limit flames, which are in different states of instability. In one of the previous studies, the flammability limits under microgtravity conditions were found to be between the upward and downward limits obtained in a standard flammability tube under normal gravity conditions. It was found in those experiments that there are two limits under microgravity conditions: one indicated by visible flame propagation and another indicated by an increase of pressure without observed flame propagation. These limits were found to be far behind the limit for downward-propagating flame at 1 g (6.3% C3H8) and close to the limit for upward-propagating flame at 1 g (9.2% C3H8). It was decided in the present work to apply a special schlieren system and instant temperature measuring system for drop tower experiments to observe combustion development during propagation of the flame front. A small cubic closed vessel (inner side, 9 cm 9 cm 9 cm) with schlieren quality glass windows were used to study limit flames under gravity and microgravity conditions. Flame development in rich limit mixtures, not visible in previous experiments under microgravity conditions for strait photography, was identified with the use of the schlieren method and instant temperature measuring system. It was found in experiments in a small vessel that there is practically no difference in flammability limits under gravity and microgravity conditions. In this paper, the mechanism of flame propagation under these different conditions is systematically studied and compared and limit burning velocity is estimated.