Accurate simulation of the evolution of freak waves by the wave phase focusing method requires accurate linear and nonlinear properties,especially in deep-water conditions.In this paper,we analyze the ability to simul...Accurate simulation of the evolution of freak waves by the wave phase focusing method requires accurate linear and nonlinear properties,especially in deep-water conditions.In this paper,we analyze the ability to simulate deep-water focused waves of a two-layer Boussinesq-type(BT)model,which has been shown to have excellent linear and nonlinear performance.To further improve the numerical accuracy and stability,the internal wavegenerated method is introduced into the two-layer Boussinesq-type model.Firstly,the sensitivity of the numerical results to the grid resolution is analyzed to verify the convergence of the model;secondly,the focused wave propagating in two opposite directions is simulated to prove the symmetry of the numerical results and the feasibility of the internal wave-generated method;thirdly,the limiting focused wave condition is simulated to compare and analyze the wave surface and the horizontal velocity of the profile at the focusing position,which is in good agreement with the measured values.Meanwhile the simulation of focused waves in very deep waters agrees well with the measured values,which further demonstrates the capability of the two-layer BT model in simulating focused waves in deep waters.展开更多
A new simulation model for the development of gas condensate reservoirs is introduced based on the influence that phase change,non-Darcy flow,and capillary pressure have on the production of gas condensates.The model ...A new simulation model for the development of gas condensate reservoirs is introduced based on the influence that phase change,non-Darcy flow,and capillary pressure have on the production of gas condensates.The model predicts well performance,including bottom-hole pressure,oil/gas production rate,oil/gas recovery,gaseoil ratio,and the change in produced fluid composition.It also calculates dynamic characters,such as the change of pressure field and oil/gas saturation field during the development of gas condensate reservoirs.The model is applicable to different boundary conditions(both constant-pressure and sealed boundary)and different production modes(both constant-pressure and constant-volume production modes).Model validation attempted using numerical simulation results for sealed boundary conditions with constant-pressure production mode has shown a relatively good match,proving its validity.For constant-pressure boundary conditions with constant-volume production mode,four stages are defined according to the dynamic behavior of production performance in the development of gas condensate reservoirs.展开更多
Air brake systems are critical equipment for railway trains, which affects the running safety of the trains significantly. To study air braking characteristics of long freight trains, an approach for simulating air br...Air brake systems are critical equipment for railway trains, which affects the running safety of the trains significantly. To study air braking characteristics of long freight trains, an approach for simulating air brake systems based on fuid dynamics theory was proposed. The structures and working mechanisms of locomotive and wagon air brakes are introduced, and mathematical models of the pipes, brake valves, reservoirs or chambers, cylinders, etc., are presented.Besides, the dynamic motions of parts in the main valve are considered. The simulation model of the whole air brake system is then formulated, and the solving method based on the finite-difference method is used. New efficient pipe boundary conditions without iterations are developed for brake pipes and branch pipes, which can achieve higher computational efficiency. The proposed approach for simulating the air brake system is validated by comparing with published measured data. Simulation results of different train formations indicate that models that consider the dynamic behavior of brake pipes are recommended for predicting the characteristics of long trains under service braking conditions.展开更多
Stochastic finite-fault simulations are effective for simulating ground motions and are widely used in engineering to determine the impacts of ground motion and develop relevant predictive equations.In this study,the ...Stochastic finite-fault simulations are effective for simulating ground motions and are widely used in engineering to determine the impacts of ground motion and develop relevant predictive equations.In this study,the source,path,and site amplification coefficient of western Sichuan Province,China,and stochastic finite-fault simulations were used to simulate the acceleration time series,Fourier amplitude spectra,and 5%damped response spectra of 28 strong-motion stations with rupture distances within 300 km of the 2022 MS6.8 Luding earthquake.The simulation results of 14 stations at rupture distances of 45-185 km match the observation.However,the simulation results of 3 near-and 6 far-field stations at rupture distances of 12-36 km and 222-286 km,respectively,were obviously deviated from the observations.Simulation results of the near-field stations are larger than the observations at high frequencies(>6 Hz).The discrepancy likely comes from the nonlinear site effect of near-field stations,which reduced the site amplification at high frequencies.Simulation result of the far-field stations is smaller than the observation at frequencies above 1 Hz.As these stations are located close to the Longmenshan Fault Zone(LFZ),thus,we obtained a new quality factor(Q)from data of historical events and stations located around LFZ.Using the new Q value,the discrepancies of the high-frequency simulation results of the far-field stations were corrected.This result indicated that the laterally varying Q values can be used to address the impact of strong crustal lateral heterogeneity on simulation.展开更多
Considering the complex topographic forcing and large cryosphere concentration,the present study utilized the polar-optimized WRF model(Polar WRF)to conduct downscaling simulations over the Qinghai-Tibet Plateau(TP)an...Considering the complex topographic forcing and large cryosphere concentration,the present study utilized the polar-optimized WRF model(Polar WRF)to conduct downscaling simulations over the Qinghai-Tibet Plateau(TP)and its surrounding regions.Multi-group experiments with the 10 km horizontal resolution are used to evaluate the modeling of precipitation.Firstly,on the basis of the model ground surface properties upgrade and the optimized Noah-MP,the“better-performing”configuration suite for modeling precipitation is comprehensively examined.Various model parameters such as nudging options,five cumulus parameterization schemes,two planetary boundary layer schemes,and six microphysics schemes are investigated to further refine the Polar WRF configuration.Moreover,the precipitation simulation for a full calendar year is compared with multiple reanalyses and observations.The simulations demonstrate that the Polar WRF model successfully captures the general features of precipitation over this region and is sensitive to model parameters.Based on the results,it is recommended to use grid nudging with q intensity coefficient of 0.0002,the multi-scale kain-fritsch cumulus parameterization,the Yonsei University boundary layer scheme,and the Morrison 2-mom microphysics with reduced default droplet concentration value of 100 cm-3.Overall,the model performance is better than the ERA-interim and TRMM 3b42.It is comparable to,and in some cases slightly better than,the CRA-Land,especially in the prediction for the western part of the plateau where in situ observations are limited,and the cryosphere-atmosphere interaction is more pronounced.展开更多
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.展开更多
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.展开更多
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.展开更多
The response of biological phantoms against high velocity impact is actively sought for applications in defense,space,soft robotics and sensing.Towards this end,we study the ballistic response of silicone based skin s...The response of biological phantoms against high velocity impact is actively sought for applications in defense,space,soft robotics and sensing.Towards this end,we study the ballistic response of silicone based skin simulant against fragment impact.Using a pneumatic gas gun setup,six chisel-nosed and three regular shaped(sphere,cylinder,and cube)fragments were impacted on the skin simulant.The resulting skin simulant response was studied in terms of ballistic limit velocities,energy densities,failure pattern,and the mechanics of interaction.The results indicate that the shape of the fragment affects the ballistic limit velocities.The ballistic limit velocities,energy densities of the chisel-nosed fragment simulating projectiles were relatively insensitive to the size(mass),except for the smallest(0.16 g)and largest(2.79 g)chisel-nosed fragment.For the same size(1 g),ballistic limit velocities and failure are dependent on the shape of the fragment.The skin simulant failed by combined plugging and elastic hole enlargement.Failure in the spherical fragment was dominated by the elastic hole enlargement,whereas plugging failure was dominant in all other fragments.The spherical,cylindrical,and chisel-nosed fragments created circular cavities,and the cubical fragment created a square cavity.In the case of the spherical fragment,slipping of the fragment within the skin simulant was seen.Cubical fragments created lateral cracks emanating from the corners of the square cavity.Interestingly,for all the fragments,the maximum deformation corresponding to the perforation was lower than the non-perforation indicating rate dependent,stress driven failure.The maximum deformation was also dependent on the shape of the fragment.Overall,these results provide unique insights into the mechanical response of a soft simulant against ballistic impact.Results have utility in the calibration and validation of computational models,design of personal protective equipment,and antipersonnel systems.展开更多
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.展开更多
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.展开更多
The spin-transfer-torque(STT)magnetic tunneling junction(MTJ)device is one of the prominent candidates for spintronic logic circuit and neuromorphic computing.Therefore,building a simulation framework of hybrid STT-MT...The spin-transfer-torque(STT)magnetic tunneling junction(MTJ)device is one of the prominent candidates for spintronic logic circuit and neuromorphic computing.Therefore,building a simulation framework of hybrid STT-MTJ/CMOS(complementary metal-oxide-semiconductor)circuits is of great value for designing a new kind of computing paradigm based on the spintronic devices.In this work,we develop a simulation framework of hybrid STT-MTJ/CMOS circuits based on MATLAB/Simulink,which is mainly composed of a physics-based STT-MTJ model,a controlled resistor,and a current sensor.In the proposed framework,the STT-MTJ model,based on the Landau-Lifshitz-Gilbert-Slonczewsk(LLGS)equation,is implemented using the MATLAB script.The proposed simulation framework is modularized design,with the advantage of simple-to-use and easy-to-expand.To prove the effectiveness of the proposed framework,the STT-MTJ model is benchmarked with experimental results.Furthermore,the pre-charge sense amplifier(PCSA)circuit consisting of two STT-MTJ devices is validated and the electrical coupling of two spin-torque oscillators is simulated.The results demonstrate the effectiveness of our simulation framework.展开更多
基金The National Natural Science Foundation under contract Nos 52171247,51779022,52071057,and 51709054.
文摘Accurate simulation of the evolution of freak waves by the wave phase focusing method requires accurate linear and nonlinear properties,especially in deep-water conditions.In this paper,we analyze the ability to simulate deep-water focused waves of a two-layer Boussinesq-type(BT)model,which has been shown to have excellent linear and nonlinear performance.To further improve the numerical accuracy and stability,the internal wavegenerated method is introduced into the two-layer Boussinesq-type model.Firstly,the sensitivity of the numerical results to the grid resolution is analyzed to verify the convergence of the model;secondly,the focused wave propagating in two opposite directions is simulated to prove the symmetry of the numerical results and the feasibility of the internal wave-generated method;thirdly,the limiting focused wave condition is simulated to compare and analyze the wave surface and the horizontal velocity of the profile at the focusing position,which is in good agreement with the measured values.Meanwhile the simulation of focused waves in very deep waters agrees well with the measured values,which further demonstrates the capability of the two-layer BT model in simulating focused waves in deep waters.
基金supported by a National Science and Technology Major Project(2016ZX05048003).
文摘A new simulation model for the development of gas condensate reservoirs is introduced based on the influence that phase change,non-Darcy flow,and capillary pressure have on the production of gas condensates.The model predicts well performance,including bottom-hole pressure,oil/gas production rate,oil/gas recovery,gaseoil ratio,and the change in produced fluid composition.It also calculates dynamic characters,such as the change of pressure field and oil/gas saturation field during the development of gas condensate reservoirs.The model is applicable to different boundary conditions(both constant-pressure and sealed boundary)and different production modes(both constant-pressure and constant-volume production modes).Model validation attempted using numerical simulation results for sealed boundary conditions with constant-pressure production mode has shown a relatively good match,proving its validity.For constant-pressure boundary conditions with constant-volume production mode,four stages are defined according to the dynamic behavior of production performance in the development of gas condensate reservoirs.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51825504, 51735012, and 52072317)
文摘Air brake systems are critical equipment for railway trains, which affects the running safety of the trains significantly. To study air braking characteristics of long freight trains, an approach for simulating air brake systems based on fuid dynamics theory was proposed. The structures and working mechanisms of locomotive and wagon air brakes are introduced, and mathematical models of the pipes, brake valves, reservoirs or chambers, cylinders, etc., are presented.Besides, the dynamic motions of parts in the main valve are considered. The simulation model of the whole air brake system is then formulated, and the solving method based on the finite-difference method is used. New efficient pipe boundary conditions without iterations are developed for brake pipes and branch pipes, which can achieve higher computational efficiency. The proposed approach for simulating the air brake system is validated by comparing with published measured data. Simulation results of different train formations indicate that models that consider the dynamic behavior of brake pipes are recommended for predicting the characteristics of long trains under service braking conditions.
基金supported by the Special Fund of the Institute of Geophysics,China Earthquake Administration(Nos.DQJB2B28 and DQJB22X09)the National Natural Science Foundation of China(No.52192675).
文摘Stochastic finite-fault simulations are effective for simulating ground motions and are widely used in engineering to determine the impacts of ground motion and develop relevant predictive equations.In this study,the source,path,and site amplification coefficient of western Sichuan Province,China,and stochastic finite-fault simulations were used to simulate the acceleration time series,Fourier amplitude spectra,and 5%damped response spectra of 28 strong-motion stations with rupture distances within 300 km of the 2022 MS6.8 Luding earthquake.The simulation results of 14 stations at rupture distances of 45-185 km match the observation.However,the simulation results of 3 near-and 6 far-field stations at rupture distances of 12-36 km and 222-286 km,respectively,were obviously deviated from the observations.Simulation results of the near-field stations are larger than the observations at high frequencies(>6 Hz).The discrepancy likely comes from the nonlinear site effect of near-field stations,which reduced the site amplification at high frequencies.Simulation result of the far-field stations is smaller than the observation at frequencies above 1 Hz.As these stations are located close to the Longmenshan Fault Zone(LFZ),thus,we obtained a new quality factor(Q)from data of historical events and stations located around LFZ.Using the new Q value,the discrepancies of the high-frequency simulation results of the far-field stations were corrected.This result indicated that the laterally varying Q values can be used to address the impact of strong crustal lateral heterogeneity on simulation.
基金National Natural Science Foundation of China(91937000,91637208)。
文摘Considering the complex topographic forcing and large cryosphere concentration,the present study utilized the polar-optimized WRF model(Polar WRF)to conduct downscaling simulations over the Qinghai-Tibet Plateau(TP)and its surrounding regions.Multi-group experiments with the 10 km horizontal resolution are used to evaluate the modeling of precipitation.Firstly,on the basis of the model ground surface properties upgrade and the optimized Noah-MP,the“better-performing”configuration suite for modeling precipitation is comprehensively examined.Various model parameters such as nudging options,five cumulus parameterization schemes,two planetary boundary layer schemes,and six microphysics schemes are investigated to further refine the Polar WRF configuration.Moreover,the precipitation simulation for a full calendar year is compared with multiple reanalyses and observations.The simulations demonstrate that the Polar WRF model successfully captures the general features of precipitation over this region and is sensitive to model parameters.Based on the results,it is recommended to use grid nudging with q intensity coefficient of 0.0002,the multi-scale kain-fritsch cumulus parameterization,the Yonsei University boundary layer scheme,and the Morrison 2-mom microphysics with reduced default droplet concentration value of 100 cm-3.Overall,the model performance is better than the ERA-interim and TRMM 3b42.It is comparable to,and in some cases slightly better than,the CRA-Land,especially in the prediction for the western part of the plateau where in situ observations are limited,and the cryosphere-atmosphere interaction is more pronounced.
基金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.
基金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.
基金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 response of biological phantoms against high velocity impact is actively sought for applications in defense,space,soft robotics and sensing.Towards this end,we study the ballistic response of silicone based skin simulant against fragment impact.Using a pneumatic gas gun setup,six chisel-nosed and three regular shaped(sphere,cylinder,and cube)fragments were impacted on the skin simulant.The resulting skin simulant response was studied in terms of ballistic limit velocities,energy densities,failure pattern,and the mechanics of interaction.The results indicate that the shape of the fragment affects the ballistic limit velocities.The ballistic limit velocities,energy densities of the chisel-nosed fragment simulating projectiles were relatively insensitive to the size(mass),except for the smallest(0.16 g)and largest(2.79 g)chisel-nosed fragment.For the same size(1 g),ballistic limit velocities and failure are dependent on the shape of the fragment.The skin simulant failed by combined plugging and elastic hole enlargement.Failure in the spherical fragment was dominated by the elastic hole enlargement,whereas plugging failure was dominant in all other fragments.The spherical,cylindrical,and chisel-nosed fragments created circular cavities,and the cubical fragment created a square cavity.In the case of the spherical fragment,slipping of the fragment within the skin simulant was seen.Cubical fragments created lateral cracks emanating from the corners of the square cavity.Interestingly,for all the fragments,the maximum deformation corresponding to the perforation was lower than the non-perforation indicating rate dependent,stress driven failure.The maximum deformation was also dependent on the shape of the fragment.Overall,these results provide unique insights into the mechanical response of a soft simulant against ballistic impact.Results have utility in the calibration and validation of computational models,design of personal protective equipment,and antipersonnel systems.
基金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.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.62004223)the Science and Technology Innovation Program of Hunan Province,China(Grant No.2022RC1094)+1 种基金the Open Research Fund Program of the State Key Laboratory of Low-Dimensional Quantum Physics,China(Grant No.KF202012)the Hunan Provincial Science Innovation Project for Postgraduate,China(Grant No.CX20210086).
文摘The spin-transfer-torque(STT)magnetic tunneling junction(MTJ)device is one of the prominent candidates for spintronic logic circuit and neuromorphic computing.Therefore,building a simulation framework of hybrid STT-MTJ/CMOS(complementary metal-oxide-semiconductor)circuits is of great value for designing a new kind of computing paradigm based on the spintronic devices.In this work,we develop a simulation framework of hybrid STT-MTJ/CMOS circuits based on MATLAB/Simulink,which is mainly composed of a physics-based STT-MTJ model,a controlled resistor,and a current sensor.In the proposed framework,the STT-MTJ model,based on the Landau-Lifshitz-Gilbert-Slonczewsk(LLGS)equation,is implemented using the MATLAB script.The proposed simulation framework is modularized design,with the advantage of simple-to-use and easy-to-expand.To prove the effectiveness of the proposed framework,the STT-MTJ model is benchmarked with experimental results.Furthermore,the pre-charge sense amplifier(PCSA)circuit consisting of two STT-MTJ devices is validated and the electrical coupling of two spin-torque oscillators is simulated.The results demonstrate the effectiveness of our simulation framework.