In this work, a two-dimensional fluid model has been employed to study the characteristics of Ar/O2 radio frequency（RF） inductively coupled plasma discharges. The emphasis of this work has been put on the influence ...In this work, a two-dimensional fluid model has been employed to study the characteristics of Ar/O2 radio frequency（RF） inductively coupled plasma discharges. The emphasis of this work has been put on the influence of the external parameters（i.e., the RF power, the pressure, and the Ar/O2 gas ratio） on the plasma properties. The numerical results show that the RF power has a significant influence on the amplitude of the plasma density rather than on the spatial distribution.However, the pressure and the Ar/O2 gas ratio affect not only the amplitude of the plasma density, but also the spatial uniformity. Finally, the comparison between the simulation results and the experimental data has been made at different gas pressures and oxygen contents, and a good agreement has been achieved.展开更多
Background In the smoothed particle hydrodynamics(SPH)fluid simulation method,the smoothing length affects not only the process of neighbor search but also the calculation accuracy of the pressure solver.Therefore,it ...Background In the smoothed particle hydrodynamics(SPH)fluid simulation method,the smoothing length affects not only the process of neighbor search but also the calculation accuracy of the pressure solver.Therefore,it plays a crucial role in ensuring the accuracy and stability of SPH.Methods In this study,an adaptive SPH fluid simulation method with a variable smoothing length is designed.In this method,the smoothing length is adaptively adjusted according to the ratio of the particle density to the weighted average of the density of the neighboring particles.Additionally,a neighbor search scheme and kernel function scheme are designed to solve the asymmetry problems caused by the variable smoothing length.Results The simulation efficiency of the proposed algorithm is comparable to that of some classical methods,and the variance of the number of neighboring particles is reduced.Thus,the visual effect is more similar to the corresponding physical reality.Conclusions The precision of the interpolation calculation performed in the SPH algorithm is improved using the adaptive-smoothing length scheme;thus,the stability of the algorithm is enhanced,and a larger timestep is possible.展开更多
To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature...To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature.This high temperature can modulate the plasma characteristics to a large degree.We thus study the effect of dielectric window temperature on plasma parameters in two different ICP structures based on COMSOL software.The distributions of various plasma species are examined at different dielectric window temperatures.The concentration of neutral gas is found to be largely modulated at high dielectric window temperature,which further affects the electron collision probability with neutrals and the electron temperature.However,the electron density profiles are barely affected by the dielectric window temperature,which is mainly concentrated at the center of the reactor due to the fixed power input and pressure.展开更多
Planar radio frequency inductively coupled plasmas（ICP） are employed for low-voltage ion implantation processes,with capacitive pulse biasing of the substrate for modulation of the ion energy. In this work, a two-di...Planar radio frequency inductively coupled plasmas（ICP） are employed for low-voltage ion implantation processes,with capacitive pulse biasing of the substrate for modulation of the ion energy. In this work, a two-dimensional（2D） selfconsistent fluid model has been employed to investigate the influence of the pulsed bias power on the nitrogen plasmas for various bias voltages and pulse frequencies. The results indicate that the plasma density as well as the inductive power density increase significantly when the bias voltage varies from 0 V to-4000 V, due to the heating of the capacitive field caused by the bias power. The N＋fraction increases rapidly to a maximum at the beginning of the power-on time, and then it decreases and reaches the steady state at the end of the glow period. Moreover, it increases with the bias voltage during the power-on time, whereas the N2-＋ fraction exhibits a reverse behavior. When the pulse frequency increases to 25 kHz and40 kHz, the plasma steady state cannot be obtained, and a rapid decrease of the ion density at the substrate surface at the beginning of the glow period is observed.展开更多
t A fluid model is employed to investigate the effect of radio frequency bias on the behavior of an argon inductively coupled plasma(ICP).In particular,the effects of ICP source power,single-frequency bias power,and d...t A fluid model is employed to investigate the effect of radio frequency bias on the behavior of an argon inductively coupled plasma(ICP).In particular,the effects of ICP source power,single-frequency bias power,and dual-frequency bias power on the characteristics of ICP are simulated at a fixed pressure of 30 mTorr(1 Torr=1.33322×102 Pa).When the bias frequency is fixed at 27.12 MHz,the two-dimensional(2D)plasma density profile is significantly affected by the bias power at low ICP source power(e.g.,50 W),whereas it is weakly affected by the bias power at higher ICP source power(e.g.,100 W).When dual-frequency(27.12 MHz/2.26 MHz)bias is applied and the sum of bias powers is fixed at 500 W,a pronounced increase in the maximum argon ion density is observed with the increase of the bias power ratio in the absence of ICP source power.As the ratio of 27.12-MHz/2.26-MHz bias power decreases from 500 W/0 W to 0 W/500 W with the ICP source power fixed at 50 W,the plasma density profiles smoothly shifts from edge-high to center-high,and the effect of bias power on the plasma distribution becomes weaker with the bias power ratio decreasing.Besides,the axial ion flux at the substrate surface is characterized by a maximum at the edge of the substrate.When the ICP source power is higher,the 2D plasma density profiles,as well as the spatiotemporal and radial distributions of ion flux at the substrate surface are characterized by a peak in the reactor center,and the distributions of plasma parameters are negligibly affected by the dual-frequency bias power ratio.展开更多
In classical smoothed particle hydrodynamics(SPH)fluid simulation approaches,the smoothing length of Lagrangian particles is typically constant.One major disadvantage is the lack of adaptiveness,which may compromise a...In classical smoothed particle hydrodynamics(SPH)fluid simulation approaches,the smoothing length of Lagrangian particles is typically constant.One major disadvantage is the lack of adaptiveness,which may compromise accuracy in fluid regions such as splashes and surfaces.Attempts to address this problem used variable smoothing lengths.Yet the existing methods are computationally complex and non-efficient,because the smoothing length is typically calculated using iterative optimization.Here,we propose an efficient non-iterative SPH fluid simulation method with variable smoothing length(VSLSPH).VSLSPH correlates the smoothing length to the density change,and adaptively adjusts the smoothing length of particles with high accuracy and low computational cost,enabling large time steps.Our experimental results demonstrate the advantages of the VSLSPH approach in terms of its simulation accuracy and efficiency.展开更多
Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production eff...Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production efficiency of shale oil resources.Effectively stimulating oil shale reservoirs remains a challenging and active research topic.This investigation employed shale specimens obtained from the Longmaxi Formation.Scanning electron microscopy,fluid injection experiments,and fluid-structure interaction simulations were used to comprehensively analyze structural changes and fluid flow behavior under high temperatures from microscopic to macroscopic scales.Experimental results indicate that the temperature has little effect on the structure and permeability of shale before 300℃.However,there are two threshold temperatures within the range of 300 to 600℃that have significant effects on the structure and permeability of oil shale.The first threshold temperature is between 300 and 400℃,which causes the oil shale porosity,pore-fracture ratio,and permeability begin to increase.This is manifested by the decrease in micropores and mesopores,the increase in macropores,and the formation of a large number of isolated pores and fissures within the shale.The permeability increases but not significantly.The second threshold temperature is between 500 and 600℃,which increases the permeability of oil shale significantly.During this stage,micropores and mesopores are further reduced,and macropores are significantly enlarged.A large number of connected and penetrated pores and fissures are formed.More numerous and thicker streamlines appear inside the oil shale.The experimental results demonstrate that high temperatures significantly alter the microstructure and permeability of oil shale.At the same time,the experimental results can provide a reference for the research of in-situ heating techniques in oil shale reservoir transformation.展开更多
Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
Ingestion of microplastics by various organisms has been widely evidenced.Chemicals associated with microplastics(MPs)may be released to digestive tracts upon ingestion.However,the effect of aging and temperature on t...Ingestion of microplastics by various organisms has been widely evidenced.Chemicals associated with microplastics(MPs)may be released to digestive tracts upon ingestion.However,the effect of aging and temperature on the chemical desorption for MPs remains poorly understood.The exposure of polyethylene(PE)particles to UV radiation in dry air,tap water,and sea water was conducted to mimic the aging process of MPs in different environments.Polychlorinated biphenyls(PCBs),as a typical hydrophobic organic contaminant,were preloaded in these aged and pristine PE.The desorption was performed by exposing preloaded PE particles in simulated gastric and gut fluids at 25℃and 40℃.After UV aging,the average diameter of PE particles decreased rapidly with aging time,indicating continuously fragmentation of PE under UV exposure.The desorption of PCBs from PE particles under different conditions varied from 7%to 40%,and that from aged PE in gut fluid at 37℃was significantly higher than those under other conditions(P<0.05).Furthermore,a clear declining trend was observed as lg K_(ow)(octanol-water partition coefficient)value increased.The aging process,hydrophobicity of chemicals,and incubation temperature were important factors on the desorption of PCBs from PE.The present study helps understand the desorption of PCBs from microplastics and the potential risks of microplastics ingestion by organisms.展开更多
In order to capture stable and realistic microscopic features of fluid surface, a surface tension and adhesion method based on implicit incompressible SPH （smoothed particle hydrodynamics） is presented in this paper...In order to capture stable and realistic microscopic features of fluid surface, a surface tension and adhesion method based on implicit incompressible SPH （smoothed particle hydrodynamics） is presented in this paper. It gives a steady and fast tension model and can solve the problem of not considering adhesion. Molecular cohesion and surface minimization are considered for surface tension, and adhesion is added to show the microscopic characteristics of the surface. To simulate surface tension and adhesion stably and efficiently, the surface tension and adhesion model is integrated to an implicit incompressible SPH method. The experimental results show that the method can better simulate surface features in a variety of scenarios compared with previous methods and meanwhile ensure stability and efficiency.展开更多
This article presents a novel and flexible bubble modelling technique for multi-fluid simulations using a volume fraction representation. By combining the volume fraction data obtained from a primary multi-fluid simul...This article presents a novel and flexible bubble modelling technique for multi-fluid simulations using a volume fraction representation. By combining the volume fraction data obtained from a primary multi-fluid simulation with simple and efficient secondary bubble simulation, a range of real-world bubble phenomena are captured with a high degree of physical realism, including large bubble deformation,sub-cell bubble motion, bubble stacking over the liquid surface, bubble volume change, dissolving of bubbles,etc. Without any change in the primary multi-fluid simulator, our bubble modelling approach is applicable to any multi-fluid simulator based on the volume fraction representation.展开更多
The settling flux of biodeposition affects the environmental quality of cage culture areas and determines their environmental carrying capacity.Simple and effective simulation of the settling flux of biodeposition is ...The settling flux of biodeposition affects the environmental quality of cage culture areas and determines their environmental carrying capacity.Simple and effective simulation of the settling flux of biodeposition is extremely important for determining the spatial distribution of biodeposition.Theoretically,biodeposition in cage culture areas without specific emission rules can be simplified as point source pollution.Fluent is a fluid simulation software that can simulate the dispersion of particulate matter simply and efficiently.Based on the simplification of pollution sources and bays,the settling flux of biodeposition can be easily and effectively simulated by Fluent fluid software.In the present work,the feasibility of this method was evaluated by simulation of the settling flux of biodeposition in Maniao Bay,Hainan Province,China,and 20 sampling sites were selected for determining the settling fluxes.At sampling sites P1,P2,P3,P4,P5,Z1,Z2,Z3,Z4,A1,A2,A3,A4,B1,B2,C1,C2,C3 and C4,the measured settling fluxes of biodeposition were 26.02,15.78,10.77,58.16,6.57,72.17,12.37,12.11,106.64,150.96,22.59,11.41,18.03,7.90,19.23,7.06,11.84,5.19 and 2.57 g d^(−1)m^(−2),respectively.The simulated settling fluxes of biodeposition at the corresponding sites were 16.03,23.98,8.87,46.90,4.52,104.77,16.03,8.35,180.83,213.06,39.10,17.47,20.98,9.78,23.25,7.84,15.90,6.06 and 1.65 g d^(−1)m^(−2),respectively.There was a positive correlation between the simulated settling fluxes and measured ones(R=0.94,P=2.22×10^(−9)<0.05),which implies that the spatial differentiation of biodeposition flux was well simulated.Moreover,the posterior difference ratio of the simulation was 0.38,and the small error probability was 0.94,which means that the simulated results reached an acceptable level from the perspective of relative error.Thus,if nonpoint source pollution is simplified to point source pollution and open waters are simplified based on similarity theory,the setting flux of biodeposition in the open waters can be simply and effectively simulated by the fluid simulation software Fluent.展开更多
The multigrid method has been widely used in computational fluid dynamics(CFD)numerical calculations because of its strong convergence.To achieve real-time simulation of a fluid in computer graphics(CG),the operation ...The multigrid method has been widely used in computational fluid dynamics(CFD)numerical calculations because of its strong convergence.To achieve real-time simulation of a fluid in computer graphics(CG),the operation efficiency is also a significant factor to consider except for operational accuracy.For this problem,we introduced two multigrid cycling schemes,V-Cycle and full multigrid(FMG).Moreover,we have proposed a simple geometric multigrid method(GMG),and compared with the existing wide application of algebraic multigrid(AMG).All the calculations are the solution of parallel computing of GPU in this paper.The results showed that our approaches have improved the algorithm’s computational speed and convergence time,which prominently enhanced the efficiency of the fluid simulation.展开更多
The field of fluid simulation is developing rapidly,and data-driven methods provide many frameworks and techniques for fluid simulation.This paper presents a survey of data-driven methods used in fluid simulation in c...The field of fluid simulation is developing rapidly,and data-driven methods provide many frameworks and techniques for fluid simulation.This paper presents a survey of data-driven methods used in fluid simulation in computer graphics in recent years.First,we provide a brief introduction of physical based fluid simulation methods based on their spatial discretization,including Lagrangian,Eulerian,and hybrid methods.The characteristics of these underlying structures and their inherent connection with data driven methodologies are then analyzed.Subsequently,we review studies pertaining to a wide range of applications,including data-driven solvers,detail enhancement,animation synthesis,fluid control,and differentiable simulation.Finally,we discuss some related issues and potential directions in data-driven fluid simulation.We conclude that the fluid simulation combined with data-driven methods has some advantages,such as higher simulation efficiency,rich details and different pattern styles,compared with traditional methods under the same parameters.It can be seen that the data-driven fluid simulation is feasible and has broad prospects.展开更多
We investigated the effect of supply air rate and temperature on formaldehyde emission characteristics in an environment chamber.A three-dimensional computational fluid dynamics(CFD) chamber model for simulating forma...We investigated the effect of supply air rate and temperature on formaldehyde emission characteristics in an environment chamber.A three-dimensional computational fluid dynamics(CFD) chamber model for simulating formaldehyde emission in twelve different cases was developed for obtaining formaldehyde concentration by the area-weighted average method.Laboratory experiments were conducted in an environment chamber to validate the simulation results of twelve different cases and the formaldehyde concentration was measured by continuous sampling.The results show that there was good agreement between the model prediction and the experimental values within 4.3 difference for each case.The CFD simulation results varied in the range from 0.21 mg/m3 to 0.94 mg/m3,and the measuring results in the range from 0.17 mg/m3 to 0.87 mg/m3.The variation trend of formaldehyde concentration with supply air rate and temperature variation for CFD simulation and experiment measuring was consistent.With the existence of steady formaldehyde emission sources,formaldehyde concentration generally increased with the increase of temperature,and it decreased with the increase of air supply rate.We also provided some reasonable suggestions to reduce formaldehyde concentration and to improve indoor air quality for newly decorated rooms.展开更多
Computational fluid dynamics（CFD） simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds num...Computational fluid dynamics（CFD） simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.展开更多
According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in...According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.展开更多
Using finite element technique of the plane-strain problem in solid-liquid two-phase medium, we Studied the char acteristics of 'field precursors' and 'focus precursors' of subsurface fluid and their s...Using finite element technique of the plane-strain problem in solid-liquid two-phase medium, we Studied the char acteristics of 'field precursors' and 'focus precursors' of subsurface fluid and their spatial-temporal evolution in case of dip-slip earthquake. The results show that: ① the change of ground fluid is slow and the anomaly is not prominent in the early period which is of elastic accumulation and non-linear; ② dilatancy emerges and anomalyfocus mainly in the source region in the moderate period which is hardening and of local dilatancy. In the period the focus precursors emerge earlier than the field precursors; ③ anomalies spreed continuously in the source area and new regions with big anomaly emerge out of the source region in the middle-short period which is of large scale dilatancy.展开更多
The suspension characteristics of ultrafine powder slurry in the stirred vessel were simulated by using computational fluid dynamics.The results show that the Rushton disk turbine impeller is more conducive to maintai...The suspension characteristics of ultrafine powder slurry in the stirred vessel were simulated by using computational fluid dynamics.The results show that the Rushton disk turbine impeller is more conducive to maintaining suspended homogeneity and circulation of slurry compared with the pitch blade turbine pumping up impeller and the pitch blade turbine pumping down impeller.And the increase in stirring speed enhances turbulent fluctuation and anisotropic velocity of the fluid at the cost of more power consumption,which improves dispersibility and suspensibility of the particles.Meanwhile,the change of impeller clearance has a weak influence on the flow pattern,and the impeller clearance of 0.32T(T is the diameter of the bottom of the reactor)can achieve better dispersivity and suspensibility of the particles with lower power consumption and larger axial velocity.The experiments of surface coating modification of ultrafine titanium dioxide(TiO2)were carried out under the same conditions for those of the simulation system.The surface film morphology and photocatalytic properties of the modified TiO2 were analyzed,and the obtained data are well consistent with the simulation results.展开更多
The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely dis...The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.展开更多
基金Project supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China(Grant No.2011ZX02403-001)the National Natural Science Foundation of China(Grant No.11205025)
文摘In this work, a two-dimensional fluid model has been employed to study the characteristics of Ar/O2 radio frequency（RF） inductively coupled plasma discharges. The emphasis of this work has been put on the influence of the external parameters（i.e., the RF power, the pressure, and the Ar/O2 gas ratio） on the plasma properties. The numerical results show that the RF power has a significant influence on the amplitude of the plasma density rather than on the spatial distribution.However, the pressure and the Ar/O2 gas ratio affect not only the amplitude of the plasma density, but also the spatial uniformity. Finally, the comparison between the simulation results and the experimental data has been made at different gas pressures and oxygen contents, and a good agreement has been achieved.
基金National Natural Science Foundation of China(61976052,61876043)National Key R&D Program of China(2017YFB1002701,2017YFB1201203).
文摘Background In the smoothed particle hydrodynamics(SPH)fluid simulation method,the smoothing length affects not only the process of neighbor search but also the calculation accuracy of the pressure solver.Therefore,it plays a crucial role in ensuring the accuracy and stability of SPH.Methods In this study,an adaptive SPH fluid simulation method with a variable smoothing length is designed.In this method,the smoothing length is adaptively adjusted according to the ratio of the particle density to the weighted average of the density of the neighboring particles.Additionally,a neighbor search scheme and kernel function scheme are designed to solve the asymmetry problems caused by the variable smoothing length.Results The simulation efficiency of the proposed algorithm is comparable to that of some classical methods,and the variance of the number of neighboring particles is reduced.Thus,the visual effect is more similar to the corresponding physical reality.Conclusions The precision of the interpolation calculation performed in the SPH algorithm is improved using the adaptive-smoothing length scheme;thus,the stability of the algorithm is enhanced,and a larger timestep is possible.
基金financially supported by National Natural Science Foundation of China(NSFC)(Nos.12105035 and 11935005)Guangdong Basic and Applied Basic Research Foundation(No.2021B1515120018)+1 种基金the Fundamental Research Funds for the Central Universities(No.DUT21TD104)the Advanced Space Propulsion Laboratory of BICE and Beijing Engineering Research Center of Efficient and Green Aerospace Propulsion Technology(No.Lab ASP-2020-01)。
文摘To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature.This high temperature can modulate the plasma characteristics to a large degree.We thus study the effect of dielectric window temperature on plasma parameters in two different ICP structures based on COMSOL software.The distributions of various plasma species are examined at different dielectric window temperatures.The concentration of neutral gas is found to be largely modulated at high dielectric window temperature,which further affects the electron collision probability with neutrals and the electron temperature.However,the electron density profiles are barely affected by the dielectric window temperature,which is mainly concentrated at the center of the reactor due to the fixed power input and pressure.
基金supported by the National Natural Science Foundation of China(Grant Nos.11175034,11335004,and 11405019)the Important National Science and Technology Specific Project of China(Grant No.2011 ZX 02403-001)
文摘Planar radio frequency inductively coupled plasmas（ICP） are employed for low-voltage ion implantation processes,with capacitive pulse biasing of the substrate for modulation of the ion energy. In this work, a two-dimensional（2D） selfconsistent fluid model has been employed to investigate the influence of the pulsed bias power on the nitrogen plasmas for various bias voltages and pulse frequencies. The results indicate that the plasma density as well as the inductive power density increase significantly when the bias voltage varies from 0 V to-4000 V, due to the heating of the capacitive field caused by the bias power. The N＋fraction increases rapidly to a maximum at the beginning of the power-on time, and then it decreases and reaches the steady state at the end of the glow period. Moreover, it increases with the bias voltage during the power-on time, whereas the N2-＋ fraction exhibits a reverse behavior. When the pulse frequency increases to 25 kHz and40 kHz, the plasma steady state cannot be obtained, and a rapid decrease of the ion density at the substrate surface at the beginning of the glow period is observed.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11875101 and 11905307).
文摘t A fluid model is employed to investigate the effect of radio frequency bias on the behavior of an argon inductively coupled plasma(ICP).In particular,the effects of ICP source power,single-frequency bias power,and dual-frequency bias power on the characteristics of ICP are simulated at a fixed pressure of 30 mTorr(1 Torr=1.33322×102 Pa).When the bias frequency is fixed at 27.12 MHz,the two-dimensional(2D)plasma density profile is significantly affected by the bias power at low ICP source power(e.g.,50 W),whereas it is weakly affected by the bias power at higher ICP source power(e.g.,100 W).When dual-frequency(27.12 MHz/2.26 MHz)bias is applied and the sum of bias powers is fixed at 500 W,a pronounced increase in the maximum argon ion density is observed with the increase of the bias power ratio in the absence of ICP source power.As the ratio of 27.12-MHz/2.26-MHz bias power decreases from 500 W/0 W to 0 W/500 W with the ICP source power fixed at 50 W,the plasma density profiles smoothly shifts from edge-high to center-high,and the effect of bias power on the plasma distribution becomes weaker with the bias power ratio decreasing.Besides,the axial ion flux at the substrate surface is characterized by a maximum at the edge of the substrate.When the ICP source power is higher,the 2D plasma density profiles,as well as the spatiotemporal and radial distributions of ion flux at the substrate surface are characterized by a peak in the reactor center,and the distributions of plasma parameters are negligibly affected by the dual-frequency bias power ratio.
基金the Key Program of National Natural Science Foundation of China,No.62237001National Natural Science Foundation for Excellent Young Scholars,No.6212200101+2 种基金National Natural Science Foundation for General Program,Nos.62176066 and 61976052Guangdong Provincial Science and Technology Innovation Strategy Fund,No.2019B121203012and Guangzhou Science and Technology Plan,No.202007040005.
文摘In classical smoothed particle hydrodynamics(SPH)fluid simulation approaches,the smoothing length of Lagrangian particles is typically constant.One major disadvantage is the lack of adaptiveness,which may compromise accuracy in fluid regions such as splashes and surfaces.Attempts to address this problem used variable smoothing lengths.Yet the existing methods are computationally complex and non-efficient,because the smoothing length is typically calculated using iterative optimization.Here,we propose an efficient non-iterative SPH fluid simulation method with variable smoothing length(VSLSPH).VSLSPH correlates the smoothing length to the density change,and adaptively adjusts the smoothing length of particles with high accuracy and low computational cost,enabling large time steps.Our experimental results demonstrate the advantages of the VSLSPH approach in terms of its simulation accuracy and efficiency.
基金supported by the Chongqing Natural Science Foundation of Chongqing,China(No.CSTB2022NSCQ-MSX0333)the Science and Technology Research Program of Chongqing Municipal Education Commission(Grant No.KJZD-K202401205)+1 种基金Chongqing Three Gorges University Graduate Research and Innovation Project Funding(No.YJSKY24045)Chongqing Engineering Research Center of Disaster Prevention&Control for Banks and Structures in Three Gorges Reservoir Area(No.SXAPGC24YB14,No.SXAPGC24YB03,No.SXAPGC24YB12)。
文摘Oil shale is characterized by a dense structure,low proportion of pores and fissures,and low permeability.Pore-fracture systems serve as crucial channels for shale oil migration,directly influencing the production efficiency of shale oil resources.Effectively stimulating oil shale reservoirs remains a challenging and active research topic.This investigation employed shale specimens obtained from the Longmaxi Formation.Scanning electron microscopy,fluid injection experiments,and fluid-structure interaction simulations were used to comprehensively analyze structural changes and fluid flow behavior under high temperatures from microscopic to macroscopic scales.Experimental results indicate that the temperature has little effect on the structure and permeability of shale before 300℃.However,there are two threshold temperatures within the range of 300 to 600℃that have significant effects on the structure and permeability of oil shale.The first threshold temperature is between 300 and 400℃,which causes the oil shale porosity,pore-fracture ratio,and permeability begin to increase.This is manifested by the decrease in micropores and mesopores,the increase in macropores,and the formation of a large number of isolated pores and fissures within the shale.The permeability increases but not significantly.The second threshold temperature is between 500 and 600℃,which increases the permeability of oil shale significantly.During this stage,micropores and mesopores are further reduced,and macropores are significantly enlarged.A large number of connected and penetrated pores and fissures are formed.More numerous and thicker streamlines appear inside the oil shale.The experimental results demonstrate that high temperatures significantly alter the microstructure and permeability of oil shale.At the same time,the experimental results can provide a reference for the research of in-situ heating techniques in oil shale reservoir transformation.
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金Supported by the National Natural Science Foundation of China(No.21936004)the Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(No.SML2021SP208)the Fundamental Research Funds for the Central Universities(No.21623118)。
文摘Ingestion of microplastics by various organisms has been widely evidenced.Chemicals associated with microplastics(MPs)may be released to digestive tracts upon ingestion.However,the effect of aging and temperature on the chemical desorption for MPs remains poorly understood.The exposure of polyethylene(PE)particles to UV radiation in dry air,tap water,and sea water was conducted to mimic the aging process of MPs in different environments.Polychlorinated biphenyls(PCBs),as a typical hydrophobic organic contaminant,were preloaded in these aged and pristine PE.The desorption was performed by exposing preloaded PE particles in simulated gastric and gut fluids at 25℃and 40℃.After UV aging,the average diameter of PE particles decreased rapidly with aging time,indicating continuously fragmentation of PE under UV exposure.The desorption of PCBs from PE particles under different conditions varied from 7%to 40%,and that from aged PE in gut fluid at 37℃was significantly higher than those under other conditions(P<0.05).Furthermore,a clear declining trend was observed as lg K_(ow)(octanol-water partition coefficient)value increased.The aging process,hydrophobicity of chemicals,and incubation temperature were important factors on the desorption of PCBs from PE.The present study helps understand the desorption of PCBs from microplastics and the potential risks of microplastics ingestion by organisms.
基金This work was supported by the National Natural Science Foundation of China under Grant Nos. 61572075 and 61702036, the National Key Research and Development Program of China under Grant Nos. 2016YFB0700502 and 2016YFB1001404, and the Fun- damental Research Funds for the Central Universities of China under Grant No. 2302017FRF-TP-17-012A1.
文摘In order to capture stable and realistic microscopic features of fluid surface, a surface tension and adhesion method based on implicit incompressible SPH （smoothed particle hydrodynamics） is presented in this paper. It gives a steady and fast tension model and can solve the problem of not considering adhesion. Molecular cohesion and surface minimization are considered for surface tension, and adhesion is added to show the microscopic characteristics of the surface. To simulate surface tension and adhesion stably and efficiently, the surface tension and adhesion model is integrated to an implicit incompressible SPH method. The experimental results show that the method can better simulate surface features in a variety of scenarios compared with previous methods and meanwhile ensure stability and efficiency.
基金supported by the National High Technology Research and Development Program of China (No. 2013AA013903)the National Natural Science Foundation of China (No. 61120106007)+3 种基金Research Grant of Beijing Higher Institution Engineering Research CenterTsinghua University Initiative Scientific Research Programsupported in part by the U.S. National Science Foundationthe University of North Carolina Arts & Sciences Foundation
文摘This article presents a novel and flexible bubble modelling technique for multi-fluid simulations using a volume fraction representation. By combining the volume fraction data obtained from a primary multi-fluid simulation with simple and efficient secondary bubble simulation, a range of real-world bubble phenomena are captured with a high degree of physical realism, including large bubble deformation,sub-cell bubble motion, bubble stacking over the liquid surface, bubble volume change, dissolving of bubbles,etc. Without any change in the primary multi-fluid simulator, our bubble modelling approach is applicable to any multi-fluid simulator based on the volume fraction representation.
基金support from the National Key Research and Development Program of China(No.2018YFD0900704)the National Natural Science Foundation of China(No.31972796).
文摘The settling flux of biodeposition affects the environmental quality of cage culture areas and determines their environmental carrying capacity.Simple and effective simulation of the settling flux of biodeposition is extremely important for determining the spatial distribution of biodeposition.Theoretically,biodeposition in cage culture areas without specific emission rules can be simplified as point source pollution.Fluent is a fluid simulation software that can simulate the dispersion of particulate matter simply and efficiently.Based on the simplification of pollution sources and bays,the settling flux of biodeposition can be easily and effectively simulated by Fluent fluid software.In the present work,the feasibility of this method was evaluated by simulation of the settling flux of biodeposition in Maniao Bay,Hainan Province,China,and 20 sampling sites were selected for determining the settling fluxes.At sampling sites P1,P2,P3,P4,P5,Z1,Z2,Z3,Z4,A1,A2,A3,A4,B1,B2,C1,C2,C3 and C4,the measured settling fluxes of biodeposition were 26.02,15.78,10.77,58.16,6.57,72.17,12.37,12.11,106.64,150.96,22.59,11.41,18.03,7.90,19.23,7.06,11.84,5.19 and 2.57 g d^(−1)m^(−2),respectively.The simulated settling fluxes of biodeposition at the corresponding sites were 16.03,23.98,8.87,46.90,4.52,104.77,16.03,8.35,180.83,213.06,39.10,17.47,20.98,9.78,23.25,7.84,15.90,6.06 and 1.65 g d^(−1)m^(−2),respectively.There was a positive correlation between the simulated settling fluxes and measured ones(R=0.94,P=2.22×10^(−9)<0.05),which implies that the spatial differentiation of biodeposition flux was well simulated.Moreover,the posterior difference ratio of the simulation was 0.38,and the small error probability was 0.94,which means that the simulated results reached an acceptable level from the perspective of relative error.Thus,if nonpoint source pollution is simplified to point source pollution and open waters are simplified based on similarity theory,the setting flux of biodeposition in the open waters can be simply and effectively simulated by the fluid simulation software Fluent.
基金the National High Technology Research and Development Program of China(“863”Program)[No.2015AA016404]the Fundamental Research Funds for the Central Universities[No.3132016310].
文摘The multigrid method has been widely used in computational fluid dynamics(CFD)numerical calculations because of its strong convergence.To achieve real-time simulation of a fluid in computer graphics(CG),the operation efficiency is also a significant factor to consider except for operational accuracy.For this problem,we introduced two multigrid cycling schemes,V-Cycle and full multigrid(FMG).Moreover,we have proposed a simple geometric multigrid method(GMG),and compared with the existing wide application of algebraic multigrid(AMG).All the calculations are the solution of parallel computing of GPU in this paper.The results showed that our approaches have improved the algorithm’s computational speed and convergence time,which prominently enhanced the efficiency of the fluid simulation.
基金the Natural Key Research and Development Program of China(2018YFB1004902)the Natural Science Foundation of China(61772329,61373085).
文摘The field of fluid simulation is developing rapidly,and data-driven methods provide many frameworks and techniques for fluid simulation.This paper presents a survey of data-driven methods used in fluid simulation in computer graphics in recent years.First,we provide a brief introduction of physical based fluid simulation methods based on their spatial discretization,including Lagrangian,Eulerian,and hybrid methods.The characteristics of these underlying structures and their inherent connection with data driven methodologies are then analyzed.Subsequently,we review studies pertaining to a wide range of applications,including data-driven solvers,detail enhancement,animation synthesis,fluid control,and differentiable simulation.Finally,we discuss some related issues and potential directions in data-driven fluid simulation.We conclude that the fluid simulation combined with data-driven methods has some advantages,such as higher simulation efficiency,rich details and different pattern styles,compared with traditional methods under the same parameters.It can be seen that the data-driven fluid simulation is feasible and has broad prospects.
基金Funded by National Science Foundation(No.50778415 and No.50878177)
文摘We investigated the effect of supply air rate and temperature on formaldehyde emission characteristics in an environment chamber.A three-dimensional computational fluid dynamics(CFD) chamber model for simulating formaldehyde emission in twelve different cases was developed for obtaining formaldehyde concentration by the area-weighted average method.Laboratory experiments were conducted in an environment chamber to validate the simulation results of twelve different cases and the formaldehyde concentration was measured by continuous sampling.The results show that there was good agreement between the model prediction and the experimental values within 4.3 difference for each case.The CFD simulation results varied in the range from 0.21 mg/m3 to 0.94 mg/m3,and the measuring results in the range from 0.17 mg/m3 to 0.87 mg/m3.The variation trend of formaldehyde concentration with supply air rate and temperature variation for CFD simulation and experiment measuring was consistent.With the existence of steady formaldehyde emission sources,formaldehyde concentration generally increased with the increase of temperature,and it decreased with the increase of air supply rate.We also provided some reasonable suggestions to reduce formaldehyde concentration and to improve indoor air quality for newly decorated rooms.
基金National Natural Science Foundation of China (No.50435030)
文摘Computational fluid dynamics（CFD） simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.
基金Project(51074027)supported by the National Natural Science Foundation of China
文摘According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.
文摘Using finite element technique of the plane-strain problem in solid-liquid two-phase medium, we Studied the char acteristics of 'field precursors' and 'focus precursors' of subsurface fluid and their spatial-temporal evolution in case of dip-slip earthquake. The results show that: ① the change of ground fluid is slow and the anomaly is not prominent in the early period which is of elastic accumulation and non-linear; ② dilatancy emerges and anomalyfocus mainly in the source region in the moderate period which is hardening and of local dilatancy. In the period the focus precursors emerge earlier than the field precursors; ③ anomalies spreed continuously in the source area and new regions with big anomaly emerge out of the source region in the middle-short period which is of large scale dilatancy.
基金financially supported by the National Natural Science Foundation of China(Nos.21838003,91834301 and 21878092)the Shanghai Scientific and Technological Innovation Project(No.18JC1410600)+2 种基金the Social Development Program of Shanghai(Nos.17DZ1200900 and 18DZ2252400)the Innovation Program of Shanghai Municipal Education Commissionthe Fundamental Research Funds for the Central Universities(No.222201718002)。
文摘The suspension characteristics of ultrafine powder slurry in the stirred vessel were simulated by using computational fluid dynamics.The results show that the Rushton disk turbine impeller is more conducive to maintaining suspended homogeneity and circulation of slurry compared with the pitch blade turbine pumping up impeller and the pitch blade turbine pumping down impeller.And the increase in stirring speed enhances turbulent fluctuation and anisotropic velocity of the fluid at the cost of more power consumption,which improves dispersibility and suspensibility of the particles.Meanwhile,the change of impeller clearance has a weak influence on the flow pattern,and the impeller clearance of 0.32T(T is the diameter of the bottom of the reactor)can achieve better dispersivity and suspensibility of the particles with lower power consumption and larger axial velocity.The experiments of surface coating modification of ultrafine titanium dioxide(TiO2)were carried out under the same conditions for those of the simulation system.The surface film morphology and photocatalytic properties of the modified TiO2 were analyzed,and the obtained data are well consistent with the simulation results.
基金ChinaCommitteeofEducation theUniver sityofArizonaandtheMetropolitanWaterDistrictofSouthernCaliforni a.
文摘The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.
