LECTURE 4: GEOMETRIC PROBLEMS
|
|
- Laura Davidson
- 5 years ago
- Views:
Transcription
1 LECTURE 4: GEOMETRIC PROBLEMS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and Related Topics Isaac Newton Institute, January 2014
2 LECTURE 4: GEOMETRIC PROBLEMS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and Related Topics Isaac Newton Institute, January 2014
3 OUTLINE Shape Differential Calculus Geometric Gradient Flows Example 1: Mean Curvature Flow Example 2: Surface Diffusion Example 3: Optimal Shape Design Example 4: Biomembranes Example 5: Director Fields on Flexible Surfaces Electrowetting with Moving Contact Line
4 Shape Differential Calculus Ω U ad : J(Ω, y(ω )) = inf Ω U ad J(Ω, y(ω)), y(ω) : Ly(Ω) = f on Ω. Goal: Construct a geometric gradient flow of J(Ω) that will lead to a local minimum. This represents a continuous descent direction and is dictated by a geometric PDE. Deformation of Domains (Hadamard, Murat, Delfour, olesio) Given (Lipschitz) vector velocity V Consider mapping X(, t) : x R d X(x, t) R d defined by dx dt = V(X(x, t)), X(x, 0) = x Ω, (Ω := Ω0) Deformed domain Ω t = {X(x, t) : x Ω 0}
5 Shape Derivatives Eulerian derivative of J(Ω) in the direction V: 1 δ ΩJ(Ω; V) = lim J(Ω t 0 t) J(Ω 0) t 1 Material Derivative: ẏ(ω; V) := lim t 0 y(ω t t) X(, t) y(ω 0) Shape Derivative: y (Ω; V) := ẏ(ω; V) {z } y(ω), V {z } material derivative space derivative Notation: V := V ν, := Ω, h := κ 1 + κ 2 (total curvature) Reynolds Theorems: J(Ω) := R φ(, Ω)dx δωj(ω; V) = R Ω Ω φ (Ω; V) + R φv ds J(Ω) := R φ(, Ω)ds δj(ω; V) = R φ (Ω; V) + R νφ + φh V ds
6 Shape Derivatives Eulerian derivative of J(Ω) in the direction V: 1 δ ΩJ(Ω; V) = lim J(Ω t 0 t) J(Ω 0) t 1 Material Derivative: ẏ(ω; V) := lim t 0 y(ω t t) X(, t) y(ω 0) Shape Derivative: y (Ω; V) := ẏ(ω; V) {z } y(ω), V {z } material derivative space derivative Notation: V := V ν, := Ω, h := κ 1 + κ 2 (total curvature) Reynolds Theorems: J(Ω) := R φ(, Ω)dx δωj(ω; V) = R Ω Ω φ (Ω; V) + R φv ds J(Ω) := R φ(, Ω)ds δj(ω; V) = R φ (Ω; V) + R νφ + φh V ds
7 Shape Derivatives Eulerian derivative of J(Ω) in the direction V: 1 δ ΩJ(Ω; V) = lim J(Ω t 0 t) J(Ω 0) t 1 Material Derivative: ẏ(ω; V) := lim t 0 y(ω t t) X(, t) y(ω 0) Shape Derivative: y (Ω; V) := ẏ(ω; V) {z } y(ω), V {z } material derivative space derivative Notation: V := V ν, := Ω, h := κ 1 + κ 2 (total curvature) Reynolds Theorems: J(Ω) := R φ(, Ω)dx δωj(ω; V) = R Ω Ω φ (Ω; V) + R φv ds J(Ω) := R φ(, Ω)ds δj(ω; V) = R φ (Ω; V) + R νφ + φh V ds
8 Examples of Shape Derivatives Example 1: Volume J(Ω) = Example 2: Area J() = Ω 1 δj(ω; V ) = 1 δj(; V ) = Example 3: Bending Energy 1 J()) := 2 h2 ds δj(; V ) = because shape calculus gives V dσ HV dσ h V 1 2 h3 V + 2hκV ν (, V) = V, h (, V) = V, νh = ν 2 = (h 2 2κ) whence δ W (, φ) = = hh + (h νh h3 )V V h h ν 2 V h 3 V.
9 Geometric Gradient Flows Derivative of J(Ω) in the direction V satisfies: δj(ω; V ) = R GV for explicit functions g, f. G(x, Ω) = g(x, Ω)H + f(x, Ω) Scalar product B(, ) on and corresponding norm. Gradient flow: B(V, W ) = R GW δj(ω; V ) = B(V, V ) V 2 < 0 Energy decreasing flow! Gradient flow in operator form: BV = G
10 Implicit Time Discretizaion Implicit Euler: X n+1 = X n + τ V n+1 (Ω n Ω n+1 ) Minimization: argmin Ω n+1 J(Ω n+1 ) + 1 2τ d2 (Ω n+1, Ω n ) Distance between Domains: d 2 (Ω n+1, Ω n ) = τv n+1 2 n+1 = τ 2 V n+1 2 n+1 = τ 2 B(V n+1, V n+1 ) Implicit Minimization: argmin W n+1 J(Ω n+1 ) + τ2 W n+1 2 n+1 Optimality Condition: (Euler-Lagrange equation) B n+1 (V n+1, W ) = δj(ω n+1 ; W ) = G n+1 W, W n+1 Energy Decrease Property: Ω n+1 = Ω n + τv n+1 J(Ω n+1 ) J(Ω n + τv n+1 ) + τ 2 V n+1 2 J(Ω n ).
11 Implicit Time Discretizaion Implicit Euler: X n+1 = X n + τ V n+1 (Ω n Ω n+1 ) Minimization: argmin Ω n+1 J(Ω n+1 ) + 1 2τ d2 (Ω n+1, Ω n ) Distance between Domains: d 2 (Ω n+1, Ω n ) = τv n+1 2 n+1 = τ 2 V n+1 2 n+1 = τ 2 B(V n+1, V n+1 ) Implicit Minimization: argmin W n+1 J(Ω n+1 ) + τ2 W n+1 2 n+1 Optimality Condition: (Euler-Lagrange equation) B n+1 (V n+1, W ) = δj(ω n+1 ; W ) = G n+1 W, W n+1 Energy Decrease Property: Ω n+1 = Ω n + τv n+1 J(Ω n+1 ) J(Ω n + τv n+1 ) + τ 2 V n+1 2 J(Ω n ).
12 Linearization: Semi-Implicit Time Discretization Treat geometry explicitly: replacement in Optimality Condition Ω n+1 Ω n, ν n+1 ν n, G n+1 Ĝn Ĝn := g(x, Ω n )H n+1 + f(x, Ω n ) [Dziuk], [Bänsch, Morin, N], [Barrett, Garcke, Nürnberg], [Dziuk, Elliott] The Riesz representative Ĝn of shape gradient handles curvature implicitly. Gradient flow equation: B n+1 (V n+1, W ) = R G n+1 W, W n+1 B n (V n+1, W ) = Ĝ n W, W n B n V n+1 + g n H n+1 = f n
13 Operator Splitting Issue: How to compute H, H Differential Geometry Identity [G. Dziuk]: X = H = H ν KEY IDEA: write the problem in terms of scalar and vector quantities H, H, V, V H = H ν, V = V ν X n X n+1 = X n + τ n V n+1 H n+1 + τ n n V n+1 = n X n H n+1 H n+1 ν n = 0 B n V n+1 + g n H n+1 = f n V n+1 V n+1 ν n = 0 (Operator Splitting)
14 Operator Splitting Issue: How to compute H, H Differential Geometry Identity [G. Dziuk]: X = H = H ν KEY IDEA: write the problem in terms of scalar and vector quantities H, H, V, V H = H ν, V = V ν X n X n+1 = X n + τ n V n+1 H n+1 + τ n n V n+1 = n X n H n+1 H n+1 ν n = 0 B n V n+1 + g n H n+1 = f n V n+1 V n+1 ν n = 0 (Operator Splitting)
15 Operator Splitting Issue: How to compute H, H Differential Geometry Identity [G. Dziuk]: X = H = H ν KEY IDEA: write the problem in terms of scalar and vector quantities H, H, V, V H = H ν, V = V ν X n X n+1 = X n + τ n V n+1 H n+1 + τ n n V n+1 = n X n H n+1 H n+1 ν n = 0 B n V n+1 + g n H n+1 = f n V n+1 V n+1 ν n = 0 (Operator Splitting)
16 Semi-Implicit FEM Geometry: ν = ν n h, = n h Discrete spaces (C 0 piecewise linears): V h () H 1 (), Vh () H 1 () Fully discrete problem: seek V n+1, H n+1 V h (), V n+1, H n+1 V h () such that D H n+1, φ E E E h τ D V n+1, φh = D X n, φh φ h V h () D E H n+1, φ h H n+1 ν, φ h = 0 φ h V h () B n (V n+1, φ h ) + g n H n+1, φ h = f n, φ h φ h V h () D V n+1, φ E h DV n+1, ν φ E h = 0 φ h V h ()
17 Mean Curvature Flow: Parametric Approach (G. Dziuk) L 2 - Gradient flow: J() = R V W = δj(; W ) V = H Energy law: R V 2 + d dt = 0 Vector form: V = V ν = H ν = H = X t X X = 0 (vector heat equation on ) Semi-implicit time discretization: X n+1 τ n n X n+1 = X n
18 Mean Curvature Flow: Level Set Approach (Osher, Sethian) Levet Set Function φ: Unit normal, mean curvature and normal velocity read ν = φ tφ, H = div ν, V = φ φ Motion by Mean Curvature: V = H tφ φ div φ φ tφ nx i,j=1 δ ij x i φ x φ j φ 2 xi x j φ = 0 Viscosity Solutions: [Evans, Spruck], [Chen, Giga, Goto] Numerical Methods: [Osher-Sethian], [Walkington], [Chambolle], etc Error estimate for Regularized Flow: [Deckelnick], [Deckelnick, Dziuk, Elliott].
19 Mean Curvature Flow: Phase Field Approach (Allen-Cahn) Allen-Cahn Equation for the Phase Variable ϕ: ɛ tϕ ɛ ɛ ϕ ɛ + 1 ɛ f(ϕɛ) = 0 f(ϕ) = ϕ(ϕ2 1) Interface Convergence: Error bounds of order ɛ 2 log ɛ 2 for Haussdorff distance between zero level set of ϕ ɛ and surface moving by mean curvature in the smooth regime [Bellettini, Paolini] Viscosity Solutions: Convergence of ϕ ɛ to viscosity solution to MCF pass singularities [Evans, Soner, Souganidis]. Double Obstacle Phase Field Model: The nonlinearity f becomes f(ϕ) = ϕ + I [,1,1] (ϕ). Interface Convergence: Estimates of order ɛ 2 for the smooth regime [N, Paolini, Verdi], and of order ɛ pass singularities [N, Verdi]. Techniques: Comparison principles, barrier arguments, non-degeneracy.
20 Surface Diffusion (w. E. Bänsch and P. Morin) Epitaxially stressed solid due to mismatch with the rigid substrate. Simplified model: solid governed by the Laplace equation Dynamics of free surface (t): Rearrangement of particles on (t) so as to minimize surface tension (plastic deformation) V = `κ + ε S Ω S where ε = y 2 ( elastic energy ), and 8 y = 0 in Ω >< y = x on D y = x + periodic on S >: νy = 0 on D
21 Surface Diffusion (Continued) Domain functional: J(Ω) := y 2 + γ ds Ω {z } {z 1 } elastic energy surface energy y = 0 Ω, νy = 0 y = g D (fixed). Shape derivative y = y (Ω; V ) satisfies y = 0 Ω, νy = div (V y), y = 0 D. Shape derivative of J(Ω) in the direction V satisfies dj(ω; V ) = R y 2 + γh V ds {z } =G=f+gH
22 Surface Diffusion (Continued) Domain functional: J(Ω) := y 2 + γ ds Ω {z } {z 1 } elastic energy surface energy y = 0 Ω, νy = 0 y = g D (fixed). Shape derivative y = y (Ω; V ) satisfies y = 0 Ω, νy = div (V y), y = 0 D. Shape derivative of J(Ω) in the direction V satisfies dj(ω; V ) = R y 2 + γh V ds {z } =G=f+gH
23 Optimal Shape Design (w. G. Dogan, M. Morin, and M. Verani) Domain functional: J(Ω) = 1 2 D `y(ω) zg 2 y = 0 Ω, νy = 1 1, y = 0 2 (fixed). Shape derivative y = y (Ω; V ) satisfies y = 0 Ω, νy = div (V y) + HV 1, y = 0 2. Dual (or adjoint) solution p accounts for L 2 -norm in J(Ω) and satisfies p = χ D(y z g) Ω, νp = 0 1, p = 0 2. Derivative of J(Ω) in the direction V satisfies dj(ω; V ) = R 1 y p + Hp V ds {z } =G=f+Hg
24 Optimal Shape Design (w. G. Dogan, M. Morin, and M. Verani) Domain functional: J(Ω) = 1 2 D `y(ω) zg 2 y = 0 Ω, νy = 1 1, y = 0 2 (fixed). Shape derivative y = y (Ω; V ) satisfies y = 0 Ω, νy = div (V y) + HV 1, y = 0 2. Dual (or adjoint) solution p accounts for L 2 -norm in J(Ω) and satisfies p = χ D(y z g) Ω, νp = 0 1, p = 0 2. Derivative of J(Ω) in the direction V satisfies dj(ω; V ) = R 1 y p + Hp V ds {z } =G=f+Hg
25 L 2 - Gradient Flow k = 3 J = k = 6 J = k = 9 J = k = 22 J = L 2 -gradient flow from a non-centered ellipse to a centered circle. The exact solution is a circle of radius one centered at the origin, and the initial configuration is a small ellipse centered at (0.5, 0.7). The L 2 -dynamics is locally stable where p > 0 (upper right) and locally unstable where p < 0 (lower left). The bottom row is a zoom of the unstable region exhibiting oscillations at the mesh level.
26 (Stable) Weighted H 1 - Gradient Flow k = 0 J = k = 20 J = k = 40 J = k = 60 J = k = 80 J = k = 100 J = k = 120 J = k = 140 J = Stable dynamics: Weighted H 1 flow from a non-centered ellipse to a centered circle.
27 Biomembranes (w. A. Bonito and M.S. Pauletti) Bending energy: W () = 1 H 2 + λ 1 + p 1 2 Ω δ W = H 1 «2 H3 + 2KH ν + λhν + pν Equivalent Form of Derivative of Willmore Energy δ W : δ W (, V ) = V H H ν 2 V H 3 V. Vector Form of Derivative of Willmore Energy δ W : V H = V : H ( x+ x T ) V : H δ W (; V) = H ν V = H ν 2 1 H 2 V 2 V : H ( x+ x T ) V : H+ 1 2 Related methods: [Rusu], [Dziuk], [Barrett, Garcke, Nürnberg] H ν V H V
28 Biomembranes (w. A. Bonito and M.S. Pauletti) Bending energy: J() = 1 H 2 + λ 1 + p 1 2 Ω δ J = H 1 «2 H3 + 2KH ν + λhν + pν Equivalent Form of Derivative of Willmore Energy δ W : δ W (, V ) = V H H ν 2 V H 3 V. Vector Form of Derivative of Willmore Energy δ W : V H = V : H ( x+ x T ) V : H δ W (; V) = H ν V = H ν 2 1 H 2 V 2 V : H ( x+ x T ) V : H+ 1 2 Related methods: [Rusu], [Dziuk], [Barrett, Garcke, Nürnberg] H ν V H V
29 Biomembranes (w. A. Bonito and M.S. Pauletti) Bending energy: J() = 1 H 2 + λ 1 + p 1 2 Ω δ J = H 1 «2 H3 + 2KH ν + λhν + pν Equivalent Form of Derivative of Willmore Energy δ W : δ W (, V ) = V H H ν 2 V H 3 V. Vector Form of Derivative of Willmore Energy δ W : V H = V : H ( x+ x T ) V : H δ W (; V) = H ν V = H ν 2 1 H 2 V 2 V : H ( x+ x T ) V : H+ 1 2 Related methods: [Rusu], [Dziuk], [Barrett, Garcke, Nürnberg] H ν V H V
30 Biomembranes (w. A. Bonito and M.S. Pauletti) Bending energy: J() = 1 H 2 + λ 1 + p 1 2 Ω δ J = H 1 «2 H3 + 2KH ν + λhν + pν Equivalent Form of Derivative of Willmore Energy δ W : δ W (, V ) = V H H ν 2 V H 3 V. Vector Form of Derivative of Willmore Energy δ W : V H = V : H ( x+ x T ) V : H δ W (; V) = H ν V = H ν 2 1 H 2 V 2 V : H ( x+ x T ) V : H+ 1 2 Related methods: [Rusu], [Dziuk], [Barrett, Garcke, Nürnberg] H ν V H V
31 Biomembranes (w. A. Bonito and M.S. Pauletti) Bending energy: J() = 1 H 2 + λ 1 + p 1 2 Ω δ J = H 1 «2 H3 + 2KH ν + λhν + pν Equivalent Form of Derivative of Willmore Energy δ W : δ W (, V ) = V H H ν 2 V H 3 V. Vector Form of Derivative of Willmore Energy δ W : V H = V : H ( x+ x T ) V : H δ W (; V) = H ν V = H ν 2 1 H 2 V 2 V : H ( x+ x T ) V : H+ 1 2 Related methods: [Rusu], [Dziuk], [Barrett, Garcke, Nürnberg] H ν V H V
32 Fluid-Membrane Interaction (w. A. Bonito and M.S. Pauletti) Refs: Coutand-Shkoller (Incompressible) Navier-Stokes Equations Ω ρ v div ( pi + µd(v) ) = b in Ω t, {z } Σ div v = 0 in Ω t, [Σ]ν = δ J on t, v = ϑ on t, v(, 0) = v 0 in Ω 0. D Ω in g Ω out Membrane Force: Bending δ J = H + 1 «2 H3 2KH ν + λhν.
33 Fluid-Membrane Interaction (w. A. Bonito and M.S. Pauletti) Refs: Coutand-Shkoller (Incompressible) Navier-Stokes Equations Ω ρ v div ( pi + µd(v) ) = b in Ω t, {z } Σ div v = 0 in Ω t, [Σ]ν = δ J on t, v = ϑ on t, v(, 0) = v 0 in Ω 0. D Ω in g Ω out Membrane Force: Bending δ J = H + 1 «2 H3 2KH ν + λhν.
34 Fluid-Membrane Interaction (w. A. Bonito and M.S. Pauletti) Refs: Coutand-Shkoller (Incompressible) Navier-Stokes Equations Ω ρ v div ( pi + µd(v) ) = b in Ω t, {z } Σ div v = 0 in Ω t, [Σ]ν = δ J on t, v = ϑ on t, v(, 0) = v 0 in Ω 0. D Ω in g Ω out Membrane Force: Bending δ J = H + 1 «2 H3 2KH ν + λhν.
35 Geometric vs Fluid Red Cell: 5x5x1 Ellipsoid play
36 Fluid Red Blood Cell: Ellipsoid 5x5x1 Streamlines
37 Comparison of Geometric with Fluid Red Cell: Physical Parameters Comparison: Physical Parameters
38 Prototype Model and Applications Basic setting: Surface R 3 and director field n : S 2 n Physics: is the surface constituted by some molecules and n is the orientation of those or other molecules Question: How do shape and orientational order interact? Basic Model:, n minimize Helfrich and Frank energy with coupling F (, n) E[, n] = 1 H n 2 + F (, n) 2 2 for in fixed topology class and n a unit-length vector-field (H = div ν).
39 Prototype Model and Applications Basic setting: Surface R 3 and director field n : S 2 n Physics: is the surface constituted by some molecules and n is the orientation of those or other molecules Question: How do shape and orientational order interact? Basic Model:, n minimize Helfrich and Frank energy with coupling F (, n) E[, n] = 1 H n 2 + F (, n) 2 2 for in fixed topology class and n a unit-length vector-field (H = div ν).
40 Nonlinear Model: First Variation of the Energy E(, n) E(, n) := 1 (div ν δdiv n) 2 + λ n 2 + µ` n f(n ν) 2 2 2ε 2 Extend n constant along the normal ν, which makes sense for biomembranes. The first variation of E(, n) with respect to is (H = div ν) D δe E δ, φ = H, φ + H ν 2, φ δ H( n : ν), φ + δ H(ν T n), φ δ (div n), φ + δ div n ν 2, φ δ 2 div n( n : ν), φ + δ 2 div n(ν T n), φ 1 H(H + δdiv n) 2, φ λ ( n) T : Dν( n), φ 2 λ H n 2, φ 1 f (n ν)n, φ 1 Hf(n ν), φ, 2 ε ε for all φ C (). The variation of E(, n) with respect to n is D δe E δn, φ = δ H+δdiv n, div m +λ n, m + 2µn, m + 1 f (n ν), m ν, ε for all m C (; R 3 ). Expression simplifies if m is perpendicular to n.
41 Nonlinear Model: First Variation of the Energy E(, n) E(, n) := 1 (div ν δdiv n) 2 + λ n 2 + µ` n f(n ν) 2 2 2ε 2 Extend n constant along the normal ν, which makes sense for biomembranes. The first variation of E(, n) with respect to is (H = div ν) D δe E δ, φ = H, φ + H ν 2, φ δ H( n : ν), φ + δ H(ν T n), φ δ (div n), φ + δ div n ν 2, φ δ 2 div n( n : ν), φ + δ 2 div n(ν T n), φ 1 H(H + δdiv n) 2, φ λ ( n) T : Dν( n), φ 2 λ H n 2, φ 1 f (n ν)n, φ 1 Hf(n ν), φ, 2 ε ε for all φ C (). The variation of E(, n) with respect to n is D δe E δn, φ = δ H+δdiv n, div m +λ n, m + 2µn, m + 1 f (n ν), m ν, ε for all m C (; R 3 ). Expression simplifies if m is perpendicular to n.
42 Relaxation Dynamics: L 2 -Gradient Flow If v denotes the normal velocity of, we have to solve the system of PDE on D δe E v, φ = δ, φ for all φ C (), D δe E tn, m = δn, m for all m C (; T ns 2 ), subject to the constraint that n(t, x) S 2 for almost every (t, x). Semi-implicit Time Discretization v j 1 `Xj X j 1 ν j 1 τ with ν j 1 : j 1 S 2 the outer unit normal to j 1, and j 1X j = H j ν j 1 Refs: Dziuk, Dziuk-Elliott, Bänsch-Morin-Nochetto, Barrett-Garcke-Nürnberg.
43 Relaxation Dynamics: L 2 -Gradient Flow If v denotes the normal velocity of, we have to solve the system of PDE on D δe E v, φ = δ, φ for all φ C (), D δe E tn, m = δn, m for all m C (; T ns 2 ), subject to the constraint that n(t, x) S 2 for almost every (t, x). Semi-implicit Time Discretization v j 1 `Xj X j 1 ν j 1 τ with ν j 1 : j 1 S 2 the outer unit normal to j 1, and j 1X j = H j ν j 1 Refs: Dziuk, Dziuk-Elliott, Bänsch-Morin-Nochetto, Barrett-Garcke-Nürnberg.
44 Simulations for Surfactants Evolution from a prolate ellipsoid to a discocyte, with (δ = 1) and without (δ = 0) surfactant, and with area and volume constraints. Snapshots of the evolution after n = 100, 200, 400 and 1200 time-steps. The upper plots show the evolution in the presence of surfactants (δ = 1) while the second row shows the Helfrich-flow (δ = 0). The coupling of the surface and the director field decelerates the evolution and leads to a less pronounced biconcave shape.
45 Simulations for Biomembranes: Two Positive Degree-One Defects Two degree +1 defects pointing outwards at the poles develop a cone-like structure near poles (stomatocyte shape) 2 Defects +1
46 Simulations for Biomembranes: Two Positive Degree-One Defects Two inward and outward pointing at the poles of degree +1 develop a heart-like shape (echinocytes) 2 Defects +1
47 Simulations for Biomembranes: Two Negative Degree-One Defects Biomembrane develops saddle-shapes at the poles, where the negative degree-one defects are located (top row) and cone shapes in neighborhoods of the four positive degree-one defects located on the equator (see the cut through the equator in the lower left picture). The surface and the director field are colored by div n.
48 Phase Field Model for Electrowetting Two-Phase Flow: Phase variable φ satisfies a Cahn-Hilliard type equation tφ + u φ = div `M(φ) µ 1 µ = γ δ f(φ) δ φ 1 2 ɛ (φ) V ρ (φ) u 2 where u is fluid velocity, M(φ) mobility, µ chemical potential, f(φ) cubic nonlinearlity, ɛ(φ) electric permitivity, V electric potential, ρ(φ) density. Variable Density Navier-Stokes: Fluid velocity and pressure (u, p) satisfy D t`ρ(φ)u div `η(φ)σ(u) + p = µ φ q (V + λq) div u = 0, where q is the electric charge. Navier boundary condition, electrostatic equation for V, transport equation for q, and boundary conditions. Formulation and Analysis: [N, Salgado, Walker] [Grün et al.], [Garcke et al.]
49 Electrowetting Simulation: Droplet Splitting 2d Configuration: Some Parameters: ρ 1 = 1, ρ 2 = 100, δ = , γ = 50, λ = 0.5, M = 10 2, V = 20, η 1 = 1, η 2 = 10,
CURVATURE DRIVEN FLOWS IN DEFORMABLE DOMAINS
CURVATURE DRIVEN FLOWS IN DEFORMABLE DOMAINS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA BCAM EHU/UPV Basque Colloquium in Mathematics and its
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationA sharp diffuse interface tracking method for approximating evolving interfaces
A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex Overview Introduction Phase field models Double well and double obstacle
More informationDynamics of Biomembranes: Effect of the Bulk Fluid
Dynamics of Biomembranes: Effect of the Bulk Fluid A. Bonito 1, R.H. Nochetto 2, and M.S. Pauletti 1 1 Department of Mathematics, Texas A&M University, College Station, Texas, USA 2 Department of Mathematics,
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 informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationA finite element level set method for anisotropic mean curvature flow with space dependent weight
A finite element level set method for anisotropic mean curvature flow with space dependent weight Klaus Deckelnick and Gerhard Dziuk Centre for Mathematical Analysis and Its Applications, School of Mathematical
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationConvergence Rates for AFEM: PDE on Parametric Surfaces
Department of Mathematics and Institute for Physical Science and Technology University of Maryland Joint work with A. Bonito and R. DeVore (Texas a&m) J.M. Cascón (Salamanca, Spain) K. Mekchay (Chulalongkorn,
More informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
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 informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationFinite Element Methods for Shape Optimization and Applications
Finite Element Methods for Shape Optimization and Applications Gunay Doǧan a, Pedro Morin b, Ricardo H. Nochetto a and Marco Verani c a Department of Mathematics, University of Maryland, College Park,
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 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 informationApplication of Elastic Flows to Cell Motility. Björn Stinner
Mathematics Institute Centre for Scientific Computing Application of Elastic Flows to Cell Motility Björn Stinner (joint work with C Venkataraman and C Elliott) FBP - Theory and Applications 2012 Cell
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationNumerics for Liquid Crystals with Variable Degree of Orientation
Numerics for Liquid Crystals with Variable Degree of Orientation Ricardo H. Nochetto, Shawn W. Walker 2 and Wujun hang Department of Mathematics, University of Maryland 2 Department of Mathematics, Louisiana
More informationIsogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows
Isogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows Babak S. Hosseini 1, Stefan Turek 1, Matthias Möller 2 1 TU Dortmund University, Institute
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More 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 informationUniversität Regensburg Mathematik
Universität Regensburg Mathematik Parametric approximation of Willmore flow and related geometric evolution equations John W. Barrett, Harald Garcke and Robert Nürnberg Preprint Nr. 22/2007 Parametric
More informationOn the Parametric Finite Element Approximation of Evolving Hypersurfaces in R 3
On the Parametric Finite Element Approximation of Evolving Hypersurfaces in R 3 John W. Barrett Harald Garcke Robert Nürnberg Abstract We present a variational formulation of motion by minus the Laplacian
More informationSurface phase separation and flow in a simple model of drops and vesicles
Surface phase separation and flow in a simple model of drops and vesicles Tutorial Lecture 4 John Lowengrub Department of Mathematics University of California at Irvine Joint with J.-J. Xu (UCI), S. Li
More informationNumerical Solutions of Geometric Partial Differential Equations. Adam Oberman McGill University
Numerical Solutions of Geometric Partial Differential Equations Adam Oberman McGill University Sample Equations and Schemes Fully Nonlinear Pucci Equation! "#+ "#* "#) "#( "#$ "#' "#& "#% "#! "!!!"#$ "
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationA level set projection model of lipid vesicles in general flows
, A level set projection model of lipid vesicles in general flows D. Salac a, M. Miksis a a Northwestern University, Engineering Sciences and Applied Mathematics, Evanston, IL, 60208 Abstract A new numerical
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationexamples of equations: what and why intrinsic view, physical origin, probability, geometry
Lecture 1 Introduction examples of equations: what and why intrinsic view, physical origin, probability, geometry Intrinsic/abstract F ( x, Du, D u, D 3 u, = 0 Recall algebraic equations such as linear
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 informationIntrinsic finite element modeling of a linear membrane shell problem
RR Intrinsic finite element modeling of a linear membrane shell problem PETER HANSBO AND MATS G. LARSON School of Engineering Jönköping University Research Report No. : ISSN -8 Intrinsic finite element
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
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 informationMinimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation.
Lecture 7 Minimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation div + u = ϕ on ) = 0 in The solution is a critical point or the minimizer
More informationFinite-time singularity formation for Euler vortex sheet
Finite-time singularity formation for Euler vortex sheet Daniel Coutand Maxwell Institute Heriot-Watt University Oxbridge PDE conference, 20-21 March 2017 Free surface Euler equations Picture n N x Ω Γ=
More informationLecture III. Curvature Elasticity of Fluid Lipid Membranes. Helfrich s Approach
Lecture III. Curvature Elasticity of Fluid Lipid Membranes. Helfrich s Approach Marina Voinova Department of Applied Physics, Chalmers University of Technology and Göteborg University, SE-412 96, Göteborg,
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationA MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION
1 A MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique The most recent results were obtained in collaboration with F. de Gournay, F. Jouve, O. Pantz, A.-M. Toader.
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 informationTHE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS
THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS Alain Miranville Université de Poitiers, France Collaborators : L. Cherfils, G. Gilardi, G.R. Goldstein, G. Schimperna,
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 informationNonsmooth Schur Newton Methods and Applications
Nonsmooth Schur Newton Methods and Applications C. Gräser, R. Kornhuber, and U. Sack IMA Annual Program Year Workshop Numerical Solutions of Partial Differential Equations: Fast Solution Techniques November
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationAFEM for Geometric PDE: The Laplace-Beltrami Operator
AFEM for Geometric PDE: The Laplace-Beltrami Operator In memory of Enrico Magenes Andrea Bonito J. Manuel Cascón Pedro Morin Ricardo H. Nochetto December 17, 2012 Abstract We present several applications
More informationPast, present and space-time
Past, present and space-time Arnold Reusken Chair for Numerical Mathematics RWTH Aachen Utrecht, 12.11.2015 Reusken (RWTH Aachen) Past, present and space-time Utrecht, 12.11.2015 1 / 20 Outline Past. Past
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationSingular Diffusion Equations With Nonuniform Driving Force. Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009
Singular Diffusion Equations With Nonuniform Driving Force Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009 1 Contents 0. Introduction 1. Typical Problems 2. Variational Characterization
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationContinuum Methods 1. John Lowengrub Dept Math UC Irvine
Continuum Methods 1 John Lowengrub Dept Math UC Irvine Membranes, Cells, Tissues complex micro-structured soft matter Micro scale Cell: ~10 micron Sub cell (genes, large proteins): nanometer Macro scale
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationActive Contours Snakes and Level Set Method
Active Contours Snakes and Level Set Method Niels Chr. Overgaard 2010-09-21 N. Chr. Overgaard Lecture 7 2010-09-21 1 / 26 Outline What is an mage? What is an Active Contour? Deforming Force: The Edge Map
More informationUniversität Regensburg Mathematik
Universität Regensburg Mathematik On the stable discretization of strongly anisotropic phase field models with applications to crystal growth John W. Barrett, Harald Garcke and Robert Nürnberg Preprint
More informationLagrangian acceleration in confined 2d turbulent flow
Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole
More informationIntroduction to immersed boundary method
Introduction to immersed boundary method Ming-Chih Lai mclai@math.nctu.edu.tw Department of Applied Mathematics Center of Mathematical Modeling and Scientific Computing National Chiao Tung University 1001,
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 IB Thursday 7 June 2007 9 to 12 PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question in Section
More informationRisk Averse Shape Optimization
Risk Averse Shape Optimization Sergio Conti 2 Martin Pach 1 Martin Rumpf 2 Rüdiger Schultz 1 1 Department of Mathematics University Duisburg-Essen 2 Rheinische Friedrich-Wilhelms-Universität Bonn Workshop
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
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 informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
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 informationA posteriori error control for the binary Mumford Shah model
A posteriori error control for the binary Mumford Shah model Benjamin Berkels 1, Alexander Effland 2, Martin Rumpf 2 1 AICES Graduate School, RWTH Aachen University, Germany 2 Institute for Numerical Simulation,
More informationPropagation d interfaces avec termes non locaux
Propagation d interfaces avec termes non locaux P. Cardaliaguet Univ. Brest Janvier 2008 Joint works with G. Barles (Tours), O. Alvarez (Rouen), O. Ley (Tours), R. Monneau (CERMICS), A. Monteillet (Brest).
More informationMotivation Power curvature flow Large exponent limit Analogues & applications. Qing Liu. Fukuoka University. Joint work with Prof.
On Large Exponent Behavior of Power Curvature Flow Arising in Image Processing Qing Liu Fukuoka University Joint work with Prof. Naoki Yamada Mathematics and Phenomena in Miyazaki 2017 University of Miyazaki
More informationQUADRATIC RATE OF CONVERGENCE FOR CURVATURE DEPENDENT SMOOTH INTERFACES: A SIMPLE PROOF 1
QUADRATIC RATE OF CONVERGENCE FOR CURVATURE DEPENDENT SMOOTH INTERFACES: A SIMPLE PROOF R. H. NOCHETTO Department of Mathematics and Institute for Physical Science and Technology University of Maryland,
More informationA VOLUME MESH FINITE ELEMENT METHOD FOR PDES ON SURFACES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2012) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 10-14, 2012 A VOLUME MESH FINIE ELEMEN MEHOD FOR
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More 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 informationMathematical Theory of Non-Newtonian Fluid
Mathematical Theory of Non-Newtonian Fluid 1. Derivation of the Incompressible Fluid Dynamics 2. Existence of Non-Newtonian Flow and its Dynamics 3. Existence in the Domain with Boundary Hyeong Ohk Bae
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas Rencontres
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 informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
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 informationPriority Programme The Combination of POD Model Reduction with Adaptive Finite Element Methods in the Context of Phase Field Models
Priority Programme 1962 The Combination of POD Model Reduction with Adaptive Finite Element Methods in the Context of Phase Field Models Carmen Gräßle, Michael Hinze Non-smooth and Complementarity-based
More informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More informationShape Optimization for Drag Minimization Using the Navier-Stokes Equation
Louisiana State University LSU Digital Commons LSU Master's Theses Graduate School 2015 Shape Optimization for Drag Minimization Using the Navier-Stokes Equation Chukwudi Paul Chukwudozie Louisiana State
More information2 phase problem for the Navier-Stokes equations in the whole space.
2 phase problem for the Navier-Stokes equations in the whole space Yoshihiro Shibata Department of Mathematics & Research Institute of Sciences and Engineerings, Waseda University, Japan Department of
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 informationSmoluchowski Navier-Stokes Systems
Smoluchowski Navier-Stokes Systems Peter Constantin Mathematics, U. of Chicago CSCAMM, April 18, 2007 Outline: 1. Navier-Stokes 2. Onsager and Smoluchowski 3. Coupled System Fluid: Navier Stokes Equation
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
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 informationExercise 1 (Formula for connection 1-forms) Using the first structure equation, show that
1 Stokes s Theorem Let D R 2 be a connected compact smooth domain, so that D is a smooth embedded circle. Given a smooth function f : D R, define fdx dy fdxdy, D where the left-hand side is the integral
More informationOn the Willmore Functional and Applications
On the Willmore Functional and Applications Yann Bernard Monash University La Trobe University January 26, 2017 Some History ca. 1740: Leonhard Euler and Daniel Bernoulli study the 1-dim. elastica ca.
More informationBuilding Solutions to Nonlinear Elliptic and Parabolic Partial Differential Equations
Building Solutions to Nonlinear Elliptic and Parabolic Partial Differential Equations Adam Oberman University of Texas, Austin http://www.math.utexas.edu/~oberman Fields Institute Colloquium January 21,
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationMATH 228: Calculus III (FALL 2016) Sample Problems for FINAL EXAM SOLUTIONS
MATH 228: Calculus III (FALL 216) Sample Problems for FINAL EXAM SOLUTIONS MATH 228 Page 2 Problem 1. (2pts) Evaluate the line integral C xy dx + (x + y) dy along the parabola y x2 from ( 1, 1) to (2,
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationFIRST VARIATION OF THE GENERAL CURVATURE-DEPENDENT SURFACE ENERGY
FIRST VARIATION OF THE GENERAL CURVATURE-DEPENDENT SURFACE ENERGY G. DOĞAN1 AND R.H. NOCHETTO 2 Abstract. We consider general surface energies, which are weighted integrals over a closed surface with a
More informationENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More information