Solving a magnetic diffusion equation by the Mixed Hybrid Finite Element method.

Size: px
Start display at page:

Download "Solving a magnetic diffusion equation by the Mixed Hybrid Finite Element method."

Transcription

1 Solving a magnetic diffusion equation by the Mixed Hybrid Finite Element method. Corinne Aymard, Christophe Fochesato CEA, DAM, DIF, F Arpajon CEA/DAM/DIF 6 juillet

2 Outline Objectives Model Discretization Results Conclusion CEA/DAM/DIF 6 juillet

3 Outline Objectives Model Discretization Results, soon Conclusion CEA/DAM/DIF 6 juillet

4 Objectives to add a MHD model in an existing multimaterial lagrangian hydrodynamical code to study the behaviour of materials under electromagnetic conditions magnetoforming zpinch type applications CEA/DAM/DIF 6 juillet

5 Model CEA/DAM/DIF 6 juillet

6 Model Hydrodynamics Euler s equations Electromagnetics Maxwell s equations CEA/DAM/DIF 6 juillet

7 Model classical hypotheses for MHD => convective/magnetic diffusion equation modelling choice: 2D configurations axisymmetrical flow with azimutal magnetic field B θ transverse plane flow with axial current B z transverse plane flow with axial magnetic field j z splitting in time: hydro / convective evolution of B / diffusive evolution of B CEA/DAM/DIF 6 juillet

8 Model: diffusive part axisymmetrical flow with azimutal magnetic field transverse plane flow with axial current transverse plane flow with axial magnetic field + global constraint by material for fixing temporal function Ec(t) CEA/DAM/DIF 6 juillet

9 Model: generic diffusion equation Diffusive system for the different configurations + initial conditions + Dirichlet (imposed q) and Neumann boundary conditions (imposed X.n) + eventually a global constraint on q properties of the diffusion equation discontinuous coefficients large ratios between different materials coefficients distorted mesh by hydrodynamical flow CEA/DAM/DIF 6 juillet

10 Discretization CEA/DAM/DIF 6 juillet

11 Discretization: Finite Element Method Weak form of the diffusive system: looking for and on with on Approximation by Mixed Hybrid Finite Element in in Ω m Г a CEA/DAM/DIF 6 juillet

12 Discretization: MHFE Discrete weak form of the diffusive system with local boundaries conditions if internal edge if Neumann boundary for any if Dirichlet boundary Degrees of freedom X.n q λ CEA/DAM/DIF 6 juillet

13 Discretization: basis functions for RT0 square reference element 1 η -1 1 ξ basis functions on the reference element -1 Piola transform CEA/DAM/DIF 6 juillet

14 Discretization: equation for the fluxes in each element with mass matrix we get CEA/DAM/DIF 6 juillet

15 Discretization: equation for q discrete divergence obtained from Green s formula in a matrix-like writing CEA/DAM/DIF 6 juillet

16 Discretization: local system local system for each element with inter-element continuity for edges not in CEA/DAM/DIF 6 juillet

17 Discretization: elimination of the fluxes local inversion of the mass matrix eliminating fluxes unknowns X m if if we get the system m edge mk k CEA/DAM/DIF 6 juillet

18 Discretization: time scheme and matrix system with implicit Euler s time scheme Diagonal system for 1st equation: easily inverted => system in Λ to be solved with the matrix is definite positive solved by a preconditioned conjugate gradient method CEA/DAM/DIF 6 juillet

19 Solving: bordering algorithm for the constraint the constraint is discretized use of the bordering algorithm for an arrow matrix with A V = a C and A W = B A (sparse) a C X = B we get : ( a a R V ) x = ( b a R W ) a R a x b X = W V x CEA/DAM/DIF 6 juillet

20 Coupling with hydrodynamics Euler s equations u Q1-Q0 approach + leapfrog in time ρ, p, e CEA/DAM/DIF 6 juillet

21 adding MHD coupling terms Euler s equations in Ω Laplace Laplace force needed at nodes Joule Joule effect needed at center CEA/DAM/DIF 6 juillet

22 Coupling terms using all the information from the MHFE discretization of the diffusion equation Laplace Force deduced from the weak form of the definition of j for example : the axisymmetrical case Joule effect CEA/DAM/DIF 6 juillet

