A class of hybrid algorithms of real-time simulation based on evaluation of non-integerstep right-hand side function are presented in this paper. And some results of the convergence and stability of the algorithms are...A class of hybrid algorithms of real-time simulation based on evaluation of non-integerstep right-hand side function are presented in this paper. And some results of the convergence and stability of the algorithms are given. Using the class of algorithms, evaluation for the right-hand side function is needed once in every integration-step. Moreover, comparing with the other methods with the same amount of work, their numerical stability regions are larger and the method errors are smaller, and the numerical experiments show that the algorithms are very effective.展开更多
In this paper, a mathematical model of real-time simulation is given, and the problem of convergence on real-time Runge-Kutta algorithms is analysed. At last a theorem on the relation between the order of compensation...In this paper, a mathematical model of real-time simulation is given, and the problem of convergence on real-time Runge-Kutta algorithms is analysed. At last a theorem on the relation between the order of compensation and the convergent order of real-time algorithm is proved.展开更多
With the development of accelerograph, strong ground motion data can be widely applied to many fields. Especially, it is an important milestone for strong motion observation to expand application fields into earthquak...With the development of accelerograph, strong ground motion data can be widely applied to many fields. Especially, it is an important milestone for strong motion observation to expand application fields into earthquake monitoring that real-time simulation of ground displacement can be obtained by strong motion records for determining three earthquake parameters. For the purpose of application, on the basis of principle of seismic response of single-degree-of-freedom (SDOF) system, this paper presents a suit of formula of simulating ground displacement records by using strong ground motion records with the help of simulator of SDOF system. The research results show that the technique is very efficient and can be widely applied to earthquake monitoring.展开更多
This paper extends the SPH method to gas simulation. The SPH (Smoothed Particles Hydrodynamics) method is the most popular method of flow simulation, which is widely used in large-scale liquid simulation. However, i...This paper extends the SPH method to gas simulation. The SPH (Smoothed Particles Hydrodynamics) method is the most popular method of flow simulation, which is widely used in large-scale liquid simulation. However, it is not found to apply to gas simulation, since those methods based on SPH can't be used in real-time simulation due to their enormous particles and huge computation. This paper proposes a method for gas simulation based on SPH with a small number of particles. Firstly, the method computes the position and density of each particle in each point-in-time, and outlines the shape of the simulated gas based on those particles. Secondly the method uses the grid technique to refine the shape with the diffusion of particle's density under the control of grid, and get more lifelike simulation result. Each grid will be assigned density according to the particles in it. The density determines the final appearance of the grid. For ensuring the natural transition of the color between adjacent grids, we give a diffuse process of density between these grids and assign appropriate values to vertexes of these grids. The experimental results show that the proposed method can give better gas simulation and meet the request of real-time.展开更多
In recent years,frequent fire disasters have led to enormous damage in China.Effective firefighting rescues can minimize the losses caused by fires.During the rescue processes,the travel time of fire trucks can be sev...In recent years,frequent fire disasters have led to enormous damage in China.Effective firefighting rescues can minimize the losses caused by fires.During the rescue processes,the travel time of fire trucks can be severely affected by traffic conditions,changing the effective coverage of fire stations.However,it is still challenging to determine the effective coverage of fire stations considering dynamic traffic conditions.This paper addresses this issue by combining the traveling time calculationmodelwith the effective coverage simulationmodel.In addition,it proposes a new index of total effective coverage area(TECA)based on the time-weighted average of the effective coverage area(ECA)to evaluate the urban fire services.It also selects China as the case study to validate the feasibility of the models,a fire station(FS-JX)in Changsha.FS-JX station and its surrounding 9,117 fire risk points are selected as the fire service supply and demand points,respectively.A total of 196 simulation scenarios throughout a consecutiveweek are analyzed.Eventually,1,933,815 sets of valid sample data are obtained.The results showed that the TECA of FS-JX is 3.27 km^(2),which is far below the standard requirement of 7.00 km^(2) due to the traffic conditions.The visualization results showed that three rivers around FS-JX interrupt the continuity of its effective coverage.The proposed method can provide data support to optimize the locations of fire stations by accurately and dynamically determining the effective coverage of fire stations.展开更多
Background Petrochemical products possess a high risk of flammability,explosivity,and toxicity,making petrochemical accidents exceedingly destructive.Therefore,disaster analysis,prediction,and real-time simulations ha...Background Petrochemical products possess a high risk of flammability,explosivity,and toxicity,making petrochemical accidents exceedingly destructive.Therefore,disaster analysis,prediction,and real-time simulations have become important means of controlling and reducing accident hazards.Methods This study proposes a complete real-time simulation solution of gas diffusion with coordinate and concentration data,which was mainly aimed at simulating the types of harmful gas leakage and diffusion accidents in the petrochemical industry.The rendering effect was more continuous and accurate through grid homogenization and trilinear interpolation.This study presents a data processing and rendering parallelization process to enhance simulation efficiency.Gas concentration and fragment transparency were combined to synthesize transparent pixels in a scene.To ensure the approximate accuracy of the rendering effect,improve the efficiency of real-time rendering,and meet the requirement of intuitive perception using concentration data,a weighted blended order-independent transparency(OIT)with enhanced alpha weight is presented,which can provide a more intuitive perception of the hierarchical information of concentration data while preserving depth information.This study compares and analyzes three OIT algorithms-depth peeling,weighted blended OIT,and weighted blended OIT with enhanced alpha weight-in terms of rendering image quality,rendering time,required memory,and hierarchical information.Results Using weighted blended OIT with an enhanced alpha weight technique,the rendering time was shortened by 53.2%compared with that of the depth peeling algorithm,and the texture memory required was significantly smaller than that of the depth peeling algorithm.The rendering results of weighted blended OIT with an enhanced alpha weight were approximately accurate compared with those of the depth peeling algorithm as the ground truth,and there was no popping when surfaces passed through one another.Simultaneously,compared with weighted blended OIT,weighted blended OIT with an enhanced alpha weight achieved an intuitive perception of the hierarchical information of concentration data.展开更多
Computational models that ensure accurate and fast responses to the variations in operating conditions,such as the cell tem-perature and relative humidity(RH),are essential monitoring tools for the real-time control o...Computational models that ensure accurate and fast responses to the variations in operating conditions,such as the cell tem-perature and relative humidity(RH),are essential monitoring tools for the real-time control of proton exchange membrane(PEM)fuel cells.To this end,fast cell-area-averaged numerical simulations are developed and verifi ed against the present experiments under various RH levels.The present simulations and measurements are found to agree well based on the cell voltage(polarization curve)and power density under variable RH conditions(RH=40%,RH=70%,and RH=100%),which verifi es the model accuracy in predicting PEM fuel cell performance.In addition,computationally feasible reduced-order models are found to deliver a fast output dataset to evaluate the charge/heat/mass transfer phenomena as well as water production and two-phase fl ow transport.Such fast and accurate evaluations of the overall fuel cell operation can be used to inform the real-time control systems that allow for the improved optimization of PEM fuel cell performance.展开更多
In this paper a class of real-time parallel modified Rosenbrock methods of numerical simulation is constructed for stiff dynamic systems on a multiprocessor system, and convergence and numerical stability of these met...In this paper a class of real-time parallel modified Rosenbrock methods of numerical simulation is constructed for stiff dynamic systems on a multiprocessor system, and convergence and numerical stability of these methods are discussed. A-stable real-time parallel formula of two-stage third-order and A(α)-stable real-time parallel formula with o ≈ 89.96° of three-stage fourth-order are particularly given. The numerical simulation experiments in parallel environment show that the class of algorithms is efficient and applicable, with greater speedup.展开更多
The electric vehicle(EV)charging station is a critical part of the infrastructure for the wide adoption of EVs.Realtime simulation of an EV station plays an essential role in testing its operation under different oper...The electric vehicle(EV)charging station is a critical part of the infrastructure for the wide adoption of EVs.Realtime simulation of an EV station plays an essential role in testing its operation under different operating modes.However,the large numbers of high-frequency power electronic switches contained in EV chargers pose great challenges for real-time simulation.This paper proposes a compact electromagnetic transient program(C-EMTP)algorithm for FPGA-based real-time simulation of an EV station with multiple high-frequency chargers.The C-EMTP algorithm transforms the traditional EMTP algorithm into two parallel sub-tasks only consisting of simple matrix operations,to fully utilize the high parallelism of FPGA.The simulation time step can be greatly reduced compared with that of the traditional EMTP algorithm,and so the simulation accuracy for high-frequency power electronics is improved.The EV chargers can be decoupled with each other and simulated in parallel.A CPU-FPGA-based realtime simulation platform is developed and the proposed simulation of the EV station is implemented.The control strategy is simulated in a CPU with 100μs time-step,while the EV station circuit topology is simulated in a single FPGA with a 250 ns time-step.In the case studies,the EV station consists of a two-level rectifier and five dual-active bridge(DAB)EV chargers.It is tested under different scenarios,and the real-time simulation results are validated using PSCAD/EMTDC.展开更多
The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameteri...The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.展开更多
Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging ...Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.展开更多
We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscal...We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.展开更多
文摘A class of hybrid algorithms of real-time simulation based on evaluation of non-integerstep right-hand side function are presented in this paper. And some results of the convergence and stability of the algorithms are given. Using the class of algorithms, evaluation for the right-hand side function is needed once in every integration-step. Moreover, comparing with the other methods with the same amount of work, their numerical stability regions are larger and the method errors are smaller, and the numerical experiments show that the algorithms are very effective.
文摘In this paper, a mathematical model of real-time simulation is given, and the problem of convergence on real-time Runge-Kutta algorithms is analysed. At last a theorem on the relation between the order of compensation and the convergent order of real-time algorithm is proved.
基金National Nature Science Foundation of China (50378086) Commonweal Foundation of the Ministry of Science and+1 种基金Technology (2003DIB2J099) key project of China EarthquakeAdministrationduring the tenth Five-year Plan.
文摘With the development of accelerograph, strong ground motion data can be widely applied to many fields. Especially, it is an important milestone for strong motion observation to expand application fields into earthquake monitoring that real-time simulation of ground displacement can be obtained by strong motion records for determining three earthquake parameters. For the purpose of application, on the basis of principle of seismic response of single-degree-of-freedom (SDOF) system, this paper presents a suit of formula of simulating ground displacement records by using strong ground motion records with the help of simulator of SDOF system. The research results show that the technique is very efficient and can be widely applied to earthquake monitoring.
基金Supported by National Natural Science Foundation of China(No.61272024)
文摘This paper extends the SPH method to gas simulation. The SPH (Smoothed Particles Hydrodynamics) method is the most popular method of flow simulation, which is widely used in large-scale liquid simulation. However, it is not found to apply to gas simulation, since those methods based on SPH can't be used in real-time simulation due to their enormous particles and huge computation. This paper proposes a method for gas simulation based on SPH with a small number of particles. Firstly, the method computes the position and density of each particle in each point-in-time, and outlines the shape of the simulated gas based on those particles. Secondly the method uses the grid technique to refine the shape with the diffusion of particle's density under the control of grid, and get more lifelike simulation result. Each grid will be assigned density according to the particles in it. The density determines the final appearance of the grid. For ensuring the natural transition of the color between adjacent grids, we give a diffuse process of density between these grids and assign appropriate values to vertexes of these grids. The experimental results show that the proposed method can give better gas simulation and meet the request of real-time.
基金support from the National Natural Science Foundation of China (No.52204202)the Hunan Provincial Natural Science Foundation of China (No.2023JJ40058)the Science and Technology Program of Hunan Provincial Departent of Transportation (No.202122).
文摘In recent years,frequent fire disasters have led to enormous damage in China.Effective firefighting rescues can minimize the losses caused by fires.During the rescue processes,the travel time of fire trucks can be severely affected by traffic conditions,changing the effective coverage of fire stations.However,it is still challenging to determine the effective coverage of fire stations considering dynamic traffic conditions.This paper addresses this issue by combining the traveling time calculationmodelwith the effective coverage simulationmodel.In addition,it proposes a new index of total effective coverage area(TECA)based on the time-weighted average of the effective coverage area(ECA)to evaluate the urban fire services.It also selects China as the case study to validate the feasibility of the models,a fire station(FS-JX)in Changsha.FS-JX station and its surrounding 9,117 fire risk points are selected as the fire service supply and demand points,respectively.A total of 196 simulation scenarios throughout a consecutiveweek are analyzed.Eventually,1,933,815 sets of valid sample data are obtained.The results showed that the TECA of FS-JX is 3.27 km^(2),which is far below the standard requirement of 7.00 km^(2) due to the traffic conditions.The visualization results showed that three rivers around FS-JX interrupt the continuity of its effective coverage.The proposed method can provide data support to optimize the locations of fire stations by accurately and dynamically determining the effective coverage of fire stations.
基金Supported by National Key R&D Program of China (2020YFB1710400)。
文摘Background Petrochemical products possess a high risk of flammability,explosivity,and toxicity,making petrochemical accidents exceedingly destructive.Therefore,disaster analysis,prediction,and real-time simulations have become important means of controlling and reducing accident hazards.Methods This study proposes a complete real-time simulation solution of gas diffusion with coordinate and concentration data,which was mainly aimed at simulating the types of harmful gas leakage and diffusion accidents in the petrochemical industry.The rendering effect was more continuous and accurate through grid homogenization and trilinear interpolation.This study presents a data processing and rendering parallelization process to enhance simulation efficiency.Gas concentration and fragment transparency were combined to synthesize transparent pixels in a scene.To ensure the approximate accuracy of the rendering effect,improve the efficiency of real-time rendering,and meet the requirement of intuitive perception using concentration data,a weighted blended order-independent transparency(OIT)with enhanced alpha weight is presented,which can provide a more intuitive perception of the hierarchical information of concentration data while preserving depth information.This study compares and analyzes three OIT algorithms-depth peeling,weighted blended OIT,and weighted blended OIT with enhanced alpha weight-in terms of rendering image quality,rendering time,required memory,and hierarchical information.Results Using weighted blended OIT with an enhanced alpha weight technique,the rendering time was shortened by 53.2%compared with that of the depth peeling algorithm,and the texture memory required was significantly smaller than that of the depth peeling algorithm.The rendering results of weighted blended OIT with an enhanced alpha weight were approximately accurate compared with those of the depth peeling algorithm as the ground truth,and there was no popping when surfaces passed through one another.Simultaneously,compared with weighted blended OIT,weighted blended OIT with an enhanced alpha weight achieved an intuitive perception of the hierarchical information of concentration data.
基金by the Natural Sciences and Engineering Research Council of Canada(NSERC)via a Discovery Grant,Canadian Urban Transit Research and Innovation Consortium(CUTRIC)(No.160028).
文摘Computational models that ensure accurate and fast responses to the variations in operating conditions,such as the cell tem-perature and relative humidity(RH),are essential monitoring tools for the real-time control of proton exchange membrane(PEM)fuel cells.To this end,fast cell-area-averaged numerical simulations are developed and verifi ed against the present experiments under various RH levels.The present simulations and measurements are found to agree well based on the cell voltage(polarization curve)and power density under variable RH conditions(RH=40%,RH=70%,and RH=100%),which verifi es the model accuracy in predicting PEM fuel cell performance.In addition,computationally feasible reduced-order models are found to deliver a fast output dataset to evaluate the charge/heat/mass transfer phenomena as well as water production and two-phase fl ow transport.Such fast and accurate evaluations of the overall fuel cell operation can be used to inform the real-time control systems that allow for the improved optimization of PEM fuel cell performance.
基金This project was supported by the National Natural Science Foundation of China (No. 19871080).
文摘In this paper a class of real-time parallel modified Rosenbrock methods of numerical simulation is constructed for stiff dynamic systems on a multiprocessor system, and convergence and numerical stability of these methods are discussed. A-stable real-time parallel formula of two-stage third-order and A(α)-stable real-time parallel formula with o ≈ 89.96° of three-stage fourth-order are particularly given. The numerical simulation experiments in parallel environment show that the class of algorithms is efficient and applicable, with greater speedup.
基金supported by China Postdoctoral Science Foundation(BX20200221,2020 M671122)National Key Research and Development Program of China(2019YFE012784)National Natural Science Foundation of China(51877133).
文摘The electric vehicle(EV)charging station is a critical part of the infrastructure for the wide adoption of EVs.Realtime simulation of an EV station plays an essential role in testing its operation under different operating modes.However,the large numbers of high-frequency power electronic switches contained in EV chargers pose great challenges for real-time simulation.This paper proposes a compact electromagnetic transient program(C-EMTP)algorithm for FPGA-based real-time simulation of an EV station with multiple high-frequency chargers.The C-EMTP algorithm transforms the traditional EMTP algorithm into two parallel sub-tasks only consisting of simple matrix operations,to fully utilize the high parallelism of FPGA.The simulation time step can be greatly reduced compared with that of the traditional EMTP algorithm,and so the simulation accuracy for high-frequency power electronics is improved.The EV chargers can be decoupled with each other and simulated in parallel.A CPU-FPGA-based realtime simulation platform is developed and the proposed simulation of the EV station is implemented.The control strategy is simulated in a CPU with 100μs time-step,while the EV station circuit topology is simulated in a single FPGA with a 250 ns time-step.In the case studies,the EV station consists of a two-level rectifier and five dual-active bridge(DAB)EV chargers.It is tested under different scenarios,and the real-time simulation results are validated using PSCAD/EMTDC.
基金supported by the National Natural Science Foundation of China(Grant Nos.42175099,42027804,42075073)the Innovative Project of Postgraduates in Jiangsu Province in 2023(Grant No.KYCX23_1319)+3 种基金supported by the National Natural Science Foundation of China(Grant No.42205080)the Natural Science Foundation of Sichuan(Grant No.2023YFS0442)the Research Fund of Civil Aviation Flight University of China(Grant No.J2022-037)supported by the National Key Scientific and Technological Infrastructure project“Earth System Science Numerical Simulator Facility”(Earth Lab)。
文摘The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.
基金supported by the National Natural Science Foundation of China(NNSFC)grants 42074202,42274196Strategic Priority Research Program of Chinese Academy of Sciences grant XDB41000000ISSI-BJ International Team Interaction between magnetic reconnection and turbulence:From the Sun to the Earth。
文摘Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.
基金support from the European Space Agency(ESA)PRODEX(PROgramme de Développement d’Expériences scientifiques)Project mission(No.PEA4000134960)Partial funding was provided by the Romanian Ministry of Research,Innovation and Digitalization under Romanian National Core Program LAPLAS VII(Contract No.30N/2023)+2 种基金the Belgian Solar-Terrestrial Centre of Excellencesupported by the project Belgian Research Action through Interdisciplinary Networks(BRAIN-BE)2.0(Grant No.B2/223/P1/PLATINUM)funded by the Belgian Office for Research(BELSPO)partially supported by a grant from the Romanian Ministry of Education and Research(CNCS-UEFISCDI,Project No.PN-III-P1-1.1TE-2021-0102)。
文摘We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.