1 / 33

Zoran Dragojlovic, Farrokh Najmabadi, Marcus Day

Prediction of Fluid Dynamics in The Inertial Confinement Fusion Chamber by Godunov Solver With Adaptive Grid Refinement. Zoran Dragojlovic, Farrokh Najmabadi, Marcus Day. Accomplishments. IFE Chamber dynamics code SPARTAN is developed:

chicom
Download Presentation

Zoran Dragojlovic, Farrokh Najmabadi, Marcus Day

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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Prediction of Fluid Dynamics in The Inertial Confinement Fusion Chamber by Godunov Solver With Adaptive Grid Refinement Zoran Dragojlovic, Farrokh Najmabadi, Marcus Day

  2. Accomplishments • IFE Chamber dynamics code SPARTAN is developed: Simulation of Physics by Algorithms based on RobustTurbulentApproximation of Navier-Stokes Equations. • SPARTAN current features: • 2-D Navier Stokes equations, viscosity and thermal conductivity included • arbitrary geometry • adaptive mesh refinement • SPARTAN tests: • Acoustic wave propagation. • Viscous channel flow. • Mach reflection. • Analysis of discretization errors to find code accuracy. • Initial conditions from BUCKY code are used for simulations. • Two Journal articles on SPARTAN are in preparation.

  3. Governing Equations of Fluid Flow 2-D Navier-Stokes equations in conservative form: Solution vector:

  4. Governing Equations of Fluid Flow Flux terms in x and y direction:

  5. Viscosity In general, m depends on (T, Z) T – temperature, Z – average ionization stage 2 extreme cases: neutral gas and fully ionized gas on a temperature range (10,000-60,000)K neutral gas: := (2.6 x 10-4 – 6.55 x 10-4) Pas fully ionized gas: := (4.9 x 10-11 – 4.4 x 10-9) Pas

  6. Thermal Conductivity Thermal conductivity depends on (T, Z), as well as viscosity. 2 extreme cases: neutral gas and fully ionized gas on a temperature range (10,000-60,000)K neutral gas: := (6.2 x 10-2 – 0.156) W/(mK) fully ionized gas: := (0.022 – 1.94) W/m-K

  7. Godunov Method Introduced in 1959, as a finite volume method with a special method of upwinding. Uses solution to a 1-D Riemann problem in order to estimate fluxes at the interface between cells. Formulation of Riemann problem: 1-D governing equation initial condition U(x,0) does not necessarily satisfy the conservation laws, breaks into fans, shocks and contact discontinuities.

  8. Godunov Method Solution procedure for Riemann problem uses Hugoniot jump relations and second law or thermodynamics to estimate pressure at the discontinuity, speed of propagation of discontinuities and the corresponding values of state variables. Numerical application: u,t wave diagram for Riemann problem uj+1 uj Uj-1 x j j-1 j+1

  9. wall beam channel chamber Adaptive Mesh Refinement • Motivation: efficient grid distribution results in reasonable CPU time. • Grid organized into levels from coarse to fine. • Solution tagging based on density and energy gradients. • Grid is refined at every time step. • Solution interpolated in space and time between the grid levels. • Referenced in: Almgren et al., 1993.

  10. max min Example of Adaptive Mesh Refinement geometry density contour plot

  11. Embedded Boundary Algorithm Arbitrary geometry imposed onto regular grid. This results in formation of cells irregular in shapes and sizes. Conservative update of irregular cells needs to be consistent and stable. fluid domain embedded boundary “forbidden” domain

  12. Embedded Boundary Algorithm Conservative update of irregular cells, consistent with split formulation of governing equations: Fnorth FB Fwest Feast Fsouth

  13. Embedded Boundary Algorithm non-conservative update preliminary update irregular cell’s neighborhood fluid domain cell for nonconservative update irregular cell

  14. Error Analysis Reflection of shock waves from the walls of cylindrical chamber was studied. Initial condition imposed by rotation of 1-D BUCKY solution about the center of the chamber. The wave was propagated for 0.001s, until it reflected from the wall and started converging back towards the center.

  15. Error Analysis Two cases considered:non-diffusive flow (m, k=0.0) and diffusive flow (m, k given by Sutherland law). Error analysis was done to determine the influence of the viscous and thermal diffusion on the accuracy of the solution. The errors were estimated by 4th order Richardson extrapolation.

  16. Error Analysis error no diffusion avg. = 3% 520x520 grid density no diffusion min=3.68e-5 kg/m3 max=1.19e-3 kg/m3 error diffusive flow avg. = 3.1% 520x520 grid density diffusive flow min=3.69e-5 kg/m3 max=2.09e-3 kg/m3

  17. Error Analysis error no diffusion avg. = 3.6% 520x520 grid X-momentum no diffusion min= -1.67kg/(m2s) max= 1.63kg/(m2s) error diffusive flow avg. = 3.2% 520x520 grid X-momentum diffusive flow min=-1.07kg/(m2s) max=0.96kg/(m2s)

  18. Error Analysis error no diffusion avg. = 1.57% 520x520 grid pressure no diffusion min=50.01 Pa max=1.262e3 Pa error diffusive flow avg. = 1.52% 520x520 grid pressure diffusive flow min=50.42 Pa max=1.24e3 Pa

  19. Error Analysis error no diffusion avg. = 1.86% 520x520 grid energy no diffusion min=130.59J max=4.492e3J error diffusive flow avg. = 1.7% 520x520 grid energy diffusive flow min=131.19J max=3.66e3 J

  20. Validation Problems cylindrical shock reflection from rectangular geometry

  21. initial disturbance Mach reflection 10.1 11.2 12.1 12.9 r[m]

  22. Viscous Flow Through the Channel

  23. IFE Chamber Dynamics Simulations Objectives • Determine the influence of the following factors on the chamber state at 100 ms: • viscosity • blast position in the chamber • heat conduction from gas to the wall. • Chamber density, pressure, temperature, and velocity distribution prior to insertion of next target are calculated.

  24. Numerical Simulations IFE Chamber Simulation • 2-D cylindrical chamber with a laser beam channel on the side. • 160 MJ NRL target • Boundary conditions: • Zero particle flux, Reflective velocity • Zero energy flux or determined by heat conduction. • Physical time: 500 ms (BUCKY initial conditions) to 100 ms.

  25. Numerical Simulations Initial Conditions • 1-D BUCKY solution for density, velocity and temperature at 500 ms imposed by rotation and interpolation. • Target blast has arbitrary location near the center of the chamber. • Solution was advanced by SPARTAN code until 100 ms were reached.

  26. max min Effect of Viscosity on Chamber State at 100 ms inviscid flow at 100 ms viscous flow at 100 ms pressure, pmean = 569.69 Pa pressure, pmean = 564.87 Pa temperature, Tmean = 5.08 104 K (rCvT)mean = 1.412 103 J/m3 temperature, Tmean = 4.7 104 K (rCvT)mean = 1.424 103 J/m3 Viscosity makes a difference due to it’s strong dependence on temperature.

  27. max min Effect of Blast Position on Chamber State at 100ms centered blast at 100 ms eccentric blast at 100 ms pressure, pmean = 564.87 Pa pressure, pmean = 564.43 Pa temperature, Tmean = 4.7 104 K temperature, Tmean = 4.74 104 K Large disturbance due to eccentricity of blast and small numerical disturbances have the same effect after 100 ms.

  28. Effect of Blast Position on Chamber State at 100 ms pressure at the wall pressure at the mirror • Mirror is normal to the beam tube. • Pressure is conservative by an order of magnitude. • Pressure on the mirror is so small that the mechanical response is negligible.

  29. max min Effect of Wall Heat Conduction on Chamber State at 100 ms insulated wall wall conduction pressure, pmean = 564.431 Pa pressure, pmean = 402.073 Pa temperature, tmean = 4.736 104 K temperature, tmean = 2.537 104 K Ionized gas or plasma makes a difference by the means of heat conduction.

  30. max min Chamber Gas Dynamics pressure temperature density

  31. Prediction of chamber condition at long time scale is the goal of chamber simulation research. • Chamber dynamics simulation program is on schedule. Program is based on: • Staged development of Spartan simulation code. • Periodic release of the code and extensive simulations while development of next-stage code is in progress. • Documentation and Release of Spartan (v1.0) • Two papers are under preparation • Exercise Spartan (v1.x) Code • Use hybrid models for viscosity and thermal conduction. • Parametric survey of chamber conditions for different initial conditions (gas constituent, pressure, temperature, etc.) • Need a series of Bucky runs as initial conditions for these cases. • We should run Bucky using Spartan results to model the following shot and see real “equilibrium” condition. • Investigate scaling effects to define simulation experiments.

  32. Several upgrades are planned for Spartan (v2.0) Numeric: • Implementation of multi-species capability: • Neutral gases, ions, and electrons to account for different thermal conductivity, viscosity, and radiative losses. Physics: • Evaluation of long-term transport of various species in the chamber (e.g., material deposition on the wall, beam tubes, mirrors) • Atomics and particulate release from the wall; • Particulates and aerosol formation and transport in the chamber. • Improved modeling of temperature/pressure evolution in the chamber: • Radiation heat transport; • Equation of state; • Turbulence models.

More Related