Development of a Dynamical Core for the Next-Generation Atmospheric Meso-Scale Numerical Model

Development of a Dynamical Core for the Next-Generation Atmospheric Meso-Scale Numerical Model

Development of an urban meteorological numerical model

In Cartesian coordinate

余偉明(Weiming Sha)

Geophysical Institute, Graduate School of Science, Tohoku University

Aoba-Ku, Sendai, 980-8578, Japan

Tel.: 022-217-5783; E-mail:

1. Introduction

Since terrain-following vertical coordinate (sigma) system (Phillips 1957; Gal-Chen and Somerville 1975) has been used extensively to accommodate orography in models for atmospheric flows, most of existing community meso-scale atmospheric numerical model in the word are using the terrain-following coordinate as the vertical coordinate. However, a problem that has received attention rather early in the development of sigma system primitive equation models is that of the noncancellation errors in the two terms of the gradient force in the momentum equation (Smagoringsky, et al. 1967). Mesinger and Janjic(1985), among others, have found that errors in computing the horizontal pressure gradient force in models using a sigma coordinate can be substantial in the vicinity of steep topography. To minimize this error, a step-mountain vertical coordinate, the so-called “eta coordinate”, is implemented in the National Centers for Environmental prediction (NCEP) Meso Eta Model (Mesinger et al. 1988) in which the topography is represented as discrete steps (step mountain). Recently, representation of topography, i.e., the “shaved cell”approach, and the numerical scheme for the equations of geophysical flows, in which the height is used as vertical coordinate, have been proposed (Adcroft et al. 1997; Marshall et al. 1997; Bonaventura 2000). It is reasonable to expect that a goal of running a regional model at a very high horizontal resolution may be attainable in the near future, and the topography and objects (e.g., buildings in urban city) could be more accurately represented. In such a situation, it seems natural to search for an alternative that will be better suited to handle the step topography and complex objects on the surface for the high-resolution models currently used as well as a future model expected to run with a finer resolution.

2.Concept of the numerical framework

In this numerical development, Finite Volume Method (FVM) in conjunction with the SIMPLER (Semi-Implicit Method for Pressure-Linked Equation Revised; Patankar 1980) algorithms is used for calculations of the unsteady, three dimensional, compressible Navier-Stokes equations on a staggered grid. Abandoning the customary terrain-following normalization, the Cartesian coordinate, in which the height is used as the vertical one, is chose. A Cartesian-grid system approach, which consists of the variable regular cells and a special treatment of the boundary cells, is proposed for expression of the arbitrarily complex geometries. Blocking-off Method (Patankar 1980) is then introduced to handle the steep topography and complex objects above the Earth’s sea-level surface, thus resulting in a ;000,enturarobust and efficient numerical scheme which allows for applications to meso-scale flows over complex orography. The spatial discretization is obtained by a finite volume technique on the staggered grid, and higher-order upwind convection scheme is employed to relate the flux at each control volume face. For the temporal integration of the equation, the fully time implicit scheme is utilized. As the fully implicit temporal discretization is used,the time step can be determined only physical criteria and accuracy considerations.

3. Some results of cases test

As a test, the model has been run on calculatingflows over cube/steep mountain by Direct Numerical Simulation (DNS), and turbulent flow in an urban city, i.e., Ootemati, Tokyo (Fig.1) by Large Eddy Simulation (LES), respectively. Detailed results will be presented

in the meeting.

4. Reference

(1) Phillips, N. A., 1957, J. Meteor. 14, 184-185.

(2) Gal-Chen, T., and R. Somerville, 1975, J. Comput. Phys., 17, 209-228.

(3) Smagorinsky, J., J.L.Holloway and G.D. Hembree, 1967, Proc. Inter. Symp. Dynamics Large Scale Atmospheric Processes, 70-134.

(4) Mesinger, F. and Z.I.Janjic, 1985, Large Scale Computations in Fluid Mechanics, Amer. Math. Soc., 81-120.

(5) Mesinger, F., S.Nickovic, D.Gavrilov and D.G.Deaven, 1988, Mon. Wea. Rev., 116, 1493-1518.

(6)Adcroft, A., C.Hill and J. Marshall, 1997, Mon., Wea. Rev., 125,2293-2315.

(7) Marshall, J., A.Adcroft, C..Hill, L.Perelman and C. Heisey, 1997, J. Geophys., Res., 102, 5753-5766.

(8) Bonaventura, L., 2000, J. Comput. Phys., 158, 186-213.

(9) Patankar, S.V., 1980, Numerical Heat and Transfer and Fluid Flow.

Fig. 1 Turbulent urban flow in Ootemati, Tokyo.