The '''variational multiscale method (VMS)''' is a technique used for deriving models and numerical methods for multiscale phenomena.<ref name="Hughes_Scovazzi_Franca_2004">{{cite book |last1=Hughes |first1=T.J.R. |last2=Scovazzi |first2=G. |last3=Franca |first3=L.P. |editor-last1=Stein |editor-first1=Erwin |editor-last2=de Borst |editor-first2=René |editor-last3=Hughes |editor-first3=Thomas J.R. |title=Encyclopedia of Computational Mechanics |publisher=John Wiley & Sons |date=2004 |pages=5–59 |chapter=Chapter 2: Multiscale and Stabilized Methods |isbn=0-470-84699-2}}</ref> The VMS framework has been mainly applied to design stabilized finite element methods in which stability of the standard Galerkin method is not ensured both in terms of singular perturbation and of compatibility conditions with the finite element spaces.<ref name="Codina_Badia_Baiges_Principe_2004">{{cite book |last1=Codina |first1=R. |last2=Badia |first2=S. |last3=Baiges |first3=J. |last4=Principe |first4=J. |editor-last1=Stein |editor-first1=Erwin |editor-last2=de Borst |editor-first2=René |editor-last3=Hughes |editor-first3=Thomas J.R. |title=Encyclopedia of Computational Mechanics Second Edition|publisher=John Wiley & Sons |date=2017 |pages=1–28 |chapter=Chapter 2: Variational Multiscale Methods in Computational Fluid Dynamics| isbn=9781119003793}}</ref>
Stabilized methods are getting increasing attention in computational fluid dynamics because they are designed to solve drawbacks typical of the standard Galerkin method: advection-dominated flows problems and problems in which an arbitrary combination of interpolation functions may yield to unstable discretized formulations.<ref name="Masud_2004">{{cite journal |last1=Masud |first1=Arif |title=Preface |journal=Computer Methods in Applied Mechanics and Engineering |date=April 2004 |volume=193 |issue=15–16 |pages=iii–iv |doi=10.1016/j.cma.2004.01.003}}</ref><ref name= "Quarteroni"/> The milestone of stabilized methods for this class of problems can be considered the Streamline Upwind Petrov-Galerkin method (SUPG), designed during 80s for convection dominated-flows for the incompressible Navier–Stokes equations by Brooks and Hughes.<ref>{{cite journal |last1=Brooks |first1=Alexander N. |last2=Hughes |first2=Thomas J.R. |title=Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier–Stokes equations |journal=Computer Methods in Applied Mechanics and Engineering |date=September 1982 |volume=32 |issue=1–3 |pages=199–259 |doi=10.1016/0045-7825(82)90071-8|bibcode=1982CMAME..32..199B }}</ref><ref>{{cite journal |last1=Masud |first1=Arif |last2=Calderer |first2=Ramon |title=A variational multiscale stabilized formulation for the incompressible Navier–Stokes equations |journal=Computational Mechanics |date=3 February 2009 |volume=44 |issue=2 |pages=145–160 |doi=10.1007/s00466-008-0362-3|bibcode=2009CompM..44..145M |s2cid=7036642 }}</ref> Variational Multiscale Method (VMS) was introduced by Hughes in 1995.<ref name="Hughes_2005"/> Broadly speaking, VMS is a technique used to get mathematical models and numerical methods which are able to catch multiscale phenomena;<ref name= "Hughes_Scovazzi_Franca_2004" /> in fact, it is usually adopted for problems with huge scale ranges, which are separated into a number of scale groups.<ref name="Rasthofer_Gravemeier_2017">{{cite journal |last1=Rasthofer |first1=Ursula |last2=Gravemeier |first2=Volker |title=Recent Developments in Variational Multiscale Methods for Large-Eddy Simulation of Turbulent Flow |journal=Archives of Computational Methods in Engineering |date=27 February 2017 |volume=25 |issue=3 |pages=647–690 |doi=10.1007/s11831-017-9209-4|s2cid=29169067 |hdl=20.500.11850/129122 |hdl-access=free }}</ref> The main idea of the method is to design a sum decomposition of the solution as <math> u = \bar u + u' </math>, where <math> \bar u </math> is denoted as coarse-scale solution and it is solved numerically, whereas <math> u' </math> represents the fine scale solution and is determined analytically eliminating it from the problem of the coarse scale equation.<ref name= "Hughes_Scovazzi_Franca_2004" />
== The abstract framework ==
===Abstract Dirichlet problem with variational formulation===
Consider an open bounded domain <math>\Omega \subset \mathbb R^d </math> with smooth boundary <math>\Gamma \subset \mathbb R^{d-1} </math>, being <math> d \geq 1 </math> the number of space dimensions. Denoting with <math> \mathcal L </math> a generic, second order, nonsymmetric differential operator, consider the following boundary value problem:<ref name="Quarteroni" /> :<math> \text{find } u: \Omega \to \mathbb R \text{ such that}: </math> :<math> \begin{cases} \mathcal L u = f & \text{ in } \Omega \\ u = g & \text{ on } \Gamma\\ \end{cases} </math>
being <math>f: \Omega \to \mathbb R </math> and <math>g: \Gamma \to \mathbb R </math> given functions. Let <math> H^1 (\Omega)</math> be the Hilbert space of square-integrable functions with square-integrable derivatives:<ref name="Quarteroni"/> :<math> H^1(\Omega)= \{ f \in L^2(\Omega): \nabla f \in L^2 (\Omega)\}. </math> Consider the trial solution space <math> \mathcal V_g </math> and the weighting function space <math> \mathcal V </math> defined as follows:<ref name="Quarteroni"/>
:<math> \mathcal V_g = \{u \in H^1 (\Omega): \, u=g \text{ on } \Gamma \}, </math>
:<math> \mathcal V = H_0^1(\Omega) = \{v \in H^1 (\Omega): \, v=0 \text{ on } \Gamma \}. </math>
The variational formulation of the boundary value problem defined above reads:<ref name="Quarteroni"/>
:<math> \text{find } u \in \mathcal V_g \text{ such that: } a(v, u) = f(v) \, \, \, \, \forall v \in \mathcal V </math>,
being <math> a(v, u) </math> the bilinear form satisfying <math> a(v, u) = (v, \mathcal L u) </math>, <math> f(v)=(v, f) </math> a bounded linear functional on <math> \mathcal V </math> and <math> (\cdot, \cdot) </math> is the <math> L^2(\Omega) </math> inner product.<ref name="Codina_Badia_Baiges_Principe_2004"/> Furthermore, the dual operator <math> \mathcal L^* </math> of <math> \mathcal L </math> is defined as that differential operator such that <math> \mathcal (v, \mathcal L u) = (\mathcal L^*v, u)\, \, \, \forall u, \, v \in \mathcal V </math>.<ref name="Hughes_2005"/>
===Variational multiscale method=== thumb|upright=2.0| One dimensional representation of <math>u</math>, <math>\bar u</math> and <math>u'</math> In VMS approach, the function spaces are decomposed through a multiscale direct sum decomposition for both <math>\mathcal V_g</math> and <math>\mathcal V </math> into coarse and fine scales subspaces as:<ref name = "Hughes_Scovazzi_Franca_2004" /> :<math> \mathcal V= \bar {\mathcal V} \oplus \mathcal V' </math> and :<math> \mathcal V_g= \bar {\mathcal V_g} \oplus \mathcal V_g'. </math>
Hence, an ''overlapping'' sum decomposition is assumed for both <math> u </math> and <math> v </math> as:
:<math> u = \bar{u} + u' \text{ and } v = \bar{v} + v'</math>,
where <math> \bar u </math> represents the ''coarse'' (resolvable) scales and <math> u' </math> the ''fine'' (subgrid) scales, with <math> \bar{u} \in \bar{\mathcal V_g} </math>, <math> {u'} \in {\mathcal V_g}' </math>, <math> \bar{v} \in \bar{\mathcal V} </math> and <math> v' \in {\mathcal V}' </math>. In particular, the following assumptions are made on these functions:<ref name = "Hughes_Scovazzi_Franca_2004" />
:<math> \begin{align} \bar u = g & & \text{ on } \Gamma & & \forall& \bar u \in \bar{\mathcal{V_g}}, \\ u' = 0 & & \text{ on } \Gamma & & \forall& u' \in {\mathcal{V_g}}', \\ \bar v = 0 & & \text{ on } \Gamma & & \forall& \bar v \in \bar{\mathcal{V}}, \\ v' = 0 & & \text{ on } \Gamma & & \forall& v' \in {\mathcal{V}}'. \end{align} </math>
With this in mind, the variational form can be rewritten as
:<math> a(\bar v + v', \bar u + u') = f (\bar v + v') </math>
and, by using bilinearity of <math> a (\cdot , \cdot) </math> and linearity of <math> f (\cdot) </math>,
:<math> a (\bar v, \bar u) + a (\bar v, u') + a (v', \bar u) + a(v', u') = f(\bar v) + f (v'). </math>
Last equation, yields to a coarse scale and a fine scale problem:
:<math> \text{find } \bar u \in \bar{\mathcal V}_g \text{ and } u' \in \mathcal V' \text{ such that: } </math> :<math> \begin{align} & & a (\bar v, \bar u) + a (\bar v, u') &= f(\bar v) & & \forall \bar v \in \bar{\mathcal V} & \, \, \, \, \text{coarse-scale problem}\\ & & a (v', \bar u) + a (v', u') &= f( v') & & \forall v' \in {\mathcal V}' & \, \, \, \, \text{fine-scale problem} \\ \end{align} </math> or, equivalently, considering that <math> a(v, u) = (v, \mathcal L u) </math> and <math> f(v)=(v, f) </math>: :<math> \text{find } \bar u \in \bar{\mathcal V}_g \text{ and } u' \in \mathcal V' \text{ such that: } </math> :<math> \begin{align} & & (\bar v, \mathcal L \bar u) + (\bar v, \mathcal L u') &= (\bar v, f) & & \forall \bar v \in \bar{\mathcal V}, \\ & & (v', \mathcal L \bar u) + (v', \mathcal L u') &= ( v', f) & & \forall v' \in {\mathcal V}'.\\ \end{align} </math> By rearranging the second problem as <math> (v', \mathcal L u') = - (v', \mathcal L \bar u - f) </math>, the corresponding Euler–Lagrange equation reads:<ref name="Hughes_2005"/> :<math> \begin{cases} \mathcal L u' = - (\mathcal L \bar u - f) & \text{ in } \Omega \\ u' = 0 & \text{ on } \Gamma \end{cases} </math> which shows that the fine scale solution <math> u' </math> depends on the strong residual of the coarse scale equation <math> \mathcal L \bar u - f </math>.<ref name="Hughes_2005">{{cite journal |last1=Hughes |first1=Thomas J.R. |title=Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods |journal=Computer Methods in Applied Mechanics and Engineering |date=November 1995 |volume=127 |issue=1–4 |pages=387–401 |doi=10.1016/0045-7825(95)00844-9|bibcode=1995CMAME.127..387H |doi-access=free }}</ref> The fine scale solution can be expressed in terms of <math> \mathcal L \bar u - f </math> through the Green's function <math>G: \Omega \times \Omega \to \mathbb R \text{ with } G=0 \text{ on } \Gamma \times \Gamma </math>: :<math> u'(y) = - \int_\Omega G(x, y) (\mathcal L \bar u - f )(x)\,d \Omega_x \, \, \, \forall y \in \Omega. </math> Let <math> \delta</math> be the Dirac delta function, by definition, the Green's function is found by solving <math> \forall y \in \Omega </math> :<math> \begin{cases} \mathcal L^* G(x, y) = \delta (x-y) & \text{ in } \Omega \\ G(x,y)=0 & \text{ on } \Gamma \end{cases} </math> Moreover, it is possible to express <math> u' </math> in terms of a new differential operator <math> \mathcal M </math> that approximates the differential operator <math> - \mathcal L^{-1} </math> as <ref name = "Hughes_Scovazzi_Franca_2004" /> : <math> u' = \mathcal M (\mathcal L \bar u - f), </math> with <math> \mathcal M \approx - \mathcal L^{-1}</math>. In order to eliminate the explicit dependence in the coarse scale equation of the sub-grid scale terms, considering the definition of the dual operator, the last expression can be substituted in the second term of the coarse scale equation:<ref name = "Hughes_Scovazzi_Franca_2004" /> :<math> (\bar v, \mathcal L u') = (\mathcal L^* \bar v, u') = (\mathcal L^* \bar v, \mathcal M (\mathcal L \bar u - f)). </math> Since <math> \mathcal M </math> is an approximation of <math> - \mathcal L^{-1} </math>, the Variational Multiscale Formulation will consist in finding an approximate solution <math> \tilde{\bar u} \approx \bar u </math> instead of <math> \bar u </math>. The coarse problem is therefore rewritten as:<ref name = "Hughes_Scovazzi_Franca_2004" /> :<math> \text{find } \tilde{\bar u} \in \mathcal{\bar V}_g: \; \; \; a (\bar v, \tilde{\bar u}) + (\mathcal L^* \bar v, \mathcal M (\mathcal L \tilde{\bar u} - f)) = (\bar v, f) \; \; \; \forall \bar{v} \in \mathcal {\bar V}, </math> being :<math> (\mathcal L^* \bar v, \mathcal M (\mathcal L \tilde{\bar u} - f)) = - \int_{\Omega} \int_{\Omega} (\mathcal L^* \bar v)(y) G(x, y) (\mathcal L \tilde{\bar u} - f)(x) \,d \Omega_x \,d\Omega_y. </math> Introducing the form <ref name="Hughes_2005"/> :<math> B(\bar v, \tilde{\bar u}, G) = a (\bar v, \tilde{\bar u}) + (\mathcal L^* \bar v, \mathcal M (\mathcal L \tilde{\bar u})) </math> and the functional :<math> L(\bar v, G)= (\bar v, f) + (\mathcal L^* \bar v, \mathcal M f) </math>, the VMS formulation of the coarse scale equation is rearranged as:<ref name="Hughes_2005"/> :<math> \text{find } \tilde{\bar u} \in \mathcal{\bar V}_g: \, B(\bar v, \tilde{\bar u}, G) = L(\bar v, G) \, \, \, \forall \bar{v} \in \mathcal {\bar V}. </math> Since commonly it is not possible to determine both <math> \mathcal M </math> and <math> G </math>, one usually adopt an approximation. In this sense, the coarse scale spaces <math>\bar{\mathcal V}_g</math> and <math>\bar{\mathcal V}</math> are chosen as finite dimensional space of functions as:<ref name = "Hughes_Scovazzi_Franca_2004" /> :<math> \bar{\mathcal V}_g \equiv \mathcal V_{g_h} : = \mathcal V_g \cap X_r^h(\Omega) </math> and :<math> \bar{\mathcal V} \equiv \mathcal V_{h} : = \mathcal V \cap X_h^r(\Omega), </math> being <math>X_r^h(\Omega)</math> the Finite Element space of Lagrangian polynomials of degree <math> r \geq 1 </math> over the mesh built in <math> \Omega </math> .<ref name="Quarteroni">{{cite book |author-link1= Alfio Quarteroni |last1=Quarteroni |first1=Alfio |title=Numerical models for differential problems |publisher=Springer |isbn=978-3-319-49316-9 |edition=Third |date=2017-10-10 }}</ref> Note that <math> \mathcal{V}_g' </math> and <math> \mathcal{V}' </math> are infinite-dimensional spaces, while <math> \mathcal{V}_{g_h} </math> and <math> \mathcal{V}_h </math> are finite-dimensional spaces.
Let <math> u_h \in \mathcal V_{g_h} </math> and <math> v_h \in \mathcal V_{h} </math> be respectively approximations of <math> \tilde{\bar u} </math> and <math> {\bar v} </math>, and let <math> \tilde G </math> and <math> \tilde{\mathcal M}</math> be respectively approximations of <math> G </math> and <math> {\mathcal M}</math>. The VMS problem with Finite Element approximation reads:<ref name="Hughes_2005"/> :<math> \text{find } u_h \in \mathcal V_{g_h}: B(v_h, u_h, \tilde G) = L( v_h, \tilde G) \, \, \, \forall {v}_h \in \mathcal { V}_h </math> or, equivalently: :<math> \text{find } u_h \in \mathcal V_{g_h}: a (v_h, u_h) + (\mathcal L^* v_h, \mathcal {\tilde{M}} (\mathcal L { u_h} - f)) = ( v_h, f) \, \, \, \forall {v}_h \in \mathcal { V}_h </math>
== VMS and stabilized methods ==
Consider an advection–diffusion problem:<ref name="Quarteroni" /> :<math> \begin{cases} -\mu \Delta u + \boldsymbol b \cdot \nabla u = f & \text{ in } \Omega \\ u=0 & \text{ on } \partial \Omega \end{cases} </math> where <math> \mu \in \mathbb R </math> is the diffusion coefficient with <math> \mu>0 </math> and <math> \boldsymbol b \in \mathbb R^d </math> is a given advection field. Let <math> \mathcal{V}= H^1_0(\Omega) </math> and <math> u \in \mathcal V </math>, <math> \boldsymbol b \in [L^2(\Omega)]^d </math>, <math> f \in L^2(\Omega) </math>.<ref name="Quarteroni"/> Let <math> \mathcal L = \mathcal L_{diff} + \mathcal L_{adv} </math>, being <math> \mathcal L_{diff} = - \mu \Delta </math> and <math> \mathcal L_{adv} = \boldsymbol b \cdot \nabla </math>.<ref name="Hughes_Scovazzi_Franca_2004"/> The variational form of the problem above reads:<ref name= "Quarteroni"/> :<math> \text{find} \, u \in \mathcal V: \; \; \; a(v, u) = (f, v) \; \; \; \forall v \in \mathcal V, </math> being :<math> a(v, u) = (\nabla v, \mu \nabla u) + (v, \boldsymbol b \cdot \nabla u). </math> Consider a Finite Element approximation in space of the problem above by introducing the space <math> \mathcal V_h = \mathcal V \cap X_h^r </math> over a grid <math> \Omega_h = \bigcup_{k=1}^{N} \Omega_k </math> made of <math> N </math> elements, with <math> u_h \in \mathcal V_h </math>.
The standard Galerkin formulation of this problem reads<ref name= "Quarteroni"/> :<math> \text{find } u_h \in \mathcal V_h: \; \; \; a(v_h, u_h) = (f, v_h) \; \; \; \forall v \in \mathcal V, </math> Consider a strongly consistent stabilization method of the problem above in a finite element framework: :<math> \text{ find } u_h \in \mathcal V_h: \, \, \, a(v_h, u_h) + \mathcal L_h (u_h, f; v_h)= (f, v_h) \, \, \, \forall v_h \in \mathcal V_h </math> for a suitable form <math> \mathcal L_h </math> that satisfies:<ref name="Quarteroni"/> :<math> \mathcal L_h (u, f; v_h) = 0 \, \, \, \forall v_h \in \mathcal V_h. </math> The form <math> \mathcal L_h </math> can be expressed as <math> (\mathbb L v_h, \tau (\mathcal L u_h - f))_{\Omega_h} </math>, being <math> \mathbb L </math> a differential operator such as:<ref name="Hughes_Scovazzi_Franca_2004"/> :<math> \mathbb L= \begin{cases} + \mathcal L & \, \, \, & \text{ Galerkin/least squares (GLS)} \\ + \mathcal L_{adv} & \, \, \, & \text{ Streamline Upwind Petrov-Galerkin (SUPG)} \\ - \mathcal L^* & \, \, \, & \text{ Multiscale} \\ \end{cases} </math> and <math> \tau </math> is the stabilization parameter. A stabilized method with <math> \mathbb L = -\mathcal L^* </math> is typically referred to '''multiscale stabilized method '''. In 1995, Thomas J.R. Hughes showed that a stabilized method of multiscale type can be viewed as a sub-grid scale model where the stabilization parameter is equal to :<math> \tau= - \tilde{\mathcal M} \approx - \mathcal M </math> or, in terms of the Green's function as :<math> \tau \delta (x-y) = \tilde G(x, y) \approx G(x,y), </math> which yields the following definition of <math> \tau </math>: :<math> \tau = \frac{1}{|\Omega_k|} \int_{\Omega_k} \int_{\Omega_k} G(x, y) \,d\Omega_x \,d\Omega_y. </math><ref name="Hughes_2005"/>
=== Stabilization Parameter Properties === For the 1-d advection diffusion problem, with an appropriate choice of basis functions and <math>\tau</math>, VMS provides a projection in the approximation space.<ref>{{cite journal |last1= Hughes|first1= T.J.|last2= Sangalli|first2= G.|title= Variational multiscale analysis: the fine-scale Green's function, projection, optimization, localization, and stabilized methods|journal= SIAM Journal on Numerical Analysis|date= 2007|volume= 45|issue= 2|pages= 539–557|publisher= SIAM |doi= 10.1137/050645646}}</ref> Further, an adjoint-based expression for <math>\tau</math> can be derived,<ref name="Garg_Stogner_CMAME_2019">{{cite journal |last1= Garg|first1= V.V.|last2= Stogner |first2= R.|title= Local enhancement of functional evaluation and adjoint error estimation for variational multiscale formulations|journal= Computer Methods in Applied Mechanics and Engineering|date= 2019|volume= 354|pages= 119–142| publisher= Elsevier|doi= 10.1016/j.cma.2019.05.023}}</ref> :<math> \tau_e = -\frac{\mathcal L(\tilde{z}, u_h)_e}{(\phi_{e}L_h(u_h)),L^*(\tilde{z}))_e} </math> where <math>\tau_e</math> is the element wise stabilization parameter, <math>\mathcal L(\tilde{z}, u_h)_e</math> is the element wise residual and the adjoint <math>\tilde{z}</math> problem solves, :<math> \mathcal a(\tilde{z}, v) + L_h(\tilde{z}, v) = \int_{\Omega_e} v \, dx </math> In fact, one can show that the <math>\tau</math> thus calculated allows one to compute the linear functional <math> \int_{\Omega} u \, dx </math> exactly.<ref name="Garg_Stogner_CMAME_2019" />
== VMS turbulence modeling for large-eddy simulations of incompressible flows ==
The idea of VMS turbulence modeling for Large Eddy Simulations(LES) of incompressible Navier–Stokes equations was introduced by Hughes et al. in 2000 and the main idea was to use - instead of classical filtered techniques - variational projections.<ref name="Hughes_Mazzei_Jansen_2000">{{cite journal |last1=Hughes |first1=Thomas J.R. |last2=Mazzei |first2=Luca |last3=Jansen |first3=Kenneth E. |title=Large Eddy Simulation and the variational multiscale method |journal=Computing and Visualization in Science |date=May 2000 |volume=3 |issue=1–2 |pages=47–59 |doi=10.1007/s007910050051|s2cid=120207183 }}</ref><ref name="Bazilevs_Calo_Cottrell_Hughes_Scovazzi_2007">{{cite journal |last1=Bazilevs |first1=Y. |last2=Calo |first2=V.M. |last3=Cottrell |first3=J.A. |last4=Hughes |first4=T.J.R. |last5=Reali |first5=A. |last6=Scovazzi |first6=G. |title=Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows |journal=Computer Methods in Applied Mechanics and Engineering |date=December 2007 |volume=197 |issue=1–4 |pages=173–201 |doi=10.1016/j.cma.2007.07.016|bibcode=2007CMAME.197..173B }}</ref>
=== Incompressible Navier–Stokes equations === Consider the incompressible Navier–Stokes equations for a Newtonian fluid of constant density <math> \rho </math> in a domain <math> \Omega \in \mathbb R^d </math> with boundary <math> \partial \Omega = \Gamma_D \cup \Gamma_N </math>, being <math> \Gamma_D </math> and <math> \Gamma_N </math> portions of the boundary where respectively a Dirichlet and a Neumann boundary condition is applied (<math> \Gamma_D \cap \Gamma_N = \emptyset </math>):<ref name = "Quarteroni" /> :<math> \begin{cases} \rho \dfrac{\partial \boldsymbol u}{\partial t} + \rho (\boldsymbol u \cdot \nabla) \boldsymbol u - \nabla \cdot \boldsymbol \sigma (\boldsymbol u, p) = \boldsymbol f & \text{ in } \Omega \times (0, T) \\ \nabla \cdot \boldsymbol u = 0 & \text{ in } \Omega \times (0, T) \\ \boldsymbol u = \boldsymbol g & \text{ on } \Gamma_D \times (0, T) \\ \sigma (\boldsymbol u, p) \boldsymbol{\hat n} = \boldsymbol h & \text{ on } \Gamma_N \times (0, T) \\ \boldsymbol u(0)= \boldsymbol u_0 & \text{ in } \Omega \times \{ 0\} \end{cases} </math> being <math> \boldsymbol u </math> the fluid velocity, <math> p </math> the fluid pressure, <math> \boldsymbol f </math> a given forcing term, <math> \boldsymbol{\hat n} </math> the outward directed unit normal vector to <math> \Gamma_N </math>, and <math> \boldsymbol \sigma (\boldsymbol u, p) </math> the viscous stress tensor defined as: :<math> \boldsymbol \sigma (\boldsymbol u, p) = -p \boldsymbol I + 2 \mu \boldsymbol \epsilon (\boldsymbol u). </math> Let <math> \mu </math> be the dynamic viscosity of the fluid, <math> \boldsymbol I </math> the second order identity tensor and <math> \boldsymbol \epsilon (\boldsymbol u) </math> the strain-rate tensor defined as: :<math> \boldsymbol \epsilon (\boldsymbol u) = \frac{1}{2} ((\nabla \boldsymbol u) + (\nabla \boldsymbol u)^T). </math> The functions <math> \boldsymbol g </math> and <math> \boldsymbol h </math> are given Dirichlet and Neumann boundary data, while <math> \boldsymbol u_0 </math> is the initial condition.<ref name= "Quarteroni"/>
=== Global space time variational formulation ===
In order to find a variational formulation of the Navier–Stokes equations, consider the following infinite-dimensional spaces:<ref name="Quarteroni"/> :<math> \mathcal V_g= \{ \boldsymbol u \in [H^1(\Omega)]^d : \boldsymbol u = \boldsymbol g \text{ on } \Gamma_D \}, </math> :<math> \mathcal V_0 = [H^1_0(\Omega)]^d=\{ \boldsymbol u \in [H^1(\Omega)]^d : \boldsymbol u = \boldsymbol 0 \text{ on } \Gamma_D \}, </math> :<math> \mathcal Q =L^2(\Omega). </math> Furthermore, let <math> \boldsymbol \mathcal V_g = \mathcal V_g \times \mathcal Q </math> and <math> \boldsymbol \mathcal V_0 = \mathcal V_0 \times \mathcal Q </math>. The weak form of the unsteady-incompressible Navier–Stokes equations reads:<ref name="Quarteroni"/> given <math> \boldsymbol u_{0}</math>, :<math> \forall t \in (0, T), \; \text{find } (\boldsymbol u, p) \in \boldsymbol \mathcal V_g \text{ such that } </math> :<math> \begin{align} \bigg( \boldsymbol v, \rho \dfrac{\partial \boldsymbol u}{\partial t}\bigg) + a (\boldsymbol v, \boldsymbol u) + c (\boldsymbol v, \boldsymbol u, \boldsymbol u) - b (\boldsymbol v, p) + b (\boldsymbol u, q) = (\boldsymbol v, \boldsymbol f) + (\boldsymbol v, \boldsymbol h)_{\Gamma_N} \; \; \forall (\boldsymbol v, q) \in \boldsymbol \mathcal V_0 \end{align} </math> where <math> (\cdot, \cdot) </math> represents the <math> L^2(\Omega) </math> inner product and <math> (\cdot, \cdot)_{\Gamma_N} </math> the <math> L^2(\Gamma_N) </math> inner product. Moreover, the bilinear forms <math> a(\cdot, \cdot) </math>, <math> b(\cdot, \cdot) </math> and the trilinear form <math> c(\cdot, \cdot, \cdot) </math> are defined as follows:<ref name="Quarteroni"/> :<math> \begin{align} a (\boldsymbol v, \boldsymbol u) = & (\nabla \boldsymbol v, \mu ((\nabla \boldsymbol u) + (\nabla \boldsymbol u)^T)), \\ b (\boldsymbol v, q) = &(\nabla \cdot \boldsymbol v, q), \\ c (\boldsymbol v, \boldsymbol u, \boldsymbol u) = &(\boldsymbol v, \rho (\boldsymbol u \cdot \nabla) \boldsymbol u). \end{align} </math>
===Finite element method for space discretization and VMS-LES modeling=== In order to discretize in space the Navier–Stokes equations, consider the function space of finite element :<math> X_r^h = \{u^h \in C^0 (\overline\Omega): u^h|_k \in \mathbb P_r, \; \forall k \in \Tau_h\} </math> of piecewise Lagrangian Polynomials of degree <math> r \geq 1 </math> over the domain <math> \Omega </math> triangulated with a mesh <math> \Tau_h</math> made of tetrahedrons of diameters <math> h_k </math>, <math> \forall k \in \Tau_h </math>. Following the approach shown above, let introduce a multiscale direct-sum decomposition of the space <math>\boldsymbol \mathcal V </math> which represents either <math>\boldsymbol \mathcal V_g</math> and <math>\boldsymbol \mathcal V_0</math>:<ref name="Forti_Dede_2015" >{{cite journal |last1=Forti |first1=Davide |last2=Dedè |first2=Luca |title=Semi-implicit BDF time discretization of the Navier–Stokes equations with VMS-LES modeling in a High Performance Computing framework |journal=Computers & Fluids |date=August 2015 |volume=117 |pages=168–182 |doi=10.1016/j.compfluid.2015.05.011}}</ref> :<math> \boldsymbol \mathcal V = \boldsymbol \mathcal V_h \oplus \boldsymbol \mathcal V', </math> being :<math> \boldsymbol \mathcal V_h = \mathcal V_{g_h} \times \mathcal Q \text{ or } \boldsymbol \mathcal V_h = \mathcal V_{0_h} \times \mathcal Q </math> the finite dimensional function space associated to the '''coarse scale''', and :<math> \boldsymbol \mathcal V' = \mathcal V_g' \times \mathcal Q \text{ or } \boldsymbol \mathcal V' = \mathcal V_0' \times \mathcal Q </math> the infinite-dimensional '''fine scale''' function space, with :<math> \mathcal V_{g_h} = \mathcal V_g \cap X_r^h </math>, :<math> \mathcal V_{0_h} = \mathcal V_0 \cap X_r^h </math> and :<math> \mathcal Q_h = \mathcal Q \cap X_r^h </math>. An overlapping sum decomposition is then defined as:<ref name="Bazilevs_Calo_Cottrell_Hughes_Scovazzi_2007"/><ref name="Forti_Dede_2015"/> :<math> \begin{align} & \boldsymbol u = \boldsymbol u^h + \boldsymbol u' \text{ and } p = p^h + p' \\ & \boldsymbol v = \boldsymbol v^h + \boldsymbol v' \;\text{ and } q = q^h + q' \end{align} </math> By using the decomposition above in the variational form of the Navier–Stokes equations, one gets a coarse and a fine scale equation; the fine scale terms appearing in the coarse scale equation are integrated by parts and the fine scale variables are modeled as:<ref name="Bazilevs_Calo_Cottrell_Hughes_Scovazzi_2007"/> :<math> \begin{align} \boldsymbol u' \approx & -\tau_M (\boldsymbol u^h) \boldsymbol r_M (\boldsymbol u^h, p^h), \\ p' \approx & -\tau_C (\boldsymbol u^h) \boldsymbol r_C (\boldsymbol u^h). \end{align} </math> In the expressions above, <math>\boldsymbol r_M (\boldsymbol u^h, p^h) </math> and <math> \boldsymbol r_C (\boldsymbol u^h) </math> are the residuals of the momentum equation and continuity equation in strong forms defined as: :<math> \begin{align} \boldsymbol r_M (\boldsymbol u^h, p^h) = & \rho \dfrac{\partial \boldsymbol u^h}{\partial t} + \rho (\boldsymbol u^h \cdot \nabla) \boldsymbol u^h - \nabla \cdot \boldsymbol \sigma (\boldsymbol u^h, p^h) - \boldsymbol f,\\ \boldsymbol r_C (\boldsymbol u^h) = & \nabla \cdot \boldsymbol u^h, \end{align} </math> while the stabilization parameters are set equal to:<ref name="Forti_Dede_2015"/> :<math> \begin{align} \tau_M (\boldsymbol u^h) = & \bigg ( \frac{\sigma^2 \rho^2}{\Delta t^2} + \frac{\rho^2 }{h_k^2} |\boldsymbol u^h|^2 + \frac{\mu^2}{h_k^4}C_r\bigg )^{-1/2}, \\ \tau_C (\boldsymbol u^h) = & \frac{h_k^2}{\tau_M (\boldsymbol u^h) }, \end{align} </math> where <math> C_r = 60 \cdot 2^{r-2}</math> is a constant depending on the polynomials's degree <math> r </math>, <math> \sigma </math> is a constant equal to the order of the backward differentiation formula (BDF) adopted as temporal integration scheme and <math> \Delta t </math> is the time step.<ref name="Forti_Dede_2015"/> The semi-discrete variational multiscale multiscale formulation (VMS-LES) of the incompressible Navier–Stokes equations, reads:<ref name="Forti_Dede_2015"/> given <math> \boldsymbol u_{0}</math>, :<math> \forall t \in (0, T), \; \text{find } \boldsymbol U^h = \{\boldsymbol u^h, p^h\} \in \boldsymbol \mathcal V_{g_h} \text{ such that } A(\boldsymbol V^h, \boldsymbol U^h ) = F(\boldsymbol V^h) \; \; \forall \boldsymbol V^h= \{\boldsymbol v^h, q^h\} \in \boldsymbol \mathcal V_{0_h}, </math> being :<math> A(\boldsymbol V^h, \boldsymbol U^h ) = A^{NS}(\boldsymbol V^h, \boldsymbol U^h ) + A^{VMS}(\boldsymbol V^h, \boldsymbol U^h ), </math> and :<math> F(\boldsymbol V^h) = (\boldsymbol v, \boldsymbol f) + (\boldsymbol v, \boldsymbol h)_{\Gamma_N}. </math> The forms <math> A^{NS}(\cdot, \cdot) </math> and <math> A^{VMS}(\cdot, \cdot) </math> are defined as:<ref name="Forti_Dede_2015"/> :<math> \begin{align} A^{NS}(\boldsymbol V^h, \boldsymbol U^h )= & \bigg( \boldsymbol v^h, \rho \dfrac{\partial \boldsymbol u^h}{\partial t}\bigg) + a (\boldsymbol v^h, \boldsymbol u^h) + c (\boldsymbol v^h, \boldsymbol u^h, \boldsymbol u^h) - b (\boldsymbol v^h, p^h) + b (\boldsymbol u^h, q^h), \\ A^{VMS}(\boldsymbol V^h, \boldsymbol U^h ) = & \underbrace{\big( \rho \boldsymbol u^h \cdot \nabla \boldsymbol v^h + \nabla q^h, \tau_M(\boldsymbol u^h) \boldsymbol r_M(\boldsymbol u^h, p^h) \big)}_{\text{SUPG}} - \underbrace{(\nabla \cdot \boldsymbol v^h, \tau_c(\boldsymbol u_h)\boldsymbol r_C(\boldsymbol u^h)) + \big( \rho \boldsymbol u^h \cdot (\nabla \boldsymbol u^h)^T, \tau_M(\boldsymbol u^h) \boldsymbol r_M (\boldsymbol u^h, p^h) \big )}_{\text{VMS}} - \underbrace{(\nabla \boldsymbol v^h, \tau_M(\boldsymbol u^h) \boldsymbol r_M(\boldsymbol u^h, p^h) \otimes \tau_M(\boldsymbol u^h) \boldsymbol r_M(\boldsymbol u^h, p^h) )}_{\text{LES}}. \end{align} </math> From the expressions above, one can see that:<ref name="Forti_Dede_2015"/> * the form <math>A^{NS}(\cdot, \cdot)</math> contains the standard terms of the Navier–Stokes equations in variational formulation; * the form <math>A^{VMS}(\cdot, \cdot)</math> contain four terms: # the first term is the classical SUPG stabilization term; # the second term represents a stabilization term additional to the SUPG one; # the third term is a stabilization term typical of the VMS modeling; # the fourth term is peculiar of the LES modeling, describing the Reynolds cross-stress.
==See also == * Navier–Stokes equations * Large eddy simulation * Finite element method * Backward differentiation formula * Computational fluid dynamics * Streamline upwind Petrov–Galerkin pressure-stabilizing Petrov–Galerkin formulation for incompressible Navier–Stokes equations
==References== <references />
Category:Mathematical modeling Category:Numerical analysis Category:Computational fluid dynamics