23 Conclusion In the modelling context, need for a generic planar 2D diffusion solver 3 configurations global constraint eventually Discretization of electromagnetic part in a Lagrangian hydrodynamical code Mixed Finite Element Method for the diffusion equation Euler s implicit scheme in time preconditionned conjugate gradient method for the linear system same discretization for coupling terms : Laplace s Force and Joule s effect bordering algorithm for the additional constraint hoping for results soon if problem with distorted quadrilateral meshes use of modified RT element by Boffi & Gastaldi 2009 : adding a quadratic function to the basis functions CEA/DAM/DIF 6 juillet

ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS

ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS VINCENT FONTAINE 1,2 AND ANIS YOUNES 2 1 Laboratoire de Physique des Bâtiments et des Systèmes, Université de La Réunion, 15 avenue

More information

Maxwell s equations. based on S-54. electric field charge density. current density

Maxwell s equations. based on S-54. electric field charge density. current density Maxwell s equations based on S-54 Our next task is to find a quantum field theory description of spin-1 particles, e.g. photons. Classical electrodynamics is governed by Maxwell s equations: electric field

More information

Math background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids

Math background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure

More information

One dimensional steady state diffusion, with and without source. Effective transfer coefficients

One dimensional steady state diffusion, with and without source. Effective transfer coefficients One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport

More information

NUMERICAL METHODS IN ASTROPHYSICS An Introduction

NUMERICAL METHODS IN ASTROPHYSICS An Introduction -1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,

More information

Numerical Analysis of Electromagnetic Fields

Numerical Analysis of Electromagnetic Fields Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts

More information

PROBLEM SET. Heliophysics Summer School. July, 2013

PROBLEM 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 information

A local-structure-preserving local discontinuous Galerkin method for the Laplace equation

A local-structure-preserving local discontinuous Galerkin method for the Laplace equation A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous

More information

Macroscopic plasma description

Macroscopic plasma description Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion

More information

The Convergence of Mimetic Discretization

The Convergence of Mimetic Discretization The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University

More information

A Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers

A 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 information

Finite Element Solver for Flux-Source Equations

Finite Element Solver for Flux-Source Equations Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University

More information

Computations with Discontinuous Basis Functions

Computations with Discontinuous Basis Functions Computations with Discontinuous Basis Functions Carl Sovinec University of Wisconsin-Madison NIMROD Team Meeting November 12, 2011 Salt Lake City, Utah Motivation The objective of this work is to make

More information

7.4 The Saddle Point Stokes Problem

7.4 The Saddle Point Stokes Problem 346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler

More information

ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow

ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow Michel Kern Inria Paris Rocquencourt Maison de la Simulation C2S@Exa Days, Inria Paris Centre, Novembre 2016 M. Kern

More information

Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement

Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic

More information

A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere

A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background

More information

Partial Differential Equations II

Partial Differential Equations II Partial Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Partial Differential Equations II 1 / 28 Almost Done! Homework

More information

Solving PDEs with Multigrid Methods p.1

Solving PDEs with Multigrid Methods p.1 Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration

More information

Partial Differential Equations. Examples of PDEs

Partial Differential Equations. Examples of PDEs Partial Differential Equations Almost all the elementary and numerous advanced parts of theoretical physics are formulated in terms of differential equations (DE). Newton s Laws Maxwell equations Schrodinger

More information

BETI for acoustic and electromagnetic scattering

BETI for acoustic and electromagnetic scattering BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary

More information

Recovery-Based A Posteriori Error Estimation

Recovery-Based A Posteriori Error Estimation Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements

More information

Maxwell s equations. electric field charge density. current density

Maxwell s equations. electric field charge density. current density Maxwell s equations based on S-54 Our next task is to find a quantum field theory description of spin-1 particles, e.g. photons. Classical electrodynamics is governed by Maxwell s equations: electric field

More information

Chapter 4. Two-Dimensional Finite Element Analysis

Chapter 4. Two-Dimensional Finite Element Analysis Chapter 4. Two-Dimensional Finite Element Analysis general two-dimensional boundary-value problems 4.1 The Boundary-Value Problem 2nd-order differential equation to consider α α β φ Ω (4.1) Laplace, Poisson

More information

NIMEQ: MHD Equilibrium Solver for NIMROD

NIMEQ: MHD Equilibrium Solver for NIMROD NIMEQ: MHD Equilibrium Solver for NIMOD E.C.Howell, C..Sovinec University of Wisconsin-Madison 5 th Annual Meeting of Division of Plasma Physics Dallas, Texas, Nov. 17-Nov. 1,8 1 Abstract A Grad-Shafranov

