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.展开更多
To quantify the relative contributions of Arctic sea ice and unforced atmospheric internal variability to the “warm Arctic, cold East Asia”(WACE) teleconnection, this study analyses three sets of large-ensemble simu...To quantify the relative contributions of Arctic sea ice and unforced atmospheric internal variability to the “warm Arctic, cold East Asia”(WACE) teleconnection, this study analyses three sets of large-ensemble simulations carried out by the Norwegian Earth System Model with a coupled atmosphere–land surface model, forced by seasonal sea ice conditions from preindustrial, present-day, and future periods. Each ensemble member within the same set uses the same forcing but with small perturbations to the atmospheric initial state. Hence, the difference between the present-day(or future) ensemble mean and the preindustrial ensemble mean provides the ice-loss-induced response, while the difference of the individual members within the present-day(or future) set is the effect of atmospheric internal variability. Results indicate that both present-day and future sea ice loss can force a negative phase of the Arctic Oscillation with a WACE pattern in winter. The magnitude of ice-induced Arctic warming is over four(ten) times larger than the ice-induced East Asian cooling in the present-day(future) experiment;the latter having a magnitude that is about 30% of the observed cooling. Sea ice loss contributes about 60%(80%) to the Arctic winter warming in the present-day(future) experiment. Atmospheric internal variability can also induce a WACE pattern with comparable magnitudes between the Arctic and East Asia. Ice-lossinduced East Asian cooling can easily be masked by atmospheric internal variability effects because random atmospheric internal variability may induce a larger magnitude warming. The observed WACE pattern occurs as a result of both Arctic sea ice loss and atmospheric internal variability, with the former dominating Arctic warming and the latter dominating East Asian cooling.展开更多
A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDO...A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDOT).The simulation results show that the coating of primary NMC particles significantly reduces the stress generation by efficiently accommodating the volume change associated with the lithium diffusion,and the coating layer plays roles both as a cushion against the volume change and a channel for the lithium transport,promoting the lithium distribution across the secondary particles more homogeneously.Besides,the lower stiffness,higher ionic conductivity,and larger thickness of the coating layer improve the stress mitigation.This paper provides a mathematical framework for calculating the chemo-mechanical responses of anisotropic electrode materials and fundamental insights into how the coating of NMC active particles mitigates stress levels.展开更多
A mesoscale convective system(MCS) occurred over the East China coastal provinces and the East China Sea on 30April 2021, producing damaging surface winds near the coastal city Nantong with observed speeds reaching 45...A mesoscale convective system(MCS) occurred over the East China coastal provinces and the East China Sea on 30April 2021, producing damaging surface winds near the coastal city Nantong with observed speeds reaching 45 m s^(–1). A simulation using the Weather Research and Forecasting model with a 1.5-km grid spacing generally reproduces the development and subsequent organization of this convective system into an MCS, with an eastward protruding bow segment over the sea. In the simulation, an east-west-oriented high wind swath is generated behind the gust front of the MCS. Descending dry rear-to-front inflows behind the bow and trailing gust front are found to feed the downdrafts in the main precipitation regions. The inflows help to establish spreading cold outflows and enhance the downdrafts through evaporative cooling. Meanwhile, front-to-rear inflows from the south are present, associated with severely rearward-tilted updrafts initially forming over the gust front. Such inflows descend behind(north of) the gust front, significantly enhancing downdrafts and near-surface winds within the cold pool. Consistently, calculated trajectories show that these parcels that contribute to the derecho originate primarily from the region ahead(south) of the east-west-oriented gust front, and dry southwesterly flows in the low-to-middle levels contribute to strong downdrafts within the MCS. Moreover, momentum budget analyses reveal that a large westward-directed horizontal pressure gradient force within the simulated cold pool produced rapid flow acceleration towards Nantong. The analyses enrich the understanding of damaging wind characteristics over coastal East China and will prove helpful to operational forecasters.展开更多
MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDI...MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDIA,AMD,Intel,and Apple GPUs.Moreover,MicroMagnetic.jl supports Monte Carlo simulations for atomistic models and implements the nudged-elastic-band method for energy barrier computations.With built-in support for double and single precision modes and a design allowing easy extensibility to add new features,MicroMagnetic.jl provides a versatile toolset for researchers in micromagnetics and atomistic simulations.展开更多
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 effect of undercooling DT and the interface energy anisotropy parameter e4 on the shape of the equiaxed dendritic tip has been investigated by using a quantitative phase-field model for solidification of binary al...The effect of undercooling DT and the interface energy anisotropy parameter e4 on the shape of the equiaxed dendritic tip has been investigated by using a quantitative phase-field model for solidification of binary alloys.It was found that the tip radius r increases and the tip shape amplitude coefficient A4 decreases with the increase of the fitting range for all cases.The dendrite tip shape selection parameter sdecreases and then stabilizes with the increase of the fitting range,and sincreases with the increase of e4.The relationship between sand e4 follows a power-law function sµea 4,and a is independent of DT but dependent on the fitting range.Numerical results demonstrate that the predicted sis consistent with the curve of microscopic solvability theory(MST)for e4<0.02,and sobtained from our phase-field simulations is sensitive to the undercooling when e4 is fixed.展开更多
Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pr...Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.展开更多
Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods...Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.展开更多
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.展开更多
The Makran accretionary wedge has the smallest subduction angle among any accretionary prism in the world. The factors controlling the spacing and morphological development of its deep thrust faults, as well as the fo...The Makran accretionary wedge has the smallest subduction angle among any accretionary prism in the world. The factors controlling the spacing and morphological development of its deep thrust faults, as well as the formation mechanism of shallow normal faults, remain unclear. Meanwhile, the factors affecting the continuity of plane faults must be comprehensively discussed. Clarifying the development characteristics and deformation mechanisms of the Makran accretionary wedge is crucial to effectively guide the exploration of gas hydrate deposits in the area. This study aims to interpret seismic data to identify typical structures in the Makran accretionary wedge, including deep imbricate thrust faults, shallow and small normal faults, wedge-shaped piggyback basins, mud diapirs with fuzzy and disorderly characteristics of reflection, décollements with a northward tilt of 1° – 2°, and large seamounts. Physical simulation-based experiments are performed to comprehensively analyze the results of the plane, section, and slices of the wedge. Results reveal that the distances between and shapes of thrust faults in the deep parts of the Makran accretionary wedge are controlled by the bottom décollement. The uplift of the thrust fault-related folds and the upwelling of the mud diapirs primarily contribute to the formation of small normal faults in the shallow part of the area. The mud diapirs originate from plastic material at the bottom, while those that have developed in the area near the trench are larger. Seamounts and mud diapirs break the continuity of fault plane distribution.展开更多
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.展开更多
The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To re...The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.展开更多
Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectil...Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.展开更多
The effects of impurities on ion temperature gradient(ITG)driven turbulence transport in tokamak core plasmas are investigated numerically via global simulations of microturbulence with carbon impurities and adiabatic...The effects of impurities on ion temperature gradient(ITG)driven turbulence transport in tokamak core plasmas are investigated numerically via global simulations of microturbulence with carbon impurities and adiabatic electrons.The simulations use an extended fluid code(ExFC)based on a four-field gyro-Landau-fluid(GLF)model.The multispecies form of the normalized GLF equations is presented,which guarantees the self-consistent evolution of both bulk ions and impurities.With parametric profiles of the cyclone base case,well-benchmarked ExFC is employed to perform simulations focusing on different impurity density profiles.For a fixed temperature profile,it is found that the turbulent heat diffusivity of bulk ions in a quasi-steady state is usually lower than that without impurities,which is contrary to the linear and quasilinear predictions.The evolutions of the temperature gradient and heat diffusivity exhibit a fast relaxation process,indicating that the destabilization of the outwardly peaked impurity profile is a transient state response.Furthermore,the impurity effects from different profiles can obviously influence the nonlinear critical temperature gradient,which is likely to be dominated by linear effects.These results suggest that the improvement in plasma confinement could be attributed to the impurities,most likely through adjusting both heat diffusivity and the critical temperature gradient.展开更多
The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Re...The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary ...As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.展开更多
基金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.
基金supported by the Chinese-Norwegian Collaboration Projects within Climate Systems jointly funded by the National Key Research and Development Program of China (Grant No.2022YFE0106800)the Research Council of Norway funded project MAPARC (Grant No.328943)+2 种基金the support from the Research Council of Norway funded project BASIC (Grant No.325440)the Horizon 2020 project APPLICATE (Grant No.727862)High-performance computing and storage resources were performed on resources provided by Sigma2 - the National Infrastructure for High-Performance Computing and Data Storage in Norway (through projects NS8121K,NN8121K,NN2345K,NS2345K,NS9560K,NS9252K,and NS9034K)。
文摘To quantify the relative contributions of Arctic sea ice and unforced atmospheric internal variability to the “warm Arctic, cold East Asia”(WACE) teleconnection, this study analyses three sets of large-ensemble simulations carried out by the Norwegian Earth System Model with a coupled atmosphere–land surface model, forced by seasonal sea ice conditions from preindustrial, present-day, and future periods. Each ensemble member within the same set uses the same forcing but with small perturbations to the atmospheric initial state. Hence, the difference between the present-day(or future) ensemble mean and the preindustrial ensemble mean provides the ice-loss-induced response, while the difference of the individual members within the present-day(or future) set is the effect of atmospheric internal variability. Results indicate that both present-day and future sea ice loss can force a negative phase of the Arctic Oscillation with a WACE pattern in winter. The magnitude of ice-induced Arctic warming is over four(ten) times larger than the ice-induced East Asian cooling in the present-day(future) experiment;the latter having a magnitude that is about 30% of the observed cooling. Sea ice loss contributes about 60%(80%) to the Arctic winter warming in the present-day(future) experiment. Atmospheric internal variability can also induce a WACE pattern with comparable magnitudes between the Arctic and East Asia. Ice-lossinduced East Asian cooling can easily be masked by atmospheric internal variability effects because random atmospheric internal variability may induce a larger magnitude warming. The observed WACE pattern occurs as a result of both Arctic sea ice loss and atmospheric internal variability, with the former dominating Arctic warming and the latter dominating East Asian cooling.
基金the National Research Foundation of Korea(Nos.2018R1A5A7023490 and 2022R1A2C1003003)。
文摘A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDOT).The simulation results show that the coating of primary NMC particles significantly reduces the stress generation by efficiently accommodating the volume change associated with the lithium diffusion,and the coating layer plays roles both as a cushion against the volume change and a channel for the lithium transport,promoting the lithium distribution across the secondary particles more homogeneously.Besides,the lower stiffness,higher ionic conductivity,and larger thickness of the coating layer improve the stress mitigation.This paper provides a mathematical framework for calculating the chemo-mechanical responses of anisotropic electrode materials and fundamental insights into how the coating of NMC active particles mitigates stress levels.
基金primarily supported by the Ministry of Science and Technology of the People's Republic of China (MOST)(Grant No. 2018YFC1507303)National Natural Science Foundation of China (Grant Nos. 419505044,41941007, and 42230607)+1 种基金by the Talent Research Start-Up Fund of Nanjing University of Aeronautics and Astronautics(Grant No. 1007-90YAH22046)supported by The High Performance Computing Platform of Nanjing University of Aeronautics and Astronautics。
文摘A mesoscale convective system(MCS) occurred over the East China coastal provinces and the East China Sea on 30April 2021, producing damaging surface winds near the coastal city Nantong with observed speeds reaching 45 m s^(–1). A simulation using the Weather Research and Forecasting model with a 1.5-km grid spacing generally reproduces the development and subsequent organization of this convective system into an MCS, with an eastward protruding bow segment over the sea. In the simulation, an east-west-oriented high wind swath is generated behind the gust front of the MCS. Descending dry rear-to-front inflows behind the bow and trailing gust front are found to feed the downdrafts in the main precipitation regions. The inflows help to establish spreading cold outflows and enhance the downdrafts through evaporative cooling. Meanwhile, front-to-rear inflows from the south are present, associated with severely rearward-tilted updrafts initially forming over the gust front. Such inflows descend behind(north of) the gust front, significantly enhancing downdrafts and near-surface winds within the cold pool. Consistently, calculated trajectories show that these parcels that contribute to the derecho originate primarily from the region ahead(south) of the east-west-oriented gust front, and dry southwesterly flows in the low-to-middle levels contribute to strong downdrafts within the MCS. Moreover, momentum budget analyses reveal that a large westward-directed horizontal pressure gradient force within the simulated cold pool produced rapid flow acceleration towards Nantong. The analyses enrich the understanding of damaging wind characteristics over coastal East China and will prove helpful to operational forecasters.
基金supported by the National Key R&D Program of China(Grant No.2022YFA1403603)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDB33030100)+2 种基金the National Natural Science Fund for Distinguished Young Scholar(Grant No.52325105)the National Natural Science Foundation of China(Grant Nos.12374098,11974021,and 12241406)the CAS Project for Young Scientists in Basic Research(Grant No.YSBR-084).
文摘MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDIA,AMD,Intel,and Apple GPUs.Moreover,MicroMagnetic.jl supports Monte Carlo simulations for atomistic models and implements the nudged-elastic-band method for energy barrier computations.With built-in support for double and single precision modes and a design allowing easy extensibility to add new features,MicroMagnetic.jl provides a versatile toolset for researchers in micromagnetics and atomistic simulations.
文摘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 National Key Research and De-velopment Program of China(Grant No.2021YFB3502600)Shenzhen Science and Technology Program(Grant No.JCYJ20220530161813029).
文摘The effect of undercooling DT and the interface energy anisotropy parameter e4 on the shape of the equiaxed dendritic tip has been investigated by using a quantitative phase-field model for solidification of binary alloys.It was found that the tip radius r increases and the tip shape amplitude coefficient A4 decreases with the increase of the fitting range for all cases.The dendrite tip shape selection parameter sdecreases and then stabilizes with the increase of the fitting range,and sincreases with the increase of e4.The relationship between sand e4 follows a power-law function sµea 4,and a is independent of DT but dependent on the fitting range.Numerical results demonstrate that the predicted sis consistent with the curve of microscopic solvability theory(MST)for e4<0.02,and sobtained from our phase-field simulations is sensitive to the undercooling when e4 is fixed.
基金supported by the National Natural Science Foundation of China(Project Nos.12272270,11972261).
文摘Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.
基金National Natural Science Foundation of China(No.11575080)Hunan Provincial Natural Science Foundation of China(No.2022JJ30482)Hunan Provincial Innovation Foundation for Postgraduate(No.QL20220206).
文摘Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.
基金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.
基金funded by the National Natural Science Foundation of China(No.42076069).
文摘The Makran accretionary wedge has the smallest subduction angle among any accretionary prism in the world. The factors controlling the spacing and morphological development of its deep thrust faults, as well as the formation mechanism of shallow normal faults, remain unclear. Meanwhile, the factors affecting the continuity of plane faults must be comprehensively discussed. Clarifying the development characteristics and deformation mechanisms of the Makran accretionary wedge is crucial to effectively guide the exploration of gas hydrate deposits in the area. This study aims to interpret seismic data to identify typical structures in the Makran accretionary wedge, including deep imbricate thrust faults, shallow and small normal faults, wedge-shaped piggyback basins, mud diapirs with fuzzy and disorderly characteristics of reflection, décollements with a northward tilt of 1° – 2°, and large seamounts. Physical simulation-based experiments are performed to comprehensively analyze the results of the plane, section, and slices of the wedge. Results reveal that the distances between and shapes of thrust faults in the deep parts of the Makran accretionary wedge are controlled by the bottom décollement. The uplift of the thrust fault-related folds and the upwelling of the mud diapirs primarily contribute to the formation of small normal faults in the shallow part of the area. The mud diapirs originate from plastic material at the bottom, while those that have developed in the area near the trench are larger. Seamounts and mud diapirs break the continuity of fault plane distribution.
基金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.
基金co-supported by the National Natural Science Foundation of China under Grants(42130101,42074007)the Fundamental Research Funds for the Central Universities under Grants 2042023kf1035。
文摘The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.
基金supported by the National Natural Science Foundation of China(No.12102256).
文摘Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.
基金supported by National Natural Science Foundation of China(Nos.U1967206 and 12275071)National Key R&D Program of China(No.2017YFE0301201)。
文摘The effects of impurities on ion temperature gradient(ITG)driven turbulence transport in tokamak core plasmas are investigated numerically via global simulations of microturbulence with carbon impurities and adiabatic electrons.The simulations use an extended fluid code(ExFC)based on a four-field gyro-Landau-fluid(GLF)model.The multispecies form of the normalized GLF equations is presented,which guarantees the self-consistent evolution of both bulk ions and impurities.With parametric profiles of the cyclone base case,well-benchmarked ExFC is employed to perform simulations focusing on different impurity density profiles.For a fixed temperature profile,it is found that the turbulent heat diffusivity of bulk ions in a quasi-steady state is usually lower than that without impurities,which is contrary to the linear and quasilinear predictions.The evolutions of the temperature gradient and heat diffusivity exhibit a fast relaxation process,indicating that the destabilization of the outwardly peaked impurity profile is a transient state response.Furthermore,the impurity effects from different profiles can obviously influence the nonlinear critical temperature gradient,which is likely to be dominated by linear effects.These results suggest that the improvement in plasma confinement could be attributed to the impurities,most likely through adjusting both heat diffusivity and the critical temperature gradient.
基金supported by the National MCF Energy R&D Program of China(No.2019YFE03060000)National Natural Science Foundation of China(Nos.12005063,12375215 and 12175034)the Collaborative Innovation Program of Hefei Science Center,CAS(No.2022HSC-CIP008).
文摘The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.
基金Project supported by the National Natural Science Foundation of China(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
基金supported by State Grid Ningxia Electric Power Co.,Ltd.under Grant 5229CG220006Natural Science Foundation of Ningxia Province under Grant 2022AAC03629.
文摘As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.
