Three-Dimensional Direct Numerical Simulation of Flow past A Near-Wall Circular Cylinder:Combined Effects of Gap Ratio and Boundary Layer Thickness on Flow Profiles and Pressure Distribution

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.

Quasi-direct numerical simulations of the flow characteristics of a thermal plasma reactor with counterflow jet

Three-dimensional quasi-direct numerical simulations have been performed to investigate a thermal plasma reactor with a counterflow jet. The effects of the momentum flux ratio and distance between the counterflow jet and the thermal plasma jet on the flow characteristics are addressed. The numerical results show that the dimensionless location of the stagnation layer is significantly affected by the momentum flux ratio, but it is not dependent on the distance.Specifically, the stagnation layer is closer to the plasma torch outlet with the increase of the momentum flux ratio. Furthermore, the flow regimes of the stagnation layer and the flow characteristics of the thermal plasma jet are closely related to the momentum flux ratio. The characteristic frequencies associated with the different regimes are identified. The deflecting oscillation flow regimes are found when the momentum flux ratio is low, which provokes axial velocity fluctuations inside the thermal plasma jet. By contrast, for cases with a high momentum flux ratio, flapping flow regimes are distinguished. The thermal plasma jets are very stable and the axial velocity fluctuations mainly exist in the stagnation layer.

Numerical Simulation of Multi-track and Multi-layer Temperature Field on Laser Direct Metal Shaping

To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS process.Based on the“el- ement birth and death”technique of finite element method,a three-dimensional multi-track and multi-layer model for the transient temperature field analysis of LDMS is developed by ANSYS Parametric Design Language(APDL)for the first time.In the fab- ricated modal,X-direction parallel reciprocating scanning paths is introduced.Using the same process parameters,the simulation results show good agreement with the microstructure features of samples which fabricated by LDMS.

Direct numerical simulation of turbulent boundary layer over an anisotropic compliant wall

Direct numerical simulation of a spatially developing turbulent boundary layer over a compliant wall with anisotropic wall material properties is performed. The Reynolds number varies from 300 to approximately 860 along the streamwise direction, based on the external flow velocity and the momentum thickness. Eight typical cases are selected for numerical investigation under the guidance of the monoharmonic analysis. The instantaneous flow fields exhibit the traveling wavy motion of the compliant wall, and the frequency-wavenumber power spectrum of wall pressure fluctuation is computed to quantify the mutual influence of the wall compliance and the turbulent flow at different wave numbers. It is shown that the Reynolds shear stress and the pressure fluctuation are generally enhanced by the wall compliance with the parameters considered in the present study. A dynamical decomposition of the skin-friction coefficient is derived, and a new term (CW) appears due to the wall-induced Reynolds shear stress. The influence of the anisotropic compliant wall motion on the turbulent boundary layer through the wall-induced negative Reynolds shear stress is discussed. To elucidate the underlying mechanism, the budget analysis of the Reynolds stresses transportation is further carried out. The impact of the wall compliance on the turbulent flow is disclosed by examining the variations of the diffusion and velocity-pressure correlation terms. It is shown that increase of the Reynolds stresses inside the flow domain is caused by enhancement of the velocity-pressure correlation term, possibly through the long-range influence of the wall compliance on the pressure field, rather than diffusion of the wall-induced Reynolds shear stress into the fluid flow.

AN ANALYSIS OF SUBGRID-RESOLVED SCALE INTERACTIONS WITH USE OF RESULTS FROM DIRECT NUMERICAL SIMULATIONS

Subgrid nonlinear interaction and energy transfer are analyzed usingdirect numerical simulations of isotropic turbulence. Influences of cutoff wave numberat different ranges of scale on the energetics and dynamics have been investigated.It is observed that subgrid-subgrid interaction dominates the turbulent dynamicswhen cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing thesubgrid energy transfer and nonlinear interaction into ’forward’ and ’backward’ groupsaccording to the sign of triadic interaction, we find that individually each grouphas very large contribution, but the net of them is much smaller, implying thattremendous cancellation happens between these two groups.

Direct numerical simulation of viscoelastic-fluid-based nanofluid turbulent channel flow with heat transfer

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.

Numerical Simulation of Multi-Directional Random Seas

Multi-directional irregular waves are simulated on the basis of the given directional spectrum using a double summation model, a single direction per frequency model and a single summation model. Their results are compared. It is shown that the single direction per frequency model proposed in this paper can generate a realistic wave field. The effects of the model parameters on the simulated results are also studied in this paper and corresponding suggestions are given.

Direct numerical simulation on relevance of fluctuating velocities and drag reduction in turbulent channel flow with spanwise space-dependent electromagnetic force

Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.

A study on periodic boundary condition in direct numerical simulation for gas–solid flow

Direct numerical simulation(DNS) of gas–solid flow at high resolution has been carried out by coupling the lattice Boltzmann method(LBM) for gas flow and the discrete element method(DEM) for solid particles. However,the body force periodic boundary condition(FPBC) commonly used to cut down the huge computational cost of such simulation has faced accuracy concerns. In this study, a novel two-region periodic boundary condition(TPBC) is presented to remedy this problem, with the flow driven in the region with body force and freely evolving in the other region. With simulation cases for simple circulating fluidized bed risers, the validity and advantages of TPBC are demonstrated with more reasonable heterogeneity of the particle distribution as compared to the corresponding case with FPBC.

A New Method to Determine the Grid Directions in Reservoir Numerical Simulation

Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.

