A Hybridizable Discontinuous Galerkin formulation for Fluid-Structure Interaction
|
|
- Adam Todd
- 5 years ago
- Views:
Transcription
1 A Hybridizable Discontinuous Galerkin formulation for Fluid-Structure Interaction Jason P. Sheldon Sco/ T. Miller Jonathan S. Pi/ Distribu8on Statement A Approved for public release DistribuMon unlimited Funding: Naval Sea Systems Command, Advanced Submarine Systems Development Office (SEA 073R) Fi#h deal.ii Users and Developers Workshop Texas A&M University, College Station, TX, USA; August 3-7,
2 Motivation Fluid- Structure Interac8on (FSI) The interacmon of a deformable structure with a surrounding or internal fluid Source
3 Problem statement Problem statement High- order accurate FSI calculamons are o#en too computamonally expensive for real- world problems Solu8on requirements High- order accurate computamons Reasonable computamonal expenses Proposed solu8on Develop the first FSI model that uses the state of the art hybridizable disconmnuous Galerkin (HDG) method 3
4 HDG best of both worlds ConMnuous Galerkin (CG) DisconMnuous Galerkin (DG) ComputaMonally efficient for low orders Hybridizable DisconMnuous Galerkin (HDG) UMlizes opmmized solvers such as stamc condensamon Scales well for high order Preserves local conservamon of mass 4
5 HDG discretized domain CG DG HDG - Global DOFs - Local DOFs Local solution Global solution trace Global boundary Example: 2D 1x2 quadratic elements Variables: u, v, e - 8 global DOFS per ( ) CG: 120 global DOFs DG: 144 global DOFs Variables: u, v, e - 8 local DOFS per (, ) μ - 2 global DOFS per ( ) HDG: 42 global DOFs + 72 local DOFs per cell 5
6 Fluid-structure interaction modeling Domain informamon and boundary condimons Fluid model Solid model Mesh model Couple Solve and post process Fluid Domain Solid Domain Mesh Domain Fluid-Solid Interface Fluid Boundary Solid Boundary 6
7 Strong form Linear elastostatics: HDG discretization Div (C [e]) = b e Sym(Grad u) = 0 (Mesh MoMon) - μ - u, e 1. Weight by weighmng funcmons ( eu, e) 2. Integrate over cell () 3. Integrate by parts > boundary terms (Γ) Weak form e : e d Grad eu : C (e)d eu [C (e)n d = Sym(e) : Grad u d + Sym(e)n (u µ) d =0 Flux continuity eµ [C (e)n d =0 eu b d Stabilization [C (e) :=C (e) S (u µ) n 7
8 ARL Weak form e : e d (Mesh MoMon) Linear elastostatics: HDG discretization Grad eu : C (e)d eu [C (e)n d = Sym(e) : Grad u d + Sym(e)n (u µ) d =0 [C (e) :=C (e) S (u µ) n eµ [C (e)n d =0 eu b d local local local global global local global global sollocal sol global = RHSlocal RHS global Using the Schur complement, any equation of the form!! A B U = F!, C D G can be rearranged into the following system of equations: D CA 1 B = G CA 1 F, AU = F B. 8
9 Linear elastostatics: verification h- refinement (Q3) p- refinement L 2 error norm u 0 e 00 e 01 L 2 error norm u 0 Q1 u 0 Q2 u 0 Q log(1/h) log(1/h) All errors converge at the expected rate of h p+1 Manufactured solumon funcmon of sines and cosines Pure Dirichlet boundary condimons 9
10 Preliminary results: Sandia geometry Upper (fixed) end of test specimen Source The 2014 Sandia Fracture Challenge (SFC2) (hole ahead of the notch) (no hole ahead of the notch) Lower (actuated) end of test specimen 10
11 Preliminary results: Sandia geometry Undeformed mesh 7777 cells Deformed solumon Showing stress Element Order Local # of DOFs Increase from Q1 Global # of DOFs Increase from Q1 Q Q % % Q % % Q % % Size of global system increases minimally, even for high orders! 11
12 Nonlinear elastodynamics: governing equations Variables: u, v, E μ - local - global Restricted to hyperelastic materials (Saint Venant-irchhoff) Weak form ee:e d Sym( e E):Grad u d + Grad eu:f C (E) d eu bh d = 1 e 2E:(Grad u) T Grad u d + Sym( E)n e (u µ) d d ev v d =0 eµ bh d =0 eu b d where b h := [F C (E)] n S (u µ)) F = I + Grad u 12
13 Preliminary results: Turek and Hron CSM3 A 0 Tip displacement at point A F g x displacement (m) Turek data HDG Q2 data unit interval (s) y displacement (m) Source Proposal for numerical benchmarking of fluid- structure interacmon between an elasmc object and laminar incompressible flow, Turek Turek data HDG Q2 data unit interval (s) 13
14 Source: HDG Navier- Stokes formulamon from MarMn ronbichler s unpublished work Navier-Stokes: governing equations Variables: v, p, L - local υ, ψ - global Weak d grad ev : v v L d + ev grad p k d + ev bh d = el : L d el : grad v d + eln (v )d =0 grad ep k v d + ep k n d =0 Flux stabilization bh := [ L + pi] n + S (v ) bh := [ L] n + S (v ) e e bh d n d =0 =0 ev f d Introduce ψ to close the system := 1 p d Remove the kernel of the gradient p d =0 14
15 Arbitrary Lagrangian-Eulerian (ALE) transformation HDG ALE transform L = grad v = Grad v f F 1 m = L f [v m ], Auxiliary HDG variables simplify ALE transformamon IntegraMon domain ALE transform f (x)dx = f f (X) J m dx, F (t) F (0) 15
16 German research association (DFG) Benchmark 2D channel flow around a cylinder Reynolds number Re=100 Source 16
17 Preliminary results: flow over cylinder DFG benchmark - > Matching HDG Navier- Stokes with DFG benchmark ALE HDG Navier- Stokes proof of concept: oscilla8ng cylinder 17
18 Fluid-Structure interaction modeling Fluid model Domain informamon and boundary condimons Solid model Mesh model Solid System Matrix Fluid Solid Coupling Mesh Solid Coupling Monolithic Coupling Solid Fluid Coupling Fluid System Matrix Mesh Fluid Coupling Solid Mesh Coupling Fluid Mesh Coupling Mesh System Matrix Solid Unkns Fluid Unkns Mesh Unkns 3 2 = Solid RHS Fluid RHS Mesh RHS Γ FS Couple Ω S µ s µ m Ω F υ s υ f Solve and post process T s n f T f n s 18
19 Proposed HDG FSI formulation Problem 1 (Fluid sub-problem) Find {L h f, vh f,ph f, h f, f h}2lh f Vf h Q h f Vh f f h such h f ev f, f J m + ev f, f J m L h f v f v h m + ev f,j m Fm T Grad p k f f f D E + Grad ev f,µ f J m L h f F T elf,j m L h f f De f, b T [n] f E elf,j m Grad v h f F 1 h f \ FSI + m f + D + elf Fm T [n f ],J m vf h f Fm T Grad ep k f,jm vf h + F T f De f, b T [n] f + b T [n] s E 8{ L e f, ev f, ep f, e f, e f}2l h f Vf h Q h f V e f f h, where p 0 f + p k f = p h f bt [n] f := FSI µ fj m L h f Fm T [n f ]+S f vf h bt [n] f := J m µf L h f + p k f + f h I F T S f := µf l f + f v h f L2 I, ev f, T b [n] f E h =(ev f,j m f f ) f, =0, m [n f ] ep k h f,j m f + h h e f, S f = he f s f,j m g Nf i FSI N, D E ef,j m Fm T h [n f ] f h h f, p 0 f,j m p 0 f m [n f ]+S f v h f Problem 3 (Mesh sub-problem) Find {u h m, F h m, µ h m}2um h F h m Mm h such that Grad eu m, C m F h m Deu m, T b E [n] m = (Grad eu m, C m (I)) heu m, C m (I) n m +(eu m, D E efm, Fm h efm, Grad um h + efm n m, u h m µ h m = efm, Deµ m, T b E [n] m + eµ m, µ h\ m µ h s = heµ m, C m (I) n m i FSI h m \ + heµ m, g FSI N N, 8{eu m, e F, eµ} 2U h m F h m f M m, where bt [n] m := C m F h m n m S m u h m µ h m, S m := µ l I. T h f =0, h f, Problem 2 (Solid sub-problem) Find {u h s, vs h, E h s, µ h h s, s }2Us h Vs h E h s Ms h V h s, such that ev s + Grad ev s, F s C s E h s + Dev b E T [n] s s =(ev s, b s ) s ees, E h s Sym( e E) s, Grad u h 1 e s 2E s, Grad u h T s Grad u h s s s s D + Sym( E) e E s n, u h s µ h s eu s eu s, vs h s s eµs, u h s eµs, µ h s De s, b T [n] s h s \ FSI + De s, b T [n] s + b T [n] f E 8{eu s, ev s, E e s, eµ s, e s }2U h Vs h E h s M f s V h s, where Solid System Matrix Fluid Solid Coupling Mesh Solid Coupling F s = I + Grad u h s bt [n] s := F s C s E h s Solid Fluid Coupling Fluid System Matrix Mesh Fluid Coupling S s := µ s l s I. Solid Mesh Coupling Fluid Mesh Coupling Mesh System + h h e s, S s s f FSI ns + S s Monolithic Coupling Solid Unkns Fluid Unkns Mesh Unkns v h s 3 h s, 2 = =0, FSI Solid RHS Fluid RHS Mesh RHS = he s, g Ns i Ns,
20 Actually implementing this in deal Mul8ple Triangula8ons (for each material) and sets of dokandlers, fe- sytems, etc Using blockmatrix system Each diagonal block gets a different physics model Fluid- solid interface (off- diagonal blocks) We have two sets of global (trace) degrees of freedom unlike everywhere else in the domain Face/point tracking between abunng trias Use weak form of coupling condimons with single point quadrature from the adjacent tria Need global- global coupling! Don t want to couple the local and global solumons 20
21 Turek-Hron FSI benchmark parameter symbol value [m] property value Channel length L 2.5 s [ kg m ] Channel height H 0.41 s 0.4 Cylinder center C (0.2, 0.2) E s [e6 kg ms ] Cylinder radius r 0.05 f [ kg m ] Flag length l 0.35 µ f [ kg ms ] 1 Flag height h 0.02 Ū [ m s ] 1 Reference point A (0.6, 0.2) Reference point B (0.15, 0.2) v f in [ m y(h y) s ] 1.5 Ū (H/2) 2 Mesh: 2900 Cells 21
22 Tip displacement plots (at A) Turek- Hron FSI benchmark 8p displacement Ploqed against CG and HDG results, both using backward Euler Mme- stepping on same mesh x displacement (m) y displacement (m) Turek data CG Q1 data CG Q2 data HDG Q1 data unit interval (s) Turek data CG Q1 data CG Q2 data HDG Q1 data unit interval (s) Q1 HDG results match oscillamon period for Q1 CG results Q1 HDG results are closer in amplitude to Q2 CG results than Q1 22
23 Summary Goal: Develop new computamonal tool, a HDG FSI model HDG allows for arbitrarily high order polynomials and their inherently parallel nature leads to significant increase in computamonal efficiency Progress: Models implemented, tested, and coupled Implemented and verified HDG elasmcity and Navier- Stokes models Considered preliminary test cases Implemented HDG FSI formulamon HDG FSI: Benchmarked against Turek and Hron (backward Euler in Mme) Q1 HDG results match oscillamon period for Q1 CG results Q1 HDG results are closer in amplitude to Q2 CG results than Q1 23
24 Future research tasks / questions Implement higher order Mme- stepping, such as Diagonally implicit Runge- uqa Perform verificamon with the method of manufactured solumons for the fully- coupled FSI problem. Implement superconvergent postprocess capabilimes for HDG FSI as done for individual components. Are the various HDG formulamons for the different physics models algebraically equivalent? Are there opmmal formulamons in terms of accuracy, computamonal efficiency, and/or stability? Can we quanmfy the effects on convergence and dissipamon of varying the stability parameter S? Plonng error vs. wall Mme, at what order element does using the HDG method for FSI become more computamonally efficient than using other methods, such as our CG FSI model? 24
25 Questions Thank you QuesMons? 25
Proposal for numerical benchmarking of fluid-structure interaction between an elastic object and laminar incompressible flow
Proposal for numerical benchmarking of fluid-structure interaction between an elastic object and laminar incompressible flow Stefan Turek and Jaroslav Hron Institute for Applied Mathematics and Numerics,
More informationThe Pennsylvania State University The Graduate School A HYBRIDIZABLE DISCONTINUOUS GALERKIN METHOD FOR MODELING FLUID STRUCTURE INTERACTION
The Pennsylvania State University The Graduate School A HYBRIDIZABLE DISCONTINUOUS GALERKIN METHOD FOR MODELING FLUID STRUCTURE INTERACTION A Dissertation in Engineering Science and Mechanics by Jason
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationA monolithic FEM solver for fluid structure
A monolithic FEM solver for fluid structure interaction p. 1/1 A monolithic FEM solver for fluid structure interaction Stefan Turek, Jaroslav Hron jaroslav.hron@mathematik.uni-dortmund.de Department of
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 informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationAn arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem Yifan Wang University of Houston, Department of Mathematics
More informationLarge Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes
Large Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes Thomas Gallinger * Alexander Kupzok Roland Wüchner Kai-Uwe Bletzinger Lehrstuhl für Statik Technische Universität München
More informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
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 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 informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
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 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 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 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 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 informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
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 informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationA Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow
A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow Michael Köster, Michael Hinze, Stefan Turek Michael Köster Institute for Applied Mathematics TU Dortmund Trier,
More informationExtrapolation methods for accelerating unsteady partitioned fluid-structure interaction simulations
Int. Jnl. of Multiphysics Volume 5 Number 4 211 287 Extrapolation methods for accelerating unsteady partitioned fluid-structure interaction simulations Stephen Sachs*, Marcel Streitenberger, Dörte C. Sternel**,
More informationA monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion
1 / 41 A monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion Chen-Yu CHIANG O. Pironneau, T.W.H. Sheu, M. Thiriet Laboratoire Jacques-Louis Lions
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 informationPREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures
PREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures Kevin G. Wang Virginia Tech Patrick Lea, Alex Main, Charbel Farhat Stanford University
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
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 informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
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 informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
More informationUnified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems.
Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems. Stefan Girke WWU Münster Institut für Numerische und Angewandte Mathematik 10th
More informationSebastian Gjertsen Master s Thesis, Spring 2017
Development of a Verified and Validated Computational Framework for Fluid-Structure Interaction: Investigating Lifting Operators and Numerical Stability Sebastian Gjertsen Master s Thesis, Spring 2017
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
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 informationDiese Arbeit wurde vorgelegt am Lehrstuhl für computergestützte Analyse technischer Systeme
Diese Arbeit wurde vorgelegt am Lehrstuhl für computergestützte Analyse technischer Systeme Dreidimensionale Isogeometrische Elastische Festkörper für Fluid-Struktur Interaktionen Three-Dimensional Isogeometric
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More informationNumerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell
J.S.Howell Research Statement 1 Numerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell Fig. 1. Plot of the magnitude of the velocity (upper left) and
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
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 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 informationNumerical Solution of PDEs: Bounds for Functional Outputs and Certificates
Numerical Solution of PDEs: Bounds for Functional Outputs and Certificates J. Peraire Massachusetts Institute of Technology, USA ASC Workshop, 22-23 August 2005, Stanford University People MIT: A.T. Patera,
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationFinite Element Model of a complex Glass Forming Process as a Tool for Control Optimization
Excerpt from the Proceedings of the COMSOL Conference 29 Milan Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization Felix Sawo and Thomas Bernard Fraunhofer Institute
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
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 information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationNumerical Benchmarking of Fluid-Structure Interaction: A comparison of different discretization and solution approaches
Numerical Benchmarking of Fluid-Structure Interaction: A comparison of different discretization and solution approaches Stefan Turek, Jaroslav Hron, Mudassar Razzaq, Hilmar Wobker, and Michael Schäfer
More informationNumerical simulation of laminar incompressible fluid-structure interaction for elastic material with point constraints
Numerical simulation of laminar incompressible fluid-structure interaction for elastic material with point constraints M. Razzaq 1, J. Hron 2 and S. Turek 3 1 Institute of Applied Mathematics, TU Dortmund,
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
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 informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
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 informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationAcoustic Streaming Driven Mixing
Acoustic Streaming Driven Mixing Nitesh Nama, Po-Hsun Huang, Francesco Costanzo, and Tony Jun Huang Department of Engineering Science and Mechanics The Pennsylvania State University, State College, PA,
More informationTectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop
.. Tectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Politecnico di Milano - MOX, Dipartimento di
More informationElectronic Transactions on Numerical Analysis Volume 32, 2008
Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationElmer. Introduction into Elmer multiphysics. Thomas Zwinger. CSC Tieteen tietotekniikan keskus Oy CSC IT Center for Science Ltd.
Elmer Introduction into Elmer multiphysics FEM package Thomas Zwinger CSC Tieteen tietotekniikan keskus Oy CSC IT Center for Science Ltd. Contents Elmer Background History users, community contacts and
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
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 informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
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 informationIsogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows
Isogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows Babak S. Hosseini 1, Stefan Turek 1, Matthias Möller 2 1 TU Dortmund University, Institute
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 informationNumerik 2 Motivation
Numerik 2 Motivation P. Bastian Universität Heidelberg Interdisziplinäres Zentrum für Wissenschaftliches Rechnen Im Neuenheimer Feld 368, D-69120 Heidelberg email: Peter.Bastian@iwr.uni-heidelberg.de April
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationComputation for the Backward Facing Step Test Case with an Open Source Code
Computation for the Backward Facing Step Test Case with an Open Source Code G.B. Deng Equipe de Modélisation Numérique Laboratoire de Mécanique des Fluides Ecole Centrale de Nantes 1 Rue de la Noë, 44321
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationMultiscale Analysis of Vibrations of Streamers
Multiscale Analysis of Vibrations of Streamers Leszek Demkowicz Joint work with S. Prudhomme, W. Rachowicz, W. Qiu and L. Chamoin Institute for Computational Engineering and Sciences (ICES) The University
More informationTime-adaptive methods for the incompressible Navier-Stokes equations
Time-adaptive methods for the incompressible Navier-Stokes equations Joachim Rang, Thorsten Grahs, Justin Wiegmann, 29.09.2016 Contents Introduction Diagonally implicit Runge Kutta methods Results with
More informationHybridizable Discontinuous Galerkin with degree adaptivity for the incompressible Navier-Stokes equations
Hybridizable Discontinuous Galerkin with degree adaptivity for the incompressible Navier-Stokes equations Giorgio Giorgiani, Sonia Fernández-Méndez, Antonio Huerta To cite this version: Giorgio Giorgiani,
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 informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
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 informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationLocally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials
Locally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials Harald Klimach, Michael Gaida, Sabine Roller harald.klimach@uni-siegen.de 26th WSSP 2017 Motivation Fluid-Dynamic
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationA Mixed Finite Element Formulation for Solving Phase Change Problems with Convection
A Mixed Finite Element Formulation for Solving Phase Change Problems with Convection Youssef Belhamadia 1, Abdoulaye S. Kane 2, and André Fortin 3 1 University of Alberta, Campus Saint-Jean and Department
More informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
More informationMéthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel
Méthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel Julien Diaz, Équipe Inria Magique 3D Séminaire de Mathématiques Appliquées, Laboratoire de Mathématiques Nicolas Oresme
More informationMultiphysics Analysis of Electromagnetic Flow Valve
Multiphysics Analysis of Electromagnetic Flow Valve Jeffrey S. Crompton *, Kyle C. Koppenhoefer, and Sergei Yushanov AltaSim Technologies, LLC *Corresponding author: 13 E. Wilson Bridge Road, Suite 14,
More informationA Variational Multiscale Stabilized Finite Element. Method for the Solution of the Euler Equations of
A Variational Multiscale Stabilized Finite Element Method for the Solution of the Euler Equations of Nonhydrostatic Stratified Flows M. Moragues S. Marras M. Vázquez G. Houzeaux Barcelona Supercomputing
More informationSolvers for large-displacement fluid structure interaction problems: segregated versus monolithic approaches
Comput Mech (2008) 43:91 101 DOI 10.1007/s00466-008-0270-6 ORIGINAL PAPER Solvers for large-displacement fluid structure interaction problems: segregated versus monolithic approaches Matthias Heil Andrew
More informationA connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes
A connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes Sarah E. Malick Sponsor: Leo G. Rebholz Abstract We prove, for Stokes, Oseen, and Boussinesq
More information