The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resour...The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.展开更多
In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge mod...In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.展开更多
Detailed formulation of the detached-eddy simulation (DES) based on shear stress transport model (SST) was deduced. Based on finite volume method, detached-eddy simulation of super- sonic turbulent flow over rearw...Detailed formulation of the detached-eddy simulation (DES) based on shear stress transport model (SST) was deduced. Based on finite volume method, detached-eddy simulation of super- sonic turbulent flow over rearward facing step was conducted. Multi-block structured non-conform grid was designed. The obtained distribution of wall pressure coefficient along the downstream of the step was in good agreement with experiment. The structure of supersonic flow over rearward facing step is complicated and contains boundary layer-shock interaction, shock-shock interaction and sec- ondary flow separation. The predictions of the corner pressure and location of the reattachment shock are more exact than the Reynolds-averaged method. Results show that detached-eddy simulation can simulate the supersonic separated flows exactly and has advantages over Reynolds-averaged method. The partition of the detached-eddy simulation can be easily achieved by using multi-block non-conform grid, which reduces the requirement of the computational resources as well.展开更多
The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameteri...The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.展开更多
Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging ...Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.展开更多
We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscal...We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.展开更多
Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing tech...Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing techniques.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)and Lunar Environment heliospheric X-ray Imager(LEXI)missions aim to obtain soft Xray images of near-Earth space thanks to their Soft X-ray Imager(SXI)instruments.While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission,the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma.To investigate the possible signatures of ion-kinetic-scale processes in soft Xray images,we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model.The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field.We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations,providing meridional and equatorial views.We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirrormode wave structures in the magnetosheath and flux transfer events(FTEs)at the magnetopause.Our results suggest that,although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers,mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images.For instance,a local increase by 30%in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12%enhancement in the line-of-sight-and time-integrated soft X-ray emissivity originating from this region.Likewise,a proton density increase by 14%in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%.These are likely conservative estimates,given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind.These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale,transient phenomena occurring on the dayside.展开更多
The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology...The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology for soil slope stability evaluation,employing Monte Carlo Simulation(MCS)and Subset Simulation(SS)with the"UPSS 3.0 Add-in"in MS-Excel.Focused on an 11.693-meter embankment with a soil slope(inclination ratio of 2H:1V),the investigation considers earthquake coefficients(kh)and pore water pressure ratios(ru)following Indian zoning requirements.The chance of slope failure showed a considerable increase as the Coefficient of Variation(COV),seismic coefficients(kh),and pore water pressure ratios(ru)experienced an escalation.The SS approach showed exceptional efficacy in calculating odds of failure that are notably low.Within computational modeling,the study optimized the worst-case scenario using ANFIS-GA,ANFIS-GWO,ANFIS-PSO,and ANFIS-BBO models.The ANFIS-PSO model exhibits exceptional accuracy(training R2=0.9011,RMSE=0.0549;testing R2=0.8968,RMSE=0.0615),emerging as the most promising.This study highlights the significance of conducting thorough risk assessments and offers practical insights into evaluating and improving the stability of soil slopes in transportation infrastructure.These findings contribute to the enhancement of safety and reliability in real-world situations.展开更多
The advantages of a flat-panel X-ray source(FPXS)make it a promising candidate for imaging applications.Accurate imaging-system modeling and projection simulation are critical for analyzing imaging performance and res...The advantages of a flat-panel X-ray source(FPXS)make it a promising candidate for imaging applications.Accurate imaging-system modeling and projection simulation are critical for analyzing imaging performance and resolving overlapping projection issues in FPXS.The conventional analytical ray-tracing approach is limited by the number of patterns and is not applicable to FPXS-projection calculations.However,the computation time of Monte Carlo(MC)simulation is independent of the size of the patterned arrays in FPXS.This study proposes two high-efficiency MC projection simulators for FPXS:a graphics processing unit(GPU)-based phase-space sampling MC(gPSMC)simulator and GPU-based fluence sampling MC(gFSMC)simulator.The two simulators comprise three components:imaging-system modeling,photon initialization,and physical-interaction simulations in the phantom.Imaging-system modeling was performed by modeling the FPXS,imaging geometry,and detector.The gPSMC simulator samples the initial photons from the phase space,whereas the gFSMC simulator performs photon initialization from the calculated energy spectrum and fluence map.The entire process of photon interaction with the geometry and arrival at the detector was simulated in parallel using multiple GPU kernels,and projections based on the two simulators were calculated.The accuracies of the two simulators were evaluated by comparing them with the conventional analytical ray-tracing approach and acquired projections,and the efficiencies were evaluated by comparing the computation time.The results of simulated and realistic experiments illustrate the accuracy and efficiency of the proposed gPSMC and gFSMC simulators in the projection calculation of various phantoms.展开更多
Dear Editor,This letter presents a biocompatible cross-shaped magnetic soft robot and investigates its deformation mode control strategy through COMSOL modeling and simulation.Magnetic soft robots offer novel avenues ...Dear Editor,This letter presents a biocompatible cross-shaped magnetic soft robot and investigates its deformation mode control strategy through COMSOL modeling and simulation.Magnetic soft robots offer novel avenues for precise treatment within intricate regions of the human body.展开更多
Objective: To explore the effectiveness of applying patient simulators combined with Internet Plus scenario simulation teaching models on intravenous (IV) infusion nursing education, and to provide scientific evidence...Objective: To explore the effectiveness of applying patient simulators combined with Internet Plus scenario simulation teaching models on intravenous (IV) infusion nursing education, and to provide scientific evidence for the implementation of advanced teaching models in future nursing education. Methods: Enrolled 60 nurses who took the IV infusion therapy training program in our hospital from January 2022 to December 2023 for research. 30 nurses who were trained in traditional teaching models from January to December 2022 were selected as the control group, and 30 nurses who were trained with simulation-based teaching models with methods including simulated patients, internet, online meetings which can be replayed and scenario simulation, etc. from January to December 2023 were selected as the experimental group. Evaluated the learning outcomes based on the Competency Inventory for Nursing Students (CINS), Problem-Solving Inventory (PSI), comprehensive learning ability, scientific research ability, and proficiency in the theoretical knowledge and practical skills of IV infusion therapy. Nursing quality, the incidence of IV infusion therapy complications and nurse satisfaction with different teaching models were also measured. Results: The scientific research ability, PSI scores, CINS scores, and comprehensive learning ability of the experimental group were better than those of the control group (P 0.05), and their assessment results of practical skills, nursing quality of IV infusion therapy during training, and satisfaction with teaching models were all better than those of the control group with statistical significance (P < 0.05). The incidence of IV infusion therapy complications in the experimental group was lower than that in the control group, indicating statistical significance (P < 0.05). Conclusions: Teaching models based on patient simulators combined with Internet Plus scenario simulation enable nursing students to learn more directly and practice at any time and in any place, and can improve their proficiency in IV infusion theoretical knowledge and skills (e.g. PICC catheterization), core competencies, problem-solving ability, comprehensive learning ability, scientific research ability and the ability to deal with complicated cases. Also, it helps provide high-quality nursing education, improve the nursing quality of IV therapy, reduce the incidence of related complications, and ensure the safety of patients with IV therapy.展开更多
The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and w...The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.展开更多
Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experim...Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.展开更多
Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear...Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.展开更多
Porous materials are widely used in the field of protection because of their excellent energy absorption characteristics.In this work,a series of polyurethane microscopic models are established and the effect of poros...Porous materials are widely used in the field of protection because of their excellent energy absorption characteristics.In this work,a series of polyurethane microscopic models are established and the effect of porosity on the shock waves is studied with classical molecular dynamics simulations.Firstly,shock Hugoniot relations for different porosities are obtained,which compare well with the experimental data.The pores collapse and form local stress wave,which results in the complex multi-wave structure of the shock wave.The microstructure analysis shows that the local stress increases and the local velocity decreases gradually during the process of pore collapse to complete compaction.Finally,it leads to stress relaxation and velocity homogenization.The shock stress peaks can be fitted with two exponential functions,and the amplitude of attenuation coefficient decreases with the increase of density.Besides,the pore collapse under shock or non-shock are discussed by the entropy increase rate of the system.The energy is dissipated mainly through the multiple interactions of the waves under shock.The energy is dissipated mainly by the friction between atoms under non-shock.展开更多
A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten st...A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.展开更多
基金supported by National Natural Science Foundation of China (Grant No.11172040)
文摘The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.
基金Project(U1534210)supported by the National Natural Science Foundation of ChinaProject(14JJ1003)supported by the Natural Science Foundation of Hunan Province,China+2 种基金Project(2015CX003)supported by the Project of Innovation-driven Plan in Central South University,ChinaProject(14JC1003)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2015T002-A)supported by the Technological Research and Development program of China Railways Cooperation
文摘In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.
文摘Detailed formulation of the detached-eddy simulation (DES) based on shear stress transport model (SST) was deduced. Based on finite volume method, detached-eddy simulation of super- sonic turbulent flow over rearward facing step was conducted. Multi-block structured non-conform grid was designed. The obtained distribution of wall pressure coefficient along the downstream of the step was in good agreement with experiment. The structure of supersonic flow over rearward facing step is complicated and contains boundary layer-shock interaction, shock-shock interaction and sec- ondary flow separation. The predictions of the corner pressure and location of the reattachment shock are more exact than the Reynolds-averaged method. Results show that detached-eddy simulation can simulate the supersonic separated flows exactly and has advantages over Reynolds-averaged method. The partition of the detached-eddy simulation can be easily achieved by using multi-block non-conform grid, which reduces the requirement of the computational resources as well.
基金supported by the National Natural Science Foundation of China(Grant Nos.42175099,42027804,42075073)the Innovative Project of Postgraduates in Jiangsu Province in 2023(Grant No.KYCX23_1319)+3 种基金supported by the National Natural Science Foundation of China(Grant No.42205080)the Natural Science Foundation of Sichuan(Grant No.2023YFS0442)the Research Fund of Civil Aviation Flight University of China(Grant No.J2022-037)supported by the National Key Scientific and Technological Infrastructure project“Earth System Science Numerical Simulator Facility”(Earth Lab)。
文摘The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.
基金supported by the National Natural Science Foundation of China(NNSFC)grants 42074202,42274196Strategic Priority Research Program of Chinese Academy of Sciences grant XDB41000000ISSI-BJ International Team Interaction between magnetic reconnection and turbulence:From the Sun to the Earth。
文摘Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.
基金support from the European Space Agency(ESA)PRODEX(PROgramme de Développement d’Expériences scientifiques)Project mission(No.PEA4000134960)Partial funding was provided by the Romanian Ministry of Research,Innovation and Digitalization under Romanian National Core Program LAPLAS VII(Contract No.30N/2023)+2 种基金the Belgian Solar-Terrestrial Centre of Excellencesupported by the project Belgian Research Action through Interdisciplinary Networks(BRAIN-BE)2.0(Grant No.B2/223/P1/PLATINUM)funded by the Belgian Office for Research(BELSPO)partially supported by a grant from the Romanian Ministry of Education and Research(CNCS-UEFISCDI,Project No.PN-III-P1-1.1TE-2021-0102)。
文摘We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.
基金the European Research Council for starting grant 200141-QuESpace,with which the Vlasiator model was developedconsolidator grant 682068-PRESTISSIMO awarded for further development of Vlasiator and its use in scientific investigations+4 种基金Academy of Finland grant numbers 338629-AERGELC’H,339756-KIMCHI,336805-FORESAIL,and 335554-ICT-SUNVACThe Academy of Finland also supported this work through the PROFI4 grant(grant number 3189131)support from the NASA grants,80NSSC20K1670 and 80MSFC20C0019the NASA GSFC FY23 IRADHIF funds。
文摘Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing techniques.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)and Lunar Environment heliospheric X-ray Imager(LEXI)missions aim to obtain soft Xray images of near-Earth space thanks to their Soft X-ray Imager(SXI)instruments.While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission,the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma.To investigate the possible signatures of ion-kinetic-scale processes in soft Xray images,we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model.The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field.We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations,providing meridional and equatorial views.We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirrormode wave structures in the magnetosheath and flux transfer events(FTEs)at the magnetopause.Our results suggest that,although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers,mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images.For instance,a local increase by 30%in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12%enhancement in the line-of-sight-and time-integrated soft X-ray emissivity originating from this region.Likewise,a proton density increase by 14%in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%.These are likely conservative estimates,given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind.These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale,transient phenomena occurring on the dayside.
文摘The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology for soil slope stability evaluation,employing Monte Carlo Simulation(MCS)and Subset Simulation(SS)with the"UPSS 3.0 Add-in"in MS-Excel.Focused on an 11.693-meter embankment with a soil slope(inclination ratio of 2H:1V),the investigation considers earthquake coefficients(kh)and pore water pressure ratios(ru)following Indian zoning requirements.The chance of slope failure showed a considerable increase as the Coefficient of Variation(COV),seismic coefficients(kh),and pore water pressure ratios(ru)experienced an escalation.The SS approach showed exceptional efficacy in calculating odds of failure that are notably low.Within computational modeling,the study optimized the worst-case scenario using ANFIS-GA,ANFIS-GWO,ANFIS-PSO,and ANFIS-BBO models.The ANFIS-PSO model exhibits exceptional accuracy(training R2=0.9011,RMSE=0.0549;testing R2=0.8968,RMSE=0.0615),emerging as the most promising.This study highlights the significance of conducting thorough risk assessments and offers practical insights into evaluating and improving the stability of soil slopes in transportation infrastructure.These findings contribute to the enhancement of safety and reliability in real-world situations.
文摘The advantages of a flat-panel X-ray source(FPXS)make it a promising candidate for imaging applications.Accurate imaging-system modeling and projection simulation are critical for analyzing imaging performance and resolving overlapping projection issues in FPXS.The conventional analytical ray-tracing approach is limited by the number of patterns and is not applicable to FPXS-projection calculations.However,the computation time of Monte Carlo(MC)simulation is independent of the size of the patterned arrays in FPXS.This study proposes two high-efficiency MC projection simulators for FPXS:a graphics processing unit(GPU)-based phase-space sampling MC(gPSMC)simulator and GPU-based fluence sampling MC(gFSMC)simulator.The two simulators comprise three components:imaging-system modeling,photon initialization,and physical-interaction simulations in the phantom.Imaging-system modeling was performed by modeling the FPXS,imaging geometry,and detector.The gPSMC simulator samples the initial photons from the phase space,whereas the gFSMC simulator performs photon initialization from the calculated energy spectrum and fluence map.The entire process of photon interaction with the geometry and arrival at the detector was simulated in parallel using multiple GPU kernels,and projections based on the two simulators were calculated.The accuracies of the two simulators were evaluated by comparing them with the conventional analytical ray-tracing approach and acquired projections,and the efficiencies were evaluated by comparing the computation time.The results of simulated and realistic experiments illustrate the accuracy and efficiency of the proposed gPSMC and gFSMC simulators in the projection calculation of various phantoms.
基金supported by NSFC(62273019,52072015,12332019,U20A20390)the 111 Project(B13003)。
文摘Dear Editor,This letter presents a biocompatible cross-shaped magnetic soft robot and investigates its deformation mode control strategy through COMSOL modeling and simulation.Magnetic soft robots offer novel avenues for precise treatment within intricate regions of the human body.
文摘Objective: To explore the effectiveness of applying patient simulators combined with Internet Plus scenario simulation teaching models on intravenous (IV) infusion nursing education, and to provide scientific evidence for the implementation of advanced teaching models in future nursing education. Methods: Enrolled 60 nurses who took the IV infusion therapy training program in our hospital from January 2022 to December 2023 for research. 30 nurses who were trained in traditional teaching models from January to December 2022 were selected as the control group, and 30 nurses who were trained with simulation-based teaching models with methods including simulated patients, internet, online meetings which can be replayed and scenario simulation, etc. from January to December 2023 were selected as the experimental group. Evaluated the learning outcomes based on the Competency Inventory for Nursing Students (CINS), Problem-Solving Inventory (PSI), comprehensive learning ability, scientific research ability, and proficiency in the theoretical knowledge and practical skills of IV infusion therapy. Nursing quality, the incidence of IV infusion therapy complications and nurse satisfaction with different teaching models were also measured. Results: The scientific research ability, PSI scores, CINS scores, and comprehensive learning ability of the experimental group were better than those of the control group (P 0.05), and their assessment results of practical skills, nursing quality of IV infusion therapy during training, and satisfaction with teaching models were all better than those of the control group with statistical significance (P < 0.05). The incidence of IV infusion therapy complications in the experimental group was lower than that in the control group, indicating statistical significance (P < 0.05). Conclusions: Teaching models based on patient simulators combined with Internet Plus scenario simulation enable nursing students to learn more directly and practice at any time and in any place, and can improve their proficiency in IV infusion theoretical knowledge and skills (e.g. PICC catheterization), core competencies, problem-solving ability, comprehensive learning ability, scientific research ability and the ability to deal with complicated cases. Also, it helps provide high-quality nursing education, improve the nursing quality of IV therapy, reduce the incidence of related complications, and ensure the safety of patients with IV therapy.
基金supported by the National Natural Science Foundation of China (Grant Nos. 92158204, 41506001 and 42076019)a Project supported by the Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (Grant No. 311021005)。
文摘The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.
基金This work was supported by the National Natural Science Foundation of China(No.11575080)the National Natural Science Foundation of Hunan Province,China(No.2022JJ30482)the Hunan Provincial Innovation Foundation for Postgraduates(No.QL20220206).
文摘Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.51775077 and 51909023)。
文摘Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.
基金financial support from National Natural Science Foundation of China(Grant No.12172325)。
文摘Porous materials are widely used in the field of protection because of their excellent energy absorption characteristics.In this work,a series of polyurethane microscopic models are established and the effect of porosity on the shock waves is studied with classical molecular dynamics simulations.Firstly,shock Hugoniot relations for different porosities are obtained,which compare well with the experimental data.The pores collapse and form local stress wave,which results in the complex multi-wave structure of the shock wave.The microstructure analysis shows that the local stress increases and the local velocity decreases gradually during the process of pore collapse to complete compaction.Finally,it leads to stress relaxation and velocity homogenization.The shock stress peaks can be fitted with two exponential functions,and the amplitude of attenuation coefficient decreases with the increase of density.Besides,the pore collapse under shock or non-shock are discussed by the entropy increase rate of the system.The energy is dissipated mainly through the multiple interactions of the waves under shock.The energy is dissipated mainly by the friction between atoms under non-shock.
基金the National Key R&D Program(No.2023YFB3709900)the National Natural Science Foundation of China(Nos.U22A20171 and 52104343)the High Steel Central(HSC)at North China University of Science and Technology and Yanshan Univ ersity,China。
文摘A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.