More information

PDE Solvers for Fluid Flow

PDE Solvers for Fluid Flow PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic

More information

Chapter 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) 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 information

Preconditioned CG-Solvers and Finite Element Grids

Preconditioned CG-Solvers and Finite Element Grids Preconditioned CG-Solvers and Finite Element Grids R. Bauer and S. Selberherr Institute for Microelectronics, Technical University of Vienna Gusshausstrasse 27-29, A-1040 Vienna, Austria Phone +43/1/58801-3854,

More information

On two fractional step finite volume and finite element schemes for reactive low Mach number flows

On two fractional step finite volume and finite element schemes for reactive low Mach number flows Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes

More information

A Stable Spectral Difference Method for Triangles

A 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 information

Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD

Physic-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 information

Fluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17

Fluid 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 information

Discrete 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 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 information

Iterative methods for positive definite linear systems with a complex shift

Iterative methods for positive definite linear systems with a complex shift Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution

More information

Report Title Sharp Interface Algorithm for Large Density Ratio Incompressible Multiphase Magnetohydrodynamic Flows ABSTRACT A numerical algorithm and

Report Title Sharp Interface Algorithm for Large Density Ratio Incompressible Multiphase Magnetohydrodynamic Flows ABSTRACT A numerical algorithm and Report Title Sharp Interface Algorithm for Large Density Ratio Incompressible Multiphase Magnetohydrodynamic Flows ABSTRACT A numerical algorithm and the corresponding paralleled implementation for the

More information

SAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD

SAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD 1 SAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD Bruno Koobus, 1,2 Laurent Hascoët, 1 Frédéric Alauzet, 3 Adrien Loseille, 3 Youssef Mesri, 1 Alain Dervieux 1 1 INRIA

More information

J. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990

J. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990 Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME

More information

A 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 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

Multipole-Based Preconditioners for Sparse Linear Systems.

Multipole-Based Preconditioners for Sparse Linear Systems. Multipole-Based Preconditioners for Sparse Linear Systems. Ananth Grama Purdue University. Supported by the National Science Foundation. Overview Summary of Contributions Generalized Stokes Problem Solenoidal

More information

Waves in plasma. Denis Gialis

Waves 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 information

Approximation of fluid-structure interaction problems with Lagrange multiplier

Approximation of fluid-structure interaction problems with Lagrange multiplier Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline

More information

Fast Iterative Solution of Saddle Point Problems

Fast 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 information

Local discontinuous Galerkin methods for elliptic problems

Local discontinuous Galerkin methods for elliptic problems COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn

More information

Equations Introduction

Equations Introduction Chapter 9 90 Introduction Partial Differential Equations The numerical treatment of partial differential equations is, by itself, a vast subject Partial differential equations are at the heart of many,

More information

Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models

Applying 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 information

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C. Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference

More information

FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES

FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like

More information

On the choice of abstract projection vectors for second level preconditioners

On 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 information

Hybrid Simulation Method ISSS-10 Banff 2011

Hybrid Simulation Method ISSS-10 Banff 2011 Hybrid Simulation Method ISSS-10 Banff 2011 David Burgess Astronomy Unit Queen Mary University of London With thanks to Dietmar Krauss-Varban Space Plasmas: From Sun to Earth Space Plasma Plasma is (mostly)

More information

MHD flow and heat transfer due to a linearly stretching sheet. with induced magnetic field: Exact solution. Tarek M. A.

MHD flow and heat transfer due to a linearly stretching sheet. with induced magnetic field: Exact solution. Tarek M. A. MHD flow and heat transfer due to a linearly stretching sheet with induced magnetic field: Exact solution Tarek M. A. El-Mistikawy Dept. Eng. Math. & Phys., Faculty of Engineering, Cairo University, Giza

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed

More information

Electromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems

Electromagnetic 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 information

Lua, Boundary Conditions and ect..

Lua, Boundary Conditions and ect.. Lua, Boundary Conditions and ect.. Some further informations about calculating in Femm M. Rad 10.03.2017 Outline Lua scripting Border conditions About torque and force calculation in Femm Lua scripts To

More information

Parallelizing large scale time domain electromagnetic inverse problem

Parallelizing 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 information

A Finite Element Model for Numerical Analysis of Sintering

