To design the optimum acceleration control schedule for the Adaptive Cycle Engine(ACE)in the full flight envelope,this paper establishes a direct simulation model of the ACE transient state.In this model,geometric par...To design the optimum acceleration control schedule for the Adaptive Cycle Engine(ACE)in the full flight envelope,this paper establishes a direct simulation model of the ACE transient state.In this model,geometric parameters are used to replace the component state parameters.The corresponding relationship between geometric parameters and component state parameters is determined by sensitivity analysis.The geometric variables are controlled when the geometric adjustment speed exceeds the limit,and at the same time the corresponding component state parameters are iterated.The gradient optimization algorism is used to optimize the ground acceleration process of ACE,and the control schedule in terms of operating point of compression components and corrected acceleration rate is used as the full-envelope acceleration control schedule based on the similarity principle.The acceleration control schedules of the triple-bypass mode and the double-bypass mode are designed in this paper.The acceleration processes under various flight conditions are simulated using the acceleration control schedules.Compared with the acceleration process with the linear geometric adjustment schedule,the acceleration performance of ACE is improved by the acceleration control schedule,with the impulse of the acceleration process of the triple-bypass mode being increased by 8.7%-12.3% and the impulse of the double-bypass mode acceleration process being increased by 11.8%-14.1%.展开更多
Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed bounda...Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.展开更多
The direct implicit particle-in-cell is a powerful kinetic method for researching plasma characteristics.However,it is time-consuming to obtain the future electromagnetic field in such a method since the field equatio...The direct implicit particle-in-cell is a powerful kinetic method for researching plasma characteristics.However,it is time-consuming to obtain the future electromagnetic field in such a method since the field equations contain time-dependent matrix coefficients.In this work,we propose to explicitly push particles and obtain the future electromagnetic field based on the information about the particles in the future.The new method retains the form of implicit particle pusher,but the future field is obtained by solving the traditional explicit equation.Several numerical experiments,including the motion of charged particle in electromagnetic field,plasma sheath,and free diffusion of plasma into vacuum,are implemented to evaluate the performance of the method.The results demonstrate that the proposed method can suppress finite-grid-instability resulting from the coarse spatial resolution in electron Debye length through the strong damping of high-frequency plasma oscillation,while accurately describe low-frequency plasma phenomena,with the price of losing the numerical stability at large time-step.We believe that this work is helpful for people to research the bounded plasma by using particle-in-cell simulations.展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
In gas-solid flows, particle-particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical ...In gas-solid flows, particle-particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles (e.g., agglomeration and fragmentation) in spite of their initial monodisperse particle distribution. The conventional direct simulation Monte Carlo (DSMC) method for particle collision tracks equally weighted simulation particles, which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions. In this study, a new differentially weighted DSMC (DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian-Lagrangian models for hydrodynamics. Three schemes (mass, momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass, momentum or energy of the whole system unchanged respectively. A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency. The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best. Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC, compared with the equally weighted DSMC. Meanwhile, computational cost can be largely reduced in contrast with direct numerical simulation.展开更多
The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional (2-D) sustainable micro-vibration is investig...The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional (2-D) sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation (OSE). The growth rate from direct numerical simulation(DNS) is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting (T-S) wave on the boundary layer.展开更多
Direct Simulation Monte Carlo(DSMC)solves the Boltzmann equation with large Knudsen number.The Boltzmann equation generally consists of three terms:the force term,the diffusion term and the collision term.While the fi...Direct Simulation Monte Carlo(DSMC)solves the Boltzmann equation with large Knudsen number.The Boltzmann equation generally consists of three terms:the force term,the diffusion term and the collision term.While the first two terms of the Boltzmann equation can be discretized by numerical methods such as the finite volume method,the third term can be approximated by DSMC,and DSMC simulates the physical behaviors of gas molecules.However,because of the low sampling efficiency of Monte Carlo Simulation in DSMC,this part usually occupies large portion of computational costs to solve the Boltzmann equation.In this paper,by Markov Chain Monte Carlo(MCMC)and multicore programming,we develop Direct Simulation Multi-Chain Markov Chain Monte Carlo(DSMC3):a fast solver to calculate the numerical solution for the Boltzmann equation.Computational results show that DSMC3 is significantly faster than the conventional method DSMC.展开更多
A new criterion is presented to detect global convergence to steady state,and to identify local transient characteristics,during rarefied gas flow simulations performed using the direct simulation Monte Carlo(DSMC)met...A new criterion is presented to detect global convergence to steady state,and to identify local transient characteristics,during rarefied gas flow simulations performed using the direct simulation Monte Carlo(DSMC)method.Unlike deterministic computational fluid dynamics(CFD)schemes,DSMC is generally subject to large statistical scatter in instantaneous flow property evaluations,which prevents the use of residual tracking procedures as are often employed in CFD simulations.However,reliable prediction of the time to reach steady state is necessary for initialization of DSMC sampling operations.Techniques currently used in DSMC to identify steady state convergence are usually insensitive to weak transient behavior in small regions of relatively low density or recirculating flow.The proposed convergence criterion is developed with the goal of properly identifying such weak transient behavior,while adding negligible computational expense and allowing simple implementation in any existing DSMC code.Benefits of the proposed technique over existing convergence detection methods are demonstrated for representative nozzle/plume expansion flow,hypersonic blunt body flow and driven cavity flow problems.展开更多
This paper reviews the authors' recent studies on compressible turbulence by using direct numerical simulation (DNS),including DNS of isotropic(decaying) turbulence, turbulent mixing-layer,turbulent boundary-laye...This paper reviews the authors' recent studies on compressible turbulence by using direct numerical simulation (DNS),including DNS of isotropic(decaying) turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.展开更多
A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian me...A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian method was adopted to simulate the particle motion based on one-way coupling. It is found that the turbulent intensity profiles attain self-similar status in the jet downstream regions. At the Stokes number of 1, particles are concentrated largely in the outer boundaries of the large-scale vortex structures with the most uneven distribution and the widest dispersion in the lateral direction. Particles at the much smaller Stokes numbers are distributed evenly in the flow field, and the lateral dispersion is also considerable. Distribution of particles at much larger Stokes numbers is more uniform and the lateral dispersion becomes small. In addition, the inflow conditions have different effects on the particle dispersion. The direct numerical simulation (DNS) results accord with the previous experiments and numerical studies.展开更多
Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can...Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can reduce turbulent flow resistance as compared with water and enhance heat transfer as compared with VBF. Direct numerical simulation(DNS) is performed in this study to explore the mechanisms of heat transfer enhancement(HTE) and flow drag reduction(DR) for the VFBN turbulent flow. The Giesekus model is used as the constitutive equation for VFBN. Our previously proposed thermal dispersion model is adopted to take into account the thermal dispersion effects of nanoparticles in the VFBN turbulent flow. The DNS results show similar behaviors for flow resistance and heat transfer to those obtained in our previous experiments. Detailed analyses are conducted for the turbulent velocity, temperature, and conformation fields obtained by DNSs for different fluid cases, and for the friction factor with viscous, turbulent, and elastic contributions and heat transfer rate with conductive, turbulent and thermal dispersion contributions of nanoparticles, respectively. The mechanisms of HTE and DR of VFBN turbulent flows are then discussed. Based on analogy theory, the ratios of Chilton–Colburn factor to friction factor for different fluid flow cases are investigated, which from another aspect show the significant enhancement in heat transfer performance for some cases of water-based nanofluid and VFBN turbulent flows.展开更多
Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
Direct numerical simulation of decaying homogeneous isotropic turbulence (DHIT) of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we dire...Direct numerical simulation of decaying homogeneous isotropic turbulence (DHIT) of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we directly investigate the influence of polymers on velocity gradient tensor including vorticity and strain. By visualizing vortex tubes and sheets, we observe a remarkable inhibition of vortex structures in an intermediate-scale field and a small-scale field but not for a large scale field in DHIT with polymers. The geometric study indicates a strong relevance among the vorticity vector, rate-of-strain tensor, and polymer conformation tensor. Joint probability density functions show that the polymer effect can increase "strain generation resistance" and "vorticity generation resistance", i.e., inhibit the generation of vortex sheets and tubes, ultimately leading to turbulence inhibition effects.展开更多
The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation (DNS) of incompressible turbulent reacting channel flows. The instantaneous DNS results show the n...The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation (DNS) of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment (SOM) combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment (ASOM) combustion model is well validated by DNS.展开更多
We report the results of protein folding (219M, C34, N36, 2KES, 2KHK) by the method of accelerated molecular dynamics (aMD) at room temperature with the implicit solvent model. Starting from the linear structures,...We report the results of protein folding (219M, C34, N36, 2KES, 2KHK) by the method of accelerated molecular dynamics (aMD) at room temperature with the implicit solvent model. Starting from the linear structures, these proteins successfully fold to the native structure in a lO0-ns aMD simulation. In contrast, they are failed under the traditional MD simulation in the same simulation time. Then we find that the lowest root mean square deviations of helix structures from the native structures are 0.36 A, 0.63 A, 0.52 A, 1.1 A and 0.78 A. What is more, native contacts, cluster and free energy analyses show that the results of the aMD method are in accordance with the experiment very well. All analyses show that the aMD can accelerate the simulation process, thus we may apply it to the field of computer aided drug designs.展开更多
In this paper,we present a direct numerical simulation(DNS) of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number(Re = 1) in a three-dimensional straight channel flow for the first time,us...In this paper,we present a direct numerical simulation(DNS) of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number(Re = 1) in a three-dimensional straight channel flow for the first time,using the Giesekus constitutive model for the fluid.In order to generate and maintain the turbulent fluid motion in the straight channel,a sinusoidal force term is added to the momentum equation,and then the elastic turbulence is numerically realized with an initialized chaotic velocity field and a stretched conformation field.Statistical and structural characteristics of the elastic turbulence therein are analyzed based on the detailed information obtained from the DNS.The fluid mixing enhancement effect of elastic turbulence is also demonstrated for the potential applications of this phenomenon.展开更多
Accurate descriptions of matrix diffusion across the fracture/matrix interface are critical to assessing contaminant migration in fractured media. The classical transfer probability method is only applicable for relat...Accurate descriptions of matrix diffusion across the fracture/matrix interface are critical to assessing contaminant migration in fractured media. The classical transfer probability method is only applicable for relatively large diffusion coefficients and small fracture spacings, due to an intrinsic assumption of an equilibrium concentration profile in the matrix blocks. Motivated and required by practical applications, we propose a direct numerical simulation (DNS) approach without any empirical assumptions. A three-step Lagrangian algorithm was developed and validated to directly track the particle dynamics across the fracture/matrix interface, where particle's diffusive displacement across the discontinuity is controlled by an analytical, one-side reflection probability. Numerical experiments show that the DNS approach is especially efficient for small diffusion coefficients and large fracture spacings, alleviating limitations of the classical modeling approach.展开更多
The new improved directional vector simulation method foranalyzing the reliability of struc- tural systems failure probabilityis researched. This paper also points out the defects of generaldirectional vector simulati...The new improved directional vector simulation method foranalyzing the reliability of struc- tural systems failure probabilityis researched. This paper also points out the defects of generaldirectional vector simulation, and gives rise to a new higheraccuracy approximate integral formula of structural systems failureprobability. A new geometric meaning of characteristic function isobtained. A new simple method of generating uniformly distributedrandom vector sample sin n-dimensional unit hyper-spherical surfaceis put forward and strictly proved. This method is easy to put intopractice. Numerical examples are given to show the applicability andeffectiveness of the suggested approach to structural systemsreliability problems.展开更多
Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on ...Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures' behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures' elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.展开更多
Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models(GCMs)but also in various higher-resoluti...Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models(GCMs)but also in various higher-resolution limited-area models such as cloud-resolving models(CRMs)and large-eddy simulation(LES)models.Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years,this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations:multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions;discrete vs“continuous”representation of hydrometeor types;turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation;theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology;and approaches for developing bulk microphysics parameterizations.Also presented are the spectral bin scheme and particle-based scheme(especially,super-droplet method)for representing explicit microphysics.Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations.Particle-resolved direct numerical simulation(DNS)models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds.Outstanding challenges and future research directions are explored as well.展开更多
基金co-supported by the National Science and Technology Major Project,China(No.J2019-I-0015-0014)the National Natural Science Foundation of China(No.52372397).
文摘To design the optimum acceleration control schedule for the Adaptive Cycle Engine(ACE)in the full flight envelope,this paper establishes a direct simulation model of the ACE transient state.In this model,geometric parameters are used to replace the component state parameters.The corresponding relationship between geometric parameters and component state parameters is determined by sensitivity analysis.The geometric variables are controlled when the geometric adjustment speed exceeds the limit,and at the same time the corresponding component state parameters are iterated.The gradient optimization algorism is used to optimize the ground acceleration process of ACE,and the control schedule in terms of operating point of compression components and corrected acceleration rate is used as the full-envelope acceleration control schedule based on the similarity principle.The acceleration control schedules of the triple-bypass mode and the double-bypass mode are designed in this paper.The acceleration processes under various flight conditions are simulated using the acceleration control schedules.Compared with the acceleration process with the linear geometric adjustment schedule,the acceleration performance of ACE is improved by the acceleration control schedule,with the impulse of the acceleration process of the triple-bypass mode being increased by 8.7%-12.3% and the impulse of the double-bypass mode acceleration process being increased by 11.8%-14.1%.
基金financially supported by the National Key R&D Program of China (Grant No.2022YFB2603000)the National Natural Science Foundation of China (Grant Nos.51779172 and 52179076)。
文摘Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.
基金Project supported by the National Key Research and Development Program of China (Grant No.2022YFE03050001)partly by the National Natural Science Foundation of China (Grant No.12175160)the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD)。
文摘The direct implicit particle-in-cell is a powerful kinetic method for researching plasma characteristics.However,it is time-consuming to obtain the future electromagnetic field in such a method since the field equations contain time-dependent matrix coefficients.In this work,we propose to explicitly push particles and obtain the future electromagnetic field based on the information about the particles in the future.The new method retains the form of implicit particle pusher,but the future field is obtained by solving the traditional explicit equation.Several numerical experiments,including the motion of charged particle in electromagnetic field,plasma sheath,and free diffusion of plasma into vacuum,are implemented to evaluate the performance of the method.The results demonstrate that the proposed method can suppress finite-grid-instability resulting from the coarse spatial resolution in electron Debye length through the strong damping of high-frequency plasma oscillation,while accurately describe low-frequency plasma phenomena,with the price of losing the numerical stability at large time-step.We believe that this work is helpful for people to research the bounded plasma by using particle-in-cell simulations.
基金Project supported by the National Natural Science Foundation of China(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
基金supported by the National Natural Science Foundation of China(51276077 and 51390494)the National Key Basic Research and Development Program(2010CB227004)
文摘In gas-solid flows, particle-particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles (e.g., agglomeration and fragmentation) in spite of their initial monodisperse particle distribution. The conventional direct simulation Monte Carlo (DSMC) method for particle collision tracks equally weighted simulation particles, which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions. In this study, a new differentially weighted DSMC (DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian-Lagrangian models for hydrodynamics. Three schemes (mass, momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass, momentum or energy of the whole system unchanged respectively. A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency. The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best. Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC, compared with the equally weighted DSMC. Meanwhile, computational cost can be largely reduced in contrast with direct numerical simulation.
基金Supported by the National Natural Science Foundation of China(10672052)the Advanced TalentStart-Up Foundation of Jiangsu University(08JDG018)~~
文摘The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional (2-D) sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation (OSE). The growth rate from direct numerical simulation(DNS) is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting (T-S) wave on the boundary layer.
文摘Direct Simulation Monte Carlo(DSMC)solves the Boltzmann equation with large Knudsen number.The Boltzmann equation generally consists of three terms:the force term,the diffusion term and the collision term.While the first two terms of the Boltzmann equation can be discretized by numerical methods such as the finite volume method,the third term can be approximated by DSMC,and DSMC simulates the physical behaviors of gas molecules.However,because of the low sampling efficiency of Monte Carlo Simulation in DSMC,this part usually occupies large portion of computational costs to solve the Boltzmann equation.In this paper,by Markov Chain Monte Carlo(MCMC)and multicore programming,we develop Direct Simulation Multi-Chain Markov Chain Monte Carlo(DSMC3):a fast solver to calculate the numerical solution for the Boltzmann equation.Computational results show that DSMC3 is significantly faster than the conventional method DSMC.
基金NASA for financial support of this work,through grant NNX08AD02A.
文摘A new criterion is presented to detect global convergence to steady state,and to identify local transient characteristics,during rarefied gas flow simulations performed using the direct simulation Monte Carlo(DSMC)method.Unlike deterministic computational fluid dynamics(CFD)schemes,DSMC is generally subject to large statistical scatter in instantaneous flow property evaluations,which prevents the use of residual tracking procedures as are often employed in CFD simulations.However,reliable prediction of the time to reach steady state is necessary for initialization of DSMC sampling operations.Techniques currently used in DSMC to identify steady state convergence are usually insensitive to weak transient behavior in small regions of relatively low density or recirculating flow.The proposed convergence criterion is developed with the goal of properly identifying such weak transient behavior,while adding negligible computational expense and allowing simple implementation in any existing DSMC code.Benefits of the proposed technique over existing convergence detection methods are demonstrated for representative nozzle/plume expansion flow,hypersonic blunt body flow and driven cavity flow problems.
基金supported by the National Basic Research Program of China(2009CB724100)the National Natural Science Foundation of China(10632050,10872205,11072248).
文摘This paper reviews the authors' recent studies on compressible turbulence by using direct numerical simulation (DNS),including DNS of isotropic(decaying) turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.
基金Supported by the Natural Science Foundation of Zhejiang Province (No. 502047 and No. M503094)National Basic Research Program of China (No. 2003CB214500).
文摘A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian method was adopted to simulate the particle motion based on one-way coupling. It is found that the turbulent intensity profiles attain self-similar status in the jet downstream regions. At the Stokes number of 1, particles are concentrated largely in the outer boundaries of the large-scale vortex structures with the most uneven distribution and the widest dispersion in the lateral direction. Particles at the much smaller Stokes numbers are distributed evenly in the flow field, and the lateral dispersion is also considerable. Distribution of particles at much larger Stokes numbers is more uniform and the lateral dispersion becomes small. In addition, the inflow conditions have different effects on the particle dispersion. The direct numerical simulation (DNS) results accord with the previous experiments and numerical studies.
基金supported by the National Natural Science Foundation of China(Grant No.51276046)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20112302110020)+1 种基金the China Postdoctoral Science Foundation(Grant No.2014M561037)the President Fund of University of Chinese Academy of Sciences,China(Grant No.Y3510213N00)
文摘Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can reduce turbulent flow resistance as compared with water and enhance heat transfer as compared with VBF. Direct numerical simulation(DNS) is performed in this study to explore the mechanisms of heat transfer enhancement(HTE) and flow drag reduction(DR) for the VFBN turbulent flow. The Giesekus model is used as the constitutive equation for VFBN. Our previously proposed thermal dispersion model is adopted to take into account the thermal dispersion effects of nanoparticles in the VFBN turbulent flow. The DNS results show similar behaviors for flow resistance and heat transfer to those obtained in our previous experiments. Detailed analyses are conducted for the turbulent velocity, temperature, and conformation fields obtained by DNSs for different fluid cases, and for the friction factor with viscous, turbulent, and elastic contributions and heat transfer rate with conductive, turbulent and thermal dispersion contributions of nanoparticles, respectively. The mechanisms of HTE and DR of VFBN turbulent flows are then discussed. Based on analogy theory, the ratios of Chilton–Colburn factor to friction factor for different fluid flow cases are investigated, which from another aspect show the significant enhancement in heat transfer performance for some cases of water-based nanofluid and VFBN turbulent flows.
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
基金supported by the National Natural Science Foundation of China (Grant No. 10872060)the Fundamental Research Funds for the Central Universities (Grant No. HIT.BRET2.2010008)
文摘Direct numerical simulation of decaying homogeneous isotropic turbulence (DHIT) of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we directly investigate the influence of polymers on velocity gradient tensor including vorticity and strain. By visualizing vortex tubes and sheets, we observe a remarkable inhibition of vortex structures in an intermediate-scale field and a small-scale field but not for a large scale field in DHIT with polymers. The geometric study indicates a strong relevance among the vorticity vector, rate-of-strain tensor, and polymer conformation tensor. Joint probability density functions show that the polymer effect can increase "strain generation resistance" and "vorticity generation resistance", i.e., inhibit the generation of vortex sheets and tubes, ultimately leading to turbulence inhibition effects.
基金Supported by the National Natural Science Foundation of China (50606026, 50736006).
文摘The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation (DNS) of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment (SOM) combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment (ASOM) combustion model is well validated by DNS.
基金Supported by the National Natural Science Foundation of China under Grant Nos 31200545,11274206 and 11574184
文摘We report the results of protein folding (219M, C34, N36, 2KES, 2KHK) by the method of accelerated molecular dynamics (aMD) at room temperature with the implicit solvent model. Starting from the linear structures, these proteins successfully fold to the native structure in a lO0-ns aMD simulation. In contrast, they are failed under the traditional MD simulation in the same simulation time. Then we find that the lowest root mean square deviations of helix structures from the native structures are 0.36 A, 0.63 A, 0.52 A, 1.1 A and 0.78 A. What is more, native contacts, cluster and free energy analyses show that the results of the aMD method are in accordance with the experiment very well. All analyses show that the aMD can accelerate the simulation process, thus we may apply it to the field of computer aided drug designs.
基金Project supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (Grant No. 51121004)the National Natural Science Foundation of China (Grant No. 51076036)+2 种基金the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20112302110020)the Fundamental Research Funds for the Central Universities,China (Grant No. HIT.BRET1.2010008)the Japan Society for the Promotion of Science (JSPS) Research Fellowship
文摘In this paper,we present a direct numerical simulation(DNS) of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number(Re = 1) in a three-dimensional straight channel flow for the first time,using the Giesekus constitutive model for the fluid.In order to generate and maintain the turbulent fluid motion in the straight channel,a sinusoidal force term is added to the momentum equation,and then the elastic turbulence is numerically realized with an initialized chaotic velocity field and a stretched conformation field.Statistical and structural characteristics of the elastic turbulence therein are analyzed based on the detailed information obtained from the DNS.The fluid mixing enhancement effect of elastic turbulence is also demonstrated for the potential applications of this phenomenon.
基金supported by the United States Department of Energythe Desert Research Institute IR&D Funds
文摘Accurate descriptions of matrix diffusion across the fracture/matrix interface are critical to assessing contaminant migration in fractured media. The classical transfer probability method is only applicable for relatively large diffusion coefficients and small fracture spacings, due to an intrinsic assumption of an equilibrium concentration profile in the matrix blocks. Motivated and required by practical applications, we propose a direct numerical simulation (DNS) approach without any empirical assumptions. A three-step Lagrangian algorithm was developed and validated to directly track the particle dynamics across the fracture/matrix interface, where particle's diffusive displacement across the discontinuity is controlled by an analytical, one-side reflection probability. Numerical experiments show that the DNS approach is especially efficient for small diffusion coefficients and large fracture spacings, alleviating limitations of the classical modeling approach.
文摘The new improved directional vector simulation method foranalyzing the reliability of struc- tural systems failure probabilityis researched. This paper also points out the defects of generaldirectional vector simulation, and gives rise to a new higheraccuracy approximate integral formula of structural systems failureprobability. A new geometric meaning of characteristic function isobtained. A new simple method of generating uniformly distributedrandom vector sample sin n-dimensional unit hyper-spherical surfaceis put forward and strictly proved. This method is easy to put intopractice. Numerical examples are given to show the applicability andeffectiveness of the suggested approach to structural systemsreliability problems.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51276046 and 51506037)the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(Grant No.51421063)+2 种基金the China Postdoctoral Science Foundation(Grant No.2016M591526)the Heilongjiang Postdoctoral Fund,China(Grant No.LBH-Z15063)the China Postdoctoral International Exchange Program
文摘Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures' behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures' elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.
基金supported by the US Department of Energy(DOE)’s Office of Science Atmospheric Systems Research(ASR)Programthe Office of Energy Efficiency and Renewable Energy(EERE)Solar Energy Technologies Office(SETO)award(33504)+3 种基金the Brookhaven National Laboratory(BNL)’s Laboratory Directed Research&Development Program(LDRD)(22-065)The Brookhaven National Laboratory is operated by the Brookhaven Science Associates,LLC(BSA),for the US Department of Energy under Contract No.DESC0012704supported by JSPS KAKENHI Grant No.26286089MEXT KAKENHI Grant No.18H04448。
文摘Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models(GCMs)but also in various higher-resolution limited-area models such as cloud-resolving models(CRMs)and large-eddy simulation(LES)models.Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years,this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations:multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions;discrete vs“continuous”representation of hydrometeor types;turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation;theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology;and approaches for developing bulk microphysics parameterizations.Also presented are the spectral bin scheme and particle-based scheme(especially,super-droplet method)for representing explicit microphysics.Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations.Particle-resolved direct numerical simulation(DNS)models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds.Outstanding challenges and future research directions are explored as well.