Tectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop
|
|
- Dana Cameron
- 5 years ago
- Views:
Transcription
1 .. Tectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Politecnico di Milano - MOX, Dipartimento di Matematica F. Brioschi Bilbao, March 11th 2013 M X MILANO MODELING AND SCIENTIFIC COMPUTING MODELLISTICA E CALCOLO SCIENTIFICO POLITECNICO Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 1/ 55
2 Aims and motivations Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 2/ 55
3 Aims and motivations Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 3/ 55
4 Continental rift Analogue models The analogue models have the following disadvantages: they lack the thermo-mechanical description of phenomena, Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 4/ 55
5 Continental rift Analogue models The analogue models have the following disadvantages: they lack the thermo-mechanical description of phenomena, they are very expensive to setup and perform, Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 4/ 55
6 Continental rift Analogue models The analogue models have the following disadvantages: they lack the thermo-mechanical description of phenomena, they are very expensive to setup and perform, they are not general, but specific for a given region, Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 4/ 55
7 Continental rift Analogue models The analogue models have the following disadvantages: they lack the thermo-mechanical description of phenomena, they are very expensive to setup and perform, they are not general, but specific for a given region, they are not replicable. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 4/ 55
8 Aims and motivations the development of the mathematical and numerical model of mid-oceanic ridges: Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 5/ 55
9 Aims and motivations the development of the mathematical and numerical model of mid-oceanic ridges: improvement of the boundary conditions to model the ridge migration; a better thermo-mechanical model (melting); Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 5/ 55
10 Aims and motivations the development of the mathematical and numerical model of mid-oceanic ridges: improvement of the boundary conditions to model the ridge migration; a better thermo-mechanical model (melting); development of the mathematical tools for the defition of a numerical sandbox useful to reproduce the continental-rift evolution; Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 5/ 55
11 Aims and motivations the development of the mathematical and numerical model of mid-oceanic ridges: improvement of the boundary conditions to model the ridge migration; a better thermo-mechanical model (melting); development of the mathematical tools for the defition of a numerical sandbox useful to reproduce the continental-rift evolution; development and application of meshfree methods and variational integrators for the simulation of the numerical sandbox. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 5/ 55
12 Mid-ocean ridges Mathematical model We can treat the litosphere and asthenosphere as higly viscous incompressible fluid div(2ηd) p + ρg = 0 div v = 0 θ + v θ = div (κ θ) t η = η 0 e Cθ The Frank-Kamenetskii linearization is use to approximate the viscosity law η = 1 2A ( µ I 2 ) n 1 ( ) h m e E+pV RT b Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 6/ 55
13 Mid-ocean ridges The initial condition for the temperature field is the mean oceanic geotherm ( ) θ θ 0 y = erf θ M θ 0 2 κτ Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 7/ 55
14 Mid-ocean ridges Results of numerical simulations in a steady-state regime. Results of numerical simulations with ridge migration in act. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 8/ 55
15 Mid-ocean ridges Boundary conditions The previous boundary conditions have to be improved, since they do not depend upon the rheology, (in particular on the viscosity); the simulation showed they can lead to incorrect results near boundaries. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 9/ 55
16 Mid-ocean ridges Boundary conditions A different strategy has been developed; these are the main hypothesis: the boundaries are far enough from the ridge, the upwelling motion is negligible compared to shearing. y v extension. 0 x v migration h y dq u(y) = u 0 + s 0 η(q) p(y) = p 0 + y 0 f(q) dq Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 10/ 55
17 Mid-ocean ridges: Mathematical model div(2ηɛ(u)) p + f = 0 divu = 0 div(k T) + u T = 0 in Ω in Ω in Ω u = ū on Γ 1 ((2ηɛ(u) pi)n = g on Γ 2 where and T = T on Γ 3 k T n = s on Γ 4 Γ 1 Γ 2 =, Γ 3 Γ 4 =, Γ 1 Γ 2 = Ω, Γ 3 Γ 4 = Ω, Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 11/ 55
18 Mid-ocean ridges: Numerical Discretization (1) A(T) B T 0 u f B 0 0 p = C(u) T g expanding the nonlinear terms at the first order we get Newton s method A(T) B T A (T) δu f A(T) B T 0 u B 0 0 δp = 0 B 0 0 p C (u) 0 C(u) δt g 0 0 C(u) T Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 12/ 55
19 Mid-ocean ridges: Numerical Discretization (1) The Newton step can be solved using the block LU decomposition, denoting [ ] A(T) B T A = B B 0 1 = [ A (T) 0 ] B 2 = [ C (T) 0 ], so we get the following linear system [ ] [ I 0 A B T 1 B 2 A 1 I 0 S where S is the Schur complement, se we need a good Stokes solver, a good preconditioner for S ] [ ] δx = δy S = C B 2 A 1 B T 1 [ ] f g Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 13/ 55
20 Mid-ocean ridges: Preconditioners the Stokes problem is solved using the Schur complement; Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 14/ 55
21 Mid-ocean ridges: Preconditioners the Stokes problem is solved using the Schur complement; the preconditioners for A and the (Stokes-)Schur complement are, respectively, the algebraic multigrid (AMG) preconditioner (Trilinos); Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 14/ 55
22 Mid-ocean ridges: Preconditioners the Stokes problem is solved using the Schur complement; the preconditioners for A and the (Stokes-)Schur complement are, respectively, the algebraic multigrid (AMG) preconditioner (Trilinos); the viscosity-scaled pressure mass-matrix. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 14/ 55
23 Mid-ocean ridges: Preconditioners the Stokes problem is solved using the Schur complement; the preconditioners for A and the (Stokes-)Schur complement are, respectively, the algebraic multigrid (AMG) preconditioner (Trilinos); the viscosity-scaled pressure mass-matrix. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 14/ 55
24 Mid-ocean ridges: Preconditioners the Stokes problem is solved using the Schur complement; the preconditioners for A and the (Stokes-)Schur complement are, respectively, the algebraic multigrid (AMG) preconditioner (Trilinos); the viscosity-scaled pressure mass-matrix. these preconditioners are almost optimal. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 14/ 55
25 Mid-ocean ridges: Preconditioners the Stokes problem is solved using the Schur complement; the preconditioners for A and the (Stokes-)Schur complement are, respectively, the algebraic multigrid (AMG) preconditioner (Trilinos); the viscosity-scaled pressure mass-matrix. these preconditioners are almost optimal. The matrix S is preconditioned with the iterative inverse of C (this last problem is preconditioned with AMG preconditioner). Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 14/ 55
26 Mid-ocean ridges: Preconditioners the Stokes problem is solved using the Schur complement; the preconditioners for A and the (Stokes-)Schur complement are, respectively, the algebraic multigrid (AMG) preconditioner (Trilinos); the viscosity-scaled pressure mass-matrix. these preconditioners are almost optimal. The matrix S is preconditioned with the iterative inverse of C (this last problem is preconditioned with AMG preconditioner). this preconditioner seems to scale well with the number of degrees of freedom (for each numerical experiments we get about iterations to solve the temperature). Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 14/ 55
27 Mid-ocean ridges: Numerical treatment of stress BC To compute the stress boundary conditions we use the fixed point iteration method, with the following scheme given a solution (u, p, T) the boundary condition g is computed from g a new solution (u, p, T) is computed Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 15/ 55
28 Mid-ocean ridges Results of numerical simulations in a steady-state regime. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 16/ 55
29 Mid-ocean ridges Results of numerical simulations with migration Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 17/ 55
30 Mid-ocean ridges Scalability results: dof Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 18/ 55
31 Mid-ocean ridges Scalability results: dof Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 19/ 55
32 Mid-ocean ridges Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 20/ 55
33 Mid-ocean ridges Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 21/ 55
34 Mid-ocean ridges Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 22/ 55
35 Geomod example Figure : Extension: sylicon on all the base Figure : Sylicon only in the central part Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 23/ 55
36 Drucker Prager: Conf. 1 (MI, SI, V) Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 24/ 55
37 Drucker Prager: Conf. 2 (MI, SI, V) Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 25/ 55
38 Von Mises: Conf. 1 (MI, SI, V) Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 26/ 55
39 Von Mises: Conf. 2 (MI, SI, V) Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 27/ 55
40 Continental rift Mathematical model The previous model is extended to include the elastic behavior of lithosphere (upper-convected Maxwell model) div T p + ρg = 0 div v = 0 T + λt = 2ηD ρ dh dt = dp + T : D + ρr div q dt from the dimensional analysis we get the following estimates for the upwelling ρ dh dt dp dt ρr ρ dh dt div q ρ dh dt Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 28/ 55
41 Optimal transportation mesh-free method Continuum mechanic problems which involve large deformations or which are formulated with respect to the reference configuration need a continuous remeshing or some other strategies to overcome the difficulties arising from the continuous mesh update. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 29/ 55
42 Optimal transportation mesh-free method Continuum mechanic problems which involve large deformations or which are formulated with respect to the reference configuration need a continuous remeshing or some other strategies to overcome the difficulties arising from the continuous mesh update. Mesh-free methods are an alternative, since they give up the mesh and use only a point set to discretize the problem. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 29/ 55
43 Optimal transportation mesh-free method Continuum mechanic problems which involve large deformations or which are formulated with respect to the reference configuration need a continuous remeshing or some other strategies to overcome the difficulties arising from the continuous mesh update. Mesh-free methods are an alternative, since they give up the mesh and use only a point set to discretize the problem. the OTM is designed for the continuum mechanic problems, since it inherits from the continuous problem the symmetries and conservation properties, avoiding some issues of other mesh-free methods. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 29/ 55
44 Local maximum-entropy shape functions Set of nodes X = {x a R d, a = 1,..., N} Non-negative shape functions p a : conv X R 0 First order consistency conditions p a (x) = 1 a p a (x)x a = x a Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 30/ 55
45 Local maximum-entropy shape functions. Local maximum-entropy problem. For x fixed, minimize βu x (p) H(p) subject to p a 0 p a = 1 a p a x a = x a. where U x (p) = a p a x x a 2 and H(p) = a p a log p a Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 31/ 55
46 Local maximum-entropy shape functions. Local maximum-entropy problem. Defined the partition function Z(x, λ) = a e β x x a 2 +λ (x x a ). For β(x) R 0 and x int conv X. Then the unique solution of the problem is p a (x) = where λ is the unique solution of 1 Z(x, λ ) e β x x a 2 +λ (x x a ),. λ = arg min λ). λ RdZ(x, Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 32/ 55
47 Local maximum-entropy shape functions Spatial smoothness β C r = p a C r Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 33/ 55
48 Local maximum-entropy shape functions Spatial smoothness β C r = p a C r Smoothness with respect β p(β) C 0 ([0, )), p(β) C ((0, )) Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 33/ 55
49 Local maximum-entropy shape functions Spatial smoothness β C r = p a C r Smoothness with respect β p(β) C 0 ([0, )), p(β) C ((0, )) Limits of shape functions as β p a converge to a Delaunay convex approximants as β Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 33/ 55
50 Local maximum-entropy shape functions Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 34/ 55
51 Convergence Let us consider the following problem Ω = ( 1, 1)x( 1, 1). γ = βh 2 u = 2 x 2 y 2, in Ω u = 0 on Ω Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 35/ 55
52 Convergence: L 2 norm Nodes Relative error γ L 2 error Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 36/ 55
53 Convergence: H 1 norm Relative error H 1 error γ Nodes Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 37/ 55
54 Optimal transportation The motion of an inviscid fluid of non-interacting particles in R d is governed by the equations { ρ + div (ρu) = 0 (ρu) + div (ρu u) = 0 This problem can be recasted as an optimal transportation problem, that admits a variational formulation minimize b 1 J(ρ, u) = a R d 2 ρ u 2 dx dt subject to ρ + div (ρu) = 0 Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 38/ 55
55 Optimal transportation The problem can be recasted in term of deformation map φ hence the velocity and the density are given by u(x, t) = φ t (φ 1 t (x)) ρ(x, t) = ρ a (φ 1 t (x)) det φ t (φ 1 t and Benamou and Brenier (1999) showed that 1 inf J(ρ, u) = inf ρ a (x) φ b (x) x 2 dx b a R d (x)) Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 39/ 55
56 Optimal transportation - discretization Given a finite set of timestep t 0 = a, t 1, t 2,..., t N = b, the semidiscrete problem is given by J = k 1 t k+1 t k 1 2 ρ k(x) φ k k+1 (x) x 2 dx, then using a quadrature rule J = k n 1 1 t k+1 t k 2 ρ k(x n,k )v n φ k k+1 (x n,k ) x n,k 2 denoting x n,k+1 = φ k k+1 (x n,k ) and m n,k = ρ k (x n,k )v n, then the solution of this problem is x n,k+1 = x n,k + (t k+1 t k ) x n,k x n,k 1 t k t k 1 m n,k = m n,0 Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 40/ 55
57 Optimal transportation - solid and fluid flows This framework can be extended for inviscid fluids and elastic solids, adding a free energy term to the functional J J(ρ, u, C) = b a R d 1 2 ρ u 2 ρu(ρ, C) dx dt The kinetic energy term is discretized in same way, but for the free energy the trapezoidal quadrature rule is used J = 1 1 t k+1 t k 2 ρ k φ k k+1 (x) x 2 dx+ k + t [ ] k+1 t k ρ k U(ρ k, C k ) dx + ρ k+1 U(ρ k+1, C k+1 ) dx 2 Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 41/ 55
58 Optimal transportation - solid and fluid flows Two points sets are introduced: nodal set {x a }, used to construct the shape functions material set {x n }, used as quadrature points So the displacement map can be written as a linear combination of shape functions φ k k+1 (x) x = a p a (x)d a, and the problem can be completely discretized J = k n + t k+1 t k 2 1 t k+1 t k 1 2 m n,k φ k k+1 (x n,k ) x n,k 2 + [m n,k U(ρ n,k, C n,k ) + m n,k+1 U(ρ n,k+1, C n,k+1 )] Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 42/ 55
59 Optimal transportation - solid and fluid flows The solution of this problem is ( d k = (t k+1 t k )M 1 k l k + t ) k+1 t k 1 f k 2 x n,k+1 = φ k k+1 (x n,k ), x a,k+1 = φ k k+1 (x a,k ) v n,k+1 = v n,k det φ k k+1 (x a,k ) where l a,k = n m n,k x n,k x n,k 1 t k t k 1 p a,k (x n,k ) M k,ab = n m n,k p a,k (x n,k )p b,k (x n,k )I f a,k = n [ρ n,k b n,k p a,k (x n,k ) + σ n,k : p a,k (x n,k )] v n,k Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 43/ 55
60 Application to the Stokes flows ( ) inf u 2 p div u dx λ (u ū) dσ (u,p) Ω Ω Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 44/ 55
61 Optimal transportation - solid and fluid flows U = λ 2 [tre]2 + µtre 2 λ = Pa µ = Pa Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 45/ 55
62 Veselov-type discretization Choosen the discrete set of times 0 = t 0 < t 1 < t 2 < < t N 1 < t N = T then the motion is described by a sequence of positions q 0, q 1, q 2,..., q N 1, q N and the discrete action is defined by the sum S q = N 1 k=0 L d (q k, q k+1 ) where L d (q k, q k+1 ) is called discrete Lagrangian, it depends on two subsequent positions and it s a reasonable approximation of the action between the configurations q k and q k+1 L d (q k, q k+1 ) tk+1 t k L(q, q) dt Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 46/ 55
63 Variational Integrators in time: Störmer-Verlet integration Using the trapezoidal quadrature rule and approximating the velocity with finite difference we get the discrete Lagrangian [ ( ) ] 1 qk+1 q 2 k L d (q k, q k+1 ) = (t k+1 t k ) V(q k) + V(q k+1 ) 2 t k+1 t k 2 This leads to the well known Störmer-Verlet method q k+1 q k t k+1 t k q k q k 1 t k t k 1 = t k+1 t k 1 F(q k ) 2 and the related velocity Verlet method q k+1 q k = p k + t k+1 t k F(q k ) t k+1 t k 2 p k+1 p k = F(q k) + F(q k+1 ) t k+1 t k 2 Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 47/ 55
64 Midpoint rule Instead, using the midpoint quadrature rule we get the following discrete Lagrangian [ ( ) 1 qk+1 q 2 ( ) ] k qk + q k+1 L d (q k, q k+1 ) = (t k+1 t k ) V 2 t k+1 t k 2 and the related method is q k+1 q k = p k + t ( ) k+1 t k qk+1 + q k F t k+1 t k 2 2 ( ) p k+1 p k qk+1 + q k = F t k+1 t k 2 Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 48/ 55
65 Other variational integrators Most of standard symplectic methods can be written as variational integrators, among them: Newmark-β methods, the first proof of symplecticity of these methods exploits the framework of variational integrators, [ ( 1 η β (q 1 ) η β ) 2 (q 0 ) L d (q 0, q 1 ) = h V(η β (q 0 ))] 2 h symplectic partitioned Runge-Kutta Galerkin methods (Lobatto IIIA-IIIB) L d (q 0, q 1 ) ext {S NI (q) : q V N, q(0) = q 0, q(h) = q 1 } Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 49/ 55
66 Wave equation Probably the simplest model is the linear wave equation, which represents a linearization of the hyperelastic problem. The related Lagrangian density in the one dimensional space is [ ( φ L = 1 ) 2 ( c φ ) ] 2 2 t x Applying the trapezoidal quadrature rule in time we get the semi-discrete Störmer-Verlet scheme φ k+1 φ k ψ dx = p k ψ t k+1 t k c 2 φ k ψ B t k+1 t k B 2 x x dx B p k+1 p k ψ dx = c 2 t k+1 t k B x ( φk + φ k+1 2 ) ψ x dx Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 50/ 55
67 Wave equation c = 1; P 1 ; t = 0.01s; I = [0, 2π]; N x = 100. uh u log 2 2 u t position momentum Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 51/ 55
68 Wave equation: linear momentum conservation (=0) p dx log 2 2π t Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 52/ 55
69 Wave equation: energy Eh E E t Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 53/ 55
70 Conclusions and further developments improvement of BC treatment; preconditioners; preliminary results on OTM. Future developments of this project introduction of the melting; development of a fluid-structure and structure-structure method based on OTM; application of OTM and VI to geodynamic problems, coupling the mechanical problem with the thermodynamic. Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 54/ 55
71 Thank you for your attention! Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Tectonic at mid-ocean ridge 55/ 55
Well-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationThe finite difference code (fully staggered grid) includes a marker-in-cell Lagrangian marker
GSA Data Repository 2018289 Ruh et al., 2018, Shale-related minibasins atop a massive olistostrome in an active accretionary wedge setting: Two-dimensional numerical modeling applied to the Iranian Makran:
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
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 informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More 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 informationNUMERICAL INVESTIGATION OF BUOY- ANCY DRIVEN FLOWS IN TIGHT LATTICE FUEL BUNDLES
Fifth FreeFem workshop on Generic Solver for PDEs: FreeFem++ and its applications NUMERICAL INVESTIGATION OF BUOY- ANCY DRIVEN FLOWS IN TIGHT LATTICE FUEL BUNDLES Paris, December 12 th, 2013 Giuseppe Pitton
More informationNumerical Simulation of Powder Flow
Numerical Simulation of Powder Flow Stefan Turek, Abderrahim Ouazzi Institut für Angewandte Mathematik, Univ. Dortmund http://www.mathematik.uni-dortmund.de/ls3 http://www.featflow.de Models for granular
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
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 informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
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 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 information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
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 informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
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 informationDETAILS ABOUT THE TECHNIQUE. We use a global mantle convection model (Bunge et al., 1997) in conjunction with a
DETAILS ABOUT THE TECHNIQUE We use a global mantle convection model (Bunge et al., 1997) in conjunction with a global model of the lithosphere (Kong and Bird, 1995) to compute plate motions consistent
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationUT-LANL DOE/ASCR ISICLES Project PARADICE: Uncertainty Quantification for Large-Scale Ice Sheet Modeling and Simulation.
UT-LANL DOE/ASCR ISICLES Project PARADICE: Uncertainty Quantification for Large-Scale Ice Sheet Modeling and Simulation Omar Ghattas Representing the rest of the UT-Austin/LANL team: Don Blankenship Tan
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
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 informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationIMPLEMENTATION OF A PARALLEL AMG SOLVER
IMPLEMENTATION OF A PARALLEL AMG SOLVER Tony Saad May 2005 http://tsaad.utsi.edu - tsaad@utsi.edu PLAN INTRODUCTION 2 min. MULTIGRID METHODS.. 3 min. PARALLEL IMPLEMENTATION PARTITIONING. 1 min. RENUMBERING...
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationAn introduction to PDE-constrained optimization
An introduction to PDE-constrained optimization Wolfgang Bangerth Department of Mathematics Texas A&M University 1 Overview Why partial differential equations? Why optimization? Examples of PDE optimization
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More 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 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 informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationA monolithic FEM solver for fluid structure
A monolithic FEM solver for fluid structure interaction p. 1/1 A monolithic FEM solver for fluid structure interaction Stefan Turek, Jaroslav Hron jaroslav.hron@mathematik.uni-dortmund.de Department of
More informationLecture 2: Constitutive Relations
Lecture 2: Constitutive Relations E. J. Hinch 1 Introduction This lecture discusses equations of motion for non-newtonian fluids. Any fluid must satisfy conservation of momentum ρ Du = p + σ + ρg (1) Dt
More informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
More informationUQ Benchmark Problems for Multiphysics Modeling
SIAMUQ UQ Challenge Benchmarks UQ Benchmark Problems for Multiphysics Modeling Maarten Arnst March 31, 2014 SIAMUQ UQ Challenge Benchmarks 1 / 25 Motivation Previous presentation at USNCCM2013 UQ Challenge
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More 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 informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More information4.1 LAWS OF MECHANICS - Review
4.1 LAWS OF MECHANICS - Review Ch4 9 SYSTEM System: Moving Fluid Definitions: System is defined as an arbitrary quantity of mass of fixed identity. Surrounding is everything external to this system. Boundary
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationComputer Fluid Dynamics E181107
Computer Fluid Dynamics E181107 2181106 Transport equations, Navier Stokes equations Remark: foils with black background could be skipped, they are aimed to the more advanced courses Rudolf Žitný, Ústav
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationDynamic Meteorology - Introduction
Dynamic Meteorology - Introduction Atmospheric dynamics the study of atmospheric motions that are associated with weather and climate We will consider the atmosphere to be a continuous fluid medium, or
More informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationDynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet
Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x
More informationThe numerical method used for experiments is based on an explicit finite element
Bialas 1 Model Supplementary Data The numerical method used for experiments is based on an explicit finite element technique similar to the Fast Lagrangian Analysis of Continua (FLAC) method (Cundall,
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationThe semi-geostrophic equations - a model for large-scale atmospheric flows
The semi-geostrophic equations - a model for large-scale atmospheric flows Beatrice Pelloni, University of Reading with M. Cullen (Met Office), D. Gilbert, T. Kuna INI - MFE Dec 2013 Introduction - Motivation
More informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More 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 informationIntroduction to numerical simulation of fluid flows
Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationA monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion
1 / 41 A monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion Chen-Yu CHIANG O. Pironneau, T.W.H. Sheu, M. Thiriet Laboratoire Jacques-Louis Lions
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationThe Generalized Interpolation Material Point Method
Compaction of a foam microstructure The Generalized Interpolation Material Point Method Tungsten Particle Impacting sandstone The Material Point Method (MPM) 1. Lagrangian material points carry all state
More informationSubduction II Fundamentals of Mantle Dynamics
Subduction II Fundamentals of Mantle Dynamics Thorsten W Becker University of Southern California Short course at Universita di Roma TRE April 18 20, 2011 Rheology Elasticity vs. viscous deformation η
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 informationPart of the advantage : Constraint forces do no virtual. work under a set of virtual displacements compatible
FORCES OF CONSTRAINT Lagrangian formalism : Generalized coordinate Minimum set of Eqns Part of the advantage : Constraint forces do no virtual work under a set of virtual displacements compatible with
More informationLecture 3. Properties of Fluids 11/01/2017. There are thermodynamic properties of fluids like:
11/01/2017 Lecture 3 Properties of Fluids There are thermodynamic properties of fluids like: Pressure, p (N/m 2 ) or [ML -1 T -2 ], Density, ρ (kg/m 3 ) or [ML -3 ], Specific weight, γ = ρg (N/m 3 ) or
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 informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationEfficient BDF time discretization of the Navier Stokes equations with VMS LES modeling in a High Performance Computing framework
Introduction VS LES/BDF Implementation Results Conclusions Swiss Numerical Analysis Day 2015, Geneva, Switzerland, 17 April 2015. Efficient BDF time discretization of the Navier Stokes equations with VS
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationAnalysis of a non-isothermal model for nematic liquid crystals
Analysis of a non-isothermal model for nematic liquid crystals E. Rocca Università degli Studi di Milano 25th IFIP TC 7 Conference 2011 - System Modeling and Optimization Berlin, September 12-16, 2011
More information