Simulation and optimization of electromagnetohydrodynamic flows
This presentation is the property of its rightful owner.
Sponsored Links
1 / 128

Simulation and Optimization of Electromagnetohydrodynamic Flows PowerPoint PPT Presentation


  • 83 Views
  • Uploaded on
  • Presentation posted in: General

Simulation and Optimization of Electromagnetohydrodynamic Flows. Brian Dennis Ph.D. Candidate Aerospace Engineering Dept. Penn State University George Dulikravich Professor Mechanical and Aerospace Dept. University of Texas at Arlington. Overview. Introduction

Download Presentation

Simulation and Optimization of Electromagnetohydrodynamic Flows

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Simulation and optimization of electromagnetohydrodynamic flows

Simulation and Optimization of Electromagnetohydrodynamic Flows

Brian Dennis

Ph.D. Candidate

Aerospace Engineering Dept.

Penn State University

George Dulikravich

Professor

Mechanical and Aerospace Dept.

University of Texas at Arlington


Overview

Overview

  • Introduction

  • Governing equations for EMHD

  • General LSFEM formulation

  • h and p-type finite element solver

  • LSFEM for EMHD

  • Code validation cases

  • Applications

  • Conclusion and recommendations


Introduction

Introduction

  • EMHD is study of incompressible flow under the influence of electric and magnetic fields

  • Various simplified analytical models exist and have been used for numerical simulation(EHD,MHD)


Introduction cont

Introduction cont.

  • A fully consistent non-linear model for general EMHD has been recently developed

  • No numerical simulations of full EMHD has been report in open literature to date

  • A few numerical simulations of simplified or inconsistent EMHD models have been report in open literature

  • A computer code for numerical simulation for 2-D planar MHD/EMHD flows has been developed using LSFEM


Introduction cont1

Introduction cont.

  • Numerical simulation is necessary for performing optimization involving EMHD flows


Applications of emhd

Applications of EMHD

  • Manufacturing(solidification,crystal growth)

  • Flow control

  • Drag reduction/propulsion

  • Pumps with no moving parts(artificial heart, liquid metal pumps)

  • Compact heat exchangers

  • Shock absorbers, active damping


Fully consistent model of emhd

Fully Consistent Model of EMHD


Conservation of linear momentum

Conservation of Linear Momentum


Conservation of energy

Conservation of Energy

Conservation of Mass


Maxwell s equations

Maxwell’s Equations


Sub models

Sub-models

  • The LSFEM was applied to two sub-models of the fully nonlinear Electro-magneto-hydrodynamic system

  • Magneto-hydrodynamics

  • Electro-magneto-hydrodynamics with reduced number of source terms


Least squares finite element method lsfem

Least-Squares Finite Element Method (LSFEM)


Advantages of the lsfem

Advantages of the LSFEM

  • Can use equal order basis functions for pressure and velocity

  • can use the first order form of PDE’s

  • can handle any type of equation and mixed types of equations

  • Can discretize convection terms without upwinding or explicit artificial dissipation

  • Stable and robust method

  • Resulting system of equations is symmetric and positive definite


Simulation and optimization of electromagnetohydrodynamic flows

  • Simple iterative techniques such as PCG and multigrid can be used to solve the system of equations

  • Inclusion of divergence constraint for magnetic flux is straight forward


Least squares finite element method lsfem1

Least-squares finite element method (LSFEM)


H and p type finite element solver

h and p-Type Finite Element Solver


H and p type methods

h and p-Type Methods

  • h-type finite element methods use large numbers of low order accurate elements

  • p-type finite element methods use small numbers of high order accurate elements

  • h/p-type finite element methods use a combination of both types together with an adaptive strategy

  • p-type methods convergence to the exact solution is more rapid than in with h-type methods, if the underlying solution is smooth


P type method

p-Type Method

  • element approximation function are compose of p-type expansions. These expansions are composed of a summation of P polynomials

  • p-type expansions can be categorized as nodal expansions or modal expansions

  • the unknown coefficients in nodal expansions are the values of the function at the nodes of the element and therefor have some physical meaning. Nodal expansions based on Lagrange polynomials are typically used for low order finite elements

  • the unknown coefficients in modal expansions are not associated with any nodes


P type method1

p-Type Method

  • Modal expansions are hierarchical, that is, that all expansion sets less than P+1 are contained within the expansion set P+1

  • Nodal expansions are not hierarchical

  • Modal expansion typically produce finite element matrices that have a lower condition number than those produced with nodal expansions

  • Expansions are usually developed in one dimension. Approximation functions for multidimensional elements are constructed through tensor products of the one dimensional expansions.

  • Typically, the modal expansions are combined with the first order Lagrange polynomials so that continuity of the solution is satisfied between neighboring elements in a mesh


