Design of ALADIN NH
This presentation is the property of its rightful owner.
Sponsored Links
1 / 24

Content of presentation PowerPoint PPT Presentation


  • 50 Views
  • Uploaded on
  • Presentation posted in: General

Design of ALADIN NH with Vertical Finite Element Discretisation Jozef Vivoda , SHM Ú Pierre B énard, METEO France Karim Yessad, METEO France Emeline Larrieu Rosina, SHMÚ. Content of presentation. Introduction VFE – basic features VFE derivative and integration operator

Download Presentation

Content of presentation

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Content of presentation

Design of ALADIN NH with Vertical Finite Element Discretisation Jozef Vivoda, SHMÚPierre Bénard, METEO FranceKarim Yessad, METEO FranceEmeline Larrieu Rosina, SHMÚ

5-6 December, 2006


Content of presentation

Content of presentation

  • Introduction

  • VFE – basic features

  • VFE derivative and integration operator

  • Discretization of linear model

  • Analysis of stability of 2TL ICI NESC scheme

  • First 2D idealized tests

  • Conclusions and future plans

5-6 December, 2006


Introduction

Introduction

Motivation:

  • VFE scheme successfully implemented into HY model (Untch and Hortal,2004)

  • To extent the VFE scheme into NH model with HY model as a limit case

  • VFE has two times higher accuracy than FD method on the same stencil (eight order for cubic elements in the case of evenly spaced mesh)

  • More accurate vertical velocities for SL scheme

    Work done so far:

  • Benard – study of accuracy of vertical integral operators and of compatibility of VFE with existing model choices

  • Vivoda – linear analysis of stability of VFE scheme

  • Vivoda, Yessad –start of implementation of VFE into ALADIN code

5-6 December, 2006


Vfe scheme basic features

VFE scheme basic features

Starting point – Untch and Hortal,2004:

  • Basis functions – cubic B-splines with compact support

  • No staggering – all variables are defined on model full levels

  • Only integration/derivation is performed in FE space, the products of variables are done in physical space

  • In SL version of HY model (ECMWF or ARPEGE) only non-local operations in the vertical are integrations. In NH version derivatives plays crucial role (structure equation contains vertical laplacian).

5-6 December, 2006


Fe derivative operator

FE derivative operator

We expand F and f in terms of chosen set of functions:

Truncation error R is orthogonalized (required to vanish in weighted integral sense) on interval (0,1):

Finally we incorporate the transforms from physical to FE space and back:

5-6 December, 2006


Fe derivative operator basis functions

L+1

L+3

top

L-1

L+4

L

L+2

surf

FE derivative operator basis functions

Example for 35 unevenly spaced levels

  • In order to cover the physical domain (eta=<0,1>) with the full set of basis functions the 4 additional functions must be determined, the two at the bottom and the two at the top (analogical approach as Untch and Hortal, 2004).

5-6 December, 2006


Fe derivative operator accuracy

FD

FD

VFE

VFE

FE derivative operator accuracy

The 4 additional assumptions:

5-6 December, 2006


Fe integration operator

FE integration operator

We expand F and f in terms of chosen set of functions:

Truncation error R is orthogonalized (required to vanish in weighted integral sense) on interval (0,1):

Finally we incorporate the transforms from physical to FE space and back:

In: N values of f

Out: N+1 values of integral

N on full levels + integral from

Model top to surface

5-6 December, 2006


Fe integration operator basis functions

FE integration operator basis functions

Untch and Hortal, 2004

-3

-1

0

1

L

L+1

L+2

L+3

L+4

5-6 December, 2006


Vfe scheme linear isothermal nh model

Vertical operators

Linear system

Prognostic var.

Mass weighted vertical integral operators (already in HY model)

Mass weighted vertical laplacian operator (NH specific)

VFE scheme – linear isothermal NH model

System describing small aplitude waves around reference state

5-6 December, 2006


Constraint c1 on vertical operators

Constraint C1 on vertical operators

  • When eliminating the linear system to obtain single variable Helmholtz solver the following two equations for (d4,D) are obtained:

Constrain C1 is 0 in continuous case:

  • In discrete case C1 is NOT automatically 0. The FD vertical discretization is designed in order to satisfy C1. If C1 is not satisfied further elimination is not feasible.

  • In order to be consistent with existing HY and NH code the C1 constrain shall be satisfied.

5-6 December, 2006


Constraint c2 on vertical operators

Constraint C2 on vertical operators

When C1 is satisfied we could eliminate D from the system. It gives the structure equation

