Question about turbulent kinetic energy

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
nmichelet
Posts: 31
Joined: Fri Oct 02, 2015 2:24 pm
Location: Cerema

Question about turbulent kinetic energy

#1 Unread post by nmichelet »

Hello everyone,
I would like to know if it is possible to define a turbulent kinetic energy profile as an input into ROMS.
Thanks

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

Re: Question about turbulent kinetic energy

#2 Unread post by kate »

To what purpose? Would this be part of a vertical mixing scheme? It is certainly possible to modify ROMS, telling it to read new fields, but then what? Is it an initial condition? Something that needs updating from outside? Something that needs to be advected?

nmichelet
Posts: 31
Joined: Fri Oct 02, 2015 2:24 pm
Location: Cerema

Re: Question about turbulent kinetic energy

#3 Unread post by nmichelet »

Thank you Kate for your answer.

The purpose of this is to define in a channel the initial condition for the velocity and the turbulent kinetic energy. I have some data at the enter of my channel so I would like to put those as input.

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

Re: Question about turbulent kinetic energy

#4 Unread post by kate »

You should look at get_state.F to see how and when it reads TKE already. It is protected within:

Code: Select all

#  if defined GLS_MIXING || defined MY25_MIXING
#endif
and also:

Code: Select all

IF (get_var(idMtke).and.Perfect3D) THEN
END IF
Clearly, for a new run, it is getting TKE some other way (initialized to zero??). Just do some snooping to figure out what to change.

Post Reply