ROMS for highly stratified fuild

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
ashkenaz
Posts: 5
Joined: Sun Apr 18, 2004 9:58 pm
Location: Ben Gurion University

ROMS for highly stratified fuild

#1 Unread post by ashkenaz »

Hi,

I'm trying to simulate the dynamics of highly stratified fuild. I ran ROMS-2.2 with the following initial conditions: for the salinity I sued 10 PSU for the upper 5m while the rest bottom water had salinity of 278PSU. These values are characteristics salinities for the Dead Sea I'm simulating. The temperature is decreasing monotonically from 30C to 23C. No forcing was applied.
With this configuration the model blows up very fast. When using more moderate salinity gradient, or when the salinity gradient was smooth, the model runs smoothly.
Any suggestions how to handle such extreme case ??

Thanks

Yossi

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

#2 Unread post by kate »

The first thing I would check are the resulting density values. The equation of state we are using has a valid range close to that of the world oceans, not including the Dead Sea. Do you have a source for a more comprehensive equation of state? Remember we are computing in situ density from salinity and potential temperature.

ashkenaz
Posts: 5
Joined: Sun Apr 18, 2004 9:58 pm
Location: Ben Gurion University

#3 Unread post by ashkenaz »

Hi Kate,

I'm using linear equation of state with known thermal and haline coefficients for the Dead Sea. After few time steps, the model produced negative (!) salinities and negative relative densities. Constraining the salinities and relative densities to be postive in rho_eso.F didn't result in stable run. The model was also unstable when using nonlinear eos.

dupontf
Posts: 6
Joined: Fri May 02, 2003 2:54 pm
Location: Quebec-Ocean

#4 Unread post by dupontf »

your time step may be limited by strong baroclinic effects (not just the external gravity phase speed). Have you tried to compute the baroclinic phase speed for such a stratified system?

Fred.

ashkenaz
Posts: 5
Joined: Sun Apr 18, 2004 9:58 pm
Location: Ben Gurion University

#5 Unread post by ashkenaz »

The baroclinic phase speed is around 3m/s while the tracers time step is 25s; this time step is sufficiently small, isn't it ?

dupontf
Posts: 6
Joined: Fri May 02, 2003 2:54 pm
Location: Quebec-Ocean

#6 Unread post by dupontf »

25 s for tracer does seem small. what is your spatial resolution and total depth?

ashkenaz
Posts: 5
Joined: Sun Apr 18, 2004 9:58 pm
Location: Ben Gurion University

#7 Unread post by ashkenaz »

The spatial resolution is 0.5km and the maximum depth is ~320m.

dupontf
Posts: 6
Joined: Fri May 02, 2003 2:54 pm
Location: Quebec-Ocean

#8 Unread post by dupontf »

well, with a max depth of 320m, the external time step should be below 4s but with a baroclinic phase speed of 3m/s, you should be able to use a tracer timestep of up to 80s. Therefore 25s should be fine. i am at loss. My guess would be to verify your baroclinic phase speed calcultion. I may verify your calculation from a given density profile.
what is you density profile? just give me strategic depths since it is mainly linear in the deep part.

ashkenaz
Posts: 5
Joined: Sun Apr 18, 2004 9:58 pm
Location: Ben Gurion University

#9 Unread post by ashkenaz »

For the simplest case, the upper 5m have density of about 1000kg/m^3 while the rest 315m have density of ~1230kg/m^3; thus the reduced gravity is about 2.5.

dobbins

steep gradients

#10 Unread post by dobbins »

Hi;

I have had bad luck with ROMs and steep vertical gradients. I have also gotten negative salinities at river mouths, and negative populations in biology when the surface becomes nutrient depleted (another strong gradient).

I think that problem is due to the vertical advection or diffusion scheme, which use vertical spline functions in their calculations. Forgive me for being vague since I had this problem several years ago, and can't remember the details. Perhaps someone with more knowledge of the numerics could suggest exactly where the problem might be.

To solve my problem, I trapped negative populations in the my biology routine and set them to zero, and distributed my river mouths over several horizontal points. Those are both hacks which don't deal with the fundamental problem, but they got my problems to work.

csherwood
Posts: 39
Joined: Fri Apr 02, 2004 4:46 pm
Location: USGS, Woods Hole, USA

Try the MPDATA advection scheme

#11 Unread post by csherwood »

I disagree with Liz...the vertical advection scheme is awesome and can run at 5 - 10 x CFL. However, the horizontal advection scheme may be the problem...try using MPDATA, which is positive definate.
Chris Sherwood, USGS
1 508 457 2269

jklinck
Posts: 34
Joined: Mon Jun 30, 2003 2:29 pm
Location: CCPO/ODU, USA

stratified problem

#12 Unread post by jklinck »

A baroclinic speed of 3 m/s seems pretty small for strong stratification. Using 320 m depth and a reduced gravity of 2.5 gives an internal wave speed of 28.3 m/s. With a grid spacing of 500 m, I get an internal time step of about 4 sec (making the courant number =1/4, which may not be small enough).

I suggest setting the step ratio to 1.0, find a stable external time step and see if the model will run. Then try larger step ratios. My guess is that the internal and external speeds will differ by a factor of 2 or so (g=9.8, g^prime = 2.5, so the speed ratio is \sqrt (g/g^prime)= 1.97).

Before redesigning advection schemes, try some smaller time steps.

John

marcelcure
Posts: 9
Joined: Fri Oct 22, 2004 8:32 pm
Location: Numerics Warehouse Ltd.

#13 Unread post by marcelcure »

We have had problems with riverine input causing numerical instability and high salinities. Initially we just used the the interfacial wave speed to calculate the Courant Number and set the baroclinic time step accordingly - but still had problems. We found that we had neglected the effect the volume flux from rivers had due to mass continuity. We assume all the river flux is contained in an upper layer moving away from the river mouth and then add this speed to the interfacial wave speed. We use this combined velocity and a Courant number of 0.1 to calculate the (now exremely small) baroclimic time step. This is why spreading the river input over a horizontal strip (rather than just one point) helps to reduce the problem. Changing horizontal advection schemes, vertical turbulence parameterization, barotropic time steps, horizontal diffusion, spreading the river input over a greater depth, artifically mixing some of the fresh flux with salt water (increasing the total flux) - none of these helped us.

If you are spinning the model up from crude initial conditions, then it helps to 'flood' the region of river mouths with fresh water in the initial conditions - or you can spin up with a much reduced time step for a while and then restart with a bigger time step.

Marcel

Post Reply