With C2 constraint:

  • In continuous case T* is an identity matrix. In FD case T* is the tridiagonal matrix.

  • For stability reasons of SI time stepping algorithm we require:

    • L* to have real and negative eigenvalues

    • T* to have positive and real eigenvalues

  • In the case of C1 unsatisfied above conditions could be insufficient, we we still require them assuming that C1 is close to zero

  • 5-6 December, 2006


    Laplacian term fe treatment

    Laplacian term FE treatment

    In linear and nonlinear model the laplacian has the same form (thanks to fact

    that we use vertical divergence related prognostic variable):

    V2:

    V1:

    • Boundary conditions the same as in FD model version

    5-6 December, 2006


    Laplacian term fe treatment1

    Imaginary part of eigenvalues

    Real part of eigenvalues

    Laplacian term FE treatment

    V2:

    V1:

    Eigenvalues of laplacian for 100 levels

    Green – eigenvalues of V1

    Red – eigenvalues of V2

    Imaginary part multiplied by 100

    5-6 December, 2006


    Stability of 2tl non extrapolating iterative centered implicit scheme

    Stability of 2TL Non-Extrapolating Iterative Centered Implicit scheme

    Predictor:

    Nth-corrector:

    For the purpose of analysis of stability we linearise the nonlinear residual R

    around resting, isothermal, hydrostatically balanced reference state:

    We analyse the temperature instability due to fact that

    The following analyses were computed for the wave with dx=5km, dt=600s, 35 vertical levels, T*=350K, T*a=100K, p*s=101325Pa (perfect pressure)

    5-6 December, 2006


    C1 constrain modification of g

    Solid – VFE Dotted - FD

    SI

    N=1

    N=2

    N=3

    C1 constrain – modification of G*

    • from C1 constrain we redefine the operator G*. It must be computed via iterative procedure.

    In linear and nonlinear model

    In linear model only

    5-6 December, 2006


    C1 constrain modification of s

    Stability – C1 satisfied for NL model as well

    Stability – C1 satisfied in linear model only

    Solid – VFE Dotted - FD

    SI

    N=1

    N=2

    N=3

    C1 constrain – modification of S*

    • Directly from C1 constrain we can define operator S*

    The time stepping is unstable, because eigenvalues of L*A2 are complex.

    5-6 December, 2006


    C1 constrain unsatisfied

    Eigenvalues of C2 constraint

    Black – no C1

    Red - S* redefinition in lin and nonlin

    Green – G* redefinition in lin and nonlin

    Imaginary part multiplied by 10

    Circle radius 0.4

    Solid – VFE Dotted - FD

    SI

    N=1

    N=2

    N=3

    C1 constrain unsatisfied

    Dt=600s

    Dt=999999s

    • Solution of 2Lx2L soler for each spectral component – noncompatibility with existing model (LxL solver), problem with efficient introduction of map factor horizontal dependence (LSIDG resp. LESIDG key)

    • Yessad – proposal of iterative procedure, build from existing pieces of code, study of convergence required

    5-6 December, 2006


    Vfe scheme first tests

    VFE scheme – first tests

    • Eta coordinate

    • NLNH flow regime

    • dx=200m

    • V=10ms-1

    • Agnesi hill, h=1000m, L=1000m

    • CY30T1

    • Tested in 2D

    • 2TL ICI NESC n=1

    • sponge

    Ref VFE, 100lev, dt=1s

    (X-term, Z-term FD)

    Ref FD, 100lev, dt=5s

    5-6 December, 2006


    Vfe scheme first tests1

    VFE scheme – first tests

    Ref VFE, 35lev, dt=5s

    X-term, Z-term-FD

    Ref FD, 35, dt=5s

    5-6 December, 2006


    Vfe scheme first tests2

    VFE scheme – first tests

    Ref VFE, 35lev, dt=5s

    X-term FD, Z-term-VFE

    Ref VFE, 35lev, dt=5s

    X-term VFE, Z-term-FD

    5-6 December, 2006


    Vfe scheme conclusions

    VFE scheme – conclusions

    • Analysis of stability in linear framework suggests that it is possible to extend HY model VFE scheme to NH model

    • From stability point of view the crucial is the definition of vertical laplacian operator (all eigenvalues must be real and negative)

    • The VFE scheme is stable in linear context when FD boundary conditions are used in VFE (we can adopt BBC and TBC from FD model), but oscillatory behaviour near model bottom is observed (this is not true for large number of levels !)

    • So far no acceptable stable solutions that satisfy C1 was found

    5-6 December, 2006


    Vfe scheme future work

    VFE scheme – future work

    • Non-oscillatory discretization of vertical laplacian

    • Extension of just described VFE discretization concept into full NL model

    • Coding of “draft” of VFE scheme in NH ALADIN with 2Lx2L solver (finished)

    • testing

    5-6 December, 2006


    Content of presentation

    Thank you for your attention !

    Ďakujem za pozornosť !

    Merci de votre attention !

    Hvala !

    5-6 December, 2006


  • Login