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.展开更多
In this paper,a new approach called the Eulerian species method was proposed for simulating the convective and/or boiling heat transfer of nanofluids.The movement of nanoparticles in nanofluids is tracked by the speci...In this paper,a new approach called the Eulerian species method was proposed for simulating the convective and/or boiling heat transfer of nanofluids.The movement of nanoparticles in nanofluids is tracked by the species transport equation,and the boiling process of nanofluids is computed by the Eulerian multiphase method coupled with the RPI boiling model.The validity of the species transport equation for simulating nanoparticles movement was verified by conducting a simulation of nanofluids convective heat transfer.Simulation results of boiling heat transfer of nanofluids were obtained by using the commercial CFD software ANSYS Fluent and compared with experimental data and results from another numerical method(Eulerian three-phase model).Good agreement with experimental data was achieved,and it was proved the Eulerian species method is better than the Eulerian three-phase model since it can give better simulation results with higher accuracy but needs fewer computation resources.展开更多
The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifical...The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.展开更多
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.展开更多
The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and w...The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.展开更多
The Very Large Area gamma-ray Space Telescope(VLAST)is a mission concept proposed to detect gamma-ray photons through both Compton scattering and electron–positron pair production mechanisms,thus enabling the detecti...The Very Large Area gamma-ray Space Telescope(VLAST)is a mission concept proposed to detect gamma-ray photons through both Compton scattering and electron–positron pair production mechanisms,thus enabling the detection of photons with energies ranging from MeV to TeV.This project aims to conduct a comprehensive survey of the gamma-ray sky from a low-Earth orbit using an anti-coincidence detector,a tracker detector that also serves as a low-energy calorimeter,and a high-energy imaging calorimeter.We developed a Monte Carlo simulation application of the detector using the GEANT4 toolkit to evaluate the instrument performance,including the effective area,angular resolution,and energy resolution,and explored specific optimizations of the detector configuration.Our simulation-based analysis indicates that the current design of the VLAST is physically feasible,with an acceptance above 10 m^(2)sr which is four times larger than that of the Fermi-LAT,an energy resolution better than 2%at 10 GeV,and an angular resolution better than 0.2◦at 10 GeV.The VLAST project promises to make significant contributions to the field of gamma-ray astronomy and enhance our understanding of the cosmos.展开更多
基金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 theResearch Starting Programof Ludong University(Grant No.221/20220045).
文摘In this paper,a new approach called the Eulerian species method was proposed for simulating the convective and/or boiling heat transfer of nanofluids.The movement of nanoparticles in nanofluids is tracked by the species transport equation,and the boiling process of nanofluids is computed by the Eulerian multiphase method coupled with the RPI boiling model.The validity of the species transport equation for simulating nanoparticles movement was verified by conducting a simulation of nanofluids convective heat transfer.Simulation results of boiling heat transfer of nanofluids were obtained by using the commercial CFD software ANSYS Fluent and compared with experimental data and results from another numerical method(Eulerian three-phase model).Good agreement with experimental data was achieved,and it was proved the Eulerian species method is better than the Eulerian three-phase model since it can give better simulation results with higher accuracy but needs fewer computation resources.
文摘The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.
基金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.
基金supported by the National Natural Science Foundation of China (Grant Nos. 92158204, 41506001 and 42076019)a Project supported by the Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (Grant No. 311021005)。
文摘The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.
基金supported by the National Key Research and Development Program of China(No.2021YFA0718404)the National Natural Science Foundation of China(Nos.12220101003,12173098,U2031149)+2 种基金the Project for Young Scientists in Basic Research of Chinese Academy of Sciences(CAS)(No.YSBR-061)the Scientific Instrument Developing Project of CAS(No.GJJSTD20210009)the Youth Innovation Promotion Association of CAS,and the Young Elite Scientists Sponsorship Program by the China Association for Science and Technology(No.YESS20220197).
文摘The Very Large Area gamma-ray Space Telescope(VLAST)is a mission concept proposed to detect gamma-ray photons through both Compton scattering and electron–positron pair production mechanisms,thus enabling the detection of photons with energies ranging from MeV to TeV.This project aims to conduct a comprehensive survey of the gamma-ray sky from a low-Earth orbit using an anti-coincidence detector,a tracker detector that also serves as a low-energy calorimeter,and a high-energy imaging calorimeter.We developed a Monte Carlo simulation application of the detector using the GEANT4 toolkit to evaluate the instrument performance,including the effective area,angular resolution,and energy resolution,and explored specific optimizations of the detector configuration.Our simulation-based analysis indicates that the current design of the VLAST is physically feasible,with an acceptance above 10 m^(2)sr which is four times larger than that of the Fermi-LAT,an energy resolution better than 2%at 10 GeV,and an angular resolution better than 0.2◦at 10 GeV.The VLAST project promises to make significant contributions to the field of gamma-ray astronomy and enhance our understanding of the cosmos.