Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Thu Oct 18, 2018 3:44 am




Post new topic Reply to topic  [ 8 posts ] 

All times are UTC

Author Message
PostPosted: Thu May 31, 2018 5:27 am 
Offline
User avatar

Joined: Mon May 05, 2003 2:41 pm
Posts: 110
Location: The University of Western Australia, Perth, Australia
Dear ROMS users,

before duplicating already done job I am wondering if someone already did it?

I am looking into the additional mixing (via k-e turbulence) with NONBREAKING(!) waves as additional term in the
shear production term (along with shear and bouy2) as a function of depth, Hs and frequency (Aijaz et. al., 2017; Babanin & Chalikov 2012)

Aijaz, S., M. Ghantous, A. V. Babanin, I. Ginis, B. Thomas, and G. Wake (2017),
Nonbreaking wave-induced mixing in upper ocean during tropical cyclones using coupled hurricane-ocean-wave modeling, J. Geophys. Res. Oceans, 122, 3939–3963, doi:10.1002/ 2016JC012219.

Babanin, A. V., and D. Chalikov (2012), Numerical investigation of turbulence generation in non-breaking potential waves, J. Geophys. Res.,
117, C00J17, doi:10.1029/2012JC007929.


Thanks
Ivica


Top
 Profile  
Reply with quote  
PostPosted: Thu May 31, 2018 6:00 am 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3515
Location: IMS/UAF, USA
Hi Ivica - Good luck! I stole some code from MOM6 for Langmuir mixing in KPP that I'm playing with. Seems to be an active area of research!


Top
 Profile  
Reply with quote  
PostPosted: Fri Jun 01, 2018 5:54 am 
Offline
User avatar

Joined: Fri Apr 30, 2004 6:41 pm
Posts: 24
Location: IRD, FRANCE
For wave‐induced vertical mixing in KPP, see Wang et al., JGR 2010, based on Qiao et al, JPO 2004. It is now part of CROCO (ROMS_AGRIF)


Top
 Profile  
Reply with quote  
PostPosted: Fri Jun 01, 2018 3:07 pm 
Offline

Joined: Wed Jan 02, 2008 3:15 pm
Posts: 78
Location: University of Copenhagen
Yes, it's done by Wang et al., JGR 2010 based on Qiao et al, JPO 2004. I added it as a one of the switches (option) for vertical turbulent mixing :arrow: (lmd_vmix) to the CPP options list (ROMS 3.7), and is running now for the Mid-Atlantic Bight application :idea: . To extract wave induced mixing parameter (Bv: m2/s) user need to run the MASNUM model.Please take a look to the attached file.

Cheers,
Farshid


Attachments:
bv.pdf [19.05 MiB]
Downloaded 44 times
Top
 Profile  
Reply with quote  
PostPosted: Fri Jun 01, 2018 4:34 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3515
Location: IMS/UAF, USA
Some of us are too lazy to run a wave model. Check out:
Quote:
! Li, Q. and B. Fox-Kemper, 2017: Assessing the effects of Langmuir !
! turbulence on the entrainment buoyancy flux in the ocean surface !
! boundary layer, JPO, DOI: 10.1175/JPO-D-17-0085.1 !
! !
! Li, Q., B. Fox-Kemper, Ø. Breivik, and Adrean Webb, 2017: !
! Statistical models of global Langmuir mixing, Ocean Modelling, !
! DOI: 10.1016/j.ocemod.2017.03.016


Top
 Profile  
Reply with quote  
PostPosted: Sat Jun 02, 2018 5:33 am 
Offline
User avatar

Joined: Fri Apr 30, 2004 6:41 pm
Posts: 24
Location: IRD, FRANCE
Hi Kate, note that nonbreaking-wave-induced mixing and Langmuir turbulence are generally considered two different processes of wave-induced mixing (related to Stokes drift)... Patrick


Top
 Profile  
Reply with quote  
PostPosted: Thu Jun 07, 2018 6:33 am 
Offline
User avatar

Joined: Mon May 05, 2003 2:41 pm
Posts: 110
Location: The University of Western Australia, Perth, Australia
Thanks for all responses!

Farhid in which version of ROMS you've implemented it?
Looking in Rutger's version and lmd_vmix.F there is only diffusivity due to the internal waves and not because of surface wind induced waves?
lmd_vmix.F
Code:
 
!  Compute interior diffusivity due to wave breaking (Gargett and
!  Holloway.


Patrick,
you have put non-breaking surface wind wave effect into kpp in CROCO? Please if you can point me into which subroutine I should look.
My intention is to put it into gls, it should go into shear production term as addition "W"
\begin{equation}
\ni (M^2 + \alpha N^2) + W
\end{equation}

If you have information about wind induced waves from spectral model (or via coupling) it should be straight forward.
W is simple function of H_S is significant wave height, \omega is peak angular frequency, k is wavenumber, z is depth:
\begin{equation}
W = 0.0014 * k (\frac{\omega H_s}{2) e^{kz})^3
\end{equation}

Thanks!

Ivica


Top
 Profile  
Reply with quote  
PostPosted: Fri Jun 08, 2018 11:26 am 
Offline

Joined: Wed Jan 02, 2008 3:15 pm
Posts: 78
Location: University of Copenhagen
Dear Ivica,

Please have a look to "lmd_vmix.F " and the section, regarding to "Sum contributions due to internal wave breaking...", in fact non-breaking wave induced coefficient (Bv(i,j,k); (m2/s)) here is used to enhance the vertical mixing for momentum (Akv(i,j,k); (m2/s)) and tracer (Akt(i,j,k); (m2/s)) respectively. :arrow:

Akv* = Akv + Bv (1)
Akt* = Akt + cBv (2)

where "c" is weighting coefficient of wave‐induced mixing (Yonggang Wang et al. 2010).

Cheers,
Farshid


Attachments:
Yonggang Wang et al 2010.pdf [12.2 MiB]
Downloaded 37 times
Qiao_et_al-2004-Geophysical_Research_Letters.pdf [1.19 MiB]
Downloaded 19 times
Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 8 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 3 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group