1 / 32

Donald F. Hawken Ph. D.

Donald F. Hawken Ph. D. Flow Simulation Work. Detonation with finite-rate chemistry. 0.05 meter 1D domain with 1500 cells 298.15 °K and 1 atmosphere ambient 0.001 meter driver at 3000° K and 30 atmospheres Spatially first-order explicit HLLC solver for fluxes

cyrah
Download Presentation

Donald F. Hawken Ph. D.

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. Donald F. Hawken Ph. D. Flow Simulation Work

  2. Detonation with finite-rate chemistry • 0.05 meter 1D domain with 1500 cells • 298.15 °K and 1 atmosphere ambient • 0.001 meter driver at 3000° K and 30 atmospheres • Spatially first-order explicit HLLC solver for fluxes • Sources use implicit 6th order Runge-Kutta algorithm • Compared to results of CFD++ code

  3. H2-O2 detonation: 6 species, 7 reactions

  4. Nozzle flow of H2-O2 products

  5. H2-Air detonation: 12 species, 27 reactions

  6. FCT solver with solution adaptive mesh • 10 atmosphere 300 °K cylinder of air expands into one atmosphere air • Domain: 10 meters x 10 meters with 40 x 40 cells • Two levels of cell refinement in response to pressure gradients • FCT solver modified to handle cell refinement

  7. Pressure contours at 4 milliseconds (Pascal)

  8. Mesh at 4 milliseconds

  9. Level Set Tracking • Location of shock-, detonation-, or contact front is tracked using scalar value, G, that is zero at front. • Value of G in a cell is the signed distance from the front • Shock or detonation front moves at speed determined by Hugoniot relation • Contact front moves at flow speed • Two states are maintained in cells containing front or near to front. • Front states change in response to edge fluxes and front fluxes • Spatially second-order explicit HLLC solver

  10. Pressure and G profiles for tracked shock (100 cells)

  11. Pressure profiles for captured shock (100 cells)

  12. Pressure: tracked shock (100 cells) and captured shock (400 cells)

  13. Pressure contours for tracked circular shock (40x40 cells)

  14. Pressure contours for captured circular shock (40x40 cells)

  15. Contours of front distance, G, for tracked circular shock

  16. High-pressure gas compresses water (with and without 3 levels of refinement) (tracked contact)

  17. High-pressure gas bubble expands in water (with and without 3 levels of refinement) (tracked contact)

  18. Simulation of explosive blast (near field to far field) • Domain: 10.9 meters x 5.1 meters with 100 x 140 cells • Detonation front is tracked with two levels of refinement until solid explosive cylinder is consumed • Spatially second-order explicit HLLC solver used for explosive products and air • Trajectory of shock at height of burst matches raw experimental data at mid field and far field • Simulation produces better estimate of front pressures near burst compared to extrapolated experimental fit to raw experimental data

  19. Close up showing simulation initial conditions (click for movie)

  20. Blast pressure movie (direct shock overtaken by ground shock)

  21. Blast fireball movie

  22. Location of shock front at height of burst

  23. Improved near-field estimate of shock front pressure at height of burst

  24. Detonation of explosive charge mixed with inert metal particles (near field to mid field) • Expanding sphere of combustion products and particles • JWL gas equation of state - combustion products • van der Waals gas equation of state - air • Dense particle equation of state • HLLC solver - combustion products and air • FCT Solver - particles • Momentum and energy exchange between gas and particle phases

  25. Shock and particle front histories compared to experiment

  26. Shock front pressure histories compared to experiment

  27. Detonation of liquid explosive mixed with reactive metal particles (near field) • HOM condensed equation of state - liquid explosive • HOM gas equation of state - explosive combustion products • van der Waals gas equation of state - air • Particle equation of state • HLLC solver - liquid explosive, combustion products, and air • FCT solver - particles • Momentum, mass, and energy exchange between gas and particle phases

  28. Chemical Reactions • Arrhenius combustion model for nitromethane • Nitromethane  H20 + (CO2-CO) • Particle-diameter controlled oxidation rate • Metal+H20 + (CO2-CO)+O2  metal oxides

  29. Pure nitromethane detonation in 25 mm glass tube

  30. Detonation of pure nitromethane cylinder (100x200 cells, 2 refinement levels) Gas density movie

  31. Nitromethane-particle detonation in 15 mm glass tube

  32. Detonation of nitromethane-particle cylinder (100x200 cells, 2 levels of refinement) Particle-phase density movie

More Related