Numerical Simulations of One-Directional Fractional Pharmacokinetics Model

In this paper,we present a three-compartment of pharmacokinetics model with irreversible rate constants.The compartment consists of arterial blood,tissues and venous blood.Fick’s principle and the law of mass action were used to develop the model based on the diffusion process.The model is modified into a fractional pharmacokinetics model with the sense of Caputo derivative.The existence and uniqueness of the model are investigated and the positivity of the model is established.The behaviour of the model is investigated by implementing numerical algorithms for the numerical solution of the system of fractional differential equations.MATLAB software is used to plot the graphs for illustrating the variation of drug concentration concerning time.Therefore,the numerical simulations of the model are presented for different values ofαwhich verified the theoretical analysis.Besides,we also observed the pattern of the simulations at the three-compartment of the model by using different values of initial conditions.

Discussion of Direct Numerical Simulation Method for Supercritical Carbon Dioxide Jet Flow

A kind of direct numerical simulation method suitable for supercritical carbon dioxide jet flow has been discussed in this paper. The form of dimensionless nonconservative compressible Navier-Stokes equations in a two-dimensional cartesian coordinate system is derived in detail. High accurate finite difference compact schemes based on non-uniform grid system are introduced to solve the equations. The simulation results of the three vortex pairing phenomenon of plane mixing layer and a compressible axisymmetric jet flow field show that the discussed numerical simulation method is feasible to calculate the supercritical carbon dioxide jet fluid. And it is found that the difficulties of splitting the convective terms in conservation Navier-Stokes equations, which are brought by the supercritical carbon dioxide fluid pressure state equation, can be avoided by solving the nonconservative compressible Navier-Stokes equations.

Numerical simulation on sand sedimentation and erosion characteristics around HDPE sheet sand barrier under different wind angles

For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In this study,using numerical simulations,we examined the behavior of High Density Polyethylene(HDPE)sheet sand barriers under different wind angles,focusing on flow field distribution,windproof efficiency,and sedimentation erosion dynamics.This study discovered that at a steady wind speed,airflow velocity varies as the angle between the airflow and the HDPE barrier changes.Specifically,a 90°angle results in the widest low-speed airflow area on the barrier’s downwind side.If the airflow is not perpendicular to the barrier,it prompts a lateral airflow movement which decreases as the angle expands.The windproof efficiency correlates directly with this angle but inversely with the wind’s speed.Notably,with a wind angle of 90°,wind speed drops by 81%.The minimum wind speed is found at 5.1H(the sand barrier height)on the barrier’s downwind side.As the angle grows,the barrier’s windproof efficiency improves,extending its protective reach.Sedimentation is most prominent on the barrier’s downwind side,as the wind angle shifts from 30°to 90°,the sand sedimentation area on the barrier’s downwind side enlarges by 14.8H.As the angle grows,sedimentation intensifies,eventually overtakes the forward erosion and enlarges the sedimentation area.

Investigations on High-Speed Flash Boiling Atomization of Fuel Based on Numerical Simulations

Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.

Direct numerical simulation of compressible turbulent flows

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.

Direct Numerical Simulation of Particle Dispersion in Gas-Solid Compressible Turbulent Jets

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.

Experimental Study and Numerical Simulation of Directionally Solidified Turbine Blade Casting

汽轮机片样品铸造物的方向性的团结过程在工作被调查。可变退却率用一致 rate.A 在一个退却过程并且与其它相比被使用为热放射转移和方向性的团结过程的微观结构模拟的数学模型基于 CA-FD 被开发方法。温度分发和微观结构被模仿并且与试验性的结果相比。迷路的谷物被预言并且与试验性的结果相比。站台的不平的温度分发是迷路的谷物的形成的主要原因。

Numerical simulation of the direct reduction of pellets in a rotary hearth furnace for zinc-containing metallurgical dust treatment

A mathematical model was established to describe the direct reduction of pellets in a rotary hearth furnace (RHF). In the model, heat transfer, mass transfer, and gas-solid chemical reactions were taken into account. The behaviors of iron metallization and dezincification were analyzed by the numerical method, which was validated by experimental data of the direct reduction of pellets in a Si-Mo furnace. The simulation results show that if the production targets of iron metallization and dezincification are up to 80% and 90%, respectively, the furnace temperature for high-temperature sections must be set higher than 1300 C. Moreover, an undersupply of secondary air by 20% will lead to a decline in iron metallization rate of discharged pellets by 10% and a decrease in dezincing rate by 13%. In addition, if the residence time of pellets in the furnace is over 20 min, its further extension will hardly lead to an obvious increase in production indexes under the same furnace temperature curve.

Validation of the RANS-SOM Combustion Model Using Direct Numerical Simulation of Incompressible Turbulent Reacting Flows

秒顺序时刻燃烧 moGel，由作者求婚了用 incRmpressible 的直接数字的模拟(DNS ) 被验证狂暴的反应 Fhannel 流动。即时 DNS 结果显示出集中和温度变化的近墙的长带结构。DNS 统计结果在关联方程给术语的预算，证明生产和驱散术语是很重要的。DNS 统计数据被用来在 RANS 秒顺序验证闭合模型时刻(SOM ) 燃烧模型。在近墙的修正应该被做的地方，模仿的散开和生产术语与在大多数流动区域的 DNS 数据一致，这被发现，除了在近墙的区域在近墙的修正应该被做的地方，并且为驱散术语的闭合模型需要进一步的改进。代数学的秒顺序时刻(ASOM ) 燃烧模型被 DNS 很好验证。