DNS analysis of incipient drop impact dynamics using an accurate level set method

A series of 2D direct numerical simulations were performed with an accurate level set method for single drop impacts.The adopted ACLS method was validated to be efficient with perfect mass conservation in both normal and oblique impacts.A square-root correction for neck bases was modified in accuracy as well as scope of applications.In addition,process of jet formation and evolution was studied to reveal internal dynamics in drop impacts.It's found that pressure gradient and vortex are coexisting and completive reasons for jet topology while the inclined angle has a significant effect on them.Mechanisms of jet formation and evolution are different in the front and back necks.With the help of PDF distribution and correction calculation,a compromise in the competition is observed.This work lays a solid foundation for further studies of dynamics in gas-liquid flows.

Multi-scale numerical simulation of fluidized beds: Model applicability assessment

In the past few decades,multi-scale numerical methods have been developed to model dense gas-solidflow in fluidized beds with different resolutions,accuracies,and efficiencies.However,ambiguity needsto be clarified in the multi-scale numerical simulation of fluidized beds:(i)the selection of the submodels,parameters,and numerical resolution;(ii)the multivariate coupling of operating conditions,bed configurations,polydispersity,and additional forces.Accordingly,a state-of-the-art review is performed to assess the applicability of multi-scale numerical methods in predicting dense gas-solid flow influidized beds at specific fluidization regimes(e.g.,bubbling fluidization region,fast fluidization regime),with a focus on the inter-particle collision models,inter-phase interaction models,collision parameters,and polydispersity effect.A mutual restriction exists between resolution and efficiency.Higherresolution methods need more computational resources and thus are suitable for smaller-scale simulations to provide a database for closure development.Lower-resolution methods require fewercomputational resources and thus underpin large-scale simulations to explore macro-scale phenomena.Model validations need to be further conducted under multiple flow conditions and comprehensivemetrics(e.g.,velocity profiles at different heights,bubbles,or cluster characteristics)for furtherimprovement of the applicability of each numerical method.

Study of biomass gasification in an industrial-scale dual circulating fluidized bed(DCFB)using the Eulerian-Lagrangian method

Dual circulating fluidized bed(DCFB)has emerged as an efficient reactor for biomass gasification due to its unique feature of high gas-solid contact efficiency and separated reactions in two reactors,yet the understanding of complex in-furnace phenomena is still lacking.In this study,biomass gasification in an industrial-scale DCFB system was numerically studied using a multiphase particle-in-cell(MP-PIC)method featuring thermochemical sub-models(e.g.,heat transfer,heterogeneous reactions,and homogeneous reactions)under the Eulerian-Lagrangian framework.After model validation,the hydrodynamics and thermochemical characteristics(i.e.,pressure,temperature,and species)in the DCFB are comprehensively investigated.The results show that size-/density-induced segregation makes solid fuels concentrate on the bed surface.Interphase momentum exchange leads to the continuous decrease of the gas pressure axially.In the gasifier and combustor,the lower heating value(LHV)of the gas products is 5.56 MJ/Nm^(3)and 0.2 MJ/Nm^(3)and the combustible gas concentration(CGC)is 65.5%and 1.86%,respectively.The temperature in the combustor is about 100 K higher than that in the gasifier.A higher solid concentration results in a smaller value of particle heat transfer coefficient(HTC).The HTCs range from 50 to 150 W/(m^(2) K)for a solid concentration larger than 0.3 in the combustor while the HTCs range from 100 to 200 W/(m^(2 )K)in the gasifier.The Reynolds number of biomass particles is two orders of magnitude larger than that of the sand particle.The numerical results shed light on the reactor design and process optimization of biomass gasification in DCFBs.

Direct Numerical Simulation of Particle Dispersion in the Flow Around a Circular Cylinder

The major objective of the paper is to use direct numerical simulation method to get the dispersion patterns of different Stokes number particles in cylinder wake at low Reynolds numbers. The gas fluid structure is simulated by a spectral-element method with third-order accuracy. Non-reflecting boundary condition is specified for the outlet downstream boundary. Lagrangian approach is used to trace particle trajectory. The simulation results show factors to be considered to characterize the particle dispersion and find different particle dispersion pattern in a circular cylinder wake and in a plane wake.

Numerical study on three-dimensional CJ detonation waves interacting with isotropic turbulence

The three-dimensional structures of a cellular detonation wave interacting with different turbulent flows were investigated using a one-step irreversible Arrhenius kinetics model. High-resolution bandwidth-optimized WENO scheme of spatial discretization and total variation diminishing temporal integration are used to solve the three dimensional chemically reactive Navier-Stokes equations. The turbulent vertical and entropic forcing effects on the three dimensional detonation wave structures and dynam- ics are analyzed, as well as the detonation effects on tur- bulent vortex structures. It has been found that the turbulence field imposed has created small scale wrinkles embedded in the detonation front, apart from the large scale features of detonation without turbulence. The deto- nation propagating velocity over the leading shock front varies from 0.8 to 1.6 times of CJ velocity and its proba- bility density function （pdf） skews towards sub-CJ velocity and peaks at about 0.9. The recorded detonation velocity always preferentially decays with time, with very rapid accelerations through triple point interactions. Its pdf also skews to sub-CJ velocity, while its overall shape agrees well with W3. The reaction zone is greatly influenced by the vortex, much more irregular and elongated for the turbulent cases. Distributed burning pockets are more likely to be found there. The turbulent kinetic energy is amplified across the detonation, and periodically oscillates downstream the detonation. The off-diagonal components of Reynolds stress also show a rapid rise across the deto- nation and present to be non-zero downstream of detona- tion. Vortex structures are compound results of the convected vortex and the generated vortex by the collision of triple points. The convection term and baroclinic gen- eration term in the transport equation of enstrophy are compared in detail.

Correlation Analysis on the SGS Velocity between Phases in an Isotropic Gasparticle Two-phase Flow with FDF Model

Correlation of the SGS (sub-grid scale) velocity between the two phases in an isotropic gas-particle two-phase flow was numerically investigated with FDF model. The results show the SGS gas velocity seen by the particles varies with the relative velocities between particle phase and gas phase. The relative velocity between the two phases produces the effect of the anisotropic turbulence on the particles. The variation of Stokes number influences the magnitude of the interaction between the two phases.

Numerical prediction of indoor airborne particle concentration in a test chamber with drift-flux model

The time-dependent variation of airborne particle concentration for different sizes in a test chamber was numerically predicted with drift-flux model. The performance of the drift-flux model for particle transport in different kinds of airflow fields was analyzed. The results show the drift-flux model can predict the transport of indoor fine particles reasonably well. When the air flow field varies slowly, the model can predict both the time-dependent variation ratio of the particle concentration and final stable concentration very well, and the difference for particles with different sizes can be also well predicted. When the air flow varies drastically, the accuracy of the model is decreased due to the neglect of the particles’ independent convective terms in the air flow.

Group Theory Analysis of Free Convective BoundaryLayer Behavior at a Stretching Surface

In the Present study, free convection and heat transfer behavior of electrically conducting fluid in the boundary layer over a vertical continuously stretching surface is investigated. The effects of free convection, magnetic field, suction/blowing at the surface and the stretching speed of the surface on the flow and heat transfer characteristics are considered. By applying one-parametric group theory to analysis of the problem, a similarity solution is found. The governing equations of continuity, momentum and energy are solved numerically by a fourth-order Runge-Kutta scheme. The numerical results, which are obtained for the flow and heat transfer characteristics, reveal the influences of the parameters.