Failsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model
|
|
- Austin Stevens
- 6 years ago
- Views:
Transcription
1 Failsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model Matthias Möller 1 joint work with: Dmitri Kuzmin 2 1 Institute of Applied Mathematics (LS3), Dortmund University of Technology, Germany, matthias.moeller@math.tu-dortmund.de 2 Chair of Applied Mathematics III, University Erlangen-Nuremberg, Germany kuzmin@am.uni-erlangen.de
2 Motivation: Z-pinch implosion Phenomenological model by Banks and Shadid compressible Euler equations + source term coupled with tracer equation for Lorentz force Mathematical challenges high-resolution scheme for time-dependent conservation laws; positivity-preservation of density and pressure; failsafe strategy
3 High-resolution schemes for t U + A U = 0 High-order scheme M C du H dt = KU H Low-order scheme M L du L dt = LU L, L = K + D Predictor: Compute low-order solution M L du L dt = LU L U L M 1 L LU L Corrector: Apply limited antidiffusion M L U = M L U L + F, F = [M L M C ] U L DU L
4 High-resolution schemes for t U + A U = 0 High-order scheme M C du H dt = KU H Low-order scheme M L du L dt = LU L, L = K + D Predictor: Compute low-order solution M L du L dt = LU L U L M 1 L LU L Corrector: Apply limited antidiffusion M L U = M L U L + F, F = [M L M C ] U L DU L Low-order scheme must satisfy physical constraints
5 Design principles of FCT schemes Perform flux correction such that mass is conserved.
6 Design principles of FCT schemes Perform flux correction such that mass is conserved. Conservative flux decomposition m i U i = m i U L i + j i F ij, F ji = F ij
7 Design principles of FCT schemes Perform flux correction such that mass is conserved. Conservative flux decomposition and limiting m i U i = m i U L i + j i α ij F ij, F ji = F ij, α ji = α ij high-order approximation (α ij = 1) to be used in smooth regions low-order approximation (α ij = 0) to be used near steep fronts
8 Design principles of FCT schemes, cont d Perform flux correction such that certain physical quantities are bounded by the local extrema of the low-order solution.
9 Design principles of FCT schemes, cont d Perform flux correction such that certain physical quantities are bounded by the local extrema of the low-order solution. u min i u L i + 1 m i j i α u ij f u ij umax i u L i
10 Fixed fraction flux limiter Set α (0) ij := 1 and repeat r = 1,..., R Mark all nodes i that violate the local FCT constraint u min i u L i + 1 m i j i α (r 1) ij fij u umax i Eliminate fixed fraction of unacceptable antidiffusion { 1 r/r if node i or j is marked α (r) ij := otherwise α (r 1) ij Low-order solution is recovered in the worst case: U i = U L i
11 Zalesak s flux limiter Consider positive/negative antidiffusive contributions separately P + i u i P i
12 Zalesak s flux limiter Q + i P + i Q i P i Consider positive/negative antidiffusive contributions separately Limit antidiffusion if it exceeds the distance to local maximum/minimum
13 Zalesak s flux limiter Q + i P + i Q i P i Consider positive/negative antidiffusive contributions separately Limit antidiffusion if it exceeds the distance to local maximum/minimum Nodal correction factors R + i = min{1, Q + i /P + i } for positive fluxes into node i R i = min{1, Q i /P i } for negative fluxes into node i
14 Zalesak s flux limiter Q + i P + i Q i P i Consider positive/negative antidiffusive contributions separately Limit antidiffusion if it exceeds the distance to local maximum/minimum Nodal correction factors R + i = min{1, Q + i /P + i R i = min{1, Q i /P i } for positive fluxes into node i } for negative fluxes into node i Limit antidiffusive flux for edge ij by the minimum of R i and R j
15 Flux limiting for systems Apply flux limiter to a set of control variables, e.g., the primitive variables density, pressure and velocity.
16 Flux limiting for systems Apply flux limiter to a set of control variables, e.g., the primitive variables density, pressure and velocity. f ρ ij f p ij f v ij FCT FCT FCT min{α ρ ij, αp ij, αv ij}
17 Flux limiting for systems Apply flux limiter to a set of control variables, e.g., the primitive variables density, pressure and velocity. f ρ ij f p ij f v ij f ρ ij α ρ ij f p ij α p ij αρ ij f v ij FCT FCT FCT or FCT FCT FCT min{α ρ ij, αp ij, αv ij} α ρ ij α p ij αρ ij α v ijα p ij αρ ij
18 Flux limiting for systems Apply flux limiter to a set of control variables, e.g., the primitive variables density, pressure and velocity. f ρ ij f p ij f v ij f ρ ij α ρ ij f p ij α p ij αρ ij f v ij FCT FCT FCT or FCT FCT FCT min{α ρ ij, αp ij, αv ij} α ρ ij α p ij αρ ij α v ijα p ij αρ ij Nodal transformation of variables V i = T (U i )U i and G ij = T (U i )F ij
19 Failsafe flux correction algorithm 1 Compute low-order solution at time t n+1 M L U L U n t = θlu L + (1 θ)lu n prediction 2 Perform flux correction by Zalesak s limiter m i U (0) i = m i U L i + j i α ij F ij correction 3 Eliminate spurious undershoots/overshoots m i U (r) i = m i Ui L + α (r) ij [α ijf ij ] j i failsafe step
20 Double Mach reflection From: P.R. Woodward and P. Colella, JCP 54, 115 (1984) U L = cos sin Γ L Γ R U R = Γ out Γ wall
21 Double Mach reflection, cont d Solution at T = 0.2 computed by low-order scheme (α ij 0) density distribution density distribution h = 1/64 t = h = 1/128 t =
22 Double Mach reflection, cont d Solution at T = 0.2 computed by fixed fraction flux limiter density distribution density distribution h = 1/64 t = h = 1/128 t =
23 Double Mach reflection, cont d Solution at T = 0.2 computed by Zalesak s flux limiter (α ij = α p ij αρ ij ) density distribution density distribution h = 1/64 t = h = 1/128 t =
24 Double Mach reflection, cont d Solution at T = 0.2 computed by Zalesak s flux limiter (α ij = α p ij αρ ij ) density distribution pressure distribution h = 1/256 t = h = 1/256 t = Numerical studies: D. Kuzmin, M. Shashkov, J. Shadid, M.M. (to appear in JCP)
25 Idealized Z-pinch implosion model by Banks and Shadid Generalized Euler system coupled with scalar tracer equation ρ ρv 0 ρv t ρe + ρv v + pi ρev + pv = f f v ρλ ρλv 0 Equation of state p = (γ 1)ρ(E 0.5 v 2 ) Non-dimensional Lorentz force ( ) 2 f = (ρλ) I(t) I max êr r eff
26 Coupled solution algorithm time stepping loop Given: (U, ρλ) n Low-order scheme (ρλ) L,(k) M L du dt = L(U)U + S(v, ρλ) M L d(ρλ) dt = L(v)(ρλ) v L,(k) (U, ρλ) L Zalesak s limiter α ij = α ρλ ij αp ij αρ ij [+ failsafe corr.] (U, ρλ) n+1
27 Idealized Z-pinch implosion From: J.W. Banks and J.N. Shadid, JCP 61, 725 (2009) ρ = 1.0 λ = 0.0 ρ = 10 6 λ = 1.0 = 0.05 ρ = 0.5 λ = 0.0 v = 0.0, p = 1.0 everywhere
28 Idealized Z-pinch implosion, cont d density distribution
29 Idealized Z-pinch implosion, cont d density distribution density distribution
30 Idealized Z-pinch implosion, cont d density distribution density distribution density distribution
31 Idealized Z-pinch implosion, cont d radius exact R(t)=1 t 4 low order Rusanov FCT lumped mass FCT consistent mass time shell radius as a function of time z density symmetry plot, t = 0.9 density symmetry plot, t = 1.05 R(t) r
32 Idealized Z-pinch implosion, cont d Initialization radius exact R(t)=1 t 4 low order Rusanov FCT lumped mass FCT consistent mass time shell radius as a function of time z density symmetry plot, t = 0.9 density symmetry plot, t = 1.05 R(t) r
33 Conclusions and outlook Linearized flux correction algorithm for time-dependent flows mass conservation, boundedness of physical quantities, failsafe strategy if density/pressure becomes negative Coupled solution algorithm for idealized Z-pinch implosions positivity and symmetry preservation on unstructured grids Todo: comparative study and extension to realistic scenarios current drive, r-z plane, RT-instabilities, AMR
34
35 Appendix
36 Extended version of Zalesak s FCT limiter Return Input: auxiliary solution u L and antidiffusive fluxes f u ij, where f u ji f u ij 1 Sums of positive/negative antidiffusive fluxes into node i P + i = j i max{0, f u ij}, P i = j i min{0, f u ij} 2 Upper/lower bounds based on the local extrema of u L Q + i = m i (u max i u L i ), Q i = m i (u min i u L i ) 3 Correction factors αij u = αu ji to satisfy the FCT constraints α u ij = min{r ij, R ji }, R ij = { min{1, Q + i /P + i } if f u ij 0 min{1, Q i /P i } if f u ij < 0
37 Node-based transformation of control variables Return Conservative variables: density, momentum, total energy [ ] U i = [ρ i, (ρv) i, (ρe) i ], F ij = f ρ ij, f ρv ij, f ρe ij, F ji = F ij Primitive variables V = T U: density, velocity, pressure ] V i = [ρ i, v i, p i ], v i = (ρv)i ρ i, p i = (γ 1) [(ρe) i (ρv)i 2 2ρ i G ij = [ f ρ ij, f v ij, f p ij] = T (Ui )F ij, T (U j )F ji = G ji G ij Raw antidiffusive fluxes for the velocity and pressure [ fij v = f ρv ij vif ρ ij ρ i, f p ij = (γ 1) vi 2 2 f ρ ij v i f ρv ij ] + f ρe ij
38 Constrained initialization Pointwise initialization U(x i ) = U 0 (x i ) { 1.0 in Ω1 ρ = 0.01 in Ω 2 u = v = 0.0, p = 1.0
39 Constrained initialization Pointwise initialization U(x i ) = U 0 (x i ) Conservative initialization Ω wu h dx = Ω wu 0 dx { 1.0 in Ω1 ρ = 0.01 in Ω 2 u = v = 0.0, p = 1.0
40 Constrained initialization Pointwise initialization U(x i ) = U 0 (x i ) Conservative initialization Ω wu h dx = Ω wu 0 dx Consistent L 2 -projection j m ijuj H = Ω ϕ iu 0 dx { 1.0 in Ω1 ρ = 0.01 in Ω 2 u = v = 0.0, p = 1.0
41 Constrained initialization Pointwise initialization U(x i ) = U 0 (x i ) Conservative initialization Ω wu h dx = Ω wu 0 dx Consistent L 2 -projection j m ijuj H = Ω ϕ iu 0 dx Mass-lumped L 2 -projection m i U L i = Ω ϕ iu 0 dx { 1.0 in Ω1 ρ = 0.01 in Ω 2 u = v = 0.0, p = 1.0
42 Constrained initialization Pointwise initialization U(x i ) = U 0 (x i ) Conservative initialization Ω wu h dx = Ω wu 0 dx Consistent L 2 -projection j m ijuj H = Ω ϕ iu 0 dx Mass-lumped L 2 -projection m i U L i = Ω ϕ iu 0 dx { 1.0 in Ω1 ρ = 0.01 in Ω 2 u = v = 0.0, p = 1.0 Limited L 2 -projection (0 α ij 1) m i U i = m i U L i + j i α ij m ij (U L i U L j )
43 Initialization for bilinear elements (a) consistent L 2-projection (b) lumped L 2-projection (c) L 2-projection, α ij = α ρ ij bilinear elements, 3 3 Gauss rule ρ ρ h 2 min(ρ h ) max(ρ h ) (a) 1.048e e (b) 1.168e e (c) 1.103e e computed by adaptive cubature formulae
44 Initialization for linear elements Conclusions (a) consistent L 2-projection (b) lumped L 2-projection (c) L 2-projection, α ij = α ρ ij linear elements, 3-point Gauss rule ρ ρ h 2 min(ρ h ) max(ρ h ) (a) 1.206e e (b) 1.357e e (c) 1.259e e computed by adaptive cubature formulae
High-resolution finite element schemes for an idealized Z-pinch implosion model
High-resolution finite element schemes for an idealized Z-pinch implosion model Matthias Möller a, Dmitri Kuzmin b, John N. Shadid c a Institute of Applied Mathematics (LS III), Dortmund University of
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationAlgebraic Flux Correction II
Algebraic Flux Correction II Compressible Flow Problems Dmitri Kuzmin, Matthias Möller, and Marcel Gurris Abstract Flux limiting for hyperbolic systems requires a careful generalization of the design principles
More informationIsogeometric Analysis for Compressible Flow Problems in Industrial Applications
Isogeometric Analysis for Compressible Flow Problems in Industrial Applications Matthias Möller Delft Institute of Applied Mathematics, Numerical Analysis Delft University of Technology I 3 MS Seminar,
More informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationExplicit and implicit FEM-FCT algorithms with flux linearization
Explicit and implicit FEM-FCT algorithms with flux linearization Dmitri Kuzmin Institute of Applied Mathematics (LS III), Dortmund University of Technology Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationA parameter-free smoothness indicator for high-resolution finite element schemes
Cent. Eur. J. Math. 11(8) 2013 1478-1488 DOI: 10.2478/s11533-013-0254-4 Central European Journal of Mathematics A parameter-free smoothness indicator for high-resolution finite element schemes Research
More informationImplicit FEM-FCT algorithms and discrete Newton methods for transient convection problems
Implicit FEM-FCT algorithms and discrete Newton methods for transient convection problems M. Möller 1,,, D. Kuzmin 1, and D. Kourounis 2 1 Institute of Applied Mathematics (LS III), University of Dortmund
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationA flux-corrected RBF-FD method for convection dominated problems in domains and on manifolds
A flux-corrected RBF-FD method for convection dominated problems in domains and on manifolds Andriy Sokolov Dmitri Kuzmin, Oleg Davydov and Stefan Turek Institut für Angewandte Mathematik (LS3) TU Dortmund
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationOn the design of general-purpose flux limiters for finite element schemes. I. Scalar convection
On the design of general-purpose flux limiters for finite element schemes. I. Scalar convection D. Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund,
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More 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 Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationThe Design of Flux-Corrected Transport (FCT) Algorithms for Structured Grids
The Design of Flux-Corrected Transport (FCT) Algorithms for Structured Grids Steven T. Zalesak Abstract A given flux-corrected transport (FCT) algorithm consists of three components: (1) a high order algorithm
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
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 informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationSome recent results on positivity-preserving discretisations for the convection-diffusion equation
Some recent results on positivity-preserving discretisations for the convection-diffusion equation Gabriel R. Barrenechea 1, Volker John 2 & Petr Knobloch 3 1 Department of Mathematics and Statistics,
More informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationMath 5520 Homework 2 Solutions
Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas
More information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationTutorials in Optimization. Richard Socher
Tutorials in Optimization Richard Socher July 20, 2008 CONTENTS 1 Contents 1 Linear Algebra: Bilinear Form - A Simple Optimization Problem 2 1.1 Definitions........................................ 2 1.2
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 informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
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 informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationBoundary conditions. Diffusion 2: Boundary conditions, long time behavior
Boundary conditions In a domain Ω one has to add boundary conditions to the heat (or diffusion) equation: 1. u(x, t) = φ for x Ω. Temperature given at the boundary. Also density given at the boundary.
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationFinite Element Multigrid Framework for Mimetic Finite Difference Discretizations
Finite Element Multigrid Framework for Mimetic Finite ifference iscretizations Xiaozhe Hu Tufts University Polytopal Element Methods in Mathematics and Engineering, October 26-28, 2015 Joint work with:
More informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
More informationVortices in Superfluid MODD-Problems
Vortices in Superfluid MODD-Problems May 5, 2017 A. Steady filament (0.75) Consider a cylindrical beaker (radius R 0 a) of superfluid helium and a straight vertical vortex filament in its center Fig. 2.
More informationQUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES
1 QUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES O. ALLAIN(*), A. DERVIEUX(**), I. ABALAKIN(***), S. CAMARRI(****), H. GUILLARD(**), B. KOOBUS(*****),
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More informationNo. 588 June Frame-invariant directional vector limiters for discontinuous Galerkin methods. H. Hajduk, D. Kuzmin, V. Aizinger ISSN:
No. 588 June 2018 Frame-invariant directional vector limiters for discontinuous Galerkin methods H. Hajduk, D. Kuzmin, V. Aizinger ISSN: 2190-1767 Frame-invariant directional vector limiters for discontinuous
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
More informationThermal Detection of Small Cracks in a Plate
1 Rose-Hulman Institute of Technology May 20, 2011 1 With some helpful discussions with Chris Earls, Cornell University School of Civil and Environmental Engineering. 1 2 3 4 Motivation The USS Independence
More informationAngular momentum preserving CFD on general grids
B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration
More informationarxiv: v1 [math.na] 11 Jun 2018
High-order residual distribution scheme for the time-dependent Euler equations of fluid dynamics arxiv:1806.03986v1 [math.na] 11 Jun 2018 R. Abgrall, P. Bacigaluppi, S. Tokareva Institute of Mathematics,
More informationThe Janus Face of Turbulent Pressure
CRC 963 Astrophysical Turbulence and Flow Instabilities with thanks to Christoph Federrath, Monash University Patrick Hennebelle, CEA/Saclay Alexei Kritsuk, UCSD yt-project.org Seminar über Astrophysik,
More informationExtremum-Preserving Limiters for MUSCL and PPM
arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationBasic hydrodynamics. David Gurarie. 1 Newtonian fluids: Euler and Navier-Stokes equations
Basic hydrodynamics David Gurarie 1 Newtonian fluids: Euler and Navier-Stokes equations The basic hydrodynamic equations in the Eulerian form consist of conservation of mass, momentum and energy. We denote
More informationStabilization of sawteeth in tokamaks with toroidal flows
PHYSICS OF PLASMAS VOLUME 9, NUMBER 7 JULY 2002 Stabilization of sawteeth in tokamaks with toroidal flows Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationPhysics 202: Spring 1999 Solution to Homework Assignment #3
Physics 202: Spring 1999 Solution to Homework Assignment #3 Questions: Q3. (a) The net electric flux through each surface shown is zero, since every electric field line entering from one end exits through
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
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 informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationProblem Set #3: 2.11, 2.15, 2.21, 2.26, 2.40, 2.42, 2.43, 2.46 (Due Thursday Feb. 27th)
Chapter Electrostatics Problem Set #3:.,.5,.,.6,.40,.4,.43,.46 (Due Thursday Feb. 7th). Coulomb s Law Coulomb showed experimentally that for two point charges the force is - proportional to each of the
More informationMathematics 22: Lecture 10
Mathematics 22: Lecture 10 Euler s Method Dan Sloughter Furman University January 22, 2008 Dan Sloughter (Furman University) Mathematics 22: Lecture 10 January 22, 2008 1 / 14 Euler s method Consider the
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
More informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
More informationThermoacoustic Devices
Thermoacoustic Devices Peter in t panhuis Sjoerd Rienstra Han Slot 24th April 2008 Down-well power generation Vortex shedding in side branch Vortex shedding creates standing wave Porous medium near closed
More informationIntroduction to the Calculus of Variations
236861 Numerical Geometry of Images Tutorial 1 Introduction to the Calculus of Variations Alex Bronstein c 2005 1 Calculus Calculus of variations 1. Function Functional f : R n R Example: f(x, y) =x 2
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationSolutions to PS 2 Physics 201
Solutions to PS Physics 1 1. ke dq E = i (1) r = i = i k eλ = i k eλ = i k eλ k e λ xdx () (x x) (x x )dx (x x ) + x dx () (x x ) x ln + x x + x x (4) x + x ln + x (5) x + x To find the field for x, we
More informationA Numerical Study of Boson Star Binaries
A Numerical Study of Boson Star Binaries Bruno C. Mundim with Matthew W. Choptuik (UBC) 12th Eastern Gravity Meeting Center for Computational Relativity and Gravitation Rochester Institute of Technology
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationProcedural Animation. D.A. Forsyth
Procedural Animation D.A. Forsyth Big points Two important types of procedural animation slice and dice data, like texture synthesis build (approximate) physical simulation Extremely powerful issues how
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationHydrodynamic Modes of Incoherent Black Holes
Hydrodynamic Modes of Incoherent Black Holes Vaios Ziogas Durham University Based on work in collaboration with A. Donos, J. Gauntlett [arxiv: 1707.xxxxx, 170x.xxxxx] 9th Crete Regional Meeting on String
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationLinear Equations and Matrix
1/60 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Gaussian Elimination 2/60 Alpha Go Linear algebra begins with a system of linear
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationIterative Methods for Ax=b
1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative
More informationCHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD
CHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD For the second chapter in the thesis, we start with surveying the mathematical background which is used directly in the Boundary Element
More informationPeter Hertel. University of Osnabrück, Germany. Lecture presented at APS, Nankai University, China.
Balance University of Osnabrück, Germany Lecture presented at APS, Nankai University, China http://www.home.uni-osnabrueck.de/phertel Spring 2012 Linear and angular momentum and First and Second Law point
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationLecturer: Bengt E W Nilsson
9 3 19 Lecturer: Bengt E W Nilsson Last time: Relativistic physics in any dimension. Light-cone coordinates, light-cone stuff. Extra dimensions compact extra dimensions (here we talked about fundamental
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationDiscretization Techniques for Surface Tension Forces in Fully Conservative Two-Phase Flow Finite Volume Methods
Discretization Techniques for Surface Tension Forces in Fully Conservative Two-Phase Flow Finite Volume Methods Matthias Waidmann 1, Stephan Gerber 2 Michael Oevermann 3, Rupert Klein 1 1 Freie Universität
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More information