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.展开更多
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.展开更多
With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor co...With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.展开更多
A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interface...A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.展开更多
Foundational concepts of trait spaces, including phenotypic plasticity and function of traits, should be expanded and better integrated with ecological theory. This article addresses two areas where plasticity theory ...Foundational concepts of trait spaces, including phenotypic plasticity and function of traits, should be expanded and better integrated with ecological theory. This article addresses two areas where plasticity theory can become further integrated with ecological, evolutionary, and developmental thinking. First is the idea that not only trait means within environments and plasticity of trait means across environments is optimized by selection, but that the entire shape of phenotype distributions such as variance or skew should be optimized within and across environments. In order for trait distribution shape to evolve into adaptations, there must be a genetic basis for and selection upon variation in distribution shapes and their plasticities. I present published and new data demonstrating genetic control and selection for higher moments of phenotype distributions; though, plasticity in these values has not yet been tested. Genetic control of phenotype distribution moments is shown for Neurospora crassa ascospore size and shape. Selection on trait distribution moments is shown for Eurosta solidaginis gall size. Second, there is a tradition in modeling plasti- city as an adaptive strategy that pits it as an alternative to ecological specialization or generaliza- tion. However, these strategies need not be considered alternatives. Rather, with environmental fluctuation within generations plasticity may produce additive or non-additive intermediate （gener- alist） phenotypes, or something new altogether. I present published and new data on the snail Physa virgata and fish Gambusia affinis that show plasticity produces partly intermediate （general- ist） and partly unique phenotypic elements in mixed and fluctuating environments. Plasticity can thus be viewed in the context of a broader trait space and as having broader ecological roles than currently is conceived.展开更多
采用新的均匀三点中心约束多矩有限体积方法(3-point Multi-moment Constrained finite-Volume scheme for Uniform Points with Center Constraints,MCV3_UPCC),发展了一个三阶正定守恒的平流模式。三点多矩有限体积方法在单网格内定...采用新的均匀三点中心约束多矩有限体积方法(3-point Multi-moment Constrained finite-Volume scheme for Uniform Points with Center Constraints,MCV3_UPCC),发展了一个三阶正定守恒的平流模式。三点多矩有限体积方法在单网格内定义等距的3个自由度,采用多矩约束条件并通过控制方程获得时间演变方程。新的三点中心约束多矩方法能在单网格内采用等距的3个点值及中心一阶、二阶导数作为约束条件进行空间4次多项式数值重构,获得3个自由度的时间演变方程;所构建的新数值方案具有三阶精度,边界通量连续性保证了其数值严格守恒。为了抑制该方法的非物理数值振荡,引入了边界保型限制器技术,它能够把数值解控制在既定物理场最小值(最小值为0时则保持数值正定)与最大值之间。数值试验表明新发展的三阶平流模式具有良好的计算精度,能够严格保持数值解的正定性和守恒性,同其他高精度平流模式相当,在实际大气模式水汽等平流输送应用中具备良好的发展潜力。展开更多
基金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 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.
基金Dr.X.L.Li at the China Meteorological Administration.This study was supported by grants from the National Natural Science Foundation of China(Nos.41575103 and 91637210).
文摘With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.
基金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 positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.
文摘Foundational concepts of trait spaces, including phenotypic plasticity and function of traits, should be expanded and better integrated with ecological theory. This article addresses two areas where plasticity theory can become further integrated with ecological, evolutionary, and developmental thinking. First is the idea that not only trait means within environments and plasticity of trait means across environments is optimized by selection, but that the entire shape of phenotype distributions such as variance or skew should be optimized within and across environments. In order for trait distribution shape to evolve into adaptations, there must be a genetic basis for and selection upon variation in distribution shapes and their plasticities. I present published and new data demonstrating genetic control and selection for higher moments of phenotype distributions; though, plasticity in these values has not yet been tested. Genetic control of phenotype distribution moments is shown for Neurospora crassa ascospore size and shape. Selection on trait distribution moments is shown for Eurosta solidaginis gall size. Second, there is a tradition in modeling plasti- city as an adaptive strategy that pits it as an alternative to ecological specialization or generaliza- tion. However, these strategies need not be considered alternatives. Rather, with environmental fluctuation within generations plasticity may produce additive or non-additive intermediate （gener- alist） phenotypes, or something new altogether. I present published and new data on the snail Physa virgata and fish Gambusia affinis that show plasticity produces partly intermediate （general- ist） and partly unique phenotypic elements in mixed and fluctuating environments. Plasticity can thus be viewed in the context of a broader trait space and as having broader ecological roles than currently is conceived.
文摘采用新的均匀三点中心约束多矩有限体积方法(3-point Multi-moment Constrained finite-Volume scheme for Uniform Points with Center Constraints,MCV3_UPCC),发展了一个三阶正定守恒的平流模式。三点多矩有限体积方法在单网格内定义等距的3个自由度,采用多矩约束条件并通过控制方程获得时间演变方程。新的三点中心约束多矩方法能在单网格内采用等距的3个点值及中心一阶、二阶导数作为约束条件进行空间4次多项式数值重构,获得3个自由度的时间演变方程;所构建的新数值方案具有三阶精度,边界通量连续性保证了其数值严格守恒。为了抑制该方法的非物理数值振荡,引入了边界保型限制器技术,它能够把数值解控制在既定物理场最小值(最小值为0时则保持数值正定)与最大值之间。数值试验表明新发展的三阶平流模式具有良好的计算精度,能够严格保持数值解的正定性和守恒性,同其他高精度平流模式相当,在实际大气模式水汽等平流输送应用中具备良好的发展潜力。
