Finite Elements for Magnetohydrodynamics and its Optimal Control
|
|
- Lawrence Flowers
- 5 years ago
- Views:
Transcription
1 Finite Elements for Magnetohydrodynamics and its Karl Kunisch Marco Discacciati (RICAM) FEM Symposium Chemnitz September 25 27, 2006
2 Overview 1 2 3
3 What is Magnetohydrodynamics? Magnetohydrodynamics (MHD) concerns the mutual interaction of electrically conducting fluids and magnetic fields
4 What is Magnetohydrodynamics? Magnetohydrodynamics (MHD) concerns the mutual interaction of electrically conducting fluids and magnetic fields Use taylored magnetic fields for... stirring of conducting fluids flow damping (during casting or solidification) electromagnetic filtration, melting, levitation
5 Application: Casting of Aluminum Illustration: B.Q. Li
6 Application: Casting of Aluminum Features convection flow due to temperature gradients undesired inflow of impurities idea: damping by magnetic fields Illustration: B.Q. Li
7 Application: Production of Aluminum Solidified bath Anods Electrolytic bath Interface Liquid aluminum Cathod Illustration: J.-F. Gerbeau
8 Application: Production of Aluminum Solidified bath Anods Electrolytic bath Interface Liquid aluminum Cathod Features two fluids, free surface, free interface electrolytic bath shallow: huge energy savings electrolytic bath deep: damping of instabilities (stray magnetic fields) Illustration: J.-F. Gerbeau
9 Application: CZ Crystal Growth Illustration: B.Q. Li
10 Application: CZ Crystal Growth Features convection-driven flow free surface, Marangoni effect, non-local radiation idea: damping or stirring by magnetic fields Illustration: B.Q. Li
11 Application: FZ Crystal Growth Illustration: B.Q. Li
12 Application: FZ Crystal Growth Features two free interfaces free surface no mechanical support for melt phase idea: heat, confine and shape the melt phase by magnetic fields Illustration: B.Q. Li
13 Summary: Applications of MHD Numerous applications in... metallurgy crystal growth
14 Summary: Applications of MHD Numerous applications in... metallurgy crystal growth Most attractive features contactless application of a volume force (Lorentz force) induction heating
15 Interaction in MHD Interaction principles charge carries moving in magnetic field induce currents
16 Interaction in MHD Interaction principles charge carries moving in magnetic field induce currents currents induce magnetic fields
17 Interaction in MHD Interaction principles charge carries moving in magnetic field induce currents currents induce magnetic fields magnetic fields exert a Lorentz force on moving charge carries
18 MHD Equations: The Stationary Case Navier-Stokes system with Lorentz force ϱ (u )u η u + p = J B u = 0 fluid velocity u on Ω pressure p on Ω current density J on Ω magnetic field B on R 3
19 MHD Equations: The Stationary Case Navier-Stokes system with Lorentz force ϱ (u )u η u + p = J B u = 0 fluid velocity u on Ω pressure p on Ω Charge conservation and Ohm s law σ 1 J + φ = u B J = 0 current density J on Ω electric potential φ on Ω magnetic field B on R 3
20 MHD Equations: The Stationary Case Navier-Stokes system with Lorentz force ϱ (u )u η u + p = J B u = 0 fluid velocity u on Ω pressure p on Ω Charge conservation and Ohm s law σ 1 J + φ = u B J = 0 current density J on Ω electric potential φ on Ω No monopoles and Ampère s law B = 0 and (µ 1 B) = J magnetic field B on R 3
21 MHD Equations: The Stationary Case Navier-Stokes system with Lorentz force ϱ (u )u η u + p = J B u = 0 fluid velocity u on Ω pressure p on Ω Charge conservation and Ohm s law σ 1 J + φ = u B J = 0 current density J on Ω electric potential φ on Ω No monopoles and Ampère s law B = 0 and (µ 1 B) = J magnetic field B on R 3
22 MHD Equations: The Stationary Case Navier-Stokes system with Lorentz force ϱ (u )u η u + p = J B u = 0 u = h Charge conservation and Ohm s law σ 1 J + φ = u B J = 0 J n = j φ = φ c No monopoles and Ampère s law B = 0 and (µ 1 B) = J
23 MHD Equations: The Stationary Case Navier-Stokes system with Lorentz force ϱ (u )u η u + p = J B u = 0 Charge conservation and Ohm s law σ 1 J + φ = u B J = 0 u = h J n = j φ = φ c Elimination of magnetic field velocity current formulation B = B(J)(x) = µ x y J(y) dy Biot Savart law 4π x y 3 R 3
24 MHD Equations: Velocity Current Formulation Coupled system ϱ (u )u η u + p = J B(J) u = 0 σ 1 J + φ = u B(J) J = 0 u = h J n = j φ = φ c
25 MHD Equations: Velocity Current Formulation Coupled system ϱ (u )u η u + p = J ( ) B(J) + B 0 u = h u = 0 σ 1 J + φ = u ( ) B(J) + B 0 J n = j J = 0 φ = φ c
26 MHD Equations: Velocity Current Formulation Coupled system Note ϱ (u )u η u + p = J ( B(J) + B 0 ) u = 0 σ 1 J + φ = u ( B(J) + B 0 ) J = 0 all quantities confined to fluid domain Ω u = h J n = j φ = φ c
27 MHD Equations: Velocity Current Formulation Coupled system Note ϱ (u )u η u + p = J ( B(J) + B 0 ) u = 0 σ 1 J + φ = u ( B(J) + B 0 ) J = 0 all quantities confined to fluid domain Ω some adjustable quantities, e.g., φ c u = h J n = j φ = φ c
28 MHD Equations: Velocity Current Formulation Coupled system Note ϱ (u )u η u + p = J ( B(J) + B 0 ) u = 0 σ 1 J + φ = u ( B(J) + B 0 ) J = 0 all quantities confined to fluid domain Ω some adjustable quantities, e.g., φ c saddle point structure u = h J n = j φ = φ c
29 MHD Equations: Analysis Nonlinear saddle point problem A 0 (u, J) + A 1 ((u, J), (u, J)) + B (p, φ) = F B(u, J) = 0
30 MHD Equations: Analysis Nonlinear saddle point problem A 0 (u, J) + A 1 ((u, J), (u, J)) + B (p, φ) = F B(u, J) = 0 Solution u H 1 (Ω) J H(div; Ω) = {J L 2 (Ω): J L 2 (Ω)} B(J) H 1 (R 3 ) p L 2 (Ω)/R φ L 2 (Ω)/R [1]; Meir, Schmidt: SIAM Journal on Numerical Analysis, 1999 [2]: Griesse, Kunisch: SIAM Journal on Control and Optimization, to appear
31 Related Work Previous and ongoing work M. Gunzburger, A.J. Meir, P. Schmidt J.-F. Gerbeau, C. Le Bris, T. Lelièvre J. Rappaz, R. Touzani many authors in engineering S. Hou, J. Peterson, A.J. Meir, S.S. Ravindran M. Hinze and co-workers M. Gunzburger, C. Trenchea
32 Overview 1 2 3
33 FEM Discretization Conforming and stable discretization (u, p) H 1 (Ω) L 2 (Ω)/R
34 FEM Discretization Conforming and stable discretization (u, p) H 1 (Ω) L 2 (Ω)/R Taylor-Hood elements
35 FEM Discretization Conforming and stable discretization (u, p) H 1 (Ω) L 2 (Ω)/R Taylor-Hood elements (J, φ) L 2 (div; Ω) L 2 (Ω)/R
36 FEM Discretization Conforming and stable discretization (u, p) H 1 (Ω) L 2 (Ω)/R Taylor-Hood elements (J, φ) L 2 (div; Ω) L 2 (Ω)/R Raviart-Thomas elements
37 FEM Discretization Conforming and stable discretization (u, p) H 1 (Ω) L 2 (Ω)/R Taylor-Hood elements (J, φ) L 2 (div; Ω) L 2 (Ω)/R Raviart-Thomas elements Discrete stability condition inf sup (q,ψ) (u,j) b((u, J), (q, ψ)) (u, J) (q, ψ) β
38 FEM Discretization Biot-Savart law If J = 0... B(J)(x) = µ 4π R 3 x y J(y) dy x y 3 B = µj on R 3 B = 0 on R 3
39 FEM Discretization Biot-Savart law B(J)(x) = µ 4π Introduction of vector potential B = A R 3 x y J(y) dy x y 3 B = µj on R 3 B = 0 on R 3 ( A) = µj on R 3
40 FEM Discretization Biot-Savart law B(J)(x) = µ 4π Introduction of vector potential B = A R 3 x y J(y) dy x y 3 B = µj on R 3 B = 0 on R 3 ( A) = µj on R 3 A = 0 on R 3 (gauging)
41 FEM Discretization Biot-Savart law B(J)(x) = µ 4π Introduction of vector potential B = A R 3 x y J(y) dy x y 3 ( A) = µj on R 3 A = 0 on R 3 (gauging)
42 FEM Discretization Biot-Savart law B(J)(x) = µ 4π Introduction of vector potential B = A R 3 x y J(y) dy x y 3 ( A) = µj on Ω A A = 0 on Ω A (gauging)
43 FEM Discretization Biot-Savart law B(J)(x) = µ 4π Introduction of vector potential B = A R 3 x y J(y) dy x y 3 ( A) = µj on Ω A A = 0 on Ω A (gauging) A n = 0 on Γ A
44 FEM Discretization Biot-Savart law B(J)(x) = µ 4π Introduction of vector potential B = A R 3 x y J(y) dy x y 3 ( A) = µj on Ω A A = 0 on Ω A (gauging) A n = 0 on Γ A Solvability condition J = 0 on Ω A
45 FEM Discretization Curl curl equation ( A) = µj on Ω A A = 0 on Ω A A n = 0 on Γ A
46 FEM Discretization Curl curl equation ( A) = µj on Ω A A = 0 on Ω A A n = 0 on Γ A Existence and uniqueness For J = 0, there exists a unique solution in H(curl; Ω A ).
47 FEM Discretization Curl curl equation ( A) = µj on Ω A A = 0 on Ω A A n = 0 on Γ A Existence and uniqueness For J = 0, there exists a unique solution in H(curl; Ω A ). Conforming discretization A H(curl; Ω A )
48 FEM Discretization Curl curl equation ( A) = µj on Ω A A = 0 on Ω A A n = 0 on Γ A Existence and uniqueness For J = 0, there exists a unique solution in H(curl; Ω A ). Conforming discretization A H(curl; Ω A ) Nédélec elements
49 FEM Discretization Condition on the boundary conditions J = 0 on Ω A
50 FEM Discretization Condition on the boundary conditions J = 0 on Ω A { J = 0 on Ω J n = 0 on Ω
51 FEM Discretization Condition on the boundary conditions J = 0 on Ω A { J = 0 on Ω J n = 0 on Ω Unless Ω = Ω A, this excludes φ = φ c!
52 FEM Discretization Condition on the boundary conditions J = 0 on Ω A { J = 0 on Ω J n = 0 on Ω Unless Ω = Ω A, this excludes φ = φ c! Use [φ] Γφ = φ c instead.
53 FEM Discretization Condition on the boundary conditions J = 0 on Ω A { J = 0 on Ω J n = 0 on Ω Unless Ω = Ω A, this excludes φ = φ c! Use [φ] Γφ = φ c instead. [Hiptmair, Sterz]
54 FEM Discretization Newton system M F δa G[J] A[u] B C[A] δu B δp H[u] C[A] D E δj = b E δφ
55 FEM Discretization Newton system M F δa G[J] A[u] B C[A] δu B δp H[u] C[A] D E δj = b E δφ
56 FEM Discretization Newton system M F δa G[J] A[u] B C[A] δu B δp H[u] C[A] D E δj = b E δφ
57 FEM Discretization Newton system M F δa G[J] A[u] B C[A] δu B δp H[u] C[A] D E δj = b E δφ
58 FEM Discretization Newton system M F δa G[J] A[u] B C[A] δu B δp H[u] C[A] D E δj = b E δφ Challenges mixture of finite element spaces preconditioning of linear systems
59 Simulation Results (joint with M. Discacciati) 9.506e e e e e 04 Problem description B 0 induced by currents in wires
60 Simulation Results (joint with M. Discacciati) 9.506e e e e e 04 Problem description B 0 induced by currents in wires current J n = ±j at top/bottom (Ω = Ω A )
61 Simulation Results (joint with M. Discacciati) Fluid velocity (from top, slice at half height), Stokes 3.690e e e e e 03 Problem description B 0 induced by currents in wires current J n = ±j at top/bottom (Ω = Ω A ) two counter-rotating flow cells ([Gerbeau], 2000)
62 Simulation Results Problem size order of Nédélec FE Taylor-Hood FE Raviart-Thomas FE Grid 1 dofs (344 tetr.) iterations Grid 2 dofs (2752 tetr.) iterations Details iterative damped splitting scheme: A and (u, p, J, φ) implementation in Ngsolve, sparse direct solver Pardiso
63 Overview 1 2 3
64 Problem formulation Minimize f (y, u) subject to e(y, u) = 0
65 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state
66 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state Necessary conditions
67 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state Necessary conditions L y (y, u, λ) = f y (y, u) + e y (y, u) λ = 0
68 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state Necessary conditions L y (y, u, λ) = f y (y, u) + e y (y, u) λ = 0 adjoint PDE system
69 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state Necessary conditions L y (y, u, λ) = f y (y, u) + e y (y, u) λ = 0 L u (y, u, λ) = f u (y, u) + e u(y, u) λ = 0 adjoint PDE system
70 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state Necessary conditions L y (y, u, λ) = f y (y, u) + e y (y, u) λ = 0 L u (y, u, λ) = f u (y, u) + e u(y, u) λ = 0 adjoint PDE system gradient equation
71 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state Necessary conditions L y (y, u, λ) = f y (y, u) + e y (y, u) λ = 0 L u (y, u, λ) = f u (y, u) + e u(y, u) λ = 0 L λ (y, u, λ) = e(y, u) = 0 adjoint PDE system gradient equation
72 Problem formulation Minimize f (y, u) subject to e(y, u) = 0 Lagrange functional L(y, u, λ) := f (y, u) + e(y, u), λ λ adjoint state Necessary conditions L y (y, u, λ) = f y (y, u) + e y (y, u) λ = 0 L u (y, u, λ) = f u (y, u) + e u(y, u) λ = 0 L λ (y, u, λ) = e(y, u) = 0 adjoint PDE system gradient equation PDE system
73 Newton s method in function space L yy L yu e y δy L y L uy L uu eu δu = L u e y e u 0 δλ e }{{} KKT matrix
74 Newton s method in function space L yy L yu e y δy L y L uy L uu eu δu = L u e y e u 0 δλ e }{{} KKT matrix Numerical challenges
75 Newton s method in function space L yy L yu e y δy L y L uy L uu eu δu = L u e y e u 0 δλ e }{{} KKT matrix Numerical challenges large system of equations ( variables)
76 Newton s method in function space L yy L yu e y δy L y L uy L uu eu δu = L u e y e u 0 δλ e }{{} KKT matrix Numerical challenges large system of equations ( variables) symmetric, indefinite, ill conditioned
77 An MHD Problem A possible problem setup
78 An MHD Problem A possible problem setup fluid region Ω
79 An MHD Problem A possible problem setup electrodes fluid region Ω
80 An MHD Problem A possible problem setup electrodes fluid region Ω Purpose and control mechanisms influence flow pattern (stir or dampen)
81 An MHD Problem A possible problem setup electrodes fluid region Ω Purpose and control mechanisms influence flow pattern (stir or dampen) using adjustable quantities
82 An MHD Problem A possible problem setup electrodes φ = { φ c control R 0 fluid region Ω Purpose and control mechanisms influence flow pattern (stir or dampen) using adjustable quantities (applied potential difference)
83 An MHD Problem A possible problem setup electrodes φ = { φ c control R 0 fluid region Ω = Ω A Purpose and control mechanisms influence flow pattern (stir or dampen) using adjustable quantities (applied potential difference)
84 An MHD Problem Problem description Minimize J = α 2 u u d 2 L 2 (Ω) + γ 2 φ c 2 s.t. MHD system
85 An MHD Problem Problem description Minimize J = α 2 u u d 2 L 2 (Ω) + γ 2 φ c 2 s.t. MHD system φ = φ c at electrode 1 φ = 0 at electrode 2 J n = 0 elsewhere
86 An MHD Problem Problem description Minimize J = α 2 u u d 2 L 2 (Ω) + γ 2 φ c 2 s.t. MHD system φ = φ c at electrode 1 φ = 0 at electrode 2 J n = 0 elsewhere Given problem data u d desired velocity field; cost parameters α 0 and γ 0
87 An MHD Problem Problem description Minimize J = α 2 u u d 2 L 2 (Ω) + γ 2 φ c 2 s.t. MHD system φ = φ c at electrode 1 φ = 0 at electrode 2 J n = 0 elsewhere Given problem data u d desired velocity field; cost parameters α 0 and γ 0 applied magnetic field B 0
88 An MHD Problem Problem description Given problem data Minimize J = α 2 u u d 2 L 2 (Ω) + γ 2 φ c 2 s.t. MHD system φ = φ c at electrode 1 φ = 0 at electrode 2 J n = 0 elsewhere u d desired velocity field; cost parameters α 0 and γ 0 applied magnetic field B 0 u = h on the boundary Ω
89 An MHD Problem Adjoint system on Ω ϱ ( u) v ϱ (u ) v η v + q (B(J) + B 0 ) K = σ 1 K B(K u + v J) + ψ (B(J) + B 0 ) v = Incompressibility and boundary conditions v = 0 on Ω v = 0 on Ω K = 0 on Ω K n = 0 or ψ = 0 on Ω Adjoint variables on Ω v adjoint velocity K adjoint current q adjoint pressure ψ adjoint potential
90 An MHD Problem Adjoint system on Ω ϱ ( u) v ϱ (u ) v η v + q (B(J) + B 0 ) K = σ 1 K B(K u + v J) + ψ (B(J) + B 0 ) v = Incompressibility and boundary conditions v = 0 on Ω v = 0 on Ω K = 0 on Ω K n = 0 or ψ = 0 on Ω Adjoint variables on Ω v adjoint velocity K adjoint current q adjoint pressure ψ adjoint potential
91 An MHD Problem Adjoint system on Ω ϱ ( u) v ϱ (u ) v η v + q (B(J) + B 0 ) K = σ 1 K B(K u + v J) + ψ (B(J) + B 0 ) v = Incompressibility and boundary conditions v = 0 on Ω v = 0 on Ω K = 0 on Ω K n = 0 or ψ = 0 on Ω Adjoint variables on Ω v adjoint velocity K adjoint current q adjoint pressure ψ adjoint potential
92 An MHD Problem Adjoint system on Ω ϱ ( u) v ϱ (u ) v η v + q (B(J) + B 0 ) K = σ 1 K B(K u + v J) + ψ (B(J) + B 0 ) v = Incompressibility and boundary conditions v = 0 on Ω v = 0 on Ω K = 0 on Ω K n = 0 or ψ = 0 on Ω Optimality condition γφ c + K n = 0 at electrodes
93 Numerics for KKT matrix G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A[u] B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A[u] B C[A] B H[u] C[A] D E E
94 Numerics for KKT matrix G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A[u] B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A[u] B C[A] B H[u] C[A] D E E
95 Numerics for KKT matrix G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A[u] B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A[u] B C[A] B H[u] C[A] D E E
96 Numerics for KKT matrix Simplifications: G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A[u] B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A[u] B C[A] B H[u] C[A] D E E
97 Numerics for KKT matrix Simplifications: Stokes flow G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A[u] B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A[u] B C[A] B H[u] C[A] D E E
98 Numerics for KKT matrix Simplifications: Stokes flow G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A B C[A] B H[u] C[A] D E E
99 Numerics for KKT matrix Simplifications: Stokes flow, B = B 0 known (low R m approx.) G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A B C[A] B H[u] C[A] D E E
100 Numerics for KKT matrix Simplifications: Stokes flow, B = B 0 known (low R m approx.) G[λ J ] H[λ u ] M G[J] H[u] G[λ J ] A B C[A] B H[λ u ] F C[A] D E E γ M F G[J] A B C[A] B H[u] C[A] D E E
101 Numerical Results (joint with M. Discacciati) Problem data
102 Numerical Results (joint with M. Discacciati) Problem data grounded electrode (φ = 0)
103 Numerical Results (joint with M. Discacciati) Problem data control electrode (φ = φ c ) grounded electrode (φ = 0)
104 Numerical Results (joint with M. Discacciati) Problem data control electrode (φ = φ c ) grounded electrode (φ = 0) B 0 = 10 4 (0, 0, x)t
105 Numerical Results (joint with M. Discacciati) Problem data control electrode (φ = φ c ) grounded electrode (φ = 0) B 0 = 10 4 (0, 0, x)t u d = swirl flow
106 Numerical Results (joint with M. Discacciati) Problem data control electrode (φ = φ c ) grounded electrode (φ = 0) B 0 = 10 4 (0, 0, x)t u d = swirl flow Material data and solution material data of liquid Al at 700 C optimal control φ c = V at γ = 0.1 current J max = A/m 2, velocity u max = m/s R m = µσul =
107 Numerical Results Problem discretization 768 tetrahedra Taylor-Hood (u, p) order 2 1 Raviart-Th. (J, φ) order degrees of freedom implementation in Ngsolve CPU time matrix setup 30 s sparse factorization (Pardiso) 36 s solution of optimal control problem < 1 s
108 Numerical Results Optimal solution (potential φ)
109 Numerical Results Optimal solution (current J)
110 Numerical Results Optimal solution (velocity u)
111 Numerical Results Optimal solution (velocity u)
112 Concluding Remarks Summary simulation and optimization problems in MHD numerous applications in metallurgy, crystal growth discretization: Taylor-Hood, Raviart-Thomas and Nédélec FE
113 Concluding Remarks Summary simulation and optimization problems in MHD numerous applications in metallurgy, crystal growth discretization: Taylor-Hood, Raviart-Thomas and Nédélec FE Future challenges time-dependent and Navier-Stokes cases adaptivity preconditioning of linear systems
114 Concluding Remarks Summary simulation and optimization problems in MHD numerous applications in metallurgy, crystal growth discretization: Taylor-Hood, Raviart-Thomas and Nédélec FE Future challenges time-dependent and Navier-Stokes cases adaptivity preconditioning of linear systems Thank you!
A Practical Optimal Control Approach to the Stationary MHD System in Velocity Current Formulation
www.oeaw.ac.at A Practical Optimal Control Approach to the Stationary MHD System in Velocity Current Formulation R. Griesse, K. Kunisch RICAM-Report 2005-02 www.ricam.oeaw.ac.at A PRACTICAL OPTIMAL CONTROL
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 informationAlgorithms for PDE-Constrained Optimization
GAMM-Mitteilungen, 31 January 2014 Algorithms for PDE-Constrained Optimization Roland Herzog 1 and Karl Kunisch 2 1 Chemnitz University of Technology, Faculty of Mathematics, Reichenhainer Straße 41, D
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 informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
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 informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More informationSHALLOW WATER MODEL FOR ALUMINIUM ELECTROLYSIS CELLS WITH VARIABLE TOP AND BOTTOM
SHALLOW WATER MODEL FOR ALUMINIUM ELECTROLYSIS CELLS WITH VARIABLE TOP AND BOTTOM Valdis Boarevics and Koulis Pericleous University of Greenwich, School of Computing and Mathematical Sciences, 30 Park
More informationMATH 676. Finite element methods in scientific computing
MATH 676 Finite element methods in scientific computing Wolfgang Bangerth, Texas A&M University Lecture 33.25: Which element to use Part 2: Saddle point problems Consider the stationary Stokes equations:
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
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 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 informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More 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 informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
More informationCalculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
More informationSolving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Solving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics Fikriye Yılmaz 1, Bülent Karasözen
More informationSimulation based optimization
SimBOpt p.1/52 Simulation based optimization Feb 2005 Eldad Haber haber@mathcs.emory.edu Emory University SimBOpt p.2/52 Outline Introduction A few words about discretization The unconstrained framework
More informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
More informationParallelizing large scale time domain electromagnetic inverse problem
Parallelizing large scale time domain electromagnetic inverse problems Eldad Haber with: D. Oldenburg & R. Shekhtman + Emory University, Atlanta, GA + The University of British Columbia, Vancouver, BC,
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 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 informationAn Implementation of the Finite Element Method for the Velocity-Current Magnetohydrodynamics Equations. K. Daniel Brauss
An Implementation of the Finite Element Method for the Velocity-Current Magnetohydrodynamics Equations by K. Daniel Brauss A dissertation submitted to the Graduate Faculty of Auburn University in partial
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationSolving the curl-div system using divergence-free or curl-free finite elements
Solving the curl-div system using divergence-free or curl-free finite elements Alberto Valli Dipartimento di Matematica, Università di Trento, Italy or: Why I say to my students that divergence-free finite
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationON TWO OPTIMAL CONTROL PROBLEMS FOR MAGNETIC FIELDS
ON TWO OPTIMAL CONTROL PROBLEMS FOR MAGNETIC FIELDS SERGE NICAISE, SIMON STINGELIN, AND FREDI TRÖLTZSCH Abstract. Two optimal control problems for instationary magnetization processes are considered in
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
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 informationFinite element exterior calculus: A new approach to the stability of finite elements
Finite element exterior calculus: A new approach to the stability of finite elements Douglas N. Arnold joint with R. Falk and R. Winther Institute for Mathematics and its Applications University of Minnesota
More informationFEM for Stokes Equations
FEM for Stokes Equations Kanglin Chen 15. Dezember 2009 Outline 1 Saddle point problem 2 Mixed FEM for Stokes equations 3 Numerical Results 2 / 21 Stokes equations Given (f, g). Find (u, p) s.t. u + p
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationA NUMERICAL STUDY OF METAL PAD ROLLING INSTABILITY IN A SIMPLIFIED HALL-HÉROULT CELL
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK A NUMERICAL STUDY OF METAL PAD ROLLING INSTABILITY
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationA MAGNETOHYDRODYNAMIC STUDY OF BEHAVIOR IN AN ELECTROLYTE FLUID USING NUMERICAL AND EXPERIMENTAL SOLUTIONS
A MAGNETOHYDRODYNAMIC STUDY OF BEHAVIOR IN AN ELECTROLYTE FLUID USING NUMERICAL AND EXPERIMENTAL SOLUTIONS L. P. Aoki, M. G. Maunsell, and H. E. Schulz Universidade de São Paulo Escola de Engenharia de
More informationWorkshop Control of Complex Fluids
ÂÓ ÒÒ Ê ÓÒ ÁÒ Ø ØÙØ ÓÖ ÓÑÔÙØ Ø ÓÒ Ð Ò ÔÔÐ Å Ø Ñ Ø Ù ØÖ Ò ÑÝ Ó Ë Ò Workshop Control of Complex Fluids October 10 13, 2005 Schedule of Talks Monday Tuesday Wednesday Thursday 9:45 10:00 Opening 10:00 11:00
More informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
More informationSimulations of Electrical Arcs: Algorithms, Physical Scales, and Coupling. Henrik Nordborg HSR University of Applied Sciences Rapperswil
Simulations of Electrical Arcs: Algorithms, Physical Scales, and Coupling Henrik Nordborg HSR University of Applied Sciences Rapperswil What is an electrical arc? 2 Technical applications of arcs and industrial
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
More informationProper Orthogonal Decomposition in PDE-Constrained Optimization
Proper Orthogonal Decomposition in PDE-Constrained Optimization K. Kunisch Department of Mathematics and Computational Science University of Graz, Austria jointly with S. Volkwein Dynamic Programming Principle
More informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationarxiv: v2 [cs.ce] 22 Oct 2016
SIMULATION OF ELECTRICAL MACHINES A FEM-BEM COUPLING SCHEME LARS KIELHORN, THOMAS RÜBERG, JÜRGEN ZECHNER arxiv:1610.05472v2 [cs.ce] 22 Oct 2016 Abstract. Electrical machines commonly consist of moving
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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
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 informationAC & DC Magnetic Levitation and Semi-Levitation Modelling
International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, March 24-26, 2003 AC & DC Magnetic Levitation and Semi-Levitation Modelling V. Bojarevics, K. Pericleous Abstract
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 informationPhysic-based Preconditioning and B-Splines finite elements method for Tokamak MHD
Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD E. Franck 1, M. Gaja 2, M. Mazza 2, A. Ratnani 2, S. Serra Capizzano 3, E. Sonnendrücker 2 ECCOMAS Congress 2016, 5-10
More informationMotivations. Outline. Finite element exterior calculus and the geometrical basis of numerical stability. References. Douglas N.
Outline Finite element exterior calculus and the geometrical basis of numerical stability Douglas N. Arnold joint with R. Falk and R. Winther School of Mathematics University of Minnesota April 009 Motivations
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 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 informationNumerical Simulation of the VAR Process with calcosoft -2D and its Validation
Numerical Simulation of the VAR Process with calcosoft -2D and its Validation G. Reiter a, V. Maronnier b, C. Sommitsch a, M. Gäumann b, W. Schützenhöfer a, R. Schneider a a Böhler Edelstahl GmbH & Co
More informationarxiv: v4 [physics.comp-ph] 21 Jan 2019
A spectral/hp element MHD solver Alexander V. Proskurin,Anatoly M. Sagalakov 2 Altai State Technical University, 65638, Russian Federation, Barnaul, Lenin prospect,46, k2@list.ru 2 Altai State University,
More informationLectures Notes Algorithms and Preconditioning in PDE-Constrained Optimization. Prof. Dr. R. Herzog
Lectures Notes Algorithms and Preconditioning in PDE-Constrained Optimization Prof. Dr. R. Herzog held in July 2010 at the Summer School on Analysis and Numerics of PDE Constrained Optimization, Lambrecht
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More informationA fully divergence-free finite element method for magneto-hydrodynamic equations
A fully divergence-free finite element method for magneto-hydrodynamic equations R. Hiptmair and L.-X. Li and S.-P. Mao and W.-Y. Zheng Research Report No. 2017-25 May 2017 Seminar für Angewandte Mathematik
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 informationEffect of EM field frequency on 3D melt convection during floating zone growth of silicon
This work is done in the University of Latvia with the support of European Regional Development Fund project Nr. 2013/0051/2DP/2.1.1.1.0/13/APIA/VIAA/009 Effect of EM field frequency on 3D melt convection
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationCurrents (1) Line charge λ (C/m) with velocity v : in time t, This constitutes a current I = λv (vector). Magnetic force on a segment of length dl is
Magnetostatics 1. Currents 2. Relativistic origin of magnetic field 3. Biot-Savart law 4. Magnetic force between currents 5. Applications of Biot-Savart law 6. Ampere s law in differential form 7. Magnetic
More informationEddy Current Interaction of a Magnetic Dipole With a Translating Solid Bar
d h International Scientific Colloquium Modelling for Material Processing Riga, September 16-17, 21 Eddy Current Interaction of a Magnetic Dipole With a Translating Solid Bar M. Kirpo, T. Boeck, A. Thess
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationAposteriorierrorestimatesinFEEC for the de Rham complex
AposteriorierrorestimatesinFEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094 and a
More informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More informationEfficient Solvers for the Navier Stokes Equations in Rotation Form
Efficient Solvers for the Navier Stokes Equations in Rotation Form Computer Research Institute Seminar Purdue University March 4, 2005 Michele Benzi Emory University Atlanta, GA Thanks to: NSF (MPS/Computational
More information- Part 4 - Multicore and Manycore Technology: Chances and Challenges. Vincent Heuveline
- Part 4 - Multicore and Manycore Technology: Chances and Challenges Vincent Heuveline 1 Numerical Simulation of Tropical Cyclones Goal oriented adaptivity for tropical cyclones ~10⁴km ~1500km ~100km 2
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationOn the local well-posedness of compressible viscous flows with bounded density
On the local well-posedness of compressible viscous flows with bounded density Marius Paicu University of Bordeaux joint work with Raphaël Danchin and Francesco Fanelli Mathflows 2018, Porquerolles September
More informationNumerical Study of a DC Electromagnetic Liquid Metal Pump: Limits of the Model Nedeltcho Kandev
Numerical Study of a DC Electromagnetic Liquid Metal Pump: Limits of the Model Nedeltcho Kandev Excerpt from the Proceedings of the 2012 COMSOL Conference in Boston Introduction This work presents the
More informationElmer. Introduction into Elmer multiphysics. Thomas Zwinger. CSC Tieteen tietotekniikan keskus Oy CSC IT Center for Science Ltd.
Elmer Introduction into Elmer multiphysics FEM package Thomas Zwinger CSC Tieteen tietotekniikan keskus Oy CSC IT Center for Science Ltd. Contents Elmer Background History users, community contacts and
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 informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationFinite Element Modeling of Electromagnetic Systems
Finite Element Modeling of Electromagnetic Systems Mathematical and numerical tools Unit of Applied and Computational Electromagnetics (ACE) Dept. of Electrical Engineering - University of Liège - Belgium
More informationarxiv: v1 [math-ph] 25 Apr 2013
Higher-order compatible discretization on hexahedrals Jasper Kreeft and Marc Gerritsma arxiv:1304.7018v1 [math-ph] 5 Apr 013 Abstract We derive a compatible discretization method that relies heavily on
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 informationRobust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems
Robust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems M. Aaqib Afaq Institute for Applied Mathematics and Numerics (LSIII) TU Dortmund 13 July 2017
More informationOn the choice of abstract projection vectors for second level preconditioners
On the choice of abstract projection vectors for second level preconditioners C. Vuik 1, J.M. Tang 1, and R. Nabben 2 1 Delft University of Technology 2 Technische Universität Berlin Institut für Mathematik
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 informationNumerical approximation of output functionals for Maxwell equations
Numerical approximation of output functionals for Maxwell equations Ferenc Izsák ELTE, Budapest University of Twente, Enschede 11 September 2004 MAXWELL EQUATIONS Assumption: electric field ( electromagnetic
More informationA Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface
A Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface Damanik, H., Mierka, O., Ouazzi, A., Turek, S. (Lausanne, August 2013) Page 1 Motivation Polymer melts:
More informationMagnetostatics and the vector potential
Magnetostatics and the vector potential December 8, 2015 1 The divergence of the magnetic field Starting with the general form of the Biot-Savart law, B (x 0 ) we take the divergence of both sides with
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationA Numerical Study of Natural Convection in the Horizontal Bridgman Configuration under an External Magnetic Field
Proceedings of the Fourth International Conference on Thermal Engineering: Theory and Applications January 12-14, 29, Abu Dhabi, UAE A Numerical Study of Natural Convection in the Horizontal Bridgman Configuration
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationSonderforschungsbereich 609 Elektromagnetische Strömungsbeeinflussung in Metallurgie, Kristallzüchtung und Elektrochemie
SFB 609 Sonderforschungsbereich 609 Elektromagnetische Strömungsbeeinflussung in Metallurgie, Kristallzüchtung und Elektrochemie M. Hinze, S. Ziegenbalg Optimal control of the free boundary in a two-phase
More informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
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 informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
More informationAn advanced ILU preconditioner for the incompressible Navier-Stokes equations
An advanced ILU preconditioner for the incompressible Navier-Stokes equations M. ur Rehman C. Vuik A. Segal Delft Institute of Applied Mathematics, TU delft The Netherlands Computational Methods with Applications,
More information