A parallel iterative approach for the Stokes Darcy coupling
|
|
- Madison Hodge
- 5 years ago
- Views:
Transcription
1 A parallel iterative approach for the Stokes Darcy coupling Marco Discacciati*, Alfio Quarteroni**, Alberto Valli*** * Institute of Analysis and Scientific Computing, EPFL ** Institute of Analysis and Scientific Computing, EPFL, and MOX, Polytechnic of Milan *** Department of Mathematics, University of Trento A parallel iterative approach for the Stokes Darcy coupling p.1/20
2 Stokes Darcy coupling The coupling of the Stokes system and the Darcy equation has been recently considered as a model for the filtration of fluids through porous media (percolation of water from a basin through the ground). A parallel iterative approach for the Stokes Darcy coupling p.2/20
3 Stokes Darcy coupling The coupling of the Stokes system and the Darcy equation has been recently considered as a model for the filtration of fluids through porous media (percolation of water from a basin through the ground). Let us start describing the coupled problem. We denote by Ω S the fluid region, and by Ω D the ground region. Moreover, Γ := Ω S Ω D will be the interface between Ω S and Ω D. A parallel iterative approach for the Stokes Darcy coupling p.2/20
4 Stokes Darcy coupling (cont.) The Stokes system is given by { ν u + p = f in Ω S (1) divu = 0 in Ω S, where u is the velocity field, p is the pressure, ν > 0 is the kinematic viscosity and f is a given force field. A parallel iterative approach for the Stokes Darcy coupling p.3/20
5 Stokes Darcy coupling (cont.) The Stokes system is given by { ν u + p = f in Ω S (1) divu = 0 in Ω S, where u is the velocity field, p is the pressure, ν > 0 is the kinematic viscosity and f is a given force field. The Darcy equation is given by ( ) K (2) div N ϕ = 0 in Ω D, where ϕ is the piezometric head, K is the hydraulic conductivity tensor and N > 0 is the volumetric porosity. A parallel iterative approach for the Stokes Darcy coupling p.3/20
6 Stokes Darcy coupling (cont.) For simplicity, as boundary conditions let us assume that: (3) u = 0 on Ω S \ Γ ϕ = 0 on Ω D \ Γ. A parallel iterative approach for the Stokes Darcy coupling p.4/20
7 Stokes Darcy coupling (cont.) For simplicity, as boundary conditions let us assume that: (3) u = 0 on Ω S \ Γ ϕ = 0 on Ω D \ Γ. The interface conditions [matching of the normal velocity and of the normal stress] are: (4) u n = K N ϕ n on Γ T(u,p) n = gϕn on Γ where g is the gravity acceleration, and the fluid stress tensor is given by T ij (u,p) := ν(d i u j + D j u i ) pδ ij. Here n denotes the unit normal vector on Γ, pointing from Ω S into Ω D. A parallel iterative approach for the Stokes Darcy coupling p.4/20
8 Well-posedness and iterative solution algorithms It can be proved that the coupled problem (1) (4) has a unique solution [Discacciati and Quarteroni, ENUMATH 2001; Layton, Schieweck and Yotov, SINUM 2003]. A parallel iterative approach for the Stokes Darcy coupling p.5/20
9 Well-posedness and iterative solution algorithms It can be proved that the coupled problem (1) (4) has a unique solution [Discacciati and Quarteroni, ENUMATH 2001; Layton, Schieweck and Yotov, SINUM 2003]. We are interested in devising an efficient solution algorithm that uses as building blocks a Stokes solver and an elliptic solver. To do this, it is natural to introduce iterative algorithms. A parallel iterative approach for the Stokes Darcy coupling p.5/20
10 Well-posedness and iterative solution algorithms It can be proved that the coupled problem (1) (4) has a unique solution [Discacciati and Quarteroni, ENUMATH 2001; Layton, Schieweck and Yotov, SINUM 2003]. We are interested in devising an efficient solution algorithm that uses as building blocks a Stokes solver and an elliptic solver. To do this, it is natural to introduce iterative algorithms. Let us present some of the iterative algorithms that have been proposed for the solution of (1) (4). A parallel iterative approach for the Stokes Darcy coupling p.5/20
11 Dirichlet/Neumann algorithm Dirichlet/Neumann: given ϕ 0 on Γ, for m 0 A parallel iterative approach for the Stokes Darcy coupling p.6/20
12 Dirichlet/Neumann algorithm Dirichlet/Neumann: given ϕ 0 on Γ, for m 0 solve the Stokes problem with T(u m,p m ) n = gϕ m n on Γ A parallel iterative approach for the Stokes Darcy coupling p.6/20
13 Dirichlet/Neumann algorithm Dirichlet/Neumann: given ϕ 0 on Γ, for m 0 solve the Stokes problem with T(u m,p m ) n = gϕ m n on Γ solve the Darcy problem with K N ϕm+1/2 n = u m n on Γ A parallel iterative approach for the Stokes Darcy coupling p.6/20
14 Dirichlet/Neumann algorithm Dirichlet/Neumann: given ϕ 0 on Γ, for m 0 solve the Stokes problem with T(u m,p m ) n = gϕ m n on Γ solve the Darcy problem with pose on Γ K N ϕm+1/2 n = u m n on Γ ϕ m+1 := θϕ m+1/2 + (1 θ)ϕ m, where θ > 0 is an acceleration parameter. A parallel iterative approach for the Stokes Darcy coupling p.6/20
15 Dirichlet/Neumann algorithm (cont.) [The name Dirichlet/Neumann is somehow arbitrary: indeed, we are solving two Neumann problems. However, for the velocity field (u ΩS, ϕ ΩD ) the step in Ω S is a Neumann step, whereas the step in Ω D is a Dirichlet step (for the normal component...).] A parallel iterative approach for the Stokes Darcy coupling p.7/20
16 Dirichlet/Neumann algorithm (cont.) [The name Dirichlet/Neumann is somehow arbitrary: indeed, we are solving two Neumann problems. However, for the velocity field (u ΩS, ϕ ΩD ) the step in Ω S is a Neumann step, whereas the step in Ω D is a Dirichlet step (for the normal component...).] The algorithm is convergent (for a suitable choice of θ) [Discacciati and Quarteroni, Comput. Visual. Sci. 2004]. A parallel iterative approach for the Stokes Darcy coupling p.7/20
17 Dirichlet/Neumann algorithm (cont.) [The name Dirichlet/Neumann is somehow arbitrary: indeed, we are solving two Neumann problems. However, for the velocity field (u ΩS, ϕ ΩD ) the step in Ω S is a Neumann step, whereas the step in Ω D is a Dirichlet step (for the normal component...).] The algorithm is convergent (for a suitable choice of θ) [Discacciati and Quarteroni, Comput. Visual. Sci. 2004]. However: for a finite element approximation, the convergence is independent of the mesh parameter h, but depends heavily on the viscosity ν and the conductivity K A parallel iterative approach for the Stokes Darcy coupling p.7/20
18 Robin/Robin algorithm 1 Sequential Robin/Robin: given η 0 on Γ, for m 0 A parallel iterative approach for the Stokes Darcy coupling p.8/20
19 Robin/Robin algorithm 1 Sequential Robin/Robin: given η 0 on Γ, for m 0 solve the Darcy problem with γ D K N ϕm+1 n + gϕ m+1 = η m on Γ A parallel iterative approach for the Stokes Darcy coupling p.8/20
20 Robin/Robin algorithm 1 Sequential Robin/Robin: given η 0 on Γ, for m 0 solve the Darcy problem with γ D K N ϕm+1 n + gϕ m+1 = η m on Γ solve the Stokes problem with n T(u m+1,p m+1 ) n + γ S u m+1 n = gϕ m+1 γ S K N ϕ m+1 n on Γ τ T(u m+1,p m+1 ) n = 0 on Γ A parallel iterative approach for the Stokes Darcy coupling p.8/20
21 Robin/Robin algorithm 1 Sequential Robin/Robin: given η 0 on Γ, for m 0 solve the Darcy problem with γ D K N ϕm+1 n + gϕ m+1 = η m on Γ solve the Stokes problem with pose on Γ n T(u m+1,p m+1 ) n + γ S u m+1 n = gϕ m+1 γ S K N ϕ m+1 n on Γ τ T(u m+1,p m+1 ) n = 0 on Γ η m+1 := n T(u m+1,p m+1 ) n + γ D u m+1 n. A parallel iterative approach for the Stokes Darcy coupling p.8/20
22 Robin/Robin algorithm 1 (cont.) Here γ S > 0 and γ D > 0 are suitable acceleration parameters, and τ is a unit tangent vector on Γ. A parallel iterative approach for the Stokes Darcy coupling p.9/20
23 Robin/Robin algorithm 1 (cont.) Here γ S > 0 and γ D > 0 are suitable acceleration parameters, and τ is a unit tangent vector on Γ. The algorithm is proved to be convergent (at least for γ S = γ D ) [Discacciati, PhD Thesis, EPFL 2004]. A parallel iterative approach for the Stokes Darcy coupling p.9/20
24 Robin/Robin algorithm 1 (cont.) Here γ S > 0 and γ D > 0 are suitable acceleration parameters, and τ is a unit tangent vector on Γ. The algorithm is proved to be convergent (at least for γ S = γ D ) [Discacciati, PhD Thesis, EPFL 2004]. Moreover: A parallel iterative approach for the Stokes Darcy coupling p.9/20
25 Robin/Robin algorithm 1 (cont.) Here γ S > 0 and γ D > 0 are suitable acceleration parameters, and τ is a unit tangent vector on Γ. The algorithm is proved to be convergent (at least for γ S = γ D ) [Discacciati, PhD Thesis, EPFL 2004]. Moreover: it can be interpreted as an alternating direction algorithm (which is useful for tuning the parameters γ S and γ D A parallel iterative approach for the Stokes Darcy coupling p.9/20
26 Robin/Robin algorithm 1 (cont.) Here γ S > 0 and γ D > 0 are suitable acceleration parameters, and τ is a unit tangent vector on Γ. The algorithm is proved to be convergent (at least for γ S = γ D ) [Discacciati, PhD Thesis, EPFL 2004]. Moreover: it can be interpreted as an alternating direction algorithm (which is useful for tuning the parameters γ S and γ D for a finite element approximation, the convergence is independent of the mesh parameter h, and, for suitable choices of γ S and γ D, in the numerical computations it looks also independent of the viscosity ν and the conductivity K A parallel iterative approach for the Stokes Darcy coupling p.9/20
27 Robin/Robin algorithm 2 The preceding Robin/Robin algorithm is sequential. Our aim now is to find an efficient parallel algorithm. A parallel iterative approach for the Stokes Darcy coupling p.10/20
28 Robin/Robin algorithm 2 The preceding Robin/Robin algorithm is sequential. Our aim now is to find an efficient parallel algorithm. The leading idea comes from the Neumann/Neumann algorithm (but it cannot be applied straightforwardly: lack of regularity!). A parallel iterative approach for the Stokes Darcy coupling p.10/20
29 Robin/Robin algorithm 2 The preceding Robin/Robin algorithm is sequential. Our aim now is to find an efficient parallel algorithm. The leading idea comes from the Neumann/Neumann algorithm (but it cannot be applied straightforwardly: lack of regularity!). Recalling that the Neumann/Neumann algorithm for the Laplace operator is indeed given by a couple of Dirichlet solvers, followed by a couple of (homogeneous) Neumann solvers, which play the role of correctors, we propose a double parallel Robin/Robin algorithm. A parallel iterative approach for the Stokes Darcy coupling p.10/20
30 Robin/Robin algorithm 2 (cont.) Parallel Robin/Robin: given η 0 on Γ, for m 0 A parallel iterative approach for the Stokes Darcy coupling p.11/20
31 Robin/Robin algorithm 2 (cont.) Parallel Robin/Robin: given η 0 on Γ, for m 0 solve in parallel the Stokes problem and the Darcy problem with (5) n T(u m+1,p m+1 ) n γ 1 u m+1 n = η m = gϕ m+1 + γ 1 K N ϕ m+1 n on Γ τ T(u m+1,p m+1 ) n = 0 on Γ A parallel iterative approach for the Stokes Darcy coupling p.11/20
32 Robin/Robin algorithm 2 (cont.) Parallel Robin/Robin: given η 0 on Γ, for m 0 solve in parallel the Stokes problem and the Darcy problem with (5) n T(u m+1,p m+1 ) n γ 1 u m+1 n = η m = gϕ m+1 + γ 1 K N ϕ m+1 n on Γ τ T(u m+1,p m+1 ) n = 0 on Γ A parallel iterative approach for the Stokes Darcy coupling p.11/20
33 Robin/Robin algorithm 2 (cont.) Parallel Robin/Robin: given η 0 on Γ, for m 0 solve in parallel the Stokes problem and the Darcy problem with (5) n T(u m+1,p m+1 ) n γ 1 u m+1 n = η m = gϕ m+1 + γ 1 K N ϕ m+1 n on Γ τ T(u m+1,p m+1 ) n = 0 on Γ pose (6) σ m+1 := u m+1 n + K N ϕm+1 n A parallel iterative approach for the Stokes Darcy coupling p.11/20
34 Robin/Robin algorithm 2 (cont.) solve in parallel the homogeneous (f = 0) Stokes problem and the Darcy problem with (7) n T(ω m+1,π m+1 ) n + γ 2 ω m+1 n = γ 2 σ m+1 = gχ m+1 γ 2 K N χ m+1 n on Γ τ T(ω m+1,π m+1 ) n = 0 on Γ A parallel iterative approach for the Stokes Darcy coupling p.12/20
35 Robin/Robin algorithm 2 (cont.) solve in parallel the homogeneous (f = 0) Stokes problem and the Darcy problem with (7) + n T(ω m+1,π m+1 ) n + γ 2 ω m+1 n = γ 2 σ m+1 = + gχ m+1 γ 2 K N χ m+1 n on Γ τ T(ω m+1,π m+1 ) n = 0 on Γ A parallel iterative approach for the Stokes Darcy coupling p.12/20
36 Robin/Robin algorithm 2 (cont.) solve in parallel the homogeneous (f = 0) Stokes problem and the Darcy problem with (7) n T(ω m+1,π m+1 ) n + γ 2 ω m+1 n = γ 2 σ m+1 = gχ m+1 γ 2 K N χ m+1 n on Γ τ T(ω m+1,π m+1 ) n = 0 on Γ pose on Γ (8) η m+1 := η m θ ( n T(ω m+1,π m+1 ) n + gχ m+1), where γ 1 > 0, γ 2 > 0 and θ > 0 are acceleration parameters. A parallel iterative approach for the Stokes Darcy coupling p.12/20
37 Consistency of the Robin/Robin algorithm 2 Let us show that the algorithm is consistent, namely, at convergence we have the solution. A parallel iterative approach for the Stokes Darcy coupling p.13/20
38 Consistency of the Robin/Robin algorithm 2 Let us show that the algorithm is consistent, namely, at convergence we have the solution. At convergence, from (8) it holds T(ω,π ) n = gχ n, therefore normal stresses of the corrections are matching. A parallel iterative approach for the Stokes Darcy coupling p.13/20
39 Consistency of the Robin/Robin algorithm 2 Let us show that the algorithm is consistent, namely, at convergence we have the solution. At convergence, from (8) it holds T(ω,π ) n = gχ n, therefore normal stresses of the corrections are matching. On the other hand, from (7) normal velocities of the corrections are jumping ω n + K N χ n = 2g γ 2 χ. A parallel iterative approach for the Stokes Darcy coupling p.13/20
40 Consistency of the Robin/Robin algorithm 2 (cont.) But this jump gives an additional positive term in the energy of the problem, namely, one obtains 2ν D i ωj 2 g + Ω S ij Ω D N K χ χ + Γ 2g 2 γ 2 χ 2 = 0. A parallel iterative approach for the Stokes Darcy coupling p.14/20
41 Consistency of the Robin/Robin algorithm 2 (cont.) But this jump gives an additional positive term in the energy of the problem, namely, one obtains 2ν D i ωj 2 g + Ω S ij Ω D N K χ χ + Γ 2g 2 γ 2 χ 2 = 0. Hence, ω = 0 and χ = 0, and consequently σ = 0. A parallel iterative approach for the Stokes Darcy coupling p.14/20
42 Consistency of the Robin/Robin algorithm 2 (cont.) Finally: from (6), σ is the jump of the normal velocities, therefore they are matching: u n = K N ϕ n A parallel iterative approach for the Stokes Darcy coupling p.15/20
43 Consistency of the Robin/Robin algorithm 2 (cont.) Finally: from (6), σ is the jump of the normal velocities, therefore they are matching: u n = K N ϕ n from the Robin interface condition (5), if the normal velocities are matching also the normal stresses are matching: T(u,p ) n = gϕ n A parallel iterative approach for the Stokes Darcy coupling p.15/20
44 Consistency of the Robin/Robin algorithm 2 (cont.) Finally: from (6), σ is the jump of the normal velocities, therefore they are matching: u n = K N ϕ n from the Robin interface condition (5), if the normal velocities are matching also the normal stresses are matching: T(u,p ) n = gϕ n We have found the right solution! A parallel iterative approach for the Stokes Darcy coupling p.15/20
45 Convergence of the Robin/Robin algorithm 2 The parallel Robin/Robin algorithm can be rewritten as a preconditioned Richardson scheme: η m+1 = η m + θ(k S + K D )[ψ (H S + H D )η m ], for a suitable right hand side ψ. A parallel iterative approach for the Stokes Darcy coupling p.16/20
46 Convergence of the Robin/Robin algorithm 2 The parallel Robin/Robin algorithm can be rewritten as a preconditioned Richardson scheme: η m+1 = η m + θ(k S + K D )[ψ (H S + H D )η m ], for a suitable right hand side ψ. The operators K S, K D, H S and H D are the main building blocks of the algorithm, and are defined as follows. A parallel iterative approach for the Stokes Darcy coupling p.16/20
47 Convergence of the Robin/Robin algorithm 2 (cont.) K S σ := n T(ω,π) n, where ω is a solution with the Robin datum n T(ω,π) n + γ 2 ω n = γ 2 σ on Γ τ T(ω,π) n = 0 on Γ. K S is a Robin-to-Neumann operator A parallel iterative approach for the Stokes Darcy coupling p.17/20
48 Convergence of the Robin/Robin algorithm 2 (cont.) K S σ := n T(ω,π) n, where ω is a solution with the Robin datum n T(ω,π) n + γ 2 ω n = γ 2 σ on Γ τ T(ω,π) n = 0 on Γ. K S is a Robin-to-Neumann operator K D σ := gχ, where χ is a solution with the Robin datum gχ γ 2 K N χ n = γ 2 σ on Γ. K D is a Robin-to-Dirichlet operator A parallel iterative approach for the Stokes Darcy coupling p.17/20
49 Convergence of the Robin/Robin algorithm 2 (cont.) H S η := u n, where u is a solution (for f = 0) with the Robin datum n T(u,p) n γ 1 u n = η on Γ τ T(u,p) n = 0 on Γ. H S is a Robin-to-Dirichlet operator A parallel iterative approach for the Stokes Darcy coupling p.18/20
50 Convergence of the Robin/Robin algorithm 2 (cont.) H S η := u n, where u is a solution (for f = 0) with the Robin datum n T(u,p) n γ 1 u n = η on Γ τ T(u,p) n = 0 on Γ. H S is a Robin-to-Dirichlet operator H D η := K N ϕ n, where ϕ is a solution with the Robin datum K gϕ + γ 1 N ϕ n = η on Γ. H D is a Robin-to-Neumann operator A parallel iterative approach for the Stokes Darcy coupling p.18/20
51 Convergence of the Robin/Robin algorithm 2 (cont.) These operators have the following properties (in the space L 2 (Γ)): A parallel iterative approach for the Stokes Darcy coupling p.19/20
52 Convergence of the Robin/Robin algorithm 2 (cont.) These operators have the following properties (in the space L 2 (Γ)): K S is symmetric, continuous and positive definite A parallel iterative approach for the Stokes Darcy coupling p.19/20
53 Convergence of the Robin/Robin algorithm 2 (cont.) These operators have the following properties (in the space L 2 (Γ)): K S is symmetric, continuous and positive definite K D is symmetric, continuous and non-negative A parallel iterative approach for the Stokes Darcy coupling p.19/20
54 Convergence of the Robin/Robin algorithm 2 (cont.) These operators have the following properties (in the space L 2 (Γ)): K S is symmetric, continuous and positive definite K D is symmetric, continuous and non-negative H S is symmetric, continuous and non-negative (for γ 1 small enough) A parallel iterative approach for the Stokes Darcy coupling p.19/20
55 Convergence of the Robin/Robin algorithm 2 (cont.) These operators have the following properties (in the space L 2 (Γ)): K S is symmetric, continuous and positive definite K D is symmetric, continuous and non-negative H S is symmetric, continuous and non-negative (for γ 1 small enough) H D is symmetric, continuous and positive definite A parallel iterative approach for the Stokes Darcy coupling p.19/20
56 Convergence of the Robin/Robin algorithm 2 (cont.) These operators have the following properties (in the space L 2 (Γ)): K S is symmetric, continuous and positive definite K D is symmetric, continuous and non-negative H S is symmetric, continuous and non-negative (for γ 1 small enough) H D is symmetric, continuous and positive definite Therefore, both the operator (H S + H D ) and the preconditioner (K S + K D ) 1 are symmetric, continuous and positive definite: convergence is achieved! (for a suitable choice of the parameter θ) A parallel iterative approach for the Stokes Darcy coupling p.19/20
57 Comments It is well-known that the Dirichlet-to-Neumann operator is symmetric, continuous and positive definite from the energy trace space H 1/2 (Γ) into ts dual; here we have seen that: the Robin-to-Neumann operator is symmetric, continuous and positive definite in L 2 (Γ), whereas: the Robin-to-Dirichlet operator is symmetric, continuous and non-negative in L 2 (Γ) A parallel iterative approach for the Stokes Darcy coupling p.20/20
58 Comments It is well-known that the Dirichlet-to-Neumann operator is symmetric, continuous and positive definite from the energy trace space H 1/2 (Γ) into ts dual; here we have seen that: the Robin-to-Neumann operator is symmetric, continuous and positive definite in L 2 (Γ), whereas: the Robin-to-Dirichlet operator is symmetric, continuous and non-negative in L 2 (Γ) the above results also hold for the finite element numerical approximation, uniformly with respect to the mesh parameter h A parallel iterative approach for the Stokes Darcy coupling p.20/20
Application of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution
Application of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution P. BURDA a,, J. NOVOTNÝ b,, J. ŠÍSTE a, a Department of Mathematics Czech University of Technology
More informationControl Volume Finite Difference On Adaptive Meshes
Control Volume Finite Difference On Adaptive Meshes Sanjay Kumar Khattri, Gunnar E. Fladmark, Helge K. Dahle Department of Mathematics, University Bergen, Norway. sanjay@mi.uib.no Summary. In this work
More informationPrecise FEM solution of corner singularity using adjusted mesh applied to 2D flow
Precise FEM solution of corner singularity using adjusted mesh applied to 2D flow Jakub Šístek, Pavel Burda, Jaroslav Novotný Department of echnical Mathematics, Czech echnical University in Prague, Faculty
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationENERGY-224 Reservoir Simulation Project Report. Ala Alzayer
ENERGY-224 Reservoir Simulation Project Report Ala Alzayer Autumn Quarter December 3, 2014 Contents 1 Objective 2 2 Governing Equations 2 3 Methodolgy 3 3.1 BlockMesh.........................................
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
More informationStabilized Finite Element Method for 3D Navier-Stokes Equations with Physical Boundary Conditions
Stabilized Finite Element Method for 3D Navier-Stokes Equations with Physical Boundary Conditions Mohamed Amara, Daniela Capatina, David Trujillo Université de Pau et des Pays de l Adour Laboratoire de
More informationTHE MORTAR FINITE ELEMENT METHOD IN 2D: IMPLEMENTATION IN MATLAB
THE MORTAR FINITE ELEMENT METHOD IN D: IMPLEMENTATION IN MATLAB J. Daněk, H. Kutáková Department of Mathematics, University of West Bohemia, Pilsen MECAS ESI s.r.o., Pilsen Abstract The paper is focused
More informationEulerian Techniques for Fluid-Structure Interactions - Part II: Applications
Published in Lecture Notes in Computational Science and Engineering Vol. 103, Proceedings of ENUMATH 2013, pp. 755-762, Springer, 2014 Eulerian Techniques for Fluid-Structure Interactions - Part II: Applications
More informationNonlinear Potential Flow Solver Development in OpenFOAM
Nonlinear Potential Flow Solver Development in OpenFOAM A. Mehmood Plymouth University, UK April 19,2016 A. Mehmood Table of Contents 1 Motivation 2 Solution Methodology Mathematical Formulation Sequence
More informationGradient Free Design of Microfluidic Structures on a GPU Cluster
Gradient Free Design of Microfluidic Structures on a GPU Cluster Austen Duffy - Florida State University SIAM Conference on Computational Science and Engineering March 2, 2011 Acknowledgements This work
More informationAdaptive Isogeometric Analysis by Local h-refinement with T-splines
Adaptive Isogeometric Analysis by Local h-refinement with T-splines Michael Dörfel 1, Bert Jüttler 2, Bernd Simeon 1 1 TU Munich, Germany 2 JKU Linz, Austria SIMAI, Minisymposium M13 Outline Preliminaries:
More informationPeriod #10: Multi dimensional Fluid Flow in Soils (II)
Period #10: Multi dimensional Fluid Flow in Soils (II) A. Review Our objective is to solve multi dimensional fluid flow problems in soils. Last time, mass conservation and Darcy s Law were used to derive
More informationImagery for 3D geometry design: application to fluid flows.
Imagery for 3D geometry design: application to fluid flows. C. Galusinski, C. Nguyen IMATH, Université du Sud Toulon Var, Supported by ANR Carpeinter May 14, 2010 Toolbox Ginzburg-Landau. Skeleton 3D extension
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationNumerical Experiments
77 Chapter 4 Numerical Experiments 4.1 Error estimators and adaptive refinement Due to singularities the convergence of finite element solutions on uniform grids can be arbitrarily low. Adaptivity based
More informationFinite element method - tutorial no. 1
Martin NESLÁDEK Faculty of mechanical engineering, CTU in Prague 11th October 2017 1 / 22 Introduction to the tutorials E-mail: martin.nesladek@fs.cvut.cz Room no. 622 (6th floor - Dept. of mechanics,
More informationCHAPTER 4 CFD AND FEA ANALYSIS OF DEEP DRAWING PROCESS
54 CHAPTER 4 CFD AND FEA ANALYSIS OF DEEP DRAWING PROCESS 4.1 INTRODUCTION In Fluid assisted deep drawing process the punch moves in the fluid chamber, the pressure is generated in the fluid. This fluid
More informationGroundwater hydraulics numerical modelling
Groundwater hydraulics numerical modelling Martin Šanda, B673 martin.sanda@fsv.cvut.cz Martina Sobotková, B609 martina.sobotkova@fsv.cvut.cz CTU in Prague, FCE, Dept. of Irrigation, Drainage and landscape
More informationOptimizing Bio-Inspired Flow Channel Design on Bipolar Plates of PEM Fuel Cells
Excerpt from the Proceedings of the COMSOL Conference 2010 Boston Optimizing Bio-Inspired Flow Channel Design on Bipolar Plates of PEM Fuel Cells James A. Peitzmeier *1, Steven Kapturowski 2 and Xia Wang
More informationDetermination of Free Surface in Steady-State Seepage through a Dam with Toe Drain by the Boundary Element Method
Determination of Free Surface in Steady-State Seepage through a Dam with Toe Drain by the Boundary Element Method Somchart Chantasiriwan Faculty of Engineering, Thammasat University, Rangsit Campus, Pathum
More informationDynamic Computational Modeling of the Glass Container Forming Process
Dynamic Computational Modeling of the Glass Container Forming Process Matthew Hyre 1, Ryan Taylor, and Morgan Harris Virginia Military Institute, Lexington, Virginia, USA Abstract Recent advances in numerical
More information3D Helmholtz Krylov Solver Preconditioned by a Shifted Laplace Multigrid Method on Multi-GPUs
3D Helmholtz Krylov Solver Preconditioned by a Shifted Laplace Multigrid Method on Multi-GPUs H. Knibbe, C. W. Oosterlee, C. Vuik Abstract We are focusing on an iterative solver for the three-dimensional
More informationPRACE Workshop, Worksheet 2
PRACE Workshop, Worksheet 2 Stockholm, December 3, 2013. 0 Download files http://csc.kth.se/ rvda/prace files ws2.tar.gz. 1 Introduction In this exercise, you will have the opportunity to work with a real
More informationA Contact Angle Model for the Parallel Free Surface Lattice Boltzmann Method in walberla Stefan Donath (stefan.donath@informatik.uni-erlangen.de) Computer Science 10 (System Simulation) University of Erlangen-Nuremberg
More informationPART 1: FLUID-STRUCTURE INTERACTION Simulation of particle filtration processes in deformable media
Int. Jnl. of Multiphysics Volume 2 Number 2 2008 179 PART 1: FLUID-STRUCTURE INTERACTION Simulation of particle filtration processes in deformable media Marianne Mataln 1, Gernot Boiger 1, Bernhard Gschaider
More informationIntroduction to geodynamic modelling Introduction to DOUAR
Introduction to geodynamic modelling Introduction to DOUAR David Whipp and Lars Kaislaniemi Department of Geosciences and Geography, Univ. Helsinki 1 Goals of this lecture Introduce DOUAR, a 3D thermomechanical
More informationScalable Preconditioners for Navier Stokes Equations in Hemodynamic Simulations
1/34 Scalable Preconditioners for Navier Stokes Equations in Hemodynamic Simulations Simone Deparis, Gwenol Grandperrin, Alfio Quarteroni MATHICSE - Chair of Modelling and Scientific Computing (CMCS) EPFL,
More informationTwo-Phase flows on massively parallel multi-gpu clusters
Two-Phase flows on massively parallel multi-gpu clusters Peter Zaspel Michael Griebel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität Bonn Workshop Programming of Heterogeneous
More information(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more
Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon
More informationImprovement of Reduced Order Modeling based on POD
Author manuscript, published in "The Fifth International Conference on Computational Fluid Dynamics (28)" Improvement of Reduced Order Modeling based on POD M. Bergmann 1, C.-H. Bruneau 2, and A. Iollo
More informationPerformances and Tuning for Designing a Fast Parallel Hemodynamic Simulator. Bilel Hadri
Performances and Tuning for Designing a Fast Parallel Hemodynamic Simulator Bilel Hadri University of Tennessee Innovative Computing Laboratory Collaboration: Dr Marc Garbey, University of Houston, Department
More informationCOMPUTATIONAL METHODS FOR ENVIRONMENTAL FLUID MECHANICS
COMPUTATIONAL METHODS FOR ENVIRONMENTAL FLUID MECHANICS Tayfun Tezduyar tezduyar@rice.edu Team for Advanced Flow Simulation and Modeling (T*AFSM) Mechanical Engineering and Materials Science Rice University
More informationA finite element-based Level-Set method for fluid - elastic solid. interaction with surface tension
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2011; 00:1 34 Published online in Wiley InterScience (www.interscience.wiley.com). A finite element-based Level-Set
More informationA brief description of the particle finite element method (PFEM2). Extensions to free surface
A brief description of the particle finite element method (PFEM2). Extensions to free surface flows. Juan M. Gimenez, L.M. González, CIMEC Universidad Nacional del Litoral (UNL) Santa Fe, Argentina Universidad
More informationThe effect of irregular interfaces on the BDDC method for the Navier-Stokes equations
153 The effect of irregular interfaces on the BDDC method for the Navier-Stokes equations Martin Hanek 1, Jakub Šístek 2,3 and Pavel Burda 1 1 Introduction The Balancing Domain Decomposition based on Constraints
More informationA Novel Approach to High Speed Collision
A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.
More informationTrue 3D CAE visualization of filling imbalance in geometry-balanced runners
True 3D CAE visualization of filling imbalance in geometry-balanced runners C.C. Chien, * C.C. Chiang, W. H. Yang, Vito Tsai and David C.Hsu CoreTech System Co.,Ltd., HsinChu, Taiwan, ROC Abstract The
More informationCenter for Computational Science
Center for Computational Science Toward GPU-accelerated meshfree fluids simulation using the fast multipole method Lorena A Barba Boston University Department of Mechanical Engineering with: Felipe Cruz,
More informationNumerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods
Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods Decheng Wan 1 and Stefan Turek 2 Institute of Applied Mathematics LS III, University of Dortmund,
More informationIMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL
More informationChapter 13. Boundary Value Problems for Partial Differential Equations* Linz 2002/ page
Chapter 13 Boundary Value Problems for Partial Differential Equations* E lliptic equations constitute the third category of partial differential equations. As a prototype, we take the Poisson equation
More informationFinite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2
Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2 David W. Trott and Matthias K. Gobbert Department of Mathematics and Statistics, University of Maryland, Baltimore County,
More informationAn investigation of spatial and temporal weighting schemes for use in unstructured mesh control volume finite element methods
ANZIAM J. 44 (E) ppc759 C779, 23 C759 An investigation of spatial and temporal weighting schemes for use in unstructured mesh control volume finite element methods S. L. Truscott I. W. Turner (Received
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationIsotropic Porous Media Tutorial
STAR-CCM+ User Guide 3927 Isotropic Porous Media Tutorial This tutorial models flow through the catalyst geometry described in the introductory section. In the porous region, the theoretical pressure drop
More informationSolution of 2D Euler Equations and Application to Airfoil Design
WDS'6 Proceedings of Contributed Papers, Part I, 47 52, 26. ISBN 8-86732-84-3 MATFYZPRESS Solution of 2D Euler Equations and Application to Airfoil Design J. Šimák Charles University, Faculty of Mathematics
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More informationIMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING
Image Anal Stereol 2001;20:143-148 Original Research Paper IMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING DAVID GARCIA 1, GUY COURBEBAISSE 2 AND MICHEL JOURLIN 3 1 European Polymer Institute (PEP),
More informationBackward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn
Backward facing step Homework Department of Fluid Mechanics Budapest University of Technology and Economics Budapest, 2010 autumn Updated: October 26, 2010 CONTENTS i Contents 1 Introduction 1 2 The problem
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationFINITE POINTSET METHOD FOR 2D DAM-BREAK PROBLEM WITH GPU-ACCELERATION. M. Panchatcharam 1, S. Sundar 2
International Journal of Applied Mathematics Volume 25 No. 4 2012, 547-557 FINITE POINTSET METHOD FOR 2D DAM-BREAK PROBLEM WITH GPU-ACCELERATION M. Panchatcharam 1, S. Sundar 2 1,2 Department of Mathematics
More informationThe role of mathematics in glaciology: Ice sheet modeling and inverse methods
The role of mathematics in glaciology: Ice sheet modeling and inverse methods Martin Truffer 1 David Maxwell 2 Ed Bueler 2 1 Geophysical Institute 2 Dept. Math. and Stat. University of Alaska Fairbanks
More informationComputer Project 3. AA Computational Fluid Dyanmics University of Washington. Mishaal Aleem March 17, 2015
Computer Project 3 AA 543 - Computational Fluid Dyanmics University of Washington Mishaal Aleem March 17, 2015 Contents Introduction........................................... 1 3.1 Grid Generator.......................................
More informationThe THex Algorithm and a Simple Darcy Solver on Hexahedral Meshes
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 8C (27) 93 92 International Conference on Computational Science, ICCS 27, 2-4 June 27, Zurich, Switzerland The THex Algorithm
More informationNumerical Implementation of Overlapping Balancing Domain Decomposition Methods on Unstructured Meshes
Numerical Implementation of Overlapping Balancing Domain Decomposition Methods on Unstructured Meshes Jung-Han Kimn 1 and Blaise Bourdin 2 1 Department of Mathematics and The Center for Computation and
More informationEstimation of a Piecewise Constant Function Using Reparameterized Level-Set Functions
Estimation of a iecewise Constant Function Using Reparameterized Level-Set Functions nga Berre 1,2, Martha Lien 2,1, and Trond Mannseth 2,1 1 Department of Mathematics, University of Bergen, Johs. Brunsgt.
More informationOn Robust Parallel Preconditioning for Incompressible Flow Problems
On Robust Parallel Preconditioning for Incompressible Flow Problems Timo Heister, Gert Lube, and Gerd Rapin Abstract We consider time-dependent flow problems discretized with higher order finite element
More informationInviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics
Contents: Navier-Stokes equations Inviscid flows Boundary layers Transition, Reynolds averaging Mixing-length models of turbulence Turbulent kinetic energy equation One- and Two-equation models Flow management
More informationNUMERICAL METHOD FOR FREE SURFACE VISCOUS FLOWS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 3, September 2008 NUMERICAL METHOD FOR FREE SURFACE VISCOUS FLOWS TIBERIU IOANA AND TITUS PETRILA Abstract. In this paper we present a new algorithm
More informationConforming Vector Interpolation Functions for Polyhedral Meshes
Conforming Vector Interpolation Functions for Polyhedral Meshes Andrew Gillette joint work with Chandrajit Bajaj and Alexander Rand Department of Mathematics Institute of Computational Engineering and
More informationFEM techniques for interfacial flows
FEM techniques for interfacial flows How to avoid the explicit reconstruction of interfaces Stefan Turek, Shu-Ren Hysing (ture@featflow.de) Institute for Applied Mathematics University of Dortmund Int.
More informationFigure 6.1: Truss topology optimization diagram.
6 Implementation 6.1 Outline This chapter shows the implementation details to optimize the truss, obtained in the ground structure approach, according to the formulation presented in previous chapters.
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationConvergence Results for Non-Conforming hp Methods: The Mortar Finite Element Method
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03042-9 Convergence Results for Non-Conforming hp Methods: The Mortar Finite Element Method Padmanabhan Seshaiyer and Manil Suri 1. Introduction
More informationA novel model for biolm growth and its resolution by using the hybrid immersed interface-level set method
A novel model for biolm growth and its resolution by using the hybrid immersed interface-level set method Juan A. Asenjo & Carlos Conca Applied Mathematics Group - Basic Sciences Department - University
More informationFree Surface Flow Simulations
Free Surface Flow Simulations Hrvoje Jasak h.jasak@wikki.co.uk Wikki Ltd. United Kingdom 11/Jan/2005 Free Surface Flow Simulations p.1/26 Outline Objective Present two numerical modelling approaches for
More informationMultigrid Solvers in CFD. David Emerson. Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK
Multigrid Solvers in CFD David Emerson Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK david.emerson@stfc.ac.uk 1 Outline Multigrid: general comments Incompressible
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationPredicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity
Predicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity November 8th, 2006 Outline Motivation Motivation Motivation for Modelling Breast Deformation Mesh Generation
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationCHAPTER 3. Elementary Fluid Dynamics
CHAPTER 3. Elementary Fluid Dynamics - Understanding the physics of fluid in motion - Derivation of the Bernoulli equation from Newton s second law Basic Assumptions of fluid stream, unless a specific
More informationNUMERICAL SIMULATION OF SHALLOW WATERS EFFECTS ON SAILING SHIP "MIRCEA" HULL
NUMERICAL SIMULATION OF SHALLOW WATERS EFFECTS ON SAILING SHIP "MIRCEA" HULL Petru Sergiu ȘERBAN 1 1 PhD Student, Department of Navigation and Naval Transport, Mircea cel Batran Naval Academy, Constanța,
More informationMultigrid Third-Order Least-Squares Solution of Cauchy-Riemann Equations on Unstructured Triangular Grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids ; : 6 Prepared using fldauth.cls [Version: /9/8 v.] Multigrid Third-Order Least-Squares Solution of Cauchy-Riemann Equations
More informationNumerical studies for Flow Around a Sphere regarding different flow regimes caused by various Reynolds numbers
Numerical studies for Flow Around a Sphere regarding different flow regimes caused by various Reynolds numbers R. Jendrny, H. Damanik, O. Mierka, S. Turek Institute of Applied Mathematics (LS III), TU
More informationAERODYNAMIC OPTIMIZATION OF REAR AND FRONT FLAPS ON A CAR
Student: Giannoni Alberto Professor: Ing. Jan Pralits, Advanced Fluid Dynamics Course Co-Professor: Ing. Matteo Colli AERODYNAMIC OPTIMIZATION OF REAR AND FRONT FLAPS ON A CAR ABSTRACT In this work we
More informationChapter 1 - Basic Equations
2.20 Marine Hydrodynamics, Fall 2017 Lecture 2 Copyright c 2017 MIT - Department of Mechanical Engineering, All rights reserved. 2.20 Marine Hydrodynamics Lecture 2 Chapter 1 - Basic Equations 1.1 Description
More informationMATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES
MATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES Goals of research Dr. Uldis Raitums, Dr. Kārlis Birģelis To develop and investigate mathematical properties of algorithms
More informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
More informationThe numerical simulation of complex PDE problems. A numerical simulation project The finite element method for solving a boundary-value problem in R 2
Universidad de Chile The numerical simulation of complex PDE problems Facultad de Ciencias Físicas y Matemáticas P. Frey, M. De Buhan Year 2008 MA691 & CC60X A numerical simulation project The finite element
More informationThe Method of Approximate Particular Solutions for Solving Elliptic Problems with Variable Coefficients
The Method of Approximate Particular Solutions for Solving Elliptic Problems with Variable Coefficients C.S. Chen, C.M. Fan, P.H. Wen Abstract A new version of the method of approximate particular solutions
More informationA method for simulation of viscous, non-linear, free-surface flows
A method for simulation of viscous, non-linear, free-surface flows Ben R. Hodges Robert L. Street Yan Zang Environmental Fluid Mechanics Laboratory, Department of Civil Engineering, Stanford University,
More informationWater seepage through a dam: A finite element approach
Water seepage through a dam: A finite element approach Ashley B. Pitcher Student Number 250098269 AM466b Final Project (Undergraduate) University of Western Ontario April 15, 2005 Abstract We consider
More informationTryItNow! Step by Step Walkthrough: Spoiler Support
TryItNow! Step by Step Walkthrough: Spoiler Support 1 2015 ANSYS, Inc. March 28, 2016 TryItNow! Step by Step Walkthrough: Spoiler Support ANSYS designed this TryItNow! experience to give you quick access
More informationVery fast simulation of nonlinear water waves in very large numerical wave tanks on affordable graphics cards
Very fast simulation of nonlinear water waves in very large numerical wave tanks on affordable graphics cards By Allan P. Engsig-Karup, Morten Gorm Madsen and Stefan L. Glimberg DTU Informatics Workshop
More informationAn Eulerian Approach to Fluid-Structure Interaction and Goal-Oriented Mesh Adaptation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2002/09/18 v1.01] An Eulerian Approach to Fluid-Structure Interaction and Goal-Oriented Mesh Adaptation
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationTerminal Falling Velocity of a Sand Grain
Terminal Falling Velocity of a Sand Grain Introduction The first stop for polluted water entering a water work is normally a large tank, where large particles are left to settle. More generally, gravity
More informationChapter 5: Introduction to Differential Analysis of Fluid Motion
Chapter 5: Introduction to Differential 5-1 Conservation of Mass 5-2 Stream Function for Two-Dimensional 5-3 Incompressible Flow 5-4 Motion of a Fluid Particle (Kinematics) 5-5 Momentum Equation 5-6 Computational
More informationFLUID MECHANICS TESTS
FLUID MECHANICS TESTS Attention: there might be more correct answers to the questions. Chapter 1: Kinematics and the continuity equation T.2.1.1A flow is steady if a, the velocity direction of a fluid
More informationComparison of Central and Upwind Flux Averaging in Overlapping Finite Volume Methods for Simulation of Super-Critical Flow with Shock Waves
Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 7, 6 (pp55665) Comparison of and Flux Averaging in Overlapping Finite Volume Methods for Simulation of
More informationNUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV)
University of West Bohemia» Department of Power System Engineering NUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV) Publication was supported by project: Budování excelentního
More informationA brief introduction to fluidstructure. O. Souček
A brief introduction to fluidstructure interactions O. Souček Fluid-structure interactions Important class computational models Civil engineering Biomechanics Industry Geophysics From http://www.ihs.uni-stuttgart.de
More informationNumerical Wave Tank Modeling of Hydrodynamics of Permeable Barriers
ICHE 2014, Hamburg - Lehfeldt & Kopmann (eds) - 2014 Bundesanstalt für Wasserbau ISBN 978-3-939230-32-8 Numerical Wave Tank Modeling of Hydrodynamics of Permeable Barriers K. Rajendra & R. Balaji Indian
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More information