FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
|
|
- Britney Lamb
- 5 years ago
- Views:
Transcription
1 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 if incompressible fluids whose stress-strain relation is given through an implicit function. In particular, the stresspower law is studied and compared to the classical power law. We propose several formulations of the problem, their stable approximations and particular examples of finite element spaces. The method is demonstrated by numerical results. Key words. element method incompressible fluid, non-newtonian fluid, implicit constitutive relation, finite AMS subject classifications. 35Q3, 65N3, 76D7. Introduction. It is well-known that for many fluids the Newton s postulate of linear relation between the traceless part S of the Cauchy stress T and the symmetric velocity gradient D is invalid. One of mostly used class of non-newtonian models, taking into account the dependence of the viscosity on the shear rate, determines S as a possibly nonlinear function of D. Despite the generality, there are some fluid materials for which this explicit dependence of the stress is not true, e.g. viscoplastic fluids. As pointed out in Málek et al. [5], the concept is also unnatural from the viewpoint of causality, namely that the force is expressed in terms of the kinematic quantities that reflect the deformation. They propose a solution through the stresspower laws that express D in terms of powers of S. A general framework of Rajagopal and Srinivasa [7] unifies the above mentioned approaches allowing T and D to be related in an implicit way. To be more specific, we consider steady flows of incompressible fluids whose motion is described by the equations div S + p = ρf, div v =, G(S, D) = in Ω, (.) where Ω R d, d {, 3}, is a bounded domain with Lipschitz boundary, ρ, v, f is the density, the velocity and the body force, respectively, D = D(v) := ( v + v ) and G : R d d R d d R d d is a continuous tensor function. For example, G can take the following forms: Fluids with shear rate dependent viscosity (power-law fluids) G(S, D) = S µ( + λ D ) r D; (.) This work was supported by the GA ČR Grant /9/97 and RVO: Institute of Mathematics AS CR, Žitná 5, Praha, Czech Republic (stebel@math.cas.cz). 9
2 9 J. STEBEL Stress-power-law fluids G(S, D) = µ ( + λµ S ) n S D, n = r (r ). (.3) The scaling has been chosen in such a way that the asymptotic behaviour of both models is the same for given µ, λ > and r (, ). Observe that when G takes the form (.) or (.3), the implicit relation (.) 3 automatically implies the incompressibility constraint tr Dv = div v =. The system has to be completed by an appropriate boundary condition, for instance by the no-slip condition v = on Ω. The difficulty of the general approach stems from the fact that the stress cannot be eliminated from the governing equations of motion. Hence the traditional velocitypressure formulation does not apply. The aim of this paper is to propose several suitable formulations involving the stress and their finite element approximation. Mathematical analysis of incompressible fluids with general implicit relation (.) 3 is far from being established, the author is aware only of the references [, 5].. Weak formulations. We will present several formulations of problem (.) involving Cauchy stress or its deviatoric part as the unknown. In all cases the stress is assumed to be only integrable while the derivatives are left on the velocity. We refer to [3] for a different approach where the stress function space is a subspace of H(Ω; div). Throughout the paper the following notation will be used. } L q {ψ := L q (Ω); ψ =, W,q := { ϕ W,q (Ω); Tr ϕ = on Ω }, L q := { ξ L q (Ω; R d d sym); tr ξ L q }, Ω L q := Lq (Ω; R d d sym,tr= ). The first formulation reads: Problem (A). Find (S, v, p) L r W,r L r such that (S, Dϕ) (p, div ϕ) = (f, ϕ) ϕ W,r, (ψ, div v) = ψ L r, (G(S, Dv), ξ) = ξ L r. If one considers the full Cauchy stress as an unknown then it is possible to eliminate the pressure as it can be seen from the next formulation: Problem (B). Find (T, v) L r W,r such that (T, Dϕ) = (f, ϕ) ϕ W,r, ( G(T δ, Dv), ξ ) = ξ L r. Here T δ := T d (tr T)I stands for the deviatoric part of T. In the last case we relax the constitutive relation by adding the symmetric velocity gradient to the list of unknowns, which yields:
3 IMPLICIT CONSTITUTIVE RELATIONS 93 Problem (C). Find (T, v, D) L r W,r L r such that (T, Dϕ) = (f, ϕ) ϕ W,r, ( G(T δ, D), η ) = η L r, (D Dv, ξ) = ξ L r. Let us note that for the Stokes problem, i.e. G(S, D) = S D, r =, (.4) the existence of a unique weak solution is well known. To keep ideas clear, we will restrict to this case in the next section. 3. Numerical analysis of problems (A)-(C). In this section we investigate the general conditions on the stability of the finite element approximation of problems (A)-(C) under the assumption (.4). Let T h be a partition of Ω into simplices with the norm h := max K Th diam K. The discrete formulations are based on the Galerkin method, i.e. we seek the approximate solutions in finite dimensional spaces L h L, L h L, W h W,, and L h L, which are to be determined. In order to emphasize the structure of these problems, we define the operators A : L (L ), B : L (W, B : L (W, ) : AT, ξ := ( T δ, ξ δ), B T, ϕ := (T, Dϕ), B p, ϕ := (p, div ϕ), and define: Problem (A) h. Find (S h, v h, p h ) L h W h L h such that A B S h in (L h ), B p h = in (L h ), B B v h f in (W h ). Problem (B) h. Find (T h, v h ) L h W h such that [ ] [ ] [ ] A B Th in (Lh ) =, B v h f in (W h ). Problem (C) h. Find (D h, v h, T h ) L h W h L h such that A A D h in (L h ), B v h = f in (W h ), A B T h in (L h ). ) and Due to the saddle-point structure of these problems the choice of the finite dimensional spaces is restricted by the requirement to satisfy appropriate inf-sup conditions. At this point we recall some results of the theory of abstract saddle-point problems. We start by the classical statement (see e.g. [, 3]). Theorem 3.. Let U and P be reflexive Banach spaces, A : U U be continuous, B : P U continuous and linear. Let Z := ker B. Then the following are equivalent.
4 94 J. STEBEL. (Existence of solutions) For all f U and g P there exists (u, p) U P such that (a) A(u) Bp = f in U, B u = g in P, (b) u U/Z c g P.. (a) For all f U and u g U there exists u U such that u u g Z and (b) There exists c > such that Note that (3.) is equivalent to A(u) = f in Z, Bp U c p P. (3.) inf sup Bp, u c, p P \{} u U u U p P usually called the inf-sup condition or the Babuška-Brezzi condition. If the problem consists of more constraints, e.g. P = P P, the verification of the inf-sup condition from the above theorem may become troublesome. For this reason we mention the following result for the so-called twofold saddle-point problems. Theorem 3.. Let U, P, P be reflexive Banach spaces, B : U P, B : P P be continuous and linear. Define B : U P P by B(u, p ) := B u + B p. Then the following are equivalent.. There exist constants C, c > such that B p (U P) c p P p P and p P C u U (u, p ) ker B.. There exists c > such that B p P c p P p P and B p U c p P p ker B. 3. There exists c > such that B p (U P) c p P p P and B p P c p P p P. In the following subsections we will apply these abstract results on the particular problems and suggest examples of finite element spaces. 3.. Approximation of problem (A). In this subsection we establish the well-posedness of problem (A) h under requirements on the finite element spaces. Theorem 3.3. Let L h, L h, W h be finite dimensional subspaces of L, L, W,, respectively, which satisfy the following conditions: (i) There exists c > such that B p, ϕ sup c p ϕ W h ϕ p L h ;, (ii) {Dϕ; ϕ W h } L h. Then for every f L (Ω; R d ) problem (A) h has a unique solution (S h, p h, v h ) L h L h W h. Further there exists a constant C > independent of f such that S h + p h + v h, C f.
5 IMPLICIT CONSTITUTIVE RELATIONS 95 A particular choice of spaces that satisfy the assumptions of Theorem 3.3 is based on the P k+ /P k Taylor-Hood velocity-pressure elements. Corollary 3.4. Let k. Then L h := {S L ; S K P k (K) d d K T h }, L h := {p L C(Ω); p K P k (K) K T h }, W h := {v W, ; v K P k+ (K) d K T h } satisfy the hypothesis of Theorem Approximation of problem (B). Sufficient conditions for the well-posedness of problem (B) h are given in the following theorem. Theorem 3.5. Let L h, W h be finite dimensional subspaces of L, W,, respectively, which satisfy the following conditions: (i) {Dϕ; ϕ W h } L h ; (ii) There exists c > such that B (tr T), ϕ sup c tr T ϕ W h ϕ T L h., Then for every f L (Ω; R d ) problem (B) h has a unique solution (T h, v h ) L h W h. Further there exists a constant C > independent of f such that T h + v h, C f. The assumption (ii) arises from the fact that A is not elliptic on L. The ellipticity can be fixed if the estimate tr T c T δ T L h s.t. B T = in W h holds true, which follows from (ii). We again describe several finite element combinations that satisfy the hypothesis of Theorem 3.5. The first one is based on the P +bubble/p discontinuous velocitypressure approximation. Corollary 3.6. The spaces L h := {T L ; T K P (K) d d K T h }, W h := {v W, ; v K P (K) d B d+ (K) d K T h }, where B d+ (K) denotes the one-dimensional space of bubble functions on K, satisfy the hypothesis of Theorem 3.5. Another case with a discontinuous pressure space is the P k+ /P k discontinuous Scott-Vogelius element, which in general leads to unstable approximations, however under some additional assumptions on the topology of T h the stability has been established [8, 6, 9]. Corollary 3.7. Let T h be a barycentrically refined triangulation and k d. Then L h := {T L ; T K P k (K) d d K T h }, W h := {v W, ; v K P k+ (K) d K T h }
6 96 J. STEBEL satisfy the assumptions of Theorem 3.5. The last example is based on the edge stabilization of tr T and div v. denote the set of edges in T h and f the jump of f across a given edge. Corollary 3.8. Let k and Let E h L h := {T L ; T K P k (K) d d K T h }, W h := {v W, ; v K P k+ (K) d K T h }. Define the operators C : L (L ) and C : W, (W, ) by the formulae C T, ξ := E E h ( tr T, tr ξ ) E, C v, ϕ := E E h ( div v, div ϕ ) E. Then for every f L (Ω; R d ) and γ > the problem [ ] [ ] [ ] A + γc B Th in (Lh ) = B γc f in (W h ) v h has a unique solution (T h, v h ) L h W h. Further there exists a constant C > independent of f and γ such that T h + v h, + γ ) ( tr T h,e + div v h,e C f. E E h 3.3. Approximation of problem (C). In the case of problem (C) h the situation is similar as in the previous section. Theorem 3.9. Let L h, W h, L h be finite dimensional subspaces of L, W,, L, respectively, which satisfy the following conditions: (i) {Dϕ; ϕ W h } L h ; (ii) {T δ ; T L h } L h ; (iii) There exists c > such that B (tr T), ϕ sup c tr T ϕ W h ϕ T L h., Then for every f L (Ω; R d ) problem (C) h has a unique solution (D h, v h, T h ) L h W h L h. Further there exists a constant C > independent of f such that D h + v h, + T h C f. For the particular choice of the finite dimensional spaces, one can use the examples of L h, W h from Section 3. together with L h := {T δ ; T L h }. 4. Numerical results. The discrete problems (A) h (C) h lead to a system of nonlinear algebraic equations R(X) =, (4.)
7 IMPLICIT CONSTITUTIVE RELATIONS 97 where X R N contains the degrees of freedom. Each component of R is of the form R i (X) = r K (X, ϕ i ) + r E (X, ϕ i ), i =,... N, K T h K E E h where ϕ i is the appropriate test function. System (4.) is linearized by the Newton- Raphson method and the resulting system is solved by the sparse direct solver UMF- PACK. The numerical method has been implemented in an in-house C++ code, where the entries of the linearized matrix R are obtained from the definition of r K and r E by means of a simple automatic differentiation. The described approximations have been tested on several model problems. The constitutive relations under considerations were either the power law or the stresspower law. 4.. Stress-power law in a stenosed channel. The first example models the flow of a stress-power-law fluid through a stenosed channel, whose geometry is depicted in Figure 4.. For fixed material parameters ρ = µ = λ =, different values E Γ W Γ I.. Γ O Γ W /3 /3 Fig. 4.. Geometry of the stenosed channel. of r have been used. The following boundary conditions are prescribed: v = (x (. x ), ) on Γ I, v = on Γ W, Tn n = p + Sn n = on Γ O, v n = on Γ O. We used a computational mesh consisting of 657 triangles and approximation of problems (A)-(C) with the FE spaces described in Table 4.. In the cases (B) and (C) either the edge stabilization or the barycentric refinement were used. Apparently, (A) (B) (C) S P disc T P disc T P disc v P v P v P p P D P disc Table 4. Tested finite element spaces. all discrete problems yielded results with negligible differences. Some characteristics of the results are illustrated in Figures
8 98 J. STEBEL.5.5 n=-.38 n=-. n=.5 n=.5 n= n=-.38 n=-. n= n=.5 n= Fig. 4.. Velocity in the middle cross-section (left), along the channel (right) n=-. n= n=. n= n=-. n= n=.5 n= Fig Pressure (left) and norm of Dv (right) along the channel. 4.. Benchmark problem: Flow around cylinder. In the next example we compared the classical power-law and stress-power-law models on the benchmark problem of flow around a cylinder [4]. The geometry is depicted in Figure 4.4. The..... Fig Flow around a cylinder. following boundary conditions were prescribed: v = (.3 4x ) (.4 x ).4, on Γ I, v = on Γ W, Tn n = p + Sn n = on Γ O, v n = on Γ O. The power-law model was solved using the usual velocity-pressure formulation and P /P Taylor-Hood elements: ) ρ(v )v µ div (( + λ Dv ) r Dv + p = ρf, div v =.
9 IMPLICIT CONSTITUTIVE RELATIONS 99 For the stress-power-law model we used the formulation (B) h and P disc /P elements with jump stabilization: ρ(v )v div T = ρf, Dv = µ ( + λµ T δ ) n T δ, n = r (r ). The material parameters are ρ = λ =, µ = 3. Due to the small value of µ it is convenient to substitute T := µ T which leads to the system ρ(v )v µ div T = ρf, Dv = ( + λ T δ ) n T δ. We evaluated the following quantities: pressure drop p := p(a) p(b), drag coefficient C D := 5 Tn (, ), Γ S lift coefficient C L := 5 Tn (, ), Γ S where A = (.5,.), B = (.5,.), Γ S is the surface of the cylinder. The results obtained on a mesh consisting of 396 triangles are given in Figures For r = they are compared to the reference values computed in [4] for the Navier-Stokes system. Pressure drop over cylinder Drag coefficient power law stress power law reference power law stress power law reference Fig Flow around cylinder: Pressure drop (left), drag coefficient (right) Lift coefficient power law stress power law reference Fig Flow around cylinder: Lift coefficient. For given r the velocity fields of power-law and stress-power-law model show no differences while the pressure behaves differently. In particular, in the zone before the cylinder the pressure of the later model is always smaller, as can be seen in Figure 4.8.
10 3 J. STEBEL Velocity at channel axis: power law r=.4 r= r=3 r=4 r= Velocity at channel axis: stress power law r=.4 r= r=3 r=4 r= Fig Velocity v at x =.: power law (left), stress-power law (right) Pressure at channel axis: r=.4 power law stress power law Pressure at channel axis: r=7 power law stress power law Fig Pressure p at x =.: r =.4 (left), r = 7 (right). REFERENCES [] F. Brezzi and M. Fortin. Mixed and hybrid finite element methods. Springer-Verlag, 99. [] M. Bulíček, P. Gwiazda, J. Málek, and A. Swierczewska-Gwiazda. On steady flows of incompressible fluids with implicit power-law-like rheology. Advances in Calculus of Variations, ():9 36, 9. [3] J. Howell and N. Walkington. Inf sup conditions for twofold saddle point problems. Numerische Mathematik, 8(4): ,. [4] V. John and G. Matthies. Higher-order finite element discretizations in a benchmark problem for incompressible flows. International Journal for Numerical Methods in Fluids, 37(8): ,. [5] J. Málek, V. Průša, and K. Rajagopal. Generalizations of the Navier-Stokes fluid from a new perspective. International Journal of Engineering Science, 48:97 94,. [6] J. Qin. On the convergence of some low order mixed finite elements for incompressible fluids. PhD thesis, The Pennsylvania State University, 994. [7] K. Rajagopal and A. Srinivasa. On the thermodynamics of fluids defined by implicit constitutive relations. Zeitschrift für Angewandte Mathematik und Physik, 59(4):75 79, 8. [8] L. Scott and M. Vogelius. Norm estimates for a maximal right inverse of the divergence operator in spaces of piecewise polynomials. Modélisation mathématique et analyse numérique, 9 (): 43, 985. [9] S. Zhang. A new family of stable mixed finite elements for the 3d stokes equations. Mathematics of computation, 74(5): , 5.
Robust 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 informationGeneralized Newtonian Fluid Flow through a Porous Medium
Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the
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 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 informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
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 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 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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
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 informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
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 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 informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
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 informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationExistence of global weak solutions to implicitly constituted kinetic models of incompressible homogeneous dilute polymers
1 / 31 Existence of global weak solutions to implicitly constituted kinetic models of incompressible homogeneous dilute polymers Endre Süli Mathematical Institute, University of Oxford joint work with
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 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 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 informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationVector and scalar penalty-projection methods
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
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationAn introduction to implicit constitutive theory to describe the response of bodies
An introduction to implicit constitutive theory to describe the response of bodies Vít Průša prusv@karlin.mff.cuni.cz Mathematical Institute, Charles University in Prague 3 July 2012 Balance laws, Navier
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 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 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 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 informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
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 informationOn Mixed Methods for Signorini Problems
Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 30, 2003, Pages 45 52 ISS: 1223-6934 On Mixed Methods for Signorini Problems Faker Ben Belgacem, Yves Renard, and Leila Slimane Abstract. We
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 informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
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 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 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 informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationOn the parameter choice in grad-div stabilization for the Stokes equations
Noname manuscript No. (will be inserted by the editor) On the parameter choice in grad-div stabilization for the Stokes equations Eleanor W. Jenkins Volker John Alexander Linke Leo G. Rebholz Abstract
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 informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationVariable Exponents Spaces and Their Applications to Fluid Dynamics
Variable Exponents Spaces and Their Applications to Fluid Dynamics Martin Rapp TU Darmstadt November 7, 213 Martin Rapp (TU Darmstadt) Variable Exponent Spaces November 7, 213 1 / 14 Overview 1 Variable
More informationOn the analysis for a class of thermodynamically compatible viscoelastic fluids with stress diffusion
On the analysis for a class of thermodynamically compatible viscoelastic fluids with stress diffusion Josef Málek Nečas Center for Mathematical Modeling and Charles University, Faculty of Mathematics and
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
More informationNumerical Simulation of Powder Flow
Numerical Simulation of Powder Flow Stefan Turek, Abderrahim Ouazzi Institut für Angewandte Mathematik, Univ. Dortmund http://www.mathematik.uni-dortmund.de/ls3 http://www.featflow.de Models for granular
More informationCESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP Yiolanda Damianou, George Kaoullas,
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 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 informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationCompressible Generalized Newtonian Fluids
Nečas Center for Mathematical Modeling Compressible Generalized Newtonian Fluids Josef Málek and K. R. Rajagopal Preprint no. 2009-06 Research Team 1 Mathematical Institute of the Charles University Sokolovská
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 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 informationShape Optimization in Problems Governed by Generalised Navier Stokes Equations: Existence Analysis
Shape Optimization in Problems Governed by Generalised Navier Stokes Equations: Existence Analysis Jaroslav Haslinger 1,3,JosefMálek 2,4, Jan Stebel 1,5 1 Department of Numerical Mathematics, Faculty of
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 informationNUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES Steven
More informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationDISCONTINUOUS GALERKIN FINITE ELEMENT DISCRETIZATION FOR STEADY STOKES FLOWS WITH THRESHOLD SLIP BOUNDARY CONDITION
DISCONTINUOUS GALRKIN FINIT LMNT DISCRTIZATION FOR STADY STOKS FLOWS WITH THRSHOLD SLIP BOUNDARY CONDITION J.K. Djoko Department of Mathematics and Applied Mathematics, University of Pretoria, Private
More informationMathematical analysis of the stationary Navier-Stokes equations
Mathematical analysis of the Department of Mathematics, Sogang University, Republic of Korea The 3rd GCOE International Symposium Weaving Science Web beyond Particle Matter Hierarchy February 17-19, 2011,
More informationContents. Computational Fluid Dynamics. Fundamentals. Deformation of Materials
Contents R. Verfürth www.ruhr-uni-bochum.de/num1 Lecture Series / Bochum / Summer Term 2018 Variational Formulation of the Stokes Equations of the Stokes Equations Stationary Incompressible Navier-Stokes
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationMATH 676. Finite element methods in scientific computing
MATH 676 Finite element methods in scientific computing Wolfgang Bangerth, Texas A&M University Lecture 33.25: Which element to use Part 2: Saddle point problems Consider the stationary Stokes equations:
More informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET AND KELVIN PRINCIPLES
ON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET AND KELVIN PRINCIPLES PAVEL BOCHEV 1,3 AND MAX GUNZBURGER 2 Abstract. Least-squares finite element
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 informationA Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface
A Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface Damanik, H., Mierka, O., Ouazzi, A., Turek, S. (Lausanne, August 2013) Page 1 Motivation Polymer melts:
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationTWO-DIMENSIONAL MAGMA FLOW *
Iranian Journal of Science & Technology, Transaction A, Vol. 34, No. A2 Printed in the Islamic Republic of Iran, 2010 Shiraz University TWO-DIMENSIONAL MAGMA FLOW * A. MEHMOOD 1** AND A. ALI 2 1 Department
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
More informationOn three-dimensional flows of activated fluids
On three-dimensional flows of activated fluids Josef Málek Nečas Center for Mathematical Modeling and Mathematical institute Charles University, Faculty of Mathematics and Physics Anna Abbatiello, Tomáš
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 informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
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 informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
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 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 informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationContents. Computational Fluid Dynamics. Fundamentals. Deformation of Materials
Contents R. Verfürth Fakultät für Mathematik Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num1 Lecture Series / Milan / March 2010 of the Stokes Equations Non-Stationary Incompressible Navier-Stokes
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 informationSegment Description of Turbulence
Dynamics of PDE, Vol.4, No.3, 283-291, 2007 Segment Description of Turbulence Y. Charles Li Communicated by Y. Charles Li, received August 25, 2007. Abstract. We propose a segment description for turbulent
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
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 information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
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 informationA Numerical Investigation of Velocity-Pressure Reduced Order Models for Incompressible Flows
A Numerical Investigation of Velocity-Pressure educed Order Models for Incompressible Flows Alfonso Caiazzo a, Traian Iliescu b,1, Volker John a,c,, Swetlana Schyschlowa a a Weierstrass Institute for Applied
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 informationThe Role of the Pressure in Finite Element Methods for Incompressible Flow Problems
Weierstrass Institute for Applied Analysis and Stochastics The Role of the Pressure in Finite Element Methods for Incompressible Flow Problems Volker John Mohrenstrasse 39 10117 Berlin Germany Tel. +49
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More information