Lsfem code

LSFEM Code

  • A serial code was developed in C/C++ to solve general systems with LSFEM

  • Steady state problems only

  • Mixed triangular and quadrilateral meshes for h-type elements.

  • Quadrilateral elements for p-type elements. Element approximation functions are constructed from a modal basis derived from Jacobi polynomials

  • Support for multiple material domains such as in conjugate heat transfer problems for h-version

  • Nonlinear equations are linearized with Newton or Picard method


Lsfem code1

LSFEM Code

  • Integrals are evaluated numerically using Legendre-Gauss quadrature. A P+2 rule was used for p-type elements.

  • Static condensation is used to remove the interior degrees of freedom from p-type elements.


Sparse linear solvers

Sparse Linear Solvers

  • Direct Solver:

  • Sparse LU factorization from PETSc library

  • Iterative Solvers:

  • Jacobi preconditioned conjugate gradient method

  • Multi-p multilevel method


Lsfem for emhd

LSFEM for EMHD


Nondimensional form for emhd equations

Nondimensional Form for EMHD Equations


Nondimensional form for emhd equations1

Nondimensional Form for EMHD Equations


First order form for emhd equations

First Order Form for EMHD Equations


Iterative solution algorithm

Iterative Solution Algorithm


General first order system for 2 d

General First Order System for 2-D


General first order system for 2 d1

General First Order System for 2-D


General first order system for 2 d2

General First Order System for 2-D


Code validation

Code Validation


Verification of accuracy

Verification of Accuracy

  • Few analytic solutions for EMHD exist

  • Some analytic solutions for MHD exist

  • NSE portion of code was validated using analytic solutions for NSE and with experiment data from driven cavity flows and backward facing step

  • Heat transfer/Electric/Magnetic field portions were verified with analytic solutions


Test cases for nse

Test Cases for NSE

  • Test against analytical solutions for NSE

  • Test against driven cavity numerical benchmark solutions

  • Test against experiment data for flow over backward facing step


Mesh for analytic test cases

Mesh for analytic test cases


Comparison between computed and analytic solution for couette poiseuille flow

Comparison between computed and analytic solution for Couette-Poiseuille flow


Mesh for driven cavity benchmark

Mesh for driven cavity benchmark


Comparison with benchmark results for driven cavity at re 1000

Comparison with benchmark results for driven cavity at Re = 1000


Computed streamlines re 1000

Computed streamlinesRe = 1000


Computed pressure re 1000

Computed pressureRe = 1000


Computed vorticity re 1000

Computed vorticityRe = 1000


Comparison with benchmark results for driven cavity at re 10000

Comparison with benchmark results for driven cavity at Re = 10000


Computed streamlines re 10000

Computed streamlinesRe = 10000


Computed pressure re 10000

Computed pressureRe = 10000


Computed vorticity re 10000

Computed vorticityRe = 10000


Convergence history for driven cavity benchmark

Convergence history for driven cavity benchmark


Backward facing step problem

Backward Facing Step Problem


Mesh for backward facing step problem

Mesh for Backward Facing Step Problem


Computed streamlines

Computed Streamlines

Re = 100

Re = 400

Re = 500


Comparison of computed and measured reattachment lengths

Comparison of Computed and Measured Reattachment Lengths


Comparison of present results with results from commercial cfd codes

Comparison of Present Results with Results from Commercial CFD Codes


Mesh for buoyancy driven cavity

Mesh for Buoyancy-Driven Cavity


Computed static temperature distribution

Computed Static Temperature Distribution


Computed heat flux on vertical walls

Computed Heat Flux on Vertical Walls


Computed streamlines1

Computed Streamlines


Computed streamline in upper left corner

Computed Streamline in Upper Left Corner


Comparison of present computations with other published results

Comparison of Present Computations with Other Published Results


Test cases for mhd and emhd

Test Cases for MHD and EMHD

  • Test against analytical solutions for MHD and EMHD


Lsfem for mhd

LSFEM for MHD


First order system for mhd

First order system for MHD


First order system for 2 d mhd

First order system for 2-D MHD


First order system for 2 d mhd cont

First order system for 2-D MHD cont.


Mesh for hartmann flow test case

Mesh for Hartmann Flow Test Case


Hartmann flow

Hartmann Flow

MHD LSFEM code was compared with the analytic solution to Poisuille-Hartmann flow


Comparison of lsfem and analytic solution for hartmann flow

Comparison of LSFEM and analytic solution for Hartmann flow


Comparison of lsfem and analytic solution for hartmann flow1

