A 3D dynamic core of the non-hydrostatic model GRAPES（Global/Regional Assimilation and Prediction System） is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. T...A 3D dynamic core of the non-hydrostatic model GRAPES（Global/Regional Assimilation and Prediction System） is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. Three-dimensional Coriolis forcing is introduced to the new core, and full representation of the Coriolis forcing makes it straightforward to share code between the Yin and Yang subdomains. Similar to that in the original GRAPES model, a semi-implicit semi-Lagrangian scheme is adopted for temporal integration and advection with additional arrangement for cross-boundary transport. Under a non-centered second-order temporal and spatial discretization, the dry nonhydrostatic frame is summarized as the solution of an elliptical problem. The resulting Helmholtz equation is solved with the Generalized Conjugate Residual solver in cooperation with the classic Schwarz method. Even though the coefficients of the equation are quite different from those in the original model, the computational procedure of the new core is just the same. The bi-cubic Lagrangian interpolation serves to provide Dirichlet-type boundary conditions with data transfer between the subdomains. The dry core is evaluated with several benchmark test cases, and all the tests display reasonable numerical stability and computing performance. Persistency of the balanced flow and development of both the mountain-induced Rossby wave and Rossby–Haurwitz wave confirms the appropriate installation of the 3D Coriolis terms in the semi-implicit semi-Lagrangian dynamic core on the Yin-Yang grid.展开更多
A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrain...A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrainfollowing grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future.展开更多
A reactor core in a thorium molten salt reactor uses graphite as a moderator and reflector. The graphite core is a multi-layered arrangement of graphite bricks that are loosely connected to each other using a system o...A reactor core in a thorium molten salt reactor uses graphite as a moderator and reflector. The graphite core is a multi-layered arrangement of graphite bricks that are loosely connected to each other using a system of keys and dowels. Consequently, the graphite core is a type of discrete stack structure with highly nonlinear dynamic behavior. Hence, it is important to investigate the dynamic characteristics of the graphite core. In this study, a threedimensional single-layer graphite core model, which is a part of the thorium molten salt reactor side reflector structure, was analyzed using the explicit method in ABAQUS 2016 to study the core dynamic behavior when subjected to different excitations. The design parameters,such as the diameter of the dowel, the gap between key and keyway and the bypass flow gap between two adjacent bricks, were also considered in this model. To reduce excessive demands on available computational resources considering the effect of molten salt, the spring–dashpot model was applied to model the interaction forces between the molten salt and graphite bricks. Numerical simulation results show that the effect of molten salt is a reduction inthe peak maximal principal stress, and a larger gap between two bricks is beneficial to maintain the integrity of the graphite core under earthquake loading. The results obtained by the simulation can be used as a reference for future designs of a molten salt graphite core.展开更多
The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The ...The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The field of displacements is for- mulated using the classical broken line hypothesis and the proposed nonlinear hypothesis that generalizes the classical one. Using both hypotheses, the strains are determined as well as the stresses of each layer. The kinetic energy, the elastic strain energy, and the work of load are also determined. The system of equations of motion is derived using Hamilton＇s principle. Finally, the system of three equations is reduced to one equation of motion, in particular, the Mathieu equation. The Bubnov-Galerkin method is used to solve the system of equations of motion, and the Runge-Kutta method is used to solve the second-order differential equation. Numerical calculations are done for the chosen family of beams. The critical loads, unstable regions, angular frequencies of the beam, and the static and dynamic equilibrium paths are calculated analytically and verified numerically. The results of this study are presented in the forms of figures and tables.展开更多
Core shooting process plays a decisive role in the quality of sand cores, and core box vents distribution is one of the most important factor determining the effectiveness of core shooting process. In this paper, the ...Core shooting process plays a decisive role in the quality of sand cores, and core box vents distribution is one of the most important factor determining the effectiveness of core shooting process. In this paper, the influence of core box vents distribution on the flow dynamics of core shooting process was investigated based on in situ experimental observations with transparent core box, high-speed camera and pressure measuring system. Attention was focused on the variation of both the flow behavior of sand and pressure curves due to different vents distribution. Taking both kinetic and frictional stress into account, a kinetic-frictional constitutive model was established to describe the internal momentum transfer in the solid phase. Two-fluid model(TFM) simulation was then performed and good agreement was achieved between the experimental and simulated results on both the flow behavior of sand and the pressure curves. It was found that vents distribution has direct effect on the pressure difference of different locations in the core box, which determines the buoyancy force exerting on the sand particles and significantly influences the filling process of core sand.展开更多
The elasto-gravitational deformation response of the Earths solid parts to the perturbations of the pressure and gravity on the core-mantle boundary (CMB) and the solid inner core boundary (ICB), due to the dynamical ...The elasto-gravitational deformation response of the Earths solid parts to the perturbations of the pressure and gravity on the core-mantle boundary (CMB) and the solid inner core boundary (ICB), due to the dynamical behaviors of the fluid outer core (FOC), is discussed. The internal load Love numbers, which are formulized in a general form in this study, are employed to describe the Earths deformation. The preliminary reference Earth model (PREM) is used as an example to calculate the internal load Love numbers on the Earths surface, CMB and ICB, respectively. The characteristics of the Earths deformation variation with the depth and the perturbation periods on the boundaries of the FOC are also investigated. The numerical results indicate that the internal load Love numbers decrease quickly with the increasing degree of the spherical harmonics of the displacement and depend strongly on the perturbation frequencies, especially on the high frequencies. The results, obtained in this work, can be used to construct the boundary conditions for the core dynamics of the long-period oscillations of the Earths fluid outer core.展开更多
The influence of core sand properties on flow dynamics was investigated synchronously with various core sands, transparent core-box and high-speed camera. To confirm whether the core shooting process has significant t...The influence of core sand properties on flow dynamics was investigated synchronously with various core sands, transparent core-box and high-speed camera. To confirm whether the core shooting process has significant turbulence, the flow pattern of sand particles in the shooting head and core box was reproduced with colored core sands. By incorporating the kinetic theory of granular flow(KTGF), kinetic-frictional constitutive correlation and turbulence model, a two-fluid model(TFM) was established to study the flow dynamics of the core shooting process. Two-fluid model(TFM) simulations were then performed and a areasonable agreement was achieved between the simulation and experimental results. Based on the experimental and simulation results, the effects of turbulence, sand density, sand diameter and binder ratio were analyzed in terms of filling process, sand volume fraction(αs) and sand velocity(Vs).展开更多
To the reduct problems of decision system, the paper proposes the notion of dynamic core according to the dynamic reduct model. It describes various formal definitions of dynamic core, and discusses some properties ab...To the reduct problems of decision system, the paper proposes the notion of dynamic core according to the dynamic reduct model. It describes various formal definitions of dynamic core, and discusses some properties about dynamic core. All of these show that dynamic core possesses the essential characters of the feature core.展开更多
We give a brief review of the asymptotic theory of slender vortex filaments with emphases on (i) the choices of scalings and small parameters characterizing the physical problem,(ii) the key steps in the formulation o...We give a brief review of the asymptotic theory of slender vortex filaments with emphases on (i) the choices of scalings and small parameters characterizing the physical problem,(ii) the key steps in the formulation of the theory and (iii) the assumptions and/or restrictions on the theory of Callegari and Ting (1978).We present highlights of an extension of the 1978 asymptotic theory:the analyses for core structures with axial variation.Making use of the physical insights gained from the analyses,we present a new derivation of the evolution equations for the core structure.The new one is simpler and straightforward and shows the physics clearly.展开更多
Here,a nonhydrostatic alternative scheme(NAS)is proposed for the grey zone where the nonhydrostatic impact on the atmosphere is evident but not large enough to justify the necessity to include an implicit nonhydrostat...Here,a nonhydrostatic alternative scheme(NAS)is proposed for the grey zone where the nonhydrostatic impact on the atmosphere is evident but not large enough to justify the necessity to include an implicit nonhydrostatic solver in an atmospheric dynamical core.The NAS is designed to replace this solver,which can be incorporated into any hydrostatic models so that existing well-developed hydrostatic models can effectively serve for a longer time.Recent advances in machine learning(ML)provide a potential tool for capturing the main complicated nonlinear-nonhydrostatic relationship.In this study,an ML approach called a neural network(NN)was adopted to select leading input features and develop the NAS.The NNs were trained and evaluated with 12-day simulation results of dry baroclinic-wave tests by the Weather Research and Forecasting(WRF)model.The forward time difference of the nonhydrostatic tendency was used as the target variable,and the five selected features were the nonhydrostatic tendency at the last time step,and four hydrostatic variables at the current step including geopotential height,pressure in two different forms,and potential temperature,respectively.Finally,a practical NAS was developed with these features and trained layer by layer at a 20-km horizontal resolution,which can accurately reproduce the temporal variation and vertical distribution of the nonhydrostatic tendency.Corrected by the NN-based NAS,the improved hydrostatic solver at different horizontal resolutions can run stably for at least one month and effectively reduce most of the nonhydrostatic errors in terms of system bias,anomaly root-mean-square error,and the error of the wave spatial pattern,which proves the feasibility and superiority of this scheme.展开更多
The Global-Regional Integrated forecast System(GRIST)is the next-generation weather and climate integrated model dynamic framework developed by Chinese Academy of Meteorological Sciences.In this paper,we present sever...The Global-Regional Integrated forecast System(GRIST)is the next-generation weather and climate integrated model dynamic framework developed by Chinese Academy of Meteorological Sciences.In this paper,we present several changes made to the global nonhydrostatic dynamical(GND)core,which is part of the ongoing prototype of GRIST.The changes leveraging MPI and PnetCDF techniques were targeted at the parallelization and performance optimization to the original serial GND core.Meanwhile,some sophisticated data structures and interfaces were designed to adjust flexibly the size of boundary and halo domains according to the variable accuracy in parallel context.In addition,the I/O performance of PnetCDF decreases as the number of MPI processes increases in our experimental environment.Especially when the number exceeds 6000,it caused system-wide outages(SWO).Thus,a grouping solution was proposed to overcome that issue.Several experiments were carried out on the supercomputing platform based on Intel x86 CPUs in the National Supercomputing Center in Wuxi.The results demonstrated that the parallel GND core based on grouping solution achieves good strong scalability and improves the performance significantly,as well as avoiding the SWOs.展开更多
An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several poi...An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.展开更多
In the last several years since 2004 an important new finding has been unveiled by combined efforts due to Japanese （Satoru Tanaka）, French （Renaud Deguen, Y Albousierre and Marc Monnereau）, American and Chinese g...In the last several years since 2004 an important new finding has been unveiled by combined efforts due to Japanese （Satoru Tanaka）, French （Renaud Deguen, Y Albousierre and Marc Monnereau）, American and Chinese geophysicists （Xiaodong Song and Vernon F Cormier） who employed from unambi- guous detailed seismological evidence and explained by clear theoretical and sound laboratory arguments drawn from fluid dynamics that there exists a strong East-West hemi-spherical asymmetry on the inner- outer core boundary,展开更多
We numerically study the propagation dynamics of intense optical pulses in gas-filled hollow-core fibers（HCFs）. The spatiotemporal dynamics of the pulses show a transition from tightly confined to loosely confined c...We numerically study the propagation dynamics of intense optical pulses in gas-filled hollow-core fibers（HCFs）. The spatiotemporal dynamics of the pulses show a transition from tightly confined to loosely confined characteristics as the fiber core is increased, which manifests as a deterioration in the spatiotemporal uniformity of the beam. It is found that using the gas pressure gradient does not enhance the beam quality in large-core HCFs, while inducing a positive chirp in the pulse to lower the peak power can improve the beam quality. This indicates that the self-focusing effect in the HCFs is the main driving force for the propagation dynamics. It also suggests that pulses at longer wavelengths are more suitable for HCFs with large cores because of the lower critical power of self-focusing, which is justified by the numerical simulations. These results will benefit the generation of energetic few-cycle pulses in large-core HCFs.展开更多
Most current lattice Boltzmann （LBM） models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular ...Most current lattice Boltzmann （LBM） models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular dynamics （MD） simulation and LBM to avoid such defect. The basic idea is to first construct a molecular model based on the actual components of the rock-fluid system, then to compute the interaction force between the rock and the fluid of different densities through the MD simulation. This calculated rock-fluid interaction force, combined with the fluid-fluid force determined from the equation of state, is then used in LBM modeling. Without parameter fitting, this study presents a new systematic approach for pore-scale modeling of multi-phase flow. We have validated this ap- proach by simulating a two-phase separation process and gas-liquid-solid three-phase contact angle. Based on an actual X-ray CT image of a reservoir core, we applied our workflow to calculate the absolute permeability of the core, vapor-liquid H20 relative permeability, and capillary pressure curves.展开更多
Dynamic strain aging (DSA) is an important phenomenon in solutehardened metals and seri- ously affects the mechanical properties ofmetals. DSA is generally induced by the interaction between themoving dislocations and...Dynamic strain aging (DSA) is an important phenomenon in solutehardened metals and seri- ously affects the mechanical properties ofmetals. DSA is generally induced by the interaction between themoving dislocations and the mobile solute atoms. In this paper, onlythe interaction between moving disloca- tions and mobile solute atomsin a dislocation core area (core atmosphere) will be taken intoaccount. To es- tablish the constitutive model which can describe theDSA phenomenon, we improved the Zerilli-Armstrongdislocation-mechanics-based thermal viscoplastic constitutiverelation, and added the effect of the interaction between the movingdislocations and core atmosphere.展开更多
How to properly consider the impacts of non-hydrostatic perturbations is one of the challenging issues in developing non-hydrostatic dynamics solvers（NHDSs） for high-resolution atmospheric models. To overcome the dr...How to properly consider the impacts of non-hydrostatic perturbations is one of the challenging issues in developing non-hydrostatic dynamics solvers（NHDSs） for high-resolution atmospheric models. To overcome the drawbacks of current approaches to tackling this issue, this study analyzed the differences between hydrostatic dynamics solvers（HDSs） and their non-hydrostatic counterparts.The analysis then motivated a flexible approach to adjusting existing hydrostatic atmospheric models,especially those adopted in climate simulations for the impacts of non-hydrostatic perturbations.In this approach, the impacts of non-hydrostatic perturbations, reflecting the differences between HDSs and NHDSs, were encapsulated into a single term in the vertical momentum equation for the atmosphere. At each time step, this term was estimated by a separate sub-model, and then it was used to adjust the dynamics of the atmosphere. The adjustment was optional, and could be turned on and off flexibly by utilizing different initial conditions. The approach was illustrated using the Weather Research and Forecasting（WRF） model as an example, and was preliminarily validated by running the 3D baroclinic-wave test case in the model. Results showed that the modified dynamics solver produced simulation results that were very close to those given by the standard NHDS in the WRF model, implying that the approach was basically effective in capturing the non-hydrostatic features of the atmosphere.展开更多
In relatively coarse-resolution atmospheric models,cumulus parameterization helps account for the effect of subgridscale convection,which produces supplemental rainfall to the grid-scale precipitation and impacts the ...In relatively coarse-resolution atmospheric models,cumulus parameterization helps account for the effect of subgridscale convection,which produces supplemental rainfall to the grid-scale precipitation and impacts the diurnal cycle of precipitation.In this study,the diurnal cycle of precipitation was studied using the new simplified Arakawa-Schubert scheme in a global non-hydrostatic atmospheric model,i.e.,the Yin-Yang-grid Unified Model for the Atmosphere.Two new diagnostic closures and a convective trigger function were suggested to emphasize the job of the cloud work function corresponding to the free tropospheric large-scale forcing.Numerical results of the 0.25-degree model in 3-month batched real-case simulations revealed an improvement in the diurnal precipitation variation by using a revised trigger function with an enhanced dynamical constraint on the convective initiation and a suitable threshold of the trigger.By reducing the occurrence of convection during peak solar radiation hours,the revised scheme was shown to be effective in delaying the appearance of early-afternoon rainfall peaks over most land areas and accentuating the nocturnal peaks that were wrongly concealed by the more substantial afternoon peak.In addition,the revised scheme enhanced the simulation capability of the precipitation probability density function,such as increasing the extremely low-and high-intensity precipitation events and decreasing small and moderate rainfall events,which contributed to the reduction of precipitation bias over mid-latitude and tropical land areas.展开更多
A mathematical model based on the computational fluid dynamics method,heat and mass transfer in porous media and the unreacted shrinking core model for the oxidation reaction of an individual magnetite pellet during p...A mathematical model based on the computational fluid dynamics method,heat and mass transfer in porous media and the unreacted shrinking core model for the oxidation reaction of an individual magnetite pellet during preheating was established.The commercial software COMSOL Multiphysics was used to simulate the change in the oxidation degree of the pellet at different temperatures and oxygen concentrations,and the simulated results were compared with the exper-imental results.The model considered the influence of the exothermic heat of the reaction,and the enthalpy change was added to calculate the heat released by the oxidation.The results show that the oxidation rate on the surface of the pellet is much faster than that of the inside of the pellet.Temperature and oxygen concentration have great influence on the pellet oxidation model.Meanwhile,the exothermic calculation results show that there is a non-isothermal phenomenon inside the pellet,which leads to an increase in temperature inside the single pellet.Under the preheating condition of 873-1273 K(20%oxygen content),the heat released by the pellet oxidation reaction in a chain grate is 7.8×10^(6)-10.8×10^(6) kJ/h,which is very large and needs to be considered in the magnetite pellet oxidation modelling.展开更多
Two new two-photon absorption （TPA） molecules, named SK-G1 and NT-G1, are synthesized and the photo- physical characteristics are investigated by using linear absorption spectra, one-photon fluorescence spectra and ...Two new two-photon absorption （TPA） molecules, named SK-G1 and NT-G1, are synthesized and the photo- physical characteristics are investigated by using linear absorption spectra, one-photon fluorescence spectra and two-photon excited fluorescence spectra. Both the compounds exhibit TPA properties, and the TPA values determined by z-scan measurement are 10 GM and 39 GM for SK-G1 and NT-G1, respectively, at wavelength 80Ohm. Time-resolved spectroscopic techniques are employed to further explore the excited state dynamics of NT-G l with larger TPA cross section. The research results show that there is an ultrafast intraband energy transfer process （about 3ps） before the formation of charge transfer state with a relatively long lifetime.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 41175095)the National Key Technology R&D Program(Grant No. 2012BAC22B01)a research project of the Chinese Academy of Meteorological Sciences (Grant No. 2014Z001)
文摘A 3D dynamic core of the non-hydrostatic model GRAPES（Global/Regional Assimilation and Prediction System） is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. Three-dimensional Coriolis forcing is introduced to the new core, and full representation of the Coriolis forcing makes it straightforward to share code between the Yin and Yang subdomains. Similar to that in the original GRAPES model, a semi-implicit semi-Lagrangian scheme is adopted for temporal integration and advection with additional arrangement for cross-boundary transport. Under a non-centered second-order temporal and spatial discretization, the dry nonhydrostatic frame is summarized as the solution of an elliptical problem. The resulting Helmholtz equation is solved with the Generalized Conjugate Residual solver in cooperation with the classic Schwarz method. Even though the coefficients of the equation are quite different from those in the original model, the computational procedure of the new core is just the same. The bi-cubic Lagrangian interpolation serves to provide Dirichlet-type boundary conditions with data transfer between the subdomains. The dry core is evaluated with several benchmark test cases, and all the tests display reasonable numerical stability and computing performance. Persistency of the balanced flow and development of both the mountain-induced Rossby wave and Rossby–Haurwitz wave confirms the appropriate installation of the 3D Coriolis terms in the semi-implicit semi-Lagrangian dynamic core on the Yin-Yang grid.
基金supported by the National Key Research and Development Program of China (Grant Nos. 2017YFC1501901 and 2017YFA0603901)the Beijing Natural Science Foundation (Grant No. JQ18001)
文摘A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrainfollowing grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future.
基金supported by the“Hundred Talent Program”of the Chinese Academy of Sciences,the Ministry of Human Resources and Social Security(No.Y419016031)the Strategic Priority Research Program of Chinese Academy of Science(No.XDA02040100)
文摘A reactor core in a thorium molten salt reactor uses graphite as a moderator and reflector. The graphite core is a multi-layered arrangement of graphite bricks that are loosely connected to each other using a system of keys and dowels. Consequently, the graphite core is a type of discrete stack structure with highly nonlinear dynamic behavior. Hence, it is important to investigate the dynamic characteristics of the graphite core. In this study, a threedimensional single-layer graphite core model, which is a part of the thorium molten salt reactor side reflector structure, was analyzed using the explicit method in ABAQUS 2016 to study the core dynamic behavior when subjected to different excitations. The design parameters,such as the diameter of the dowel, the gap between key and keyway and the bypass flow gap between two adjacent bricks, were also considered in this model. To reduce excessive demands on available computational resources considering the effect of molten salt, the spring–dashpot model was applied to model the interaction forces between the molten salt and graphite bricks. Numerical simulation results show that the effect of molten salt is a reduction inthe peak maximal principal stress, and a larger gap between two bricks is beneficial to maintain the integrity of the graphite core under earthquake loading. The results obtained by the simulation can be used as a reference for future designs of a molten salt graphite core.
基金Project supported by the Ministry of Science and Higher Education of Poland(Nos.04/43/DSPB/0085and 02/21/DSPB/3464)
文摘The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The field of displacements is for- mulated using the classical broken line hypothesis and the proposed nonlinear hypothesis that generalizes the classical one. Using both hypotheses, the strains are determined as well as the stresses of each layer. The kinetic energy, the elastic strain energy, and the work of load are also determined. The system of equations of motion is derived using Hamilton＇s principle. Finally, the system of three equations is reduced to one equation of motion, in particular, the Mathieu equation. The Bubnov-Galerkin method is used to solve the system of equations of motion, and the Runge-Kutta method is used to solve the second-order differential equation. Numerical calculations are done for the chosen family of beams. The critical loads, unstable regions, angular frequencies of the beam, and the static and dynamic equilibrium paths are calculated analytically and verified numerically. The results of this study are presented in the forms of figures and tables.
基金supported by the Innovation Platform for Through Process Modeling and Simulation of Advanced Materials Processing Technologies(No.2012ZX04012011)the National Science Foundation of China(No.51575304)
文摘Core shooting process plays a decisive role in the quality of sand cores, and core box vents distribution is one of the most important factor determining the effectiveness of core shooting process. In this paper, the influence of core box vents distribution on the flow dynamics of core shooting process was investigated based on in situ experimental observations with transparent core box, high-speed camera and pressure measuring system. Attention was focused on the variation of both the flow behavior of sand and pressure curves due to different vents distribution. Taking both kinetic and frictional stress into account, a kinetic-frictional constitutive model was established to describe the internal momentum transfer in the solid phase. Two-fluid model(TFM) simulation was then performed and good agreement was achieved between the experimental and simulated results on both the flow behavior of sand and the pressure curves. It was found that vents distribution has direct effect on the pressure difference of different locations in the core box, which determines the buoyancy force exerting on the sand particles and significantly influences the filling process of core sand.
文摘The elasto-gravitational deformation response of the Earths solid parts to the perturbations of the pressure and gravity on the core-mantle boundary (CMB) and the solid inner core boundary (ICB), due to the dynamical behaviors of the fluid outer core (FOC), is discussed. The internal load Love numbers, which are formulized in a general form in this study, are employed to describe the Earths deformation. The preliminary reference Earth model (PREM) is used as an example to calculate the internal load Love numbers on the Earths surface, CMB and ICB, respectively. The characteristics of the Earths deformation variation with the depth and the perturbation periods on the boundaries of the FOC are also investigated. The numerical results indicate that the internal load Love numbers decrease quickly with the increasing degree of the spherical harmonics of the displacement and depend strongly on the perturbation frequencies, especially on the high frequencies. The results, obtained in this work, can be used to construct the boundary conditions for the core dynamics of the long-period oscillations of the Earths fluid outer core.
基金supported by the National Science Foundation of China(Grant Number 51575304)the National Science and Technology Major Project of the Ministry of Science and Technology of China(Grant Number 2012ZX04012011)
文摘The influence of core sand properties on flow dynamics was investigated synchronously with various core sands, transparent core-box and high-speed camera. To confirm whether the core shooting process has significant turbulence, the flow pattern of sand particles in the shooting head and core box was reproduced with colored core sands. By incorporating the kinetic theory of granular flow(KTGF), kinetic-frictional constitutive correlation and turbulence model, a two-fluid model(TFM) was established to study the flow dynamics of the core shooting process. Two-fluid model(TFM) simulations were then performed and a areasonable agreement was achieved between the simulation and experimental results. Based on the experimental and simulation results, the effects of turbulence, sand density, sand diameter and binder ratio were analyzed in terms of filling process, sand volume fraction(αs) and sand velocity(Vs).
文摘To the reduct problems of decision system, the paper proposes the notion of dynamic core according to the dynamic reduct model. It describes various formal definitions of dynamic core, and discusses some properties about dynamic core. All of these show that dynamic core possesses the essential characters of the feature core.
文摘We give a brief review of the asymptotic theory of slender vortex filaments with emphases on (i) the choices of scalings and small parameters characterizing the physical problem,(ii) the key steps in the formulation of the theory and (iii) the assumptions and/or restrictions on the theory of Callegari and Ting (1978).We present highlights of an extension of the 1978 asymptotic theory:the analyses for core structures with axial variation.Making use of the physical insights gained from the analyses,we present a new derivation of the evolution equations for the core structure.The new one is simpler and straightforward and shows the physics clearly.
基金supported by the National Science Foundation of China(Grant No.42230606)。
文摘Here,a nonhydrostatic alternative scheme(NAS)is proposed for the grey zone where the nonhydrostatic impact on the atmosphere is evident but not large enough to justify the necessity to include an implicit nonhydrostatic solver in an atmospheric dynamical core.The NAS is designed to replace this solver,which can be incorporated into any hydrostatic models so that existing well-developed hydrostatic models can effectively serve for a longer time.Recent advances in machine learning(ML)provide a potential tool for capturing the main complicated nonlinear-nonhydrostatic relationship.In this study,an ML approach called a neural network(NN)was adopted to select leading input features and develop the NAS.The NNs were trained and evaluated with 12-day simulation results of dry baroclinic-wave tests by the Weather Research and Forecasting(WRF)model.The forward time difference of the nonhydrostatic tendency was used as the target variable,and the five selected features were the nonhydrostatic tendency at the last time step,and four hydrostatic variables at the current step including geopotential height,pressure in two different forms,and potential temperature,respectively.Finally,a practical NAS was developed with these features and trained layer by layer at a 20-km horizontal resolution,which can accurately reproduce the temporal variation and vertical distribution of the nonhydrostatic tendency.Corrected by the NN-based NAS,the improved hydrostatic solver at different horizontal resolutions can run stably for at least one month and effectively reduce most of the nonhydrostatic errors in terms of system bias,anomaly root-mean-square error,and the error of the wave spatial pattern,which proves the feasibility and superiority of this scheme.
基金This work was supported by the National Key Research and Development Program of China under Grant No.2017YFC1502203.
文摘The Global-Regional Integrated forecast System(GRIST)is the next-generation weather and climate integrated model dynamic framework developed by Chinese Academy of Meteorological Sciences.In this paper,we present several changes made to the global nonhydrostatic dynamical(GND)core,which is part of the ongoing prototype of GRIST.The changes leveraging MPI and PnetCDF techniques were targeted at the parallelization and performance optimization to the original serial GND core.Meanwhile,some sophisticated data structures and interfaces were designed to adjust flexibly the size of boundary and halo domains according to the variable accuracy in parallel context.In addition,the I/O performance of PnetCDF decreases as the number of MPI processes increases in our experimental environment.Especially when the number exceeds 6000,it caused system-wide outages(SWO).Thus,a grouping solution was proposed to overcome that issue.Several experiments were carried out on the supercomputing platform based on Intel x86 CPUs in the National Supercomputing Center in Wuxi.The results demonstrated that the parallel GND core based on grouping solution achieves good strong scalability and improves the performance significantly,as well as avoiding the SWOs.
基金supported by The National Key Research and Development Program of China(Grants Nos.2017YFA0603901 and 2017YFC1501901)The National Natural Science Foundation of China(Grant No.41522504)。
文摘An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.
基金the Institute of Geodesy and Geophysics, CAS in Wuhanthe National Natural Science Foundation of China
文摘In the last several years since 2004 an important new finding has been unveiled by combined efforts due to Japanese （Satoru Tanaka）, French （Renaud Deguen, Y Albousierre and Marc Monnereau）, American and Chinese geophysicists （Xiaodong Song and Vernon F Cormier） who employed from unambi- guous detailed seismological evidence and explained by clear theoretical and sound laboratory arguments drawn from fluid dynamics that there exists a strong East-West hemi-spherical asymmetry on the inner- outer core boundary,
文摘We numerically study the propagation dynamics of intense optical pulses in gas-filled hollow-core fibers（HCFs）. The spatiotemporal dynamics of the pulses show a transition from tightly confined to loosely confined characteristics as the fiber core is increased, which manifests as a deterioration in the spatiotemporal uniformity of the beam. It is found that using the gas pressure gradient does not enhance the beam quality in large-core HCFs, while inducing a positive chirp in the pulse to lower the peak power can improve the beam quality. This indicates that the self-focusing effect in the HCFs is the main driving force for the propagation dynamics. It also suggests that pulses at longer wavelengths are more suitable for HCFs with large cores because of the lower critical power of self-focusing, which is justified by the numerical simulations. These results will benefit the generation of energetic few-cycle pulses in large-core HCFs.
文摘Most current lattice Boltzmann （LBM） models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular dynamics （MD） simulation and LBM to avoid such defect. The basic idea is to first construct a molecular model based on the actual components of the rock-fluid system, then to compute the interaction force between the rock and the fluid of different densities through the MD simulation. This calculated rock-fluid interaction force, combined with the fluid-fluid force determined from the equation of state, is then used in LBM modeling. Without parameter fitting, this study presents a new systematic approach for pore-scale modeling of multi-phase flow. We have validated this ap- proach by simulating a two-phase separation process and gas-liquid-solid three-phase contact angle. Based on an actual X-ray CT image of a reservoir core, we applied our workflow to calculate the absolute permeability of the core, vapor-liquid H20 relative permeability, and capillary pressure curves.
基金the Chinese Academy of Sciences and the High Technical Project.
文摘Dynamic strain aging (DSA) is an important phenomenon in solutehardened metals and seri- ously affects the mechanical properties ofmetals. DSA is generally induced by the interaction between themoving dislocations and the mobile solute atoms. In this paper, onlythe interaction between moving disloca- tions and mobile solute atomsin a dislocation core area (core atmosphere) will be taken intoaccount. To es- tablish the constitutive model which can describe theDSA phenomenon, we improved the Zerilli-Armstrongdislocation-mechanics-based thermal viscoplastic constitutiverelation, and added the effect of the interaction between the movingdislocations and core atmosphere.
基金supported by the National Key Research and Development Program[grant number 2016YFB0200805]the National Natural Science Foundation of China[grant number 41622503]
文摘How to properly consider the impacts of non-hydrostatic perturbations is one of the challenging issues in developing non-hydrostatic dynamics solvers（NHDSs） for high-resolution atmospheric models. To overcome the drawbacks of current approaches to tackling this issue, this study analyzed the differences between hydrostatic dynamics solvers（HDSs） and their non-hydrostatic counterparts.The analysis then motivated a flexible approach to adjusting existing hydrostatic atmospheric models,especially those adopted in climate simulations for the impacts of non-hydrostatic perturbations.In this approach, the impacts of non-hydrostatic perturbations, reflecting the differences between HDSs and NHDSs, were encapsulated into a single term in the vertical momentum equation for the atmosphere. At each time step, this term was estimated by a separate sub-model, and then it was used to adjust the dynamics of the atmosphere. The adjustment was optional, and could be turned on and off flexibly by utilizing different initial conditions. The approach was illustrated using the Weather Research and Forecasting（WRF） model as an example, and was preliminarily validated by running the 3D baroclinic-wave test case in the model. Results showed that the modified dynamics solver produced simulation results that were very close to those given by the standard NHDS in the WRF model, implying that the approach was basically effective in capturing the non-hydrostatic features of the atmosphere.
基金supported by the National Natural Science Foundation of China(Grant Nos.42375153,42075151).
文摘In relatively coarse-resolution atmospheric models,cumulus parameterization helps account for the effect of subgridscale convection,which produces supplemental rainfall to the grid-scale precipitation and impacts the diurnal cycle of precipitation.In this study,the diurnal cycle of precipitation was studied using the new simplified Arakawa-Schubert scheme in a global non-hydrostatic atmospheric model,i.e.,the Yin-Yang-grid Unified Model for the Atmosphere.Two new diagnostic closures and a convective trigger function were suggested to emphasize the job of the cloud work function corresponding to the free tropospheric large-scale forcing.Numerical results of the 0.25-degree model in 3-month batched real-case simulations revealed an improvement in the diurnal precipitation variation by using a revised trigger function with an enhanced dynamical constraint on the convective initiation and a suitable threshold of the trigger.By reducing the occurrence of convection during peak solar radiation hours,the revised scheme was shown to be effective in delaying the appearance of early-afternoon rainfall peaks over most land areas and accentuating the nocturnal peaks that were wrongly concealed by the more substantial afternoon peak.In addition,the revised scheme enhanced the simulation capability of the precipitation probability density function,such as increasing the extremely low-and high-intensity precipitation events and decreasing small and moderate rainfall events,which contributed to the reduction of precipitation bias over mid-latitude and tropical land areas.
基金the National Natural Science Foundation of China(51675245).
文摘A mathematical model based on the computational fluid dynamics method,heat and mass transfer in porous media and the unreacted shrinking core model for the oxidation reaction of an individual magnetite pellet during preheating was established.The commercial software COMSOL Multiphysics was used to simulate the change in the oxidation degree of the pellet at different temperatures and oxygen concentrations,and the simulated results were compared with the exper-imental results.The model considered the influence of the exothermic heat of the reaction,and the enthalpy change was added to calculate the heat released by the oxidation.The results show that the oxidation rate on the surface of the pellet is much faster than that of the inside of the pellet.Temperature and oxygen concentration have great influence on the pellet oxidation model.Meanwhile,the exothermic calculation results show that there is a non-isothermal phenomenon inside the pellet,which leads to an increase in temperature inside the single pellet.Under the preheating condition of 873-1273 K(20%oxygen content),the heat released by the pellet oxidation reaction in a chain grate is 7.8×10^(6)-10.8×10^(6) kJ/h,which is very large and needs to be considered in the magnetite pellet oxidation modelling.
基金Supported the Leading Academic Discipline Programme （T0104）, the National Natural Science Foundation of China under Grant Nos 10374020 and 10674031, and the Seventh Graduate Students Innovation Foundation of Fudan University.
文摘Two new two-photon absorption （TPA） molecules, named SK-G1 and NT-G1, are synthesized and the photo- physical characteristics are investigated by using linear absorption spectra, one-photon fluorescence spectra and two-photon excited fluorescence spectra. Both the compounds exhibit TPA properties, and the TPA values determined by z-scan measurement are 10 GM and 39 GM for SK-G1 and NT-G1, respectively, at wavelength 80Ohm. Time-resolved spectroscopic techniques are employed to further explore the excited state dynamics of NT-G l with larger TPA cross section. The research results show that there is an ultrafast intraband energy transfer process （about 3ps） before the formation of charge transfer state with a relatively long lifetime.
