Magnetized Target Fusion: Insights from Mathematical Modelling
|
|
- Teresa Ashlie Harmon
- 5 years ago
- Views:
Transcription
1 : Insights from Mathematical Modelling, UBC, PhD Candidate - Applied Math Brian Wetton, UBC, Supervisor Michael Ward, UBC, Committee Member Sandra Barsky, General Fusion, Committee Member November 28, 2014
2 Context A local Canadian research company is designing a reactor to produce fusion energy, which has never been achieved before. Success with this could hold immense promise for clean and sustainable energy sources in the future. In this talk, we will look into the mathematical analysis used to yield insights into the design work. Key components: Reactor model Numerical finite volume approach Formal asymptotic approach
3 Thanks to... Thanks to the following people for helpful discussions, collaboration, and support: Sandra Barsky and Aaron Froese, General Fusion George Bluman, Michael Ward, and Brian Wetton, UBC Math Department Randy LeVeque, University of Washington Applied Mathematics Department My colleagues and fellow graduate students
4 Fusion Fusing atomic nuclei yield new nuclei plus energy Occurs naturally in sun: pressure 300 G atm, temperature 10 7 K Potential clean energy source Tritium-deuterium plasma reaction: 3 1 H +2 1 H 4 2 He +n MeV of energy Depends on temperature, thermal distribution, etc.
5 Challenges Immense temperature and pressure must be sustained Thermal and radiative losses Lawson criterion 1 for energy yield: density temperature time cm 3 KeV s
6 General Fusion design General Fusion founded by Michel Laberge in 2002 Design magnetized target fusion reactor Isolate spheromak plasma with magnetic field Crush plasma in imploding metal cavity
7 General Fusion design How feasible is such a design?
8 Lead-Lithium Euler equations for mass and momentum conservation With density ρ, velocity v, and pressure P: ρ t + (ρv) = 0 (ρv) t + (ρv v) + P = 0 P = P(ρ) (mass) (momentum) (equation of state) Nonlinear coupled system of PDEs
9 Lead-Lithium Quadratic fit to high-pressure Lead experiments 2
10 Piston and plasma pressure Gaussian piston pressure: P(t) = P atmospheric + (P impact P atmospheric )e t2 /t 2 0 θ(t) Plasma has gas and magnetic pressure Ideal gas of volume V : Magnetic pressure 3 : P gas V γ, γ 1 P magnetic R 4 major
11 Assumptions and simplifications Spherical symmetry Isentropic (reversible) conditions so γ = 5/3 Initially P gas 0.1P magnetic Initial plasma temperature 100 ev System starts in equilibrium No mixing implies free boundary problem at plasma-metal and metal-piston interfaces: d dt r boundary(t) = v(r boundary (t), t) Neglect thermal and radiative energy losses
12 Overall model
13 Overall model In r L (t) < r < r R (t), t > 0, dimensionless system has form: ρ t + 1 r 2 (r 2 ρv) r = 0, (ρv) t + p r + 1 r 2 (r 2 ρv 2 ) r = 0 p = p(ρ), p(r L (t), t) = p L (r L (t)), v(r, ) = 0, dr L,R = v(r L,R (t), t) dt p(r R (t), t) = f (t) p(r, ) = p L (r L ( )) r L ( ) = χɛ 1/2, r R ( ) = 1
14 First-order finite volumes Given u t + (f (u)) x = 0 denote ū(x i, t) = 1 xi +h/2 h x i h/2 u(x, t)dx dū(x i,t) dt = 1 xi +h/2 h x i h/2 u t(x, t)dx = f (u(x i+1/2,t)) f (u(x i 1/2,t)) h Choose F L (u i 1, u i ) = f (u i 1 ) or f (u i ) depending on f Stable time-step k < h/ sup f, u j+1 i = u j i k h (F L (u j i, uj i+1 ) F L (u j i 1, uj i ))
15 Second-order flux-limited finite volume idea Smoothness: u t + f (u)u x = 0 Use u tt = f (u)u t u x f (u)u tx = (f (u) 2 ) x u x + f (u) 2 u xx Lax-Wendroff: u(x, t + k) = u(x, t) + ku t (x, t) + k 2 2 u tt(x, t) + O(k 3 ) Corresponding high-resolution flux F H Flux limited: take flux F = φf L + (1 φ)f H Limiter φ largest where derivative u (x) largest
16 Convergence Numerical solution u num (x, t) and exact solution u ex (x, t): E = u num (x, t) u ex (x, t) dx 0, as N Near shocks, E = O(h) Near smooth regions E = O(h n ) for n th -order method
17 Coordinate transformation Transform from moving to fixed computational domain Let τ = t and define (τ) = r R (τ) r L (τ), Γ(τ) = v R (τ) v L (τ) Set y = r r L(τ) (τ), y [0, 1] p(r R (t), t) = p(1, τ), etc. New conservation laws e.g. mass: ρ τ + { 1 [ (v L + Γy)ρ + ρv]} y = 2 r L + y Γ ρ
18 Finite volume overview Pseudocode 4 : - eigen analysis for system - constant extrapolation to ghost-points - use projections for limiters - homogeneous Riemann problem then use source - update time and boundary data Verify L 1 -convergence Table: Convergence of Numerical Scheme at fixed time N Velocity Error Density Error Mass Error Rate
19 Pulse profiles
20 Abridged sensitivity analysis Table: Min radius R min, Lawson triple product Π L, impact pressure P impact, initial plasma radius R plasma,0, initial sphere radius R lead,0, time scale t 0. System R min (cm) Π L (10 15 kev s cm 3 ) Baseline R plasma, P impact R lead, t P impact R lead, t
21 Big picture Modifications: - Gaussian impulse pressure without heaviside - linearized P(ρ) equation of state - only magnetic pressure Relevant scales: - asymptotic parameter ɛ 1 - big sound speed: bɛ 1/2 - impulse pressure: O(1) - small initial plasma radius: χɛ 1/2 - very small impulse time: O(ɛ) - very, very small initial pressure: µɛ 3/2 Matched asymptotics 5
22 Asymptotic regimes
23 Pulse formation and focusing Formation: - ρ(y, τ) 1 + ɛρ 1 + ɛ 3/2 ρ 2, v(y, τ) ɛ 1/2 v 0 + ɛv 1 - Riemann Invariants 6 - plane wave solutions in (, 0] [0, ): ρ 1 (y, τ) = 1 b 2 e (τ+y/b)2, - need divergent ρ 2, v 1 for matching Focusing: - amplitude growth: v 0 (y, τ) = 1 2 b e (τ+y/b) ρ 1 + ɛρ 1 + ɛ 3/2 ρ 2, v ɛ 1/2 v 0 + ɛv 1 (outer) ρ 1 + ɛ 1/2 ρ 1 + ɛρ 2, v v 0 + ɛ 1/2 v 1 (inner) - linear acoustic limit 7 : ρ 1,ˆt + v 0,σ + 2 σ v 0 = 0, v 0,ˆt + b2 ρ 1,σ = 0
24 Pulse reflection Complete reflection of ρ 1 and v 0, negligible compression Next order: v 1,ˆt + b2 ρ 2,σ = (ρ 1 v 0 )ˆt (v 2 0 ) σ 2 σ v 2 0 Careful integration and boundary condition trickery gives v 1 (σ, ) = 2π b 2 χσ 2 Residual velocity field v ɛ 1/2 v 1 compresses plasma! Rayleigh energy argument 8 or more dynamics
25 Plasma compression Integration of PDEs 9 with reformulation gives free boundary ODE Initially - ρ 1 + ɛ 2 ˆρ, v ɛ 1/2ˆv - O(ɛ 1/2 ) space scale, O(1) time scale - find ˆv(σ L ) = ˆv L = 2π ; inner wall at σ L b 2 χ 3/2 σ 3/2 L Finally - ρ 1 + ɛ 1/2 ˆρ, v ɛ 1/4ˆv - O(ɛ) space scale, O(ɛ 5/4 ) time scale - find ˆv(z L ) = ˆv L = 2Az L 2µ ; inner wall at z zl 2 L Matching gives A Turnaround point when ˆv L = 0
26 Minimum radius Minimum radius: r min b4 χ 3 µ ɛ π Agreement with numerics as ɛ 0 with b = 1.05, χ = 0.937, µ = π: ɛ r min,num r min,asy
27 Key insights Dimensional minimum radius R min C4 s P plasma,0 R 7 plasma,0 ϱ3 0 πp 4 impact R4 lead,0 t2 0 = 1.6 cm Symbol Meaning Symbol Meaning C s lead sound speed P plasma,0 initial pressure R plasma,0 initial plasma radius ϱ 0 lead density P impact piston pressure R lead,0 initial lead radius t 0 impulse time scale Qualitative consistency with numerics
28 Key insights consistent compression profile wave-like behaviour: sound speed dominates almost all input energy reflected: E input 8π 3 b ɛ3/2, E compression 4π2 b 4 χ 3 ɛ5/2
29 Results and future work Results: - numerics and asymptotics consistent - sensitivity to parameters, energy yield may be within reach - larger outer sphere radius and impact pressure Future directions: - incorporate more plasma physics - consider angular instabilities - more precise assessment of design 1 Lawson, Some Criteria for a Power..., Proc. Phys Society B, Rothman et al., Measurement of principal isentropes..., J. Phys D, Woodruff et al., Adiabatic compression of a doublet field..., JOFE, LeVeque, Finite Volume Methods for Hyperbolic Problems, Hinch, Perturbation Methods, Whitham, Linear and Nonlinear Waves, Landau and Lifschitz, Fluid Mechanics, Lord Rayleigh, On the Pressure Developed..., Phil. Mag. Series 6, Brennan, Cavitation and Bubble Dynamics, 1995.
Mathematical Modelling of a Magnetized Target Fusion Reactor
Mathematical Modelling of a Magnetized Target Fusion Reactor Michael PIC Assistant Adjunct Professor, Mathematics UCLA April 20, 2016 Context A Canadian research company General Fusion (2002-Present, Michel
More informationFusion: Intro Fusion: Numerics and Asymptotics Fusion: Summary Superconductor: Problem Superconductor: Results
Investigation into the Feasibility and Operation of a Magnetized Target Fusion Reactor, and Qualitative Predictions of Magnetic Field Profile Perturbations Induced by Surface Roughness in Type II Superconductors,
More informationOctober 22, :00-16:00. web: math.ubc.ca/ MLRTLM. Math: the art in problem solving. Michael Lindstrom.
October 22, 2014 15:00-16:00 email: MLRTLM@math.ubc.ca web: math.ubc.ca/ MLRTLM Math is everywhere quantifiable descriptions glacier motion, spreading of diseases, market inflation, digitally recorded
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationASYMPTOTIC ESTIMATION FOR MINIMAL PLASMA RADIUS IN A SPHERICALLY SYMMETRIC MAGNETIZED TARGET FUSION REACTOR MODEL
ASYMPTOTIC ESTIMATION FOR MINIMAL PLASMA RADIUS IN A SPHERICALLY SYMMETRIC MAGNETIZED TARGET FUSION REACTOR MODEL MICHAEL LINDSTROM Astract. Nonlinear conservation laws with moving oundaries often arise
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More informationD-D FUSION NEUTRONS FROM A STRONG SPHERICAL SHOCK WAVE FOCUSED ON A DEUTERIUM BUBBLE IN WATER. Dr. Michel Laberge General Fusion Inc.
D-D FUSION NEUTRONS FROM A STRONG SPHERICAL SHOCK WAVE FOCUSED ON A DEUTERIUM BUBBLE IN WATER Dr. Michel Laberge General Fusion Inc. SONOFUSION Sonofusion is making some noise A bit short in energy, ~mj
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
More informationWhere are we with laser fusion?
Where are we with laser fusion? R. Betti Laboratory for Laser Energetics Fusion Science Center Dept. Mechanical Engineering and Physics & Astronomy University of Rochester HEDSA HEDP Summer School August
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationOn The Collapse of a Gas Cavity by an Imploding Molten Lead Shell and Richtmyer-Meshkov Instability
On The Collapse of a Gas Cavity by an Imploding Molten Lead Shell and Richtmyer-Meshkov Instability Victoria Suponitsky, Sandra Barsky, and Aaron Froese General Fusion Inc., 8-368 Bonneville Place, Burnaby,
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 informationFluid Dynamics from Kinetic Equations
Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationWebster s horn model on Bernoulli flow
Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationTarget Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu
Muon Collider/Neutrino Factory Collaboration Meeting May 26 28, CERN, Geneva U.S. Department of Energy Target Simulations Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Center for Data Intensive
More informationDifferential criterion of a bubble collapse in viscous liquids
PHYSICAL REVIEW E VOLUME 60, NUMBER 1 JULY 1999 Differential criterion of a bubble collapse in viscous liquids Vladislav A. Bogoyavlenskiy* Low Temperature Physics Department, Moscow State University,
More informationNotes on fusion reactions and power balance of a thermonuclear plasma!
SA, 3/2017 Chapter 5 Notes on fusion reactions and power balance of a thermonuclear plasma! Stefano Atzeni See S. Atzeni and J. Meyer-ter-Vehn, The Physics of Inertial Fusion, Oxford University Press (2004,
More informationIssues in Neoclassical Tearing Mode Theory
Issues in Neoclassical Tearing Mode Theory Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX Tearing Mode Stability in Tokamaks According to standard (single-fluid)
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationGraduate Written Examination Fall 2014 Part I
Graduate Written Examination Fall 2014 Part I University of Minnesota School of Physics and Astronomy Aug. 19, 2014 Examination Instructions Part 1 of this exam consists of 10 problems of equal weight.
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 informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationIntroduction to Fusion Physics
Introduction to Fusion Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 Energy from nuclear fusion Reduction
More informationTHE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS
Journal of Hyperbolic Differential Equations Vol., No. 4 005 909 917 c World Scientific Publishing Company THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS VOLKER ELLING, and TAI-PING LIU, Dept.
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
More informationConservation Laws of Surfactant Transport Equations
Conservation Laws of Surfactant Transport Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Winter 2011 CMS Meeting Dec. 10, 2011 A. Cheviakov
More informationPhysics Dec Time Independent Solutions of the Diffusion Equation
Physics 301 10-Dec-2004 33-1 Time Independent Solutions of the Diffusion Equation In some cases we ll be interested in the time independent solution of the diffusion equation Why would be interested in
More informationClawpack Tutorial Part I
Clawpack Tutorial Part I Randall J. LeVeque Applied Mathematics University of Washington Conservation Laws Package www.clawpack.org (pdf s will be posted and green links can be clicked) Some collaborators
More informationGA A22689 SLOW LINER FUSION
GA A22689 SLOW LINER FUSION by AUGUST 1997 DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any
More informationCharacteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:
AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,
More informationINTRODUCTION TO MAGNETIC NUCLEAR FUSION
INTRODUCTION TO MAGNETIC NUCLEAR FUSION S.E. Sharapov Euratom/CCFE Fusion Association, Culham Science Centre, Abingdon, Oxfordshire OX14 3DB, UK With acknowledgments to B.Alper for use of his transparencies
More informationModels of Ultrasound Contrast Agents. Duncan Alexander Sutherland
Models of Ultrasound Contrast Agents Duncan Alexander Sutherland An essay submitted in partial fulfillment of the requirements for the degree of B.Sc. Honours) Applied Mathematics University of Sydney
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationModelling a Generalized Dialysis Device
Modelling a Generalized Dialysis Device Brian Wetton Mathematics Department University of British Columbia www.math.ubc.ca/~wetton UPC, November 11, 2015 Overview Advertising UBC Institute of Applied Mathematics
More informationTheory and simulations of hydrodynamic instabilities in inertial fusion
Theory and simulations of hydrodynamic instabilities in inertial fusion R. Betti Fusion Science Center, Laboratory for Laser Energetics, University of Rochester IPAM/UCLA Long Program PL2012 - March 12
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
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 informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationTwo Phase Transport in Porous Media
Two Phase Transport in Porous Media Lloyd Bridge Iain Moyles Brian Wetton Mathematics Department University of British Columbia www.math.ubc.ca/ wetton CRM CAMP Seminar, October 19, 2011 Overview Two phase
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationA Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More informationComputational Astrophysics 7 Hydrodynamics with source terms
Computational Astrophysics 7 Hydrodynamics with source terms Oscar Agertz Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term
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 informationMath 575-Lecture 26. KdV equation. Derivation of KdV
Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not
More informationTwo-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation
Mathematics and Computers in Simulation 76 (2007) 3 7 Two-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation S. Baboolal a,, R. Bharuthram
More informationCompressible Duct Flow with Friction
Compressible Duct Flow with Friction We treat only the effect of friction, neglecting area change and heat transfer. The basic assumptions are 1. Steady one-dimensional adiabatic flow 2. Perfect gas with
More informationWaves and characteristics: Overview 5-1
Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables
More informationTokamak Fusion Basics and the MHD Equations
MHD Simulations for Fusion Applications Lecture 1 Tokamak Fusion Basics and the MHD Equations Stephen C. Jardin Princeton Plasma Physics Laboratory CEMRACS 1 Marseille, France July 19, 21 1 Fusion Powers
More information1 Introduction to PDE MATH 22C. 1. Introduction To Partial Differential Equations Recall: A function f is an input-output machine for numbers:
1 Introduction to PDE MATH 22C 1. Introduction To Partial Differential Equations Recall: A function f is an input-output machine for numbers: y = f(t) Output y 2R Input t 2R Name of function f t=independent
More informationNatalia Tronko S.V.Nazarenko S. Galtier
IPP Garching, ESF Exploratory Workshop Natalia Tronko University of York, York Plasma Institute In collaboration with S.V.Nazarenko University of Warwick S. Galtier University of Paris XI Outline Motivations:
More informationgeneralfusion Characterization of General Fusion's Plasma Devices 2015 Nimrod Summer Workshop
Characterization of General Fusion's Plasma Devices 2015 Nimrod Summer Workshop Aaron Froese, Charlson Kim, Meritt Reynolds, Sandra Barsky, Victoria Suponitsky, Stephen Howard, Russ Ivanov, Peter O'Shea,
More informationASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS
Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationPlasma Wall Interactions in Tokamak
Plasma Wall Interactions in Tokamak Dr. C Grisolia, Association Euratom/CEA sur la fusion, CEA/Cadarache Outline 1. Conditions for Fusion in Tokamaks 2. Consequences of plasma operation on in vessel materials:
More informationQuasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths
in the presence of plasma sheaths Department of Mathematical Sciences Ulsan National Institute of Science and Technology (UNIST) joint work with Masahiro Suzuki (Nagoya) and Chang-Yeol Jung (Ulsan) The
More informationA Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007
A Study of Transonic Flow and Airfoils Presented by: Huiliang Lui 3 th April 7 Contents Background Aims Theory Conservation Laws Irrotational Flow Self-Similarity Characteristics Numerical Modeling Conclusion
More informationMAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA
University of Ljubljana Faculty of Mathematics and Physics Seminar 1 b -1st year, II. cycle MAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA Author: Lino alamon Advisor: prof. dr. Tomaº Gyergyek
More informationVIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition
VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationImposed Dynamo Current Drive
Imposed Dynamo Current Drive by T. R. Jarboe, C. Akcay, C. J. Hansen, A. C. Hossack, G. J. Marklin, K. Morgan, B. A. Nelson, D. A. Sutherland, B. S. Victor University of Washington, Seattle, WA 98195,
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
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 informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationNotes #4a MAE 533, Fluid Mechanics
Notes #4a MAE 533, Fluid Mechanics S. H. Lam lam@princeton.edu http://www.princeton.edu/ lam October 23, 1998 1 The One-dimensional Continuity Equation The one-dimensional steady flow continuity equation
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationComparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models
Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Julian Koellermeier, Manuel Torrilhon October 12th, 2015 Chinese Academy of Sciences, Beijing Julian Koellermeier,
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationPhysics & Engineering Physics University of Saskatchewan. Supported by NSERC, CRC
Fusion Energy Chijin Xiao and Akira Hirose Plasma Physics laboratory Physics & Engineering Physics University of Saskatchewan Supported by NSERC, CRC Trends in Nuclear & Medical Technologies April il6-7,
More informationTurbulence and Transport The Secrets of Magnetic Confinement
Turbulence and Transport The Secrets of Magnetic Confinement Presented by Martin Greenwald MIT Plasma Science & Fusion Center IAP January 2005 FUSION REACTIONS POWER THE STARS AND PRODUCE THE ELEMENTS
More informationICF ignition and the Lawson criterion
ICF ignition and the Lawson criterion Riccardo Betti Fusion Science Center Laboratory for Laser Energetics, University of Rochester Seminar Massachusetts Institute of Technology, January 0, 010, Cambridge
More informationElectrodynamics Qualifier Examination
Electrodynamics Qualifier Examination August 15, 2007 General Instructions: In all cases, be sure to state your system of units. Show all your work, write only on one side of the designated paper, and
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationANALYSIS OF A NONLINEAR SURFACE WIND WAVES MODEL VIA LIE GROUP METHOD
Electronic Journal of Differential Equations, Vol. 206 (206), No. 228, pp. 8. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS OF A NONLINEAR SURFACE WIND WAVES MODEL
More informationof the Schnakenberg model
Pulse motion in the semi-strong limit of the Schnakenberg model Newton Institute 2005 Jens Rademacher, Weierstraß Institut Berlin joint work with Michael Ward (UBC) Angelfish 2, 6, 12 months old [Kondo,
More informationModel adaptation in hierarchies of hyperbolic systems
Model adaptation in hierarchies of hyperbolic systems Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France February 15th, 2012 DFG-CNRS Workshop Nicolas Seguin (LJLL, UPMC) 1 / 29 Outline of the
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationMHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field
MHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field Y. Matsui, T. Yokoyama, H. Hotta and T. Saito Department of Earth and Planetary Science, University of Tokyo,
More informationInitial Boundary Value Problems for Scalar and Vector Burgers Equations
Initial Boundary Value Problems for Scalar and Vector Burgers Equations By K. T. Joseph and P. L. Sachdev In this article we stu Burgers equation and vector Burgers equation with initial and boundary conditions.
More information7 The General Circulation
7 The General Circulation 7.1 The axisymmetric state At the beginning of the class, we discussed the nonlinear, inviscid, axisymmetric theory of the meridional structure of the atmosphere. The important
More informationIncompressible MHD simulations
Incompressible MHD simulations Felix Spanier 1 Lehrstuhl für Astronomie Universität Würzburg Simulation methods in astrophysics Felix Spanier (Uni Würzburg) Simulation methods in astrophysics 1 / 20 Outline
More informationThe Stochastic Piston Problem
The Stochastic Piston Problem G. Lin, C.-H. Su and G.E. Karniadakis Division of Applied Mathematics 18 George Street Brown University Providence, RI 91 Classification: Physical Sciences: Applied Mathematics
More informationContents of lecture 2b. Lectures 2a & 2b. Physical vs. computational coordinates [2] Physical vs. computational coordinates [1]
Contents of lecture b Lectures a & b P. A. Zegeling Mathematical Institute Utrecht University The Netherlands Parameter-free non-singular moving grids in D: Theory & properties Application to resistive
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationHervé Guillard, INRIA Projet Smash, B.P. 93, Sophia-Antipolis Cedex, France,
TRAVELING WAVE ANALYSIS OF TWO-PHASE DISSIPATIVE MODELS Hervé Guillard, INRIA Projet Smash, B.P. 93, 06902 Sophia-Antipolis Cedex, France, Herve.Guillard@sophia.inria.fr Joint work with : Mathieu Labois,
More information