A Finite Element Model for Numerical Analysis of Sintering A Finite Element Model for Numerical Analysis of Sintering DANIELA CÂRSTEA High-School Group of Railways, Craiova ION CÂRSTEA Department of Computer Engineering and Communication University of Craiova

More information

A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations

A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National

More information

1 Energy dissipation in astrophysical plasmas

1 Energy dissipation in astrophysical plasmas 1 1 Energy dissipation in astrophysical plasmas The following presentation should give a summary of possible mechanisms, that can give rise to temperatures in astrophysical plasmas. It will be classified

More information

On the hydrodynamic diffusion of rigid particles

On the hydrodynamic diffusion of rigid particles On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:

More information

Fluid equations, magnetohydrodynamics

Fluid 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 information

Tdyn-CFD+HT - Validation Case 9

Tdyn-CFD+HT - Validation Case 9 Two-dimensional heat conduction with heat generation Version 14.0.0 Compass Ingeniería y Sistemas http://www.compassis.com Tel.: +34 932 181 989 - Fax.: +34 933 969 746 - E: info@compassis.com - C/ Tuset

More information

Lecture Note III: Least-Squares Method

Lecture 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 information

A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate

A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate Juan Simarro and Mariano Hortal AEMET Agencia Estatal de

More information

FINITE ELEMENT ANALYSIS USING THE TANGENT STIFFNESS MATRIX FOR TRANSIENT NON-LINEAR HEAT TRANSFER IN A BODY

FINITE ELEMENT ANALYSIS USING THE TANGENT STIFFNESS MATRIX FOR TRANSIENT NON-LINEAR HEAT TRANSFER IN A BODY Heat transfer coeff Temperature in Kelvin International Journal of Engineering Research and Applications (IJERA) ISSN: 2248-9622 FINITE ELEMENT ANALYSIS USING THE TANGENT STIFFNESS MATRIX FOR TRANSIENT

More information

Class 15 : Electromagnetic Waves

Class 15 : Electromagnetic Waves Class 15 : Electromagnetic Waves Wave equations Why do electromagnetic waves arise? What are their properties? How do they transport energy from place to place? Recap (1) In a region of space containing

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

Electromagnetics in LS-DYNA

Electromagnetics in LS-DYNA Electromagnetics in LS-DYNA Fully implicit. Double precision. 2D axisymmetric solver /3D solver. SMP and MPP versions available. Automatically coupled with LS-DYNA solid and thermal solvers. FEM for conducting

More information

From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D

From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D Luc Giraud N7-IRIT, Toulouse MUMPS Day October 24, 2006, ENS-INRIA, Lyon, France Outline 1 General Framework 2 The direct

More information

Some notes about PDEs. -Bill Green Nov. 2015

Some notes about PDEs. -Bill Green Nov. 2015 Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can

More information

Direct Numerical Simulation of Two-Phase Flows with phase change

Direct Numerical Simulation of Two-Phase Flows with phase change Direct Numerical Simulation of Two-Phase Flows with phase change Sébastien Tanguy, Lucia Rueda Villegas, Catherine Colin, Grégory Huber Institut de Mécanique des Fluides de Toulouse Maître de conférence

More information

Open boundary conditions in numerical simulations of unsteady incompressible flow

Open boundary conditions in numerical simulations of unsteady incompressible flow Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can

More information

Numerical modeling of magnetic induction and heating in injection molding tools

Numerical modeling of magnetic induction and heating in injection molding tools Downloaded from orbit.dtu.dk on: Apr 6, 08 Numerical modeling of magnetic induction and heating in injection molding tools Guerrier, Patrick; Hattel, Jesper Henri Published in: Proceedings of International

More information

Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems

Adaptive 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 information

Angular momentum preserving CFD on general grids

Angular 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 information

Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions

Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates

More information

Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems

Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems z Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems Johan Helsing Lund University Talk at Integral equation methods: fast algorithms and applications, Banff,

More information

Recovery-Based a Posteriori Error Estimators for Interface Problems: Mixed and Nonconforming Elements

Recovery-Based a Posteriori Error Estimators for Interface Problems: Mixed and Nonconforming Elements Recovery-Based a Posteriori Error Estimators for Interface Problems: Mixed and Nonconforming Elements Zhiqiang Cai Shun Zhang Department of Mathematics Purdue University Finite Element Circus, Fall 2008,

More information

Numerical Solution Techniques in Mechanical and Aerospace Engineering

Numerical Solution Techniques in Mechanical and Aerospace Engineering Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE

More information

ICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals

ICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals ICES REPORT 17-8 October 017 Direct Serendipity Finite Elements on Convex Quadrilaterals by Todd Arbogast and Zhen Tao The Institute for Computational Engineering and Sciences The University of Texas at

More information

Magnetostatic Analysis of Solenoid

Magnetostatic Analysis of Solenoid Magnetostatic Analysis of Solenoid 1. Introduction 2. Model View 3. Materials 4. Load & Restraint Information 5. Coils Information 6. Force and Torque Information 7. Study Properties 8. Results Table 9.

More information

Radiation hydrodynamics of tin targets for laser-plasma EUV sources

Radiation hydrodynamics of tin targets for laser-plasma EUV sources Radiation hydrodynamics of tin targets for laser-plasma EUV sources M. M. Basko, V. G. Novikov, A. S. Grushin Keldysh Institute of Applied Mathematics, Moscow, Russia RnD-ISAN, Troitsk, Moscow, Russia

More information

A Nested Dissection Parallel Direct Solver. for Simulations of 3D DC/AC Resistivity. Measurements. Maciej Paszyński (1,2)

A Nested Dissection Parallel Direct Solver. for Simulations of 3D DC/AC Resistivity. Measurements. Maciej Paszyński (1,2) A Nested Dissection Parallel Direct Solver for Simulations of 3D DC/AC Resistivity Measurements Maciej Paszyński (1,2) David Pardo (2), Carlos Torres-Verdín (2) (1) Department of Computer Science, AGH

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme

More information

Mixed Hybrid Finite Element Method: an introduction

Mixed Hybrid Finite Element Method: an introduction Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola

More information

NIMROD Boundary Conditions

NIMROD Boundary Conditions NIMROD Boundary Conditions Carl Sovinec University of Wisconsin-Madison and the Two-Fluid and Transport Group Plasma Science and Innovation Center Annual Meeting Seattle Washington, August 13-14, 2007

More information

A Review of Preconditioning Techniques for Steady Incompressible Flow

A 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 information

Recapitulation: Questions on Chaps. 1 and 2 #A

Recapitulation: 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 information

A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model

A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program University of Washington IEEE International Conference

More information

FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment

FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe

More information

Iterative Solvers in the Finite Element Solution of Transient Heat Conduction

Iterative Solvers in the Finite Element Solution of Transient Heat Conduction Iterative Solvers in the Finite Element Solution of Transient Heat Conduction Mile R. Vuji~i} PhD student Steve G.R. Brown Senior Lecturer Materials Research Centre School of Engineering University of

More information

Haus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, ISBN:

Haus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, ISBN: MIT OpenCourseWare http://ocw.mit.edu Haus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, 1989. ISBN: 9780132490207. Please use the following

More information

Development of discontinuous Galerkin method for linear strain gradient elasticity

Development of discontinuous Galerkin method for linear strain gradient elasticity Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace

More information

Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons

Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,

More information

Standard Finite Elements and Weighted Regularization

Standard Finite Elements and Weighted Regularization Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the

More information

Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics

Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Project work at the Department of Mathematics, TUHH Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Katja Baumbach April 4, 005 Supervisor:

More information

Index. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2

Index. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2 Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604

More information

20 Poynting theorem and monochromatic waves

20 Poynting theorem and monochromatic waves 0 Poynting theorem and monochromatic waves The magnitude of Poynting vector S = E H represents the amount of power transported often called energy flux byelectromagneticfieldse and H over a unit area transverse

More information

Code No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)

More information

An Introduction to the Discontinuous Galerkin Method

An 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 information

THE CONVECTION DIFFUSION EQUATION

THE CONVECTION DIFFUSION EQUATION 3 THE CONVECTION DIFFUSION EQUATION We next consider the convection diffusion equation ɛ 2 u + w u = f, (3.) where ɛ>. This equation arises in numerous models of flows and other physical phenomena. The

More information

Aeroacoustic and Aerodynamics of Swirling Flows*

Aeroacoustic and Aerodynamics of Swirling Flows* Aeroacoustic and Aerodynamics of Swirling Flows* Hafiz M. Atassi University of Notre Dame * supported by ONR grant and OAIAC OVERVIEW OF PRESENTATION Disturbances in Swirling Flows Normal Mode Analysis

More information