[alimea]

Originally Posted by opedrofunk

Hi All,

Happy to announce the release of blastFoam - a new solver for multi-component compressible flow with application to high explosive detonation, explosive safety and protective design. A user guide as well as multiple tutorial and validation cases are included.

https://github.com/synthetik-technologies/blastfoam

The solver features:

blastFoam currently supports the following features:

• An arbitrary number of phases/EOS's
• JLW equation of state with constant, linear, and "Miller" afterburn models
• Multiple example and tutorial cases
• Automatic mesh refinement (AMR)
• Single and multi-point detonation
• High-order (1st, 2nd, 3rd and 4th order in time; 2nd and 3rd order spatial)
• HLLC, AUSM+, Kurganov, Tadmor flux schemes
• Parallel
• Compatible with all of OpenFOAM's standard mesh generation, pre- and post-processing utilities; enhanced version of setFields to refine around material boundaries when initialized

blastFoam includes the following equations of state:

• Jones Wilkens Lee (JWL) (with afterburn)
• Ideal Gas
• Stiffened Gas
• Cochran-Chan
• Tait

Best regards,
Peter Vonk
Synthetik Applied Technologies

Dear Peter,

Thanks for sharing your valuable solver. I have a couple of questions:

1- As you know, the explosion is a turbulent phenomenon, but I didn't find any thing about it in the documents, and in tutorials files. Is your solver developed for turbulent flow regime? If yes, which turbulence model is used?

2- Can we have BKW equation of state, to compare the results of explosion with ideal gas, and JWL?

3- I want to simulate the blast wave created after an explosive. Before starting my simulations, I want to validate my solution with the experiments, and numerical investigations of Sklavounos, and Rigas which is published as:

Computer simulation of shock waves transmission in obstructed
terrains, Journal of Loss Prevention in the Process Industries 17 (2004) 407–417.

At first, I developed the sonicFoam to use the explosion energy release (Er) within the domain was given through a properly adapted step function
under the form of Eq. 4.1 of the paper (I should do that for your solver). But I couldn't get good results. After that I tried to run the case with blastFaom. However, I couldn't set the case!
I have attached main case, which is simulated using sonicFoam to this message.

http://uupload.ir/view/no2_02.1-perfectgas.rar/


Best regards,
Ali