Comparison of LSFEM and analytic solution for Hartmann flow


Comparison of lsfem and analytic solution for hartmann flow2

Comparison of LSFEM and analytic solution for Hartmann flow

Computed and analytic velocity profile


Comparison of lsfem and analytic solution for hartmann flow3

Comparison of LSFEM and analytic solution for Hartmann flow

Computed and analytic induced magnetic field


Comparison of h and p type methods for hartmann flow

Comparison of h and p-type methods for Hartmann flow


Simulation and optimization of electromagnetohydrodynamic flows

Comparison of LSFEM and analytic solution for Hartmann flow with applied electric and magnetic field


Simulation and optimization of magneto hydrodynamic flows with lsfem

Simulation and Optimization of Magneto-hydrodynamic Flowswith LSFEM


Simulation and optimization of electromagnetohydrodynamic flows

Optimization of Magneto-Hydrodynamic Control of Diffuser Flows Using Micro-Genetic Algorithms and Least-Squares Finite Elements


Simulation and optimization of electromagnetohydrodynamic flows

Goal

  • Given a fixed diffuser shape, use micro-GA and LSFEM MHD analysis to design a magnetic field distribution on the diffuser wall that will increase static pressure rise


Flow solver

Flow Solver

  • LSFEM solver for 2-D steady incompressible Navier-Stokes together with Maxwell’s equations for steady magnetic field

  • Uses hybrid quadrilateral/triangular grid

  • One analysis takes around 22 min. on a single Pentium II CPU


Bc s and parameterization

BC’s and Parameterization

  • Parabolic velocity specified at inlet

  • Static pressure specified at outlet

  • no-slip conditions on wall

  • magnetic field component along wall are specified. They were parameterized with b-spline

  • perfectly conducting wall bc used on all other solid surfaces


Simulation and optimization of electromagnetohydrodynamic flows

Hybrid mesh for diffuser analysis

Given diffuser shape


Convergence history for a typical design

Convergence History for a Typical Design


Parallel genetic algorithm

Parallel Genetic Algorithm

  • GA is a naturally coarse grained parallel algorithm

  • One node maintains the population(master) and distributes jobs to the slave nodes

  • Only simple synchronous message passing is needed to implement on distibuted memory

  • Population size need not match the number of slave nodes

  • Asynchronous models are also being developed for use when function analysis computation times vary dramatically.


Parallel computer

Parallel Computer

  • based on commodity hardware components and public domain software

  • 16 dual Pentium II 400 MHz based PC’s

  • 100 Megabits/second switched ethernet

  • total of 32 processors and 8.2 GB of main memory

  • Compressible NSE solver achieved 1.5 Gflop/sec with a LU SSOR solver on a 100x100x100 structured grid


Parallel computer1

Parallel Computer

  • based on commodity hardware components and public domain software

  • 16 dual Pentium II 400 MHz based PC’s

  • 100 Megabits/second switched ethernet


Parallel computer cont

Parallel Computer cont.

  • total of 32 processors and 8.2 GB of main memory

  • Compressible NSE solver achieved 1.5 Gflop/sec with a LU SSOR solver on a 100x100x100 structured grid

  • GA optimization of a MHD diffuser completed in 30 hours. Same problem would take 14 days on a single CPU


Genetic algorithm

Genetic Algorithm

  • Population size of 15

  • 100 generations

  • 9-bit strings for each design variable

  • elitism

  • tournament selection

  • uniform crossover

  • parallel micro-GA


Ga convergence history

GA Convergence History

Run 1

Run 2


Results

Results

Two optimizations were run simultaneously with 16 processors each.

Generation 100 was reached by both in about 30 hours

run 1 achieved a pressure increase of .207 Pa

run 2 achieved a pressure increase of .228 Pa

diffuser without an applied magnetic field achieved a pressure increase of .05 Pa.


Simulation and optimization of electromagnetohydrodynamic flows

With no applied magnetic field

With optimized applied magnetic field


Simulation of magneto hydrodynamic flows with conjugate heat transfer with lsfem

Simulation of Magneto-hydrodynamic Flows with Conjugate Heat Transferwith LSFEM


Objective

Objective

  • Simulation of flow through channel with an applied magnetic field

  • Simulation of heat transfer from the flow to a solid cold wall

  • Observe the effect of applied magnetic field on flow patterns and heat transfer characteristics


Boundary conditions

Boundary Conditions

  • Inlet temperature of 2000 K

  • Specified outlet pressure

  • Specified parabolic velocity profile at inlet

  • No-slip on walls

  • Symmetry boundary condition on top

  • Temperature of 300 K on bottom wall

  • Perfectly conducting walls except in the region 7 < x < 8 where sinusoidal magnetic field components were specified. Magnitude was varied from 0 to 5 Tesla.


