An Extended Finite Element Method for a Two-Phase Stokes problem
|
|
- Patience French
- 5 years ago
- Views:
Transcription
1 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 Physics Bubble Curvature Weak formulation Pressure enrichment Pressure enrichment Q x h Ghost penalty Velocity enrichment Velocity enrichment V x h Nitsche XFEM formulation The weighted average Choice of λ Numerical examples First example Second example Conclusion 16 1
2 The theory of this work is based on the notes [Leh15] of the XFEM - lecture in the summer term 2015 held by Dr. C. Lehrenfeld. In this project we consider an approximation of a two-phase stokes problem where the interface is not resolved by the triangulation. In section 1 we discuss the problem and the physical intepretation of the occuring equations by presenting two model problems. Depending on the parameters, the solution can have weak and strong discontinuities. Therefore new approximation spaces are designed in the section 2 and 3. In the last section we present numerical examples, where an exact solution exists and analyze the convergence of the new methods. 1 Problem description We consider an open domain Ω R 2, which is separated into two subdomains Ω 1 and Ω 2 by the interface = Ω 1 Ω 2. We assume that one phase is completely surrounded by the other, i.e. Ω =. Let {T h } h>0 be a family of shape regular meshes of Ω. We only consider simplicial meshes T h. The elements of such a mesh T h are simplices, i.e. triangles in 2D, and it holds h = max {h T T T h }, where h T := diam(t ) is the local element size. In general the interface does not coincide with the boundaries of the elements of T h. We call the interface unfitted w.r.t. the mesh T h. We assume that the resolution close to the interface is sufficiently high, such that can be resolved by the triangulation, in the sense that, if T =: T, then T can be represented as the graph of a function on a planar cross-section of T. Further the resolution close to the interface should be sufficiently high, such that T h c T for all T cut by, and with a constant c > 0 independent of h and the cut position of. We consider the two-phase Stokes problem on the domain Ω, which results as div (µ i D(u)) + p = ρ i g in Ω i, i = 1, 2, (1a) div (u) = 0 in Ω i, i = 1, 2, (1b) u = 0 on, (1c) µd(u) n + pn = f on, (1d) u = 0 on Ω, (1e) with the domainwise constant density ρ, domainwise constant viscosity µ, given right hand side g [L 2 (Ω)] 2 and surface tension force f [L 2 ()] 2. D(u) denotes the symmetric gradient D(u) := u + ( u). 1.1 Physics Before we introduce two examples of a two phase Stokes problem, we first take a closer look at the equations (1a) - (1e) and the occuring material parameters ρ, µ. Equation (1a) is the so called momentum conservation equation and can be derived by balancing the momentums in and on a control volume under the assumptions of a Newtonian fluid. Another important fact for the derivation of (1a) is that we consider a creeping flow, which means that the advective inertial forces are small compared to the viscous forces. The second equation (1b) can be derived from the conservation of mass and (1e) is a given boundary condition. The new aspects compared to an one-phase stokes problem, are the equations (1c) and (1d), which will be called the interface conditions. We are looking for continuous velocity fields what leads to (1c). Equation (1d) describes the momentum balance at the interface. The stress from the domains Ω 1 and Ω 2 2
3 Ω 1 : µ 1 = 1 ρ i g = ( ) 0 0 Ω 2 : µ 2 = 100 ( ) 0 ρ i g = 1 : f = ( ) 0 0 Figure 1: A bubble in a surrounding fluid balances with the force f on the interface. A typical f would be the surface tension force. To classify the fluid in Ω i we have two material parameters: The dynamic viscosity µ i which describes the resistance due to deformation by the stress tensor and the density ρ i. To complete our equations we consider the gravitational force ρ i g as right hand side, which acts as a pushing force on the fluid domains Bubble In the first example (Figure 1) we want to simulate a bubble with a very small density ρ 1 (for example air) surrounded by a fluid with a higher density ρ 2 (e.g. water). Also the viscosities are chosen differently with µ 1 = 1 and µ 2 = 100. The force on the interface is set to zero. Due to the different viscosities we expect a jump in the pressure to satisfy the interface condition (1d). In Figure 2 and 3 the velocity and pressure (turned 180 degrees) fields are displayed. We observe, that the gravitational force drags the mass of the outer fluid down and therefore the lighter bubble rises up and the displayed velocity field is produced. Far away from the bubble the pressure is essentially the hydrostatic pressure and increases almost linear. Closer to the bubble (Ω 2 ), we see that the pressure is still increasing towards the bottom, but due to the interaction with the bubble the increase is no longer linear. Furthermore, we see that we get a discontinuity across the interface. Only at two points (cf. Figure 3) the pressure is continuous as the normal component of D(u) vanishes there and thus equation (1e) demands p = 0. Note that the pressure gradient (from bottom to top) is not a driving force in this example! Curvature In the second example (Figure 4) we demonstrate how the surface tension force affects the velocity field. The surface tension can be modeled as σκn, where κ is the mean curvature in each point on the interface and σ is a material parameter which controls the influence of the curvature. The curvature κ has different impacts depending on the substance system. For instance the resulting surface tension force is higher for an air bubble than for an oil bubble 3
4 Figure 2: Velocity plot for example µd(u) n + pn = 0 Figure 3: Pressure plot for example
5 Ω 1 : µ 1 = 1 ρ i g = ( ) 0 0 Ω 2 : µ 2 = 1 ρ i g = ( ) 0 0 : f = σκn Figure 4: Rounded square and considered paramters for example (with the same gemoetry). When we look at this example (see Figure 4 for the considered geometrical configuration and the parameters. Note: ρ and µ are globally constant), one can see that the curvature is much higher in the corners as on the straight sides (where it is actual zero). In Figure 5 we see how this affects the velocity. In the corners the flow is directed towards the center of the bubble and on the straight sides outwards. This goes hand in hand with the physical expectations, as the bubble tends to a spherical shape if the surface tension is the dominant force. 1.2 Weak formulation We multiply (1a) by v V = [H 1 (Ω)] 2 and (1b) by q Q = L 2 (Ω). Partial integration and the interface condition (1d), then lead to the weak formulation. Approximating the spaces V and Q by V h and Q h leads to the discrete problem: Find (u, p) V h Q h, such that 1 2 µ i (D(u), D(v)) Ωi (p, div(v)) Ωi = ρ i (g, v) Ωi + (f, v) v V h, (2a) i=1,2 i=1,2 (q, div(u)) Ωi = 0 q Q h, (2b) i=1,2 where (a, b) Ωi := Ω i a b dx. The choice of the finite element spaces V h and Q h with respect to (1c) and (1e) will be discussed in the following sections. b 2 Pressure enrichment 2.1 Pressure enrichment Q x h A popular choice for the velocity-pressure space V h Q h is the Taylor-Hood element, which is LBB-stable. One chooses globally continuous piecewise polynomials of degree k for the velocity, 5
6 Figure 5: Velocity plot for example and globally continuous piecewise polynomials of degree k 1 for the pressure } V h := {v [C(Ω)] 2 v T [P k (T )] 2 for all T T h, { } Q h := q C(Ω) q T P k 1 (T ) for all T T h. In this work we consider the element with k = 2. However, due to the surface tension force f at the interface, the pressure solution p will have jump discontinuities which can not be resolved by Q h very well. In particular it holds inf p h p L p h Q 2 (Ω) h1/2 p H 1 (Ω 1,2 ), (3) h where. 2 H k (Ω 1,2 ) :=. 2 H k (Ω 1 ) +. 2 H k (Ω 2 ). This estimate is sharp, even if p is arbitrarily smooth on the subdomains Ω 1, Ω 2, as long as it has a jump across the interface. The situation can not be improved by choosing piecewise polynomials of higher degree for Q h. Therefore we introduce the enriched space Q h := Q h Q x h, where Qx h is chosen as follows: Define J := {1,..., n}, where n = dim(q h ), and let (ϕ i ) i J be the nodal basis in Q h. Let J := {j J supp ϕ j } be the index set of those basis functions in Q h which have their support intersected by. Let χ 2 be the characteristic function with χ 2 (x) = 0 for x Ω 1 and χ 2 (x) = 1 for x Ω 2. Using this, for j J we introduce the so-called enrichment function Φ j (x) := χ 2 (x) χ 2 (x j ), where x j is the vertex with index j corresponding to ϕ j. We introduce the new basis functions ϕ j := Φ jϕ j, j J. Now we can define Q x h := span { ϕ j j J }. 6
7 With the velocity-pressure space V h Q h, jump discontinuities of the pressure at the interface can be resolved, and so (3) can be improved to inf p h p L p h Q 2 (Ω) hk p H k+1 (Ω 1,2 ). (4) h (a) p h Q h (b) p h Q h Figure 6: Solution without and with pressure space enrichment using triangulation and interface position as in Figure 7a. Enriching the pressure space leads to a new problem. The discrete inf-sup condition obtained by the LBB-stable Taylor-Hood element is no longer given. In particular, with the enriched space Q h, the discretization error and the conditioning of the linear system are no longer independent of the position of the interface, c.f. Figure 6/7. y y interface interface x x (a) Almost no bad cut positions (b) Many bad cut positions Figure 7: Triangulation and interface positions 2.2 Ghost penalty To fix the problems observed, we introduce the so-called ghost penalty stabilization: Let F h be the set of facets corresponding to the elements of T h. Define Ω i,h := Ω i {T T h T } 7
8 } and let Fi,h {F := F h F Ω i,h be the set of facets which have some part in Ω i,h. Then, the ghost penalty term can be defined as J(p, q) := γ h 3 ( E i p n, E i q n ) F, (5) i=1,2 F Fi,h for (p, q) Q h Q h, where E i is the canonical extension operator (on the discrete space Q h ).The stabilization parameter γ is O(1). In particular, the term J(p, q) is independent of the cut position within the elements which results in the robustness of the method, c.f. Figure 8. The problem without the ghost penalty stabilization was, that the enriched space Q h was to big to fulfill the LBB condition. With the new stabilization the LBB condition is still not valid, but by bounding the gradients close to the interface, one can show the discrete inf-sup condition of the big system. (a) without ghost penalty (b) with ghost penalty Figure 8: p h Q h using triangulation and interface position as in Figure 7b. 3 Velocity enrichment 3.1 Velocity enrichment V x h In the past section we enriched the pressure space and obtained Q h, s.t. jumps in the pressure solution at the interface can be resolved. Now, due to (1d), if µ 1 µ 2, the velocity solution u will have weak discontinuities, i.e. kinks at the interface. Those discontinuities can not be resolved very well by the approximation space V h. In particular, for a function u which has a kink at the interface, we have the approximation error estimate inf v h u L v h V 2 (Ω) h3/2 u H 2 (Ω 1,2 ), (6) h which again, is sharp, even if u is arbitrarily smooth on the subdomains Ω 1, Ω 2. Increasing the piecewise polynomial degree in V h again will not cure the situation. To achieve better results, we introduce the enriched velocity approximation space Vh := V h Vh x. The space Vx h is constructed from V h similarly as Q x h has been constructed from Q h. The only difference is that the polynomial degree of V h is one degree higher than the degree of Q h and thus we can not refer to a nodal basis for the construction of Vh x : Therefore let again J := {1,..., n}, where n = dim(v h ) and (ϕ i ) i J a basis of V h. Now define J := {j J supp ϕ j } the index 8
9 set of basis functions which are non trivial on Ω 1 and non trivial on Ω 2. Let χ 2 again be the indicator function with χ 2 (x) = 0 for x Ω 1 and χ 2 (x) = 1 for x Ω 2. Now define ϕ j := χ 2ϕ j for j J and V x h := span { ϕ j j J }. Note, that while it is irrelevant for the theoretical construction of Vh whether we choose ϕ j := χ 2 ϕ j or ϕ j := (1 χ 2)ϕ j as new basis function for j J, it does make a difference for the arising linear system and the computational structure. For details on the right choice we refer to the lecture notes [Leh15]. 3.2 Nitsche XFEM formulation We start with a few simple calculations, the result of which facilitate the introduction of the considered Nitche discretization. Consider the term (f, v) in (2a) and equation (1d), (f, v) = σ(u, p)n v ds, where σ i (u, p) := µ i D(u i ) + Ip i. Since n 1 = n 2, we get σn v = σ 1 n 1 v 1 + σ 2 n 2 v 2 = σ 1 n 1 v 2 + σ 2 n 2 v 2 + σ 1 n 1 (v 1 v 2 ) = (σ 1 σ 2 )n }{{} 1 v 2 + σ 1 n 1 (v 1 v 2 ) =: A =f but a similar computation shows σn v = σ 1 n 1 v 1 + σ 2 n 2 v 1 + σ 2 n 2 (v 2 v 1 ) = (σ 1 σ 2 )n }{{} 1 v 1 + σ 2 n 1 (v 1 v 2 ) =: B =f Now since A = B it holds for all 0 κ 1, κ 2 1 with κ 1 + κ 2 = 1, that A = κ 1 A + κ 2 B = B. This observation shows σn v = κ 1 A + κ 2 B = f (κ 1 v 2 + κ 2 v 1 ) + (κ 1 σ 1 + κ 2 σ 2 )n 1 (v 1 v 2 ) = f (κ 1 v 2 + κ 2 v 1 ) + {σ(u, p) } n 1 v, where the weighted average is defined as {σ } := κ 1 σ 1 + κ 2 σ 2. We obtain (f, v) = N c ((u, p), (v, q)) + (f, κ 1 v 2 + κ 2 v 1 ), where N c ((u, p), (v, q)) := {σ(u, p) } n 1 v ds. We add the symmetrical counterpart N c ((v, q), (u, p)) to preserve symmetry of the bilinear form. This is consistent, due to u = 0 on for the solution u. By introducing the enriched space Vh we now can resolve kinks and jumps in the velocity. As the velocity has to be continuous, we add another symmetric (stabilization) term N s (u, v) := λ {µ } h u v ds, 9
10 which penalizes jumps in the velocity. Note that N s (u, v) = 0 for the solution u. The parameter λ has to be chosen accordingly, which will be discussed in section 3.4. We define the Nitsche term N((u, p), (v, q)) := N c ((u, p), (v, q)) + N c ((v, q), (u, p)) + N s (u, v) (7) and arrive at the Nitsche XFEM-formulation (with ghost penalty stabilization): Find (u, p) Vh Q h, such that 1 2 µ i (D(u), D(v)) Ωi (p, div(v)) Ωi + N((u, p), (v, q)) i=1,2 + (q, div(u)) Ωi + J(p, q) i=1,2 = ρ i (g, v) Ωi + (f, κ 1 v 2 + κ 2 v 1 ). i=1,2 (8) for all (v, q) V h Q h. This formulation is consistent for any valid choice of κ 1, κ 2. However, the stability of the method does depend on the precise choice of κ 1, κ The weighted average We choose the Heaviside weighting for κ 1, κ 2. This means κ i T = κ T,i, where for T T h, let { 1 if Ω 1 T > Ω 2 T κ T,1 := 0 else, (9) κ T,2 := 1 κ T,2. In particular this choice of weights guarantees, that for i = 1, 2, v Vh and all T T h κ 2 T,i T h D(v) 2 ds c tr T i D(v) 2 dx (10) holds, with a constant c tr > 0 only depending on the polynomial degree k = 2, c κ and c. This is important for the stability of the method. c κ is a fixed constant such that κ 2 i c κ T i T holds for all T T h, i = 1, 2. Since we chose the Heaviside weighting, this holds with c κ := Choice of λ For a sufficiently large λ stability of the discretization can be ensured. The drawback of choosing λ unnecessarily large is an increase in the condition number. However it is sufficient to choose λ > c λ := max {1, 4c tr } (for details see [Leh15]). In Figure 9 the differences between using and not using a velocity enrichment can be observed. Without Vh x the kink can not be approximated very well. 10
11 (a) (u h, p h ) V h Q h (b) (u h, p h ) V h Q h Figure 9: Approximation of the velocity kink 4 Numerical examples The results in the following section were calculated using the software package Netgen/NGSolve and the extension ngsxfem, which was presented in the lecture. For the new methods for the two-phase stokes problem we implemented some further extensions in the ngsxfem package. To observe the expected convergence rates, we consider an example given in [KGR15] where the exact solution is given by ( ) y u(x, y) = α(r)e r2 where r = x x 2 + y 2 { µ 1 1 if r < r α(r) = µ (µ 1 1 µ 1 2 )er2 r 2 if r r { p(x, y) = x 3 c if x Ω 1 +, 0 else where Ω := ( 1, 1) 2, Ω 1 := B 2/3 is a ball with radius r = 2/3 and a given force on the interface f = cn. Note that the coefficient α(r) is continuous and either constant (if µ 1 = µ 2 ) or not constant (if the viscosities are different), which will produce a kink in the velocity. The force f = cn will only act on the pressure and will produce a discontinuity across the interface. By that we mean that the pressure can be decomposed p = p 1 + p 2 with p 2 = χ 2 c and p 1 the solution of the problem without surface tension force. 4.1 First example In the first example we choose µ 1 = µ 2 = 1 and c = 0.5 (cf. Figure 10). 11
12 (a) velocity (b) pressure Figure 10: Solutions of the first example Due to (3) the L 2 -error of the pressure only converges with order 1/2 using no enrichment. If we use the pressure space Q h, we observe a quadradic convergence (using piecewise linears for the pressure + enrichment) (cf. Figure 11) Q h V h Ghost + Ghost h 1/2 h Figure 11: pressure error p p h L 2 (Ω) As expected we lose one order of convergence in the H 1 -error and can see that without the ghost penalty, the gradient of the pressure can not be controlled (cf. Figure 12) Q h V h + Ghost Q h V h Ghost h Figure 12: pressure error p p h H 1 (Ω) 12
13 Although we have no kinks in the velocity, we see, using no enrichment V x h, that the L2 -error of the velocity only converges with order 3/2 as the L 2 pressure error converges only with order 1/2, consider interface condition (1d) (cf. Figure 13). If the pressure space is enriched, we observe the optimal order 3 independent of the choice for the velocity space Q h V h + Ghost + Ghost h 3/2 h Figure 13: velocity error u u h L 2 (Ω) In Figure 14 we see one order less than in Figure 13, as expected for the H 1 norm Q h V h Ghost + Ghost h Figure 14: velocity error u u h H 1 (Ω) 13
14 4.2 Second example In the second example we choose µ 1 = 1, µ 2 = 10 and c = 0.5 (cf. Figure 15). (a) velocity (b) pressure Figure 15: Solutions of the second example Similar as in Figure 11 we observe the convergence rate 1/2 without enrichment, but also in the case (and also with the ghost penalty) the rate is not optimal, due to the coupling between the pressure and the velocity in equation (1d). Note, that in this example we have a kink in the velocity across the interface, which means that a bad approximation for the velocity impacts the pressure approximation (cf. Figure 17 and Figure 16). Figure 16: Pressure with and without velocity enrichment 14
15 Q h V h Ghost + Ghost h 1/2 h Figure 17: pressure error p p h L 2 (Ω) The H 1 -error shows the same situation as in Figure 17 but we see that the ghost penalty is (also in the case of a velocity enrichment) important to get an optimal order as the ghost penalty penalizes large gradients in the pressure close to the interface (cf. Figure 18) Q h V h + Ghost Q h V h Ghost h Figure 18: pressure error p p h H 1 (Ω) 15
16 In Figure 19 and 20 the L 2 - and H 1 -error of the velocity are shown. As expected due to (6), if no velocity enrichment is used, we get no better convergence rate than 3/2 respectively 1/2. Including an enrichment for the velocity and the pressure, we get the optimal rates Q h V h + Ghost + Ghost h 3/2 h Figure 19: velocity error u u h L 2 (Ω) Q h V h Ghost Q h V h Ghost h 1/2 h Figure 20: velocity error u u h H 1 (Ω) 5 Conclusion With the combination of the already known ghost penalty stabilization for the pressure enrichment and the new derived Nitsche formulation for the velocity enrichment one is now able to resolve occuring jumps and/or weak discontinuities. The numerical examples show that the optimal convergence rates can be observed. An open question is a proof for the discrete inf-sup condition although we did not observe problems in our numerical examples what indicates that the condition should be valid. Another open question is how the new method can be extended to time dependent problems, as for example the method of lines is not straight forward for velocities with weak discontinuities. 16
17 References [KGR15] Matthias Kirchhart, Sven Gross, and Arnold Reusken. Analysis of an xfem discretization for stokes interface problems Institut für Geometrie und praktische Mathematik ( [Leh15] Christoph Lehrenfeld. Extended finite element methods for interface problems, lecture notes TU Wien. 17
Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport
Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport Christoph Lehrenfeld and Arnold Reusken Preprint No. 413 December 2014 Key words: Two-phase flow, mass transport,
More informationRobust preconditioning for XFEM applied to time-dependent Stokes problems
Numerical Analysis and Scientific Computing Preprint Seria Robust preconditioning for XFEM applied to time-dependent Stokes problems S. Gross T. Ludescher M.A. Olshanskii A. Reusken Preprint #40 Department
More informationAnalysis of a high order trace finite element method for PDEs on level set surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 7 Analysis of a high order trace finite element method for PDEs on level set surfaces Jörg Grande *, Christoph Lehrenfeld and Arnold Reusken * Institut für Geometrie
More informationHigh Order Unfitted Finite Element Methods for Interface Problems and PDEs on Surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 9 High Order Unfitted Finite Element Methods for Interface Problems and PDEs on Surfaces Christoph Lehrenfeld and Arnold Reusken Institut für Geometrie und Praktische
More informationNitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem
Nitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 333 November 2011 Key words: transport problem, Nitsche method,
More informationngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve (https://github.com/ngsxfem/ngsxfem)
ngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve (https://github.com/ngsxfem/ngsxfem) Christoph Lehrenfeld (with contributions from F. Heimann, J. Preuß, M. Hochsteger,...) lehrenfeld@math.uni-goettingen.de
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
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 informationDownloaded 11/21/16 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 38, No. 6, pp. A349 A3514 c 016 Society for Industrial and Applied Mathematics ROBUST PRECONDITIONING FOR XFEM APPLIED TO TIME-DEPENDENT STOKES PROBLEMS SVEN GROSS, THOMAS LUDESCHER,
More informationAnalysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems
Analysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 340 April 2012 Key words: transport problem, Nitsche method, XFEM,
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 informationNUMERICAL SIMULATION OF INCOMPRESSIBLE TWO-PHASE FLOWS WITH A BOUSSINESQ-SCRIVEN INTERFACE STRESS TENSOR
NUMERICAL SIMULATION OF INCOMPRESSIBLE TWO-PHASE FLOWS WITH A BOUSSINESQ-SCRIVEN INTERFACE STRESS TENSOR ARNOLD REUSKEN AND YUANJUN ZHANG Abstract. We consider the numerical simulation of a three-dimensional
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 informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationOptimal preconditioners for Nitsche-XFEM discretizations of interface problems
Optimal preconditioners for Nitsche-XFEM discretizations of interface problems Christoph Lehrenfeld Arnold Reusken August 13, 2014 In the past decade, a combination of unfitted finite elements (or XFEM)
More informationA STOKES INTERFACE PROBLEM: STABILITY, FINITE ELEMENT ANALYSIS AND A ROBUST SOLVER
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,
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 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 informationAN EXTENDED FINITE ELEMENT METHOD APPLIED TO LEVITATED DROPLET PROBLEMS
AN EXTENDED FINITE ELEMENT METHOD APPLIED TO LEVITATED DROPLET PROBLEMS PATRICK ESSER, JÖRG GRANDE, ARNOLD REUSKEN Abstract. We consider a standard model for incompressible two-phase flows in which a localized
More informationarxiv: v1 [math.na] 13 Aug 2014
THE NITSCHE XFEM-DG SPACE-TIME METHOD AND ITS IMPLEMENTATION IN THREE SPACE DIMENSIONS CHRISTOPH LEHRENFELD arxiv:1408.2941v1 [math.na] 13 Aug 2014 Abstract. In the recent paper [C. Lehrenfeld, A. Reusken,
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
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 informationDiplomarbeit. Ausgefu hrt am Institut fu r. Analysis und Scientific Computing. der Technischen Universita t Wien. unter der Anleitung von
Diplomarbeit Pressure Robust Discretizations for Navier Stokes Equations: Divergence-free Reconstruction for Taylor-Hood Elements and High Order Hybrid Discontinuous Galerkin Methods Ausgefu hrt am Institut
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationUNIFORM PRECONDITIONERS FOR A PARAMETER DEPENDENT SADDLE POINT PROBLEM WITH APPLICATION TO GENERALIZED STOKES INTERFACE EQUATIONS
UNIFORM PRECONDITIONERS FOR A PARAMETER DEPENDENT SADDLE POINT PROBLEM WITH APPLICATION TO GENERALIZED STOKES INTERFACE EQUATIONS MAXIM A. OLSHANSKII, JÖRG PETERS, AND ARNOLD REUSKEN Abstract. We consider
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 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 informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
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 informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
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 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 information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationIFE for Stokes interface problem
IFE for Stokes interface problem Nabil Chaabane Slimane Adjerid, Tao Lin Virginia Tech SIAM chapter February 4, 24 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2 Problem statement
More informationFinite volume method for two-phase flows using level set formulation
Finite volume method for two-phase flows using level set formulation Peter Frolkovič,a,1, Dmitry Logashenko b, Christian Wehner c, Gabriel Wittum c a Department of Mathematics and Descriptive Geometry,
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 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 informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
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 informationA FINITE ELEMENT METHOD FOR ELLIPTIC EQUATIONS ON SURFACES
A FINITE ELEMENT METHOD FOR ELLIPTIC EQUATIONS ON SURFACES MAXIM A. OLSHANSKII, ARNOLD REUSKEN, AND JÖRG GRANDE Abstract. In this paper a new finite element approach for the discretization of elliptic
More informationMORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS VIVETTE GIRAULT, DANAIL VASSILEV, AND IVAN YOTOV Abstract. We investigate mortar multiscale numerical methods for coupled Stokes and Darcy
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
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 informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationNumerische Mathematik
Numer. Math. (006) 105:159 191 DOI 10.1007/s0011-006-0031-4 Numerische Mathematik Uniform preconditioners for a parameter dependent saddle point problem with application to generalized Stokes interface
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 informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationCoupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements
Coupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements V.J. Ervin E.W. Jenkins S. Sun Department of Mathematical Sciences Clemson University, Clemson, SC, 29634-0975, USA. Abstract We study
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 informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
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 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 informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
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 informationA VOLUME MESH FINITE ELEMENT METHOD FOR PDES ON SURFACES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2012) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 10-14, 2012 A VOLUME MESH FINIE ELEMEN MEHOD FOR
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 informationDevelopment and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Development and analysis of a Lagrange-Remap sharp interface solver
More informationOn Surface Meshes Induced by Level Set Functions
On Surface Meshes Induced by Level Set Functions Maxim A. Olshanskii, Arnold Reusken, and Xianmin Xu Bericht Nr. 347 Oktober 01 Key words: surface finite elements, level set function, surface triangulation,
More informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
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 informationHIGH DEGREE IMMERSED FINITE ELEMENT SPACES BY A LEAST SQUARES METHOD
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 4-5, Pages 604 626 c 2017 Institute for Scientific Computing and Information HIGH DEGREE IMMERSED FINITE ELEMENT SPACES BY A LEAST
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
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 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] 28 Nov 2014
On the stability of approximations for the Stokes problem using different finite element spaces for each component of the velocity F. Guillén González and J.R. Rodríguez Galván arxiv:1411.7930v1 [math.na]
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationOn atomistic-to-continuum couplings without ghost forces
On atomistic-to-continuum couplings without ghost forces Dimitrios Mitsoudis ACMAC Archimedes Center for Modeling, Analysis & Computation Department of Applied Mathematics, University of Crete & Institute
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
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 informationHyperviscous regularization of the Navier Stokes equation and the motion of slender swimmers
Hyperviscous regularization of the Navier Stokes equation and the motion of slender swimmers Giulio G. Giusteri Dipartimento di Matematica e Fisica N. Tartaglia Università Cattolica del Sacro Cuore International
More informationLecture No 2 Degenerate Diffusion Free boundary problems
Lecture No 2 Degenerate Diffusion Free boundary problems Columbia University IAS summer program June, 2009 Outline We will discuss non-linear parabolic equations of slow diffusion. Our model is the porous
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
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 informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
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 informationAdditive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems
Additive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems Atle Loneland 1, Leszek Marcinkowski 2, and Talal Rahman 3 1 Introduction In this paper
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 informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
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 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 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 informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More information