from the spectral stokes solvers to the stokes eigenmodes
play

From the Spectral Stokes solvers to the Stokes eigenmodes in - PowerPoint PPT Presentation

Introduction Continuous and time-discretized Stokes Problem Stokes Solvers Families and Properties Stokes Eigenmodes in the Square and Cube, from PrDi Solver Vorticity/Vector Potential correlations for the Stokes Flows Navier-Stokes Potential


  1. Introduction Continuous and time-discretized Stokes Problem Stokes Solvers Families and Properties Stokes Eigenmodes in the Square and Cube, from PrDi Solver Vorticity/Vector Potential correlations for the Stokes Flows Navier-Stokes Potential Formulation, and Questions From the Spectral Stokes solvers · · · to the Stokes eigenmodes in square/cube, · · · until new questions in Fluid Dynamics. G´ erard LABROSSE, Universit´ e Paris-Sud 11 e de Saint-´ Emmanuel LERICHE, Universit´ Etienne December 17, 2007 e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  2. Introduction Continuous and time-discretized Stokes Problem Stokes Solvers Families and Properties Stokes Eigenmodes in the Square and Cube, from PrDi Solver Vorticity/Vector Potential correlations for the Stokes Flows Navier-Stokes Potential Formulation, and Questions Plan Introduction 1 Continuous and time-discretized Stokes Problem 2 Stokes Solvers Families and Properties 3 UZAWA GREEN, or Influence Matrix Time-Splitting Projection-Diffusion (PrDi) Stokes Eigenmodes in the Square and Cube, from PrDi Solver 4 Vorticity/Vector Potential correlations for the Stokes Flows 5 Navier-Stokes Potential Formulation, and Questions 6 e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  3. Introduction Continuous and time-discretized Stokes Problem Stokes Solvers Families and Properties Stokes Eigenmodes in the Square and Cube, from PrDi Solver Vorticity/Vector Potential correlations for the Stokes Flows Navier-Stokes Potential Formulation, and Questions Introduction ¿ Why to pay a particular attention to the Unsteady Stokes Problem (USP) ? NS spectral numerical solutions are in fact USP solutions even for DNS of turbulence (considered as reliable) necessity of consistent and cheap USP pseudo-spectral solvers e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  4. Introduction Continuous and time-discretized Stokes Problem Stokes Solvers Families and Properties Stokes Eigenmodes in the Square and Cube, from PrDi Solver Vorticity/Vector Potential correlations for the Stokes Flows Navier-Stokes Potential Formulation, and Questions Continuous Unsteady Stokes Problem Let ( � v , p ) be solutions of ∂� v ∂ t − � ∇ 2 � v + � ∇ p = � f , in Ω × t > 0 , in ¯ � ∇ · � v = 0 Ω = Ω ∪ ∂ Ω , , V ( or ∂� v v = � � ∂ n = · · · ) on ∂ Ω , , equivalent to (1) � ∇ 2 p = � ∇ · � f , (1) � ∂ � ∂ t − � ∇ 2 ∇ 2 � � v = � ∇ × � ∇ × � f , � (2) ∇ · � v = 0 . (2) e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  5. Introduction Continuous and time-discretized Stokes Problem Stokes Solvers Families and Properties Stokes Eigenmodes in the Square and Cube, from PrDi Solver Vorticity/Vector Potential correlations for the Stokes Flows Navier-Stokes Potential Formulation, and Questions time-discretized Stokes Problem (1) Let us define • ( n ) ≡ • ( t = n δ t ). The USP high-order ( J i ) in time discretized formulation (KIO, JCP 1991) is v ( n +1) − � J i − 1 v ( n − q ) γ 0 � q =0 α q � v ( n +1) + � f ( n +1) , − � ∇ 2 � ∇ p = � in Ω , δ t v ( n +1) = 0 in ¯ � ∇ · � , Ω = Ω ∪ ∂ Ω , v ( n +1) = � V ( n +1) � , on ∂ Ω , γ 0 , α q given in Table I, p. 1390, of [E.Leriche and G.Labrosse, ”High-order direct Stokes solvers with or without temporal splitting : numerical investigations of their comparative properties”. SIAM J. Scient. Computing , 22(4) (2000), pp. 1386-1410]. e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  6. Introduction Continuous and time-discretized Stokes Problem Stokes Solvers Families and Properties Stokes Eigenmodes in the Square and Cube, from PrDi Solver Vorticity/Vector Potential correlations for the Stokes Flows Navier-Stokes Potential Formulation, and Questions time-discretized Stokes Problem (2) � �� � � J i − 1 � �� � v ( n − q ) q =0 α q � � γ 0 ∇ 2 � v ( n +1) + � f ( n +1) + δ t − � � � ∇ p = , in Ω δ t ⇓ v + � ∇ p = � H � S , in Ω , (3) � in ¯ ∇ · � v = 0 , Ω = Ω ∪ ∂ Ω , (4) ( or ∂� v v = � � V ∂ n = · · · ) on ∂ Ω . (5) , v from p , and enforce (more or less) � ¿ How to uncouple � ∇ · � v = 0 ? e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  7. Introduction Continuous and time-discretized Stokes Problem UZAWA Stokes Solvers Families and Properties GREEN, or Influence Matrix Stokes Eigenmodes in the Square and Cube, from PrDi Solver Time-Splitting Vorticity/Vector Potential correlations for the Stokes Flows Projection-Diffusion (PrDi) Navier-Stokes Potential Formulation, and Questions Stokes Solvers Families and Properties ¿ � ( � v , p ) Uncoupling Option Consistency Cost ∇ · � v = 0 ? UZAWA (’68) YES EXP. YES - GREEN or Influence Matrix YES EXP. YES (Kleiser, Schumann,’80) - Time Splitting NO CHEAP at spectral (Chorin, ’68, Temam, ’69) convergence - Projection-Diffusion (PrDi) YES CHEAP at spectral (Batoul et al., ’95) convergence CHEAP = POISSON + VECTORIAL HELMHOLTZ to be solved SAME SPACE-TIME ACCURACY e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  8. Introduction Continuous and time-discretized Stokes Problem UZAWA Stokes Solvers Families and Properties GREEN, or Influence Matrix Stokes Eigenmodes in the Square and Cube, from PrDi Solver Time-Splitting Vorticity/Vector Potential correlations for the Stokes Flows Projection-Diffusion (PrDi) Navier-Stokes Potential Formulation, and Questions UZAWA (1) The system (3-5) is space discretized, � � � Int = � ( H � v ) Int + D p S Int , � D · � v = 0 , v = � ( or ( � � V D · � n ) � v = · · · ) at the boundary nodes . with henceforth � v = ( u , v , w ) and p are column vectors of (velocity, pressure) nodal values, • Int , column vectors of internal nodal values of • , H , Helmholtz discrete operator, � D , � ∇ discrete operator. e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  9. Introduction Continuous and time-discretized Stokes Problem UZAWA Stokes Solvers Families and Properties GREEN, or Influence Matrix Stokes Eigenmodes in the Square and Cube, from PrDi Solver Time-Splitting Vorticity/Vector Potential correlations for the Stokes Flows Projection-Diffusion (PrDi) Navier-Stokes Potential Formulation, and Questions UZAWA (2) Eliminating the boundary nodal � v values through the BC, the USP reads   H u u Int  + ” � D ” p = ” � S ” H v v Int Int , (6)  H w w Int � D · � v = 0 , (7) ( H u , H v , H w ) ← H , square matrix with the BC on � v plugged in, D ” ← � ” � D , rectangular matrix. ( H u , H v , H w ) are invertible, then ... e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

  10. Introduction Continuous and time-discretized Stokes Problem UZAWA Stokes Solvers Families and Properties GREEN, or Influence Matrix Stokes Eigenmodes in the Square and Cube, from PrDi Solver Time-Splitting Vorticity/Vector Potential correlations for the Stokes Flows Projection-Diffusion (PrDi) Navier-Stokes Potential Formulation, and Questions UZAWA (3) � � H − 1 · D ” p + ” � v Int = − � ” � S ” � , Int and completing � v Int with the � v boundary values leads to � v written in term of p , and, by (7), one gets a pressure equation to be solved ... Example with � v | ∂ Ω = 0 : � v = ,, � D · � D ,, · � v Int , � D ” � H − 1 · ” � H − 1 · ” � ,, � D ,, · � p = − ,, � D ,, · � S ” ⇒ Int . � �� � UZAWA operator, full 2D/3D matrix, with a kernel (spurious pressure modes), only iteratively solved ... e de Saint-´ G´ erard LABROSSE, Universit´ e Paris-Sud 11 Emmanuel LERICHE, Universit´ From the Spectral Stokes solvers · · · to the Stokes eigenmodes Etienne

Recommend


More recommend