Vector and scalar penalty-projection methods
|
|
- Jemimah Nichols
- 6 years ago
- Views:
Transcription
1 Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille with M. Jobelin [PhD, 2006], J.-C. Latché J.-P. Caltagirone and P. Fabrie Porquerolles, April
2 Motivations and objectives Work focusing on the constraint of free divergence How to deal efficiently with the free-divergence constraint with fractional-step methods (prediction-correction steps)? How to circumvent the major drawbacks of the usual projection methods including a scalar correction step for the Lagrange multiplier with solution of a Poisson-type equation? Example of fluid-type models with pressure as Lagrange multiplier solution of unsteady incompressible Navier-Stokes equations in the primitive variables (velocity and pressure) Introduction 2
3 Outlines 1 Projection methods for incompressible flows 2 Scalar penalty-projection methods 3 Vector penalty-projection methods 4 Conclusion and perspectives Introduction 3
4 Outlines 1 Projection methods for incompressible flows Non-homogeneous incompressible flows Semi-implicit method (linearly implicit) Fractional-step and projection methods 2 Scalar penalty-projection methods 3 Vector penalty-projection methods 4 Conclusion and perspectives Projection methods for incompressible flows 4
5 Navier-Stokes problem for incompressible flows Incompressible and variable density flows of Newtonian fluids Navier-Stokes equations with mixed boundary conditions : Dirichlet on Ω D and open (outflow) B.C. on Ω N ϱ + u ϱ = 0 in Ω ]0, T [ t [ ] u ϱ + (u )u τ (u) + p = f in Ω ]0, T [ t u = 0 in Ω ]0, T [ u = u D on Ω D ]0, T [ pn + τ (u) n = f N on Ω N ]0, T [ ϱ = ϱ 0, and u = u 0 in Ω {0} τ (u) = [µ( u + u T )], or µ u (for a constant viscosity) For an homogeneous fluid with constant density, we set ϱ = 1. Projection methods for incompressible flows Non-homogeneous incompressible flows 5
6 Semi-implicit method (linearly implicit) Semi-discretization in time with low or high-order BDF schemes for all n IN such that (n + 1)δt T : D ϱ n+1 + (ū,n+1 ϱ n+1) = 0 in Ω δt ( ) Dū ϱ n+1 n+1 + (ū,n+1 )ū n+1 δt τ (ū n+1 ) + p n+1 = f n+1 in Ω ū n+1 = 0 ū n+1 = u n+1 D p n+1 n + τ (ū n+1 ) n = f n+1 N in Ω on Ω D on Ω N Resolution of an elliptic problem in space (at each time step) by methods of finite elements, finite volumes, DGM... Projection methods for incompressible flows Semi-implicit method (linearly implicit) 6
7 Semi-implicit method : coupled solver Algebraic formulation of the Navier-Stokes system inf-sup stable discretization in space or with stabilization... β q δt M ϱu n+1 + A(U n, U n 1 )U n+1 + B T P n+1 = F(U n, U n 1, f n+1, f n+1 N, un+1 D ) BU n+1 = G(u n+1 D ) Projection methods for incompressible flows Semi-implicit method (linearly implicit) 7
8 Semi-implicit method : coupled solver Algebraic formulation of the Navier-Stokes system inf-sup stable discretization in space or with stabilization... β q δt M ϱu n+1 + A(U n, U n 1 )U n+1 + B T P n+1 = F(U n, U n 1, f n+1, f n+1 N, un+1 D ) BU n+1 = G(u n+1 D ) Saddle-point problem at each time step : efficient solver? Fully-coupled solver with efficient multigrid preconditioner : Kortas, PhD Cihlář and Angot, 1999 Augmented Lagrangian method with iterative Uzawa algorithm : Fortin and Glowinski, Khadra et al., IJNMF 2000 Projection methods using the Helmoltz-Hodge-Leray decomposition of L 2 (Ω) d Projection methods for incompressible flows Semi-implicit method (linearly implicit) 7
9 Fractional-step and projection-type methods Projection scheme with scalar pressure-correction [Chorin, Temam, Goda, Van Kan, 1986] See recent review : [Guermond, Minev, Shen, CMAME 2006] Prediction ( step ) Dũ ϱ n+1 n+1 + (u,n+1 )ũ n+1 τ (ũ n+1 ) + p n = f n+1 δt ũ n+1 = u n+1 D on Ω D p n n + τ (ũ n+1 ) n = f n+1 N on Ω N Projection step β q ϱ n+1 un+1 ũ n+1 + φ = 0 δt φ n = 0 on Ω D (necessary since u n = ũ n) φ = 0 on Ω N (sufficient by ( orthogonal ) projection onto H) δt u n+1 = 0 φ = β ϱn+1 q ũ n+1 Pressure-correction step : φ pressure increment p n+1 = p n + φ Projection methods for incompressible flows Fractional-step and projection methods 8
10 Fractional-step and projection-type methods Algebraic formulation for the Navier-Stokes system β q δt MϱŨn+1 + AŨn+1 + B T P n = F L ϱ Φ = β ( ) q δt BŨn+1 G P n+1 = P n + Φ M ϱ U n+1 = M ϱ Ũ n+1 + δt β q B T Φ Projection methods for incompressible flows Fractional-step and projection methods 8
11 Fractional-step and projection-type methods Major drawbacks of the incremental projection methods Time order of the splitting error? i.e. error between the numerical solutions of the implicit (or semi-implicit) method and the fractional-step method φ n = 0 on Ω D existence of an artificial pressure boundary layer in space φ = 0 on Ω N convergence in time and space spoiled for outflow boundary conditions : splitting error varying like O(δt 1/2 ) (pressure) and no more negligible (for both velocity and pressure) with respect to the time and space discretization error Pressure-correction step strongly dependent on density and viscosity for non-homogeneous flows convergence very poor for large ratios of ϱ Projection methods for incompressible flows Fractional-step and projection methods 8
12 Numerical tests Green-Taylor vortices : Navier-Stokes with Dirichlet B.C. Velocity error (discrete L 2 (0, T ; L 2 (Ω) d ) norm) versus time step δt 1e-2 norme L2 erreur vitesse 1e-3 1e-4 1e-5 1e-6 semi implicite 1e-7 1e-4 1e-3 1e-2 1e-1 Pas de temps Time convergence in O(δt 2 ) Stagnation threshold = space discretization error in O(h 2 ) Projection methods for incompressible flows Fractional-step and projection methods 9
13 Numerical tests Green-Taylor vortices : Navier-Stokes with Dirichlet B.C. Velocity error (discrete L 2 (0, T ; L 2 (Ω) d ) norm) versus time step δt 1e-2 norme L2 erreur vitesse 1e-3 1e-4 1e-5 1e-6 semi implicite incrementale 1e-7 1e-4 1e-3 1e-2 1e-1 Pas de temps Time convergence in O(δt 2 ) Stagnation threshold = space discretization error in O(h 2 ) Projection methods for incompressible flows Fractional-step and projection methods 9
14 Outlines 1 Projection methods for incompressible flows 2 Scalar penalty-projection methods Penalty-projection methods Numerical experiments Analysis of the scalar penalty-projection method for the Stokes problem 3 Vector penalty-projection methods 4 Conclusion and perspectives Scalar penalty-projection methods 10
15 Penalty-projection methods [Jobelin et al., JCP 2006] : prediction step with augmented Lagrangian for r > 0 and consistent projection step by scalar pressure-correction [Shen, 1992] : r = 1/δt 2 with a different correction of pressure [Caltagirone and Breil, 1999] : r > 0 with a singular projection operator... Scalar penalty-projection methods Penalty-projection methods 11
16 Penalty-projection methods Penalty-prediction ( step : augmentation) parameter r 0 Dũ ϱ n+1 n+1 + (u,n+1 )ũ n+1 τ (ũ n+1 ) δt r ( ũ n+1) + p n = f n+1 ũ n+1 = u n+1 D on Ω D p n n + τ (ũ n+1 ) n = f n+1 N on Ω N Projection step β q ϱ n+1 un+1 ũ n+1 + φ = 0 δt φ n = 0 on Ω D (necessary since u n = ũ n) φ = 0 on Ω N (sufficient by( orthogonal ) projection onto H) δt u n+1 = 0 φ = β ϱn+1 q ũ n+1 Pressure-correction step : φ consistent pressure increment p n+1 = p n r ũ n+1 + φ = p n+1 + φ Scalar penalty-projection methods Penalty-projection methods 11
17 Penalty-projection methods Algebraic formulation for the Navier-Stokes system β ( ) q δt MϱŨn+1 + rb T M 1 pl BŨ n+1 G + B T P n +AŨ n+1 = F L ϱ Φ = β ( q δt ) BŨ n+1 G P n+1 = P n + rm 1 pl ( BŨ n+1 G M ϱ U n+1 = M ϱ Ũ n+1 + δt β q B T Φ ) + Φ Preconditioning the prediction step by one iteration of augmented Lagrangian and consistent scalar projection r = 0 : incremental projection method Scalar penalty-projection methods Penalty-projection methods 12
18 Penalty-projection methods Algebraic formulation for the Navier-Stokes system β ( ) q δt MϱŨn+1 + B T (rm 1 pl BŨ n+1 G + P n ) }{{} P n+1 +AŨ n+1 = F L ϱ Φ = β ( ) q BŨ n+1 G δt ( ) P n+1 = P n + rm 1 pl BŨ n+1 G +Φ }{{} P n+1 M ϱ U n+1 = M ϱ Ũ n+1 + δt B T Φ β q Preconditioning the prediction step by one iteration of augmented Lagrangian and consistent scalar projection r = 0 : incremental projection method Scalar penalty-projection methods Penalty-projection methods 12
19 Penalty-projection methods Direct implementation with the algebraic formulation or not... Why? Continuous term corresponding to the penalization T p = r Ω u v But T p does not generally vanish for a velocity field with zero discrete divergence! It depends on the spatial discretization... Hence Introduction of an additional error due to the space discretization Error which increases with the augmentation parameter r > 0 Scalar penalty-projection methods Penalty-projection methods 13
20 Numerical experiments Green-Taylor vortices : Navier-Stokes with Dirichlet B.C. Velocity error (discrete L 2 (0, T ; L 2 (Ω) d ) norm) versus time step δt 1e-2 norme L2 erreur vitesse 1e-3 1e-4 1e-5 1e-6 semi implicite incrementale 1e-7 1e-4 1e-3 1e-2 1e-1 Pas de temps Time convergence in O(δt 2 ) Stagnation threshold = space discretization error in O(h 2 ) Scalar penalty-projection methods Numerical experiments 14
21 Numerical experiments Green-Taylor vortices : Navier-Stokes with Dirichlet B.C. Velocity error (discrete L 2 (0, T ; L 2 (Ω) d ) norm) versus time step δt 1e-2 norme L2 erreur vitesse 1e-3 1e-4 1e-5 1e-6 semi implicite incrementale r=1 1e-7 1e-4 1e-3 1e-2 1e-1 Pas de temps Time convergence in O(δt 2 ) Stagnation threshold = space discretization error in O(h 2 ) Scalar penalty-projection methods Numerical experiments 14
22 Numerical experiments Green-Taylor vortices : Navier-Stokes with Dirichlet B.C. Velocity error (discrete L 2 (0, T ; L 2 (Ω) d ) norm) versus time step δt 1e-2 norme L2 erreur vitesse 1e-3 1e-4 1e-5 1e-6 semi implicite incrementale r=1 r=10 1e-7 1e-4 1e-3 1e-2 1e-1 Pas de temps Time convergence in O(δt 2 ) Stagnation threshold = space discretization error in O(h 2 ) Scalar penalty-projection methods Numerical experiments 14
23 Numerical experiments Green-Taylor vortices : Navier-Stokes with Dirichlet B.C. Velocity error (discrete L 2 (0, T ; L 2 (Ω) d ) norm) versus time step δt 1e-2 norme L2 erreur vitesse 1e-3 1e-4 1e-5 1e-6 semi implicite incrementale r=1 r=10 r=100 1e-7 1e-4 1e-3 1e-2 1e-1 Pas de temps Time convergence in O(δt 2 ) Stagnation threshold = space discretization error in O(h 2 ) Scalar penalty-projection methods Numerical experiments 14
24 Numerical experiments Green-Taylor vortices : Navier-Stokes with Dirichlet B.C. Pressure error (discrete L 2 (0, T ; L 2 (Ω)) norm) versus time step δt norme L2 erreur pression 1e-1 1e-2 1e-3 1e-4 1e-5 semi implicite incrementale r=1 r=10 r=100 1e-6 1e-4 1e-3 1e-2 1e-1 Pas de temps Stagnation threshold = space discretization error in O(h 2 ) Scalar penalty-projection methods Numerical experiments 15
25 Numerical experiments Artificial pressure boundary layer : Stokes with Dirichlet B.C. on a disk incremental projection penalty-projection r=1 p h p L (Ω) = p h p L (Ω) = Scalar penalty-projection methods Numerical experiments 16
26 Numerical experiments Artificial pressure boundary layer : Stokes with Dirichlet B.C. on a disk penalty-projection r=100 p h p L (Ω) = implicit scheme p h p L (Ω) = Scalar penalty-projection methods Numerical experiments 16
27 Numerical experiments Stokes with open boundary condition at a channel outflow Velocity error (discrete L 2 (0, T ; L 2 (Ω) d ) norm) versus time step δt Norme L2 erreur vitesse 1e-1 1e-2 1e-3 1e-4 1e-5 1e-6 1e-7 1e-8 semi implicite incrementale r=1 r=10 r=100 r=1000 r= e-9 1e-4 1e-3 1e-2 1e-1 Pas de temps Time convergence in O(δt 2 ) Stagnation threshold = space discretization error in O(h 2 ) Scalar penalty-projection methods Numerical experiments 17
28 Numerical experiments Stokes with open boundary condition at a channel outflow Pressure error (discrete L 2 (0, T ; L 2 (Ω)) norm) versus time step δt 1e+0 norme L2 erreur pression 1e-1 1e-2 1e-3 1e-4 1e-5 semi implicite incrementale r=1 r=10 r=100 r=1000 r= e-6 1e-4 1e-3 1e-2 1e-1 Pas de temps Stagnation threshold = space discretization error in O(h 2 ) Scalar penalty-projection methods Numerical experiments 18
29 Numerical experiments Navier-Stokes open flow around a circular cylinder Taylor-Hood P2-P1 finite elements Scheme 2nd order in time Reynolds number : Re = 100 (µ = 0.01) Unstructured mesh with 4200 elements Complete study of this configuration : see Khadra et al., IJNMF 2000 with FVM (MAC mesh) and iterative augmented Lagrangian algorithm. Splitting error varying as O( 1 r ) Scalar penalty-projection methods Numerical experiments 19
30 Numerical experiments Navier-Stokes open flow around a circular cylinder Norme L2 erreur de fractionnement 1e-1 1e-2 1e-3 1e-4 Pas de temps = 0.05 s Pas de temps = 0.02 s Pas de temps =0.01 s Pas de temps = s 1e-5 1e+0 1e+1 1e+2 1e+3 paramètre de pénalisation Splitting error varying as O( 1 r ) Scalar penalty-projection methods Numerical experiments 19
31 Theoretical analysis for Stokes-Dirichlet problem Unsteady homogeneous model problem with ϱ = 1 Ω connected and bounded open domain of IR d, Lipschitz boundary Ω, T > 0 u u + p = f dans Ω ]0, T [ t u = 0 dans Ω ]0, T [ u = 0 sur Ω ]0, T [ u(x, 0) = u 0 (x) dans Ω u denotes the velocity vector, p the pressure field. For f L 2 (0, T ; L 2 (Ω) d ) and u 0 V given, there exists a unique solution u L (0, T ; H) L 2 (0, T ; V ) and p L 2 (0, T ; L 2 0 (Ω)) See [Temam, Girault and Raviart, 1986] Scalar penalty-projection methods Analysis of the scalar penalty-projection method 20
32 Theoretical analysis for Stokes-Dirichlet problem [Shen, Guermond, 1996] : standard projection (pressure-correction form) [Guermond and Shen, 2003] : standard projection (velocity-correction) [Guermond and Shen, 2004] : rotational variant of [Timmermans et al., 1996] [Angot, Jobelin, Latché, Preprint 2006] : penalty-projection Analysis for small values of the penalty parameter r Theorem (Splitting error - fully discrete case in time and space) Energy estimates of splitting errors compared to Euler implicit scheme there exists c = c(ω, T, f, u 0, h) > 0 such that : for 1 n N, [ n δt e k 2 0 k=0 [ n δt ẽ k 2 0 k=0 ] 1 [ 2 n + k=0 ] 1 [ 2 n + k=0 δt ẽ k 2 0 ] 1 2 δt ɛ k 2 0 ] 1 2 c min(δt 2, δt3/2 r 1/2 ) c max(1, 1 r 1/2 )δt3/2. Scalar penalty-projection methods Analysis of the scalar penalty-projection method 21
33 Theoretical analysis for Stokes-Dirichlet problem Analysis for large values of the penalty parameter r Theorem (Splitting error - fully discrete case in time and space) Energy estimates of splitting errors compared to Euler implicit scheme there exists c = c(ω, T, f, u 0, h) > 0 such that : for 1 n N, [ n e n 0 + ẽ n 0 + δt ẽ k 2 0 [ n δt e k 2 0 k=0 ɛ n h [ n δt ɛ k 2 0 k=0 k=0 ] 1 [ 2 n + ] 1 2 k=0 δt ẽ k 2 0 ] 1 2 ] 1 2 c δt 1 2 r c δt r c 1 r 1/2 c 1 r. Scalar penalty-projection methods Analysis of the scalar penalty-projection method 22
34 Theoretical analysis for Stokes-Dirichlet problem [ n δt ẽ k 2 0 k=0 ] 1 2 : velocity splitting error ẽ n = ū n ũ n discrete L 2 (0, T ; L 2 (Ω) d ) norm versus time step δt 1e-1 1e-2 1e-3 1e-4 r = 0 r = 1 r = 5 r = 10 r = 50 r = 100 r = 500 r = e-5 1e-6 1e-7 1e-8 1e-4 1e-3 1e-2 1e-1 Pas de temps Scalar penalty-projection methods Analysis of the scalar penalty-projection method 23
35 Theoretical analysis for Stokes-Dirichlet problem [ n δt ẽ k 2 0 k=0 ] 1 2 : velocity splitting error ẽ n = ū n ũ n discrete L 2 (0, T ; L 2 (Ω) d ) norm versus penalty parameter r > 0 1e-2 1e-3 1e-4 Pas de temps: 0.1 Pas de temps: 0.05 Pas de temps: Pas de temps: Pas de temsp: Pas de temps: 3.125e-3 Pas de temps: 1.563e-3 Pas de temps: 7.813e-4 Pas de temps: 3.906e-4 1e-5 1e-6 1e-7 1e-8 1e+0 1e+1 1e+2 1e+3 Paramètre de pénalisation Scalar penalty-projection methods Analysis of the scalar penalty-projection method 23
36 Theoretical analysis for Stokes-Dirichlet problem [ n δt ɛ k 2 0 k=0 ] 1 2 : pressure splitting error ɛ n = p n p n discrete L 2 (0, T ; L 2 (Ω)) norm versus penalty parameter r > 0 1e-1 1e-2 1e-3 Pas de temps: 0.1 Pas de temps: 0.05 Pas de temps: Pas de temps: Pas de temps: Pas de temps: 3.125e-3 Pas de temps: 1.563e-3 Pas de temps: 7.813e-4 Pas de temps: 3.906e-4 1e-4 1e-5 1e-6 1e-7 1e+0 1e+1 1e+2 1e+3 Paramètre de pénalisation Scalar penalty-projection methods Analysis of the scalar penalty-projection method 23
37 Outlines 1 Projection methods for incompressible flows 2 Scalar penalty-projection methods 3 Vector penalty-projection methods New class of vector penalty-projection methods Analysis of the vector penalty-projection method for the Stokes problem 4 Conclusion and perspectives Vector penalty-projection methods 24
38 New class of vector penalty-projection methods Work in progress with Caltagirone and Fabrie A two-parameter family with vector projection step Penalty-prediction ( step : augmentation) parameter r 0 Dũ ϱ n+1 n+1 + (u,n+1 )ũ n+1 τ (ũ n+1 ) δt r ( ũ n+1) + p n = f n+1 ũ n+1 = u n+1 D on Ω D, and p n n + τ (ũ n+1 ) n = f n+1 N on Ω N p n+1 = p n r ũ n+1 Vector [ projection ( step : penalty parameter ) 0 < η δt ] η ϱ n+1 βq δt ûn+1 + (u,n+1 )û n+1 τ (û n+1 ) ( û n+1) = ( ũ n+1) û n+1 = 0 on Ω D, and ( p n+1 p n )n + τ (û n+1 ) n = 0 on Ω N Correction step : u n+1 = ũ n+1 + û n+1, and p n+1 = p n+1 1 η un+1 Vector penalty-projection methods New class of vector penalty-projection methods 25
39 New class of vector penalty-projection methods Why does it work well? Very nice conditioning properties of the vector penalty-correction step as η 0 since the right-hand side becomes more and more adapted to the left-hand side operator : only a few iterations (2 or 3) of MILU0-BiCGStab solver whatever the mesh step h! Original boundary conditions not spoiled through a scalar projection step Penalty-correction step quasi-independent on the density or viscosity when η 0 for non-homogeneous incompressible flows It is only an approximate projection scheme since u n+1 0 But both the velocity divergence and the splitting error can be made negligible with respect to the time and space discretization errors when η 0 cheap method for small values of r < 1 and η (zero machine) Vector penalty-projection methods New class of vector penalty-projection methods 26
40 Analysis for Stokes-Dirichlet problem with r = 0 Error estimates of the one-parameter family : r = 0 and 0 < η δt Theorem (Splitting error in the semi-discrete case) Energy estimates of splitting errors w.r. to Euler implicit scheme there exists C = C(Ω, T, f, u 0 ) > 0 such that : for 1 n N, [ n e n 0 + δt e k 2 0 [ n δt e k 2 0 k=0 k=0 ] 1 2 ] 1 2 C η δt C ηδt [ n η ɛ n 0 + δt ɛ k 2 0 k=0 ] 1 2 C η u n 0 = e n 0 C max ( 1, η δt ) ηδt. Vector penalty-projection methods Analysis of the vector penalty-projection method 27
41 Outlines 1 Projection methods for incompressible flows 2 Scalar penalty-projection methods 3 Vector penalty-projection methods 4 Conclusion and perspectives Conclusion and perspectives 28
42 Conclusion Scalar penality-projection methods for incompressible flows Reduce the splitting error up to make it negligible Suppress pressure boundary layers for moderate values of r > 0 Optimal convergence for r > 0 with open boundary conditions Generalization to dilatable and low Mach number flows : see [Jobelin et al., Preprint 2006] Require efficient preconditioning for large values of r multi-level preconditioner for FVM with MAC mesh : see [Kortas, PhD 1997]. Conclusion and perspectives 29
43 Perspectives Vector penality-projection methods for incompressible and non-homogeneous flows Small values of r < 1 are sufficient to get a good pressure field Approximate projection with a vector correction step all the cheaper than η 0 Same convergence properties as before Vector correction step all the less dependent on density or viscosity than η 0 Other numerical experiments (in progress) Error estimates : two-parameter class, Navier-Stokes, outflow B.C., variable density flows... Conclusion and perspectives 30
Vector Penalty-Projection Methods for the Solution of Unsteady Incompressible Flows
Author manuscript, published in "5th International Symposium on Finite Volumes for Complex Applications, Aussois : France 2008)" Vector Penalty-Projection Methods for the Solution of Unsteady Incompressible
More informationVector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy
Vector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy Philippe Angot, Rima Cheaytou To cite this version: Philippe Angot, Rima Cheaytou. Vector penalty-projection
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationEfficient Solvers for the Navier Stokes Equations in Rotation Form
Efficient Solvers for the Navier Stokes Equations in Rotation Form Computer Research Institute Seminar Purdue University March 4, 2005 Michele Benzi Emory University Atlanta, GA Thanks to: NSF (MPS/Computational
More informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationIn honor of Prof. Hokee Minn
In honor of Prof. Hokee Minn It is my greatest honor to be here today to deliver the nd Minn Hokee memorial lecture. I would like to thank the Prof. Ha for the nomination and the selecting committee for
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationLagrangian acceleration in confined 2d turbulent flow
Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole
More informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationFast discrete Helmholtz-Hodge decompositions in bounded domains
Fast discrete Helmholtz-Hodge decompositions in bounded domains Philippe Angot, Jean-Paul Caltagirone, Pierre Fabrie To cite this version: Philippe Angot, Jean-Paul Caltagirone, Pierre Fabrie. Fast discrete
More informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS
Submitted to Math. Comp. ON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS J.L. GUERMOND 1 AND JIE SHEN 2 Abstract. In this paper we study the rotational form of the pressure-correction
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More information.u= 0 ρ( u t +(u. )u)= ρ g p+.[µ( u+ t u)]
THETIS is a numerical simulation tool developed by University of Bordeaux. It is a versatile code to solve different problems: fluid flows, heat transfers, scalar transports or porous mediums. The potential
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationPostprint.
http://www.diva-portal.org Postprint This is the accepted version of a paper published in Mathematical Modelling and Analysis. This paper has been peer-reviewed but does not include the final publisher
More informationNumerical Methods for Incompressible Viscous Flow
Numerical Methods for Incompressible Viscous Flow Hans Petter Langtangen Kent-Andre Mardal Dept. of Scientific Computing, Simula Research Laboratory and Dept. of Informatics, University of Oslo Ragnar
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More information7.4 The Saddle Point Stokes Problem
346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationTowards pressure-robust mixed methods for the incompressible Navier Stokes equations
Mohrenstrasse 39 10117 Berlin Germany Tel. +49 30 20372 0 www.wias-berlin.de 2016-05-05, Pittsburgh Weierstrass Institute for Applied Analysis and Stochastics Towards pressure-robust mixed methods for
More informationRobust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems
Robust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems M. Aaqib Afaq Institute for Applied Mathematics and Numerics (LSIII) TU Dortmund 13 July 2017
More informationHIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationDiscretizations & Efficient Linear Solvers for Problems Related to Fluid Flow
Clemson University TigerPrints All Dissertations Dissertations 8-2017 Discretizations & Efficient Linear Solvers for Problems Related to Fluid Flow Ryan R. Grove Clemson University, rgrove@g.clemson.edu
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationPerformance Analysis of Parallel Alternating Directions Algorithm for Time Dependent Problems
Performance Analysis of Parallel Alternating Directions Algorithm for Time Dependent Problems Ivan Lirkov 1, Marcin Paprzycki 2, and Maria Ganzha 2 1 Institute of Information and Communication Technologies,
More informationFast Solvers for Unsteady Incompressible Flow
ICFD25 Workshop p. 1/39 Fast Solvers for Unsteady Incompressible Flow David Silvester University of Manchester http://www.maths.manchester.ac.uk/~djs/ ICFD25 Workshop p. 2/39 Outline Semi-Discrete Navier-Stokes
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationDISCRETE HODGE HELMHOLTZ
Monografías Matemáticas García de Galdeano 39, 1 10 (2014) DISCRETE HODGE HELMHOLTZ DECOMPOSITION Etienne Ahusborde, Mejdi Azaïez, Jean-Paul Caltagirone, Mark Gerritsma and Antoine Lemoine Abstract. This
More informationA Review of Solution Techniques for Unsteady Incompressible Flow
Zeist 2009 p. 1/57 A Review of Solution Techniques for Unsteady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/57 PDEs Review : 1966 1999 Update : 2000
More informationAccelerating incompressible fluid flow simulations on hybrid CPU/GPU systems
Accelerating incompressible fluid flow simulations on hybrid CPU/GPU systems Yushan Wang 1, Marc Baboulin 1,2, Karl Rupp 3,4, Yann Fraigneau 1,5, Olivier Le Maître 1,5 1 Université Paris-Sud, France 2
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationAlgebraic Multigrid Methods for the Oseen Problem
Algebraic Multigrid Methods for the Oseen Problem Markus Wabro Joint work with: Walter Zulehner, Linz www.numa.uni-linz.ac.at This work has been supported by the Austrian Science Foundation Fonds zur Förderung
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS
MATHEMATICS OF COMPUTATION Volume 73, Number 248, Pages 1719 1737 S 0025-5718(03)01621-1 Article electronically published on December 19, 2003 ON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationDNS of the Taylor-Green vortex at Re=1600
DNS of the Taylor-Green vortex at Re=1600 Koen Hillewaert, Cenaero Corentin Carton de Wiart, NASA Ames koen.hillewaert@cenaero.be, corentin.carton@cenaero.be Introduction This problem is aimed at testing
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationAdvances in the mathematical theory of the finite element immersed boundary method
Advances in the mathematical theory of the finite element immersed boundary method Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 12, 2014
More informationOn preconditioned Uzawa-type iterations for a saddle point problem with inequality constraints
On preconditioned Uzawa-type iterations for a saddle point problem with inequality constraints Carsten Gräser and Ralf Kornhuber FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 43, No., pp. 239 258 c 2005 Society for Industrial and Applied Mathematics ERROR ANALYSIS OF PRESSURE-CORRECTION SCHEMES FOR THE TIME-DEPENDENT STOKES EQUATIONS WITH OPEN BOUNDARY
More informationDirect numerical simulation and control of the flow over a backward facing ramp
Direct numerical control of the flow 2 1 IMB Université de Bordeaux, France 2 IMB Institut polytechnique de Bordeaux,France Inria Bordeaux Sud Ouest MC2 project 3rd GDR Contrôle des décollements Symposium
More informationThe Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations
The Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations George Biros Lexing Ying and Denis Zorin Courant Institute of Mathematical Sciences, New York University, NY 10012,
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationA variation on the infsup condition
Cinquième foire européenne aux éléments finis Centre International de Rencontres Mathématiques Marseille Luminy, 17-19 mai 2007 A variation on the infsup condition François Dubois CNAM Paris and University
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationAnalysis and Computation of Navier-Stokes equation in bounded domains
1 Analysis and Computation of Navier-Stokes equation in bounded domains Jian-Guo Liu Department of Physics and Department of Mathematics Duke University with Jie Liu (National University of Singapore)
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationFEM for Stokes Equations
FEM for Stokes Equations Kanglin Chen 15. Dezember 2009 Outline 1 Saddle point problem 2 Mixed FEM for Stokes equations 3 Numerical Results 2 / 21 Stokes equations Given (f, g). Find (u, p) s.t. u + p
More informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
More informationAn advanced ILU preconditioner for the incompressible Navier-Stokes equations
An advanced ILU preconditioner for the incompressible Navier-Stokes equations M. ur Rehman C. Vuik A. Segal Delft Institute of Applied Mathematics, TU delft The Netherlands Computational Methods with Applications,
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationA Projection FEM for Variable Density Incompressible Flows
Journal of Computational Physics 65, 67 88 (000) doi:0.006/jcph.000.6609, available online at http://www.idealibrary.com on A Projection FEM for Variable Density Incompressible Flows J.-L. Guermond and
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationSimulation of low Mach number flows
Simulation of low Mach number flows Improvement of accuracy and convergence of the TAU code using time derivative preconditioning Ralf Heinrich, Braunschweig, 9th February 008 Simulation of low Mach number
More informationConferences Fictitious domain methods for numerical realization of unilateral problems. J.H. (jointly with T. Kozubek and R. Kučera) p.
p. /2 Conferences 28 Fictitious domain methods for numerical realization of unilateral problems J.H. (jointly with T. Kozubek and R. Kučera) Lyon 23-24.6. 28 p. 2/2 Outline Fictitious domain methods -
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationDivergence-free wavelet projection method for incompressible viscous flow on the square
Divergence-free wavelet proection method for incompressible viscous flow on the square Souleymane Kadri Harouna, Valérie Perrier To cite this version: Souleymane Kadri Harouna, Valérie Perrier. Divergence-free
More informationProjection Methods for Rotating Flow
Projection Methods for Rotating Flow Daniel Arndt Gert Lube Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics IACM - ECCOMAS 2014 Computational Modeling of Turbulent and
More informationSchur Complements on Hilbert Spaces and Saddle Point Systems
Schur Complements on Hilbert Spaces and Saddle Point Systems Constantin Bacuta Mathematical Sciences, University of Delaware, 5 Ewing Hall 976 Abstract For any continuous bilinear form defined on a pair
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationSpectral Vanishing Viscosity Method for Large-Eddy Simulation of Turbulent Flows
Journal of Scientific Computing ( 2006) DOI: 10.1007/s10915-005-9029-9 Spectral Vanishing Viscosity Method for Large-Eddy Simulation of Turbulent Flows Richard Pasquetti 1 Received October 6, 2004; accepted
More information