A Neural-network-based Alternative Scheme to Include Nonhydrostatic Processes in an Atmospheric Dynamical Core

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.

An Improved Dynamic Core for a Non-hydrostatic Model System on the Yin-Yang Grid

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.

Improved Diurnal Cycle of Precipitation on Land in a Global Non-Hydrostatic Model Using a Revised NSAS Deep Convective Scheme

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.

Numerical study of the dynamic characteristics of a single-layer graphite core in a thorium molten salt reactor

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.

Influence of core box vents distribution on flow dynamics of core shooting process based on experiment and numerical simulation

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.

Mathematical modeling for dynamic stability of sandwich beam with variable mechanical properties of 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 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.

Earth's deformation due to the dynamical perturbations of the fluid outer core

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.

A 3D Nonhydrostatic Compressible Atmospheric Dynamic Core by Multi-moment Constrained Finite Volume Method

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.

Influence of core sand properties on flow dynamics of core shooting process based on experiment and multiphase simulation

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).

The Feature Core of Dynamic Reduct

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.

On the Dynamics of Vortex Filaments and Their Core Structure

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.

Parallelization and I/O Performance Optimization of a Global Nonhydrostatic Dynamical Core Using MPI

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 Nonhydrostatic Atmospheric Dynamical Core Using a Multi-Moment Constrained Finite Volume Method

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.

Spatiotemporal propagation dynamics of intense optical pulses in loosely confined gas-filled hollow-core fibers

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.

From molecular dynamics to lattice Boltzmann:a new approach for pore-scale modeling of multi-phase flow

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.

A DISLOCATION-MECHANICS-BASED CONSTITUTIVE MODEL FOR DYNAMIC STRAIN AGING

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.

Two-Photon Absorption and Excited State Dynamics of Two Organic Molecules Containing Donor/Acceptor Moieties

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.

核心稳定性训练降低落地动作前交叉韧带损伤的风险

背景:研究表明,动态姿势控制能力较差可能导致运动时产生异常的动作模式,进而可能增大下肢关节及前交叉韧带的损伤风险,而躯体核心的稳定是良好动态姿势控制能力的基础。目的:探讨核心稳定性训练对大学生动态姿势控制能力及落地动作损伤风险的影响,并对比不同性别训练效果差异。方法:招募35名(男=19,女=16)大学生为试验对象,对其进行6周核心稳定性训练,分析训练前后Y平衡测试、躯干伸肌耐力测试、躯干屈肌耐力、侧桥肌群耐力测试和落地错误评分系统(LESS)测试结果的差异。结果与结论:①6周核心稳定性训练可提高男女大学生的躯干伸肌耐力(P<0.001)、屈肌耐力(P<0.001)和侧桥肌群耐力(P<0.001);②核心稳定性训练可提高男女大学生Y平衡测试前向测试距离(P=0.026)、后内向距离(P<0.001)、后外向距离(P=0.005)和综合得分(P<0.001);③经过6周核心稳定性训练,男女大学生LESS评分均显著下降(P<0.001),同时可增大触地时刻膝关节(P<0.001)、髋关节屈曲角(P<0.001),减小触地时刻膝关节外翻角(P<0.001),并可增大最大屈膝角(P<0.001)和减小最大膝外翻角(P<0.001);④结论:核心稳定性训练可提高动态姿势控制能力并改善落地动作模式,提示可能有助于降低前交叉韧带损伤风险;核心稳定性训练在提高躯干屈肌耐力、侧桥肌群耐力、改善动态姿势控制能力以及降低前交叉韧带损伤风险方面无性别差异。

p53活性四聚体全原子分子动力学分析

p53是一种肿瘤抑制蛋白,对阻碍癌症发展、维持遗传完整性起着至关重要的作用.在细胞核内,4个p53分子通过高度协同的方式、通过DNA结合域与DNA结合,形成稳定的四聚体活性结构,并转录激活或抑制其靶向基因.然而,大多数肿瘤细胞中存在大量p53的突变,其中绝大部分突变发生在p53的DNA结合域,而p53的DNA结合域又是p53形成四聚体活性结构、调控下游靶基因转录的重要区域.本文通过全原子分子动力学模拟,研究了野生型p53四聚体内分子间的相互作用机制.结果表明,位于DNA两侧的对称二聚体是一个稳定的二聚体,在与DNA结合前后都能维持稳定的结构.位于DNA同侧的两个单体依靠两个接触面提供的蛋白-蛋白相互作用和DNA的骨架作用使四聚体活性结构保持稳定,这些相互作用为四聚体的形成机制提供了重要支撑.该工作厘清了p53四聚体在动力学过程中的内部相互作用机制和关键残基,揭示了四聚化过程中各个相互作用界面的关键位点,对于理解p53的抑癌机制、探索有效治癌策略、发展治癌药物具有重要意义.