Changeset 11422 for NEMO/branches/2019/fix_vvl_ticket1791/doc/latex/NEMO/subfiles/chap_conservation.tex
- Timestamp:
- 2019-08-08T15:40:47+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/fix_vvl_ticket1791/doc
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/fix_vvl_ticket1791/doc
- Property svn:ignore deleted
-
NEMO/branches/2019/fix_vvl_ticket1791/doc/latex
- Property svn:ignore
-
old new 1 *.aux 2 *.bbl 3 *.blg 4 *.dvi 5 *.fdb* 6 *.fls 7 *.idx 8 *.ilg 9 *.ind 10 *.log 11 *.maf 12 *.mtc* 13 *.out 14 *.pdf 15 *.toc 16 _minted-* 1 figures
-
- Property svn:ignore
-
NEMO/branches/2019/fix_vvl_ticket1791/doc/latex/NEMO
- Property svn:ignore deleted
-
NEMO/branches/2019/fix_vvl_ticket1791/doc/latex/NEMO/subfiles/chap_conservation.tex
r10442 r11422 21 21 horizontal kinetic energy and/or potential enstrophy of horizontally non-divergent flow, 22 22 and variance of temperature and salinity will be conserved in the absence of dissipative effects and forcing. 23 \citet{ Arakawa1966} has first pointed out the advantage of this approach.23 \citet{arakawa_JCP66} has first pointed out the advantage of this approach. 24 24 He showed that if integral constraints on energy are maintained, 25 25 the computation will be free of the troublesome "non linear" instability originally pointed out by 26 \citet{ Phillips1959}.26 \citet{phillips_TAMS59}. 27 27 A consistent formulation of the energetic properties is also extremely important in carrying out 28 28 long-term numerical simulations for an oceanographic model. 29 Such a formulation avoids systematic errors that accumulate with time \citep{ Bryan1997}.29 Such a formulation avoids systematic errors that accumulate with time \citep{bryan_JCP97}. 30 30 31 31 The general philosophy of OPA which has led to the discrete formulation presented in {\S}II.2 and II.3 is to … … 39 39 Note that in some very specific cases as passive tracer studies, the positivity of the advective scheme is required. 40 40 In that case, and in that case only, the advective scheme used for passive tracer is a flux correction scheme 41 \citep{Marti1992 , Levy1996, Levy1998}.41 \citep{Marti1992?, Levy1996?, Levy1998?}. 42 42 43 43 % ------------------------------------------------------------------------------------------------------------- … … 65 65 % \label{eq:vor_vorticity} 66 66 \int_D {{\textbf {k}}\cdot \frac{1}{e_3 }\nabla \times \left( {\varsigma 67 \;{\ rm {\bf k}}\times {\textbf {U}}_h } \right)\;dv} =067 \;{\mathrm {\mathbf k}}\times {\textbf {U}}_h } \right)\;dv} =0 68 68 \] 69 69 … … 189 189 \[ 190 190 % \label{eq:dynldf_dyn} 191 \int\limits_D {\frac{1}{e_3 }{\ rm {\bf k}}\cdot \nabla \times \left[ {\nabla191 \int\limits_D {\frac{1}{e_3 }{\mathrm {\mathbf k}}\cdot \nabla \times \left[ {\nabla 192 192 _h \left( {A^{lm}\;\chi } \right)-\nabla _h \times \left( {A^{lm}\;\zeta 193 \;{\ rm {\bf k}}} \right)} \right]\;dv} =0193 \;{\mathrm {\mathbf k}}} \right)} \right]\;dv} =0 194 194 \] 195 195 … … 197 197 % \label{eq:dynldf_div} 198 198 \int\limits_D {\nabla _h \cdot \left[ {\nabla _h \left( {A^{lm}\;\chi } 199 \right)-\nabla _h \times \left( {A^{lm}\;\zeta \;{\ rm {\bf k}}} \right)}199 \right)-\nabla _h \times \left( {A^{lm}\;\zeta \;{\mathrm {\mathbf k}}} \right)} 200 200 \right]\;dv} =0 201 201 \] … … 203 203 \[ 204 204 % \label{eq:dynldf_curl} 205 \int_D {{\ rm {\bf U}}_h \cdot \left[ {\nabla _h \left( {A^{lm}\;\chi }206 \right)-\nabla _h \times \left( {A^{lm}\;\zeta \;{\ rm {\bf k}}} \right)}205 \int_D {{\mathrm {\mathbf U}}_h \cdot \left[ {\nabla _h \left( {A^{lm}\;\chi } 206 \right)-\nabla _h \times \left( {A^{lm}\;\zeta \;{\mathrm {\mathbf k}}} \right)} 207 207 \right]\;dv} \leqslant 0 208 208 \] … … 210 210 \[ 211 211 % \label{eq:dynldf_curl2} 212 \mbox{if}\quad A^{lm}=cste\quad \quad \int_D {\zeta \;{\ rm {\bf k}}\cdot212 \mbox{if}\quad A^{lm}=cste\quad \quad \int_D {\zeta \;{\mathrm {\mathbf k}}\cdot 213 213 \nabla \times \left[ {\nabla _h \left( {A^{lm}\;\chi } \right)-\nabla _h 214 \times \left( {A^{lm}\;\zeta \;{\ rm {\bf k}}} \right)} \right]\;dv}214 \times \left( {A^{lm}\;\zeta \;{\mathrm {\mathbf k}}} \right)} \right]\;dv} 215 215 \leqslant 0 216 216 \] … … 220 220 \mbox{if}\quad A^{lm}=cste\quad \quad \int_D {\chi \;\nabla _h \cdot \left[ 221 221 {\nabla _h \left( {A^{lm}\;\chi } \right)-\nabla _h \times \left( 222 {A^{lm}\;\zeta \;{\ rm {\bf k}}} \right)} \right]\;dv} \leqslant 0222 {A^{lm}\;\zeta \;{\mathrm {\mathbf k}}} \right)} \right]\;dv} \leqslant 0 223 223 \] 224 224 … … 260 260 % \label{eq:dynzdf_vor} 261 261 \begin{aligned} 262 & \int_D {\frac{1}{e_3 }{\ rm {\bf k}}\cdot \nabla \times \left( {\frac{1}{e_3263 }\frac{\partial }{\partial k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\ rm264 {\ bf U}}_h }{\partial k}} \right)} \right)\;dv} =0 \\265 & \int_D {\zeta \,{\ rm {\bf k}}\cdot \nabla \times \left( {\frac{1}{e_3266 }\frac{\partial }{\partial k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\ rm267 {\ bf U}}_h }{\partial k}} \right)} \right)\;dv} \leq 0 \\262 & \int_D {\frac{1}{e_3 }{\mathrm {\mathbf k}}\cdot \nabla \times \left( {\frac{1}{e_3 263 }\frac{\partial }{\partial k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\mathrm 264 {\mathbf U}}_h }{\partial k}} \right)} \right)\;dv} =0 \\ 265 & \int_D {\zeta \,{\mathrm {\mathbf k}}\cdot \nabla \times \left( {\frac{1}{e_3 266 }\frac{\partial }{\partial k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\mathrm 267 {\mathbf U}}_h }{\partial k}} \right)} \right)\;dv} \leq 0 \\ 268 268 \end{aligned} 269 269 \] … … 273 273 \begin{aligned} 274 274 &\int_D {\nabla \cdot \left( {\frac{1}{e_3 }\frac{\partial }{\partial 275 k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\ rm {\bf U}}_h }{\partial k}}275 k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\mathrm {\mathbf U}}_h }{\partial k}} 276 276 \right)} \right)\;dv} =0 \\ 277 277 & \int_D {\chi \;\nabla \cdot \left( {\frac{1}{e_3 }\frac{\partial }{\partial 278 k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\ rm {\bf U}}_h }{\partial k}}278 k}\left( {\frac{A^{vm}}{e_3 }\frac{\partial {\mathrm {\mathbf U}}_h }{\partial k}} 279 279 \right)} \right)\;dv} \leq 0 \\ 280 280 \end{aligned}
Note: See TracChangeset
for help on using the changeset viewer.