Results1

Results

  • Presence of magnetic field induces a large separation in the flow field close to the wall

  • Size and complexity are proportional to the strength of the magnetic field

  • A drop in fluid/solid interface temperature was observed in the region where the magnetic field was applied


Effect of flow on magnetic field

Effect of flow on magnetic field


Solidification with mhd

Figure 14: Velocity magnitude for no applied magnetic field.

Figure 15: Streamlines for no applied magnetic field.

Solidification with MHD


Optimized mhd solidification

Figure 16: Velocity magnitude for optimized case 1.

Figure 17: Streamlines for optimized case 1.

Optimized MHD Solidification


Optimizer mhd solidification

Figure 21: Velocity magnitude for optimized case 2.

Figure 22: Streamlines for optimized case 2.

Optimizer MHD Solidification


Solidification with optimized magnetic boundary conditions

Solidification with Optimized Magnetic Boundary Conditions

Figure 10. Isotherms and streamlines without and with an optimized magnetic field (test case 2 with six sensors)


Simulation of electro magneto hydrodynamic channel flows

Simulation of Electro-magneto-hydrodynamic Channel Flows


Objective1

Objective

  • Simulation of a steady state EMHD blood pump

  • Both electric and magnetic fields are required to produce the driving force


Governing equations

Governing Equations


Boundary conditions and geometry

Boundary conditions and geometry

  • Rectangular domain with height of 4 cm and length of 40 cm

  • Triangular mesh:

  • 7021 nodes

  • 3422 elements

  • parabolic triangles

  • Specified parabolic inlet velocity profile and temperature of 310.15 K

  • No slip on walls

  • Wall temperature was 298.15 K

  • Specified exit pressure of 1 Pa

  • Positive electrode on bottom wall, negative electrode on top with 50 volts applied across them

  • Uniform magnetic field of .05 Tesla specified in Z direction


Physical parameters for emhd blood pump

Physical parameters for EMHD blood pump

  • Density(kg m-3) = 1055.0

  • Inlet height(cm) = 4

  • Length(cm) = 40

  • Inlet temp.(K) = 310

  • Wall temp(K) = 298

  • heat conductivity(W kg-1 K-1) = .51

  • specific heat(J kg-1 K-1)= 4178

  • inlet velocity(m s-1) = .05

  • dynamic viscosity(kg m-1 s-1) = .004

  • electric conductivity (S m-1) = 1.4

  • outlet pressure (Pa) = 1


Computed electric potential

Computed Electric Potential


Computed velocity vectors

Computed Velocity Vectors


Computed static pressure

Computed Static Pressure


Computed static temperature

Computed Static Temperature


Emhd channel electrode configuration

EMHD Channel Electrode Configuration


Mesh for emhd channel problem

Mesh for EMHD channel problem


Computed electric potential distribution

Computed Electric Potential Distribution


Computed electric potential on lower wall

Computed Electric Potential on Lower Wall


Compute static pressure

Compute Static Pressure


Compute static pressure1

Compute Static Pressure


Computed axial velocity distribution

Computed Axial Velocity Distribution


Emhd channel electrode configuration1

EMHD Channel Electrode Configuration


Computed electric potential1

Computed Electric Potential


Computed electric field lines

Computed Electric Field Lines


Computed static pressure1

Computed Static Pressure


Computed streamlines2

Computed Streamlines


Conclusions

Conclusions

  • A code for the simulation of MHD/EMHD was developed based on the LSFEM

  • Code was tested against analytic solutions and experimental data for separate disciplines

  • Code was applied to several MHD problems including a MHD diffuser optimization problem

  • Code was used to simulate a EMHD pump and several channel flows with finite length electrodes


Conclusions1

Conclusions

  • LSFEM works well for NSE, MHD, and EMHD when the problem are not very nonlinear.

  • The ‘squaring’ effect inherent to the method appears to enhance any nonlinearity already present in the system. Thus, as Reynolds and Hartmann numbers increase, the required number of iterations increases rapidly


Conclusions2

Conclusions

  • For LSFEM, first order forms for the PDE’s are required in order to satisfy solution continuity requirements. This leads to a relatively large number of unknowns, especially in 3D, compared to a Galerkin type method

  • The linear algebraic systems can become ill-conditioned. Better preconditioners are needed.


Future work

Future Work

  • Better preconditioning of linear systems

  • Closer look at scaling of the equations

  • Extension to 3-D

  • Better boundary conditions

  • Unsteady


  • Login