brodiepearson / mpas-model Goto Github PK
View Code? Open in Web Editor NEWThis project forked from vanroekel/mpas-model
Repository for MPAS models and shared framework releases.
License: Other
This project forked from vanroekel/mpas-model
Repository for MPAS models and shared framework releases.
License: Other
The following ADC tendency budgets require:
For example with reogranization, sometimes the parameterized pressure terms are spread across multiple ADC tendency terms where they are mixed with other processes, or they appear implicitly (in the case of the slow Rotta term). In contrast, LES and RANS they are associated with a single pressure-strain term in each budget.
After this, perhaps we should also comment the sub-plume scale budgets:
I think there are a couple of issues in the sub-plume terms in the
The Lappen and Randall papers present the sub-plume terms for the
where
The ADC scheme looks like it is applying sub-plume terms that are different from the above formulations. The
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 459 to 462 in aff0c19
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 652 to 658 in aff0c19
If the above is correct, the easiest fix would be to
Swumd
andwttend6
and wstend6
I think there are a couple of typos in several terms that need to be fixed. These typos appeared because of the switch from atmosphere to ocean orientation, which has led to some terms being flipped correctly, and others not.
Typos I potentially identified:
The
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 652 to 658 in 1a3f760
Same as above, but with
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 708 to 714 in 1a3f760
I think there is an issue with the sign of Swumd, where it should be upwelling minus downwelling, but it has been defined the other way around. Remove the minus sign at the begining of this line:
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 459 to 462 in 1a3f760
There may be a sign error in the
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 490 to 491 in 1a3f760
I think these issues stem from conversion from atmosphere to ocean systems (intense upwelling to intense downwelling). Particularly, areaFraction in the new scheme represents downwelling area as it is <0.5 when skewness is negative due to ocean convection, versus upwelling area in the Lappen and Randall papers.
The one place where the choice of
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 459 to 462 in 1a3f760
MPAS-Model/src/core_ocean/shared/mpas_ocn_adcReconstruct.F
Lines 653 to 658 in 1a3f760
The effects of Stokes drift on upper ocean mixing need to be incorporated in the ADC scheme. These contributions can be separated into
The wave-averaged equations introduce additional forces due to the Stokes drift (
We are interested in how these forces affect the turbulence statistics (fluxes, variances etc.). The quadratic Stokes term and the Coriolis-Stokes term disappear in the fluctuating momentum budgets (assuming the Stokes term is steady in time and horizontal space). The only effect that Stokes drift has in tracer (
The pressure field within a turbulent flow is coupled to the momentum equations. By taking the divergence of the (vector) momentum budget and using incompressibility, an equation relating the turbulent pressure field to the other forces that drive the turbulent flow can be derived. This means that each force in the flow has an imprint on the turbulent pressure field. As a result the pressure-strain (
Pearson et al., 2019 parameterized the effect that Stokes drift has on the pressure-strain terms in Langmuir turbulence, using a suite of LES experiments. This parameterization was:
where
The pressure-scalar terms have not been investigated in wave-affected flows. The only attempt to parameterize them was in Harcourt, 2013 where they were assumed to be identical in form to the Eulerian shear-based pressure-scalar terms:
I suspect that the above pressure-scalar closure is not appropriate for Langmuir turbulence, as it simply replicates exactly the same functional form for both the shear-driven and Stokes-driven pressure effects, with Eulerian and Stokes velocities being interchanged between the two.
However a closer look at the the tracer flux budgets (below) shows that the rapid closure exactly cancels a fraction of the shear-production of tracer fluxes. There is no shear production of the vertical tracer flux (
We will use the simplest form of the pressure-scalar closure that allows both of these closure options, which is based on Mironov, 2001 . This closure requires two closure constants for Stokes pressure-scalar closure (and analogous constants for the rapid closure, which is already coded into the ADC scheme):
where the Harcourt 2013 closure would correspond to
This form of closure for the Stokes pressure-scalar term also allows a range of other closures. For example, Canuto et al. (2001) suggested that flow constraints and idealized turbulence experiments suggest that the rapid pressure-scalar term should be parameterized as above but with
It is not clear how the Stokes drift could affect the pressure term in the skewness (
The equations governing the turbulence statistics of a flow can be found by Reynolds averaging to find equations for turbulent statistics. The inclusion of a Stokes drift, both explicitly and in pressure-related closure assumptions, leads to the following additional terms in the budgets for turbulent fluxes and variances (note that
$ \frac{D_L \overline{w'\phi'}}{D t} =\ldots - \left[\mathbf{1}-\left(\frac{\alpha_{\phi 1}^{st} - \alpha_{\phi 2}^{st}}{2}\right)\right]\left[\overline{u'\phi'}\frac{\partial u^{st}}{\partial z} + \overline{v'\phi'}\frac{\partial v^{st}}{\partial z}\right]$
A declarative, efficient, and flexible JavaScript library for building user interfaces.
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. ๐๐๐
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google โค๏ธ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.