search

CEED / Laghos

High-order Lagrangian Hydrodynamics Miniapp
http://ceed.exascaleproject.org/miniapps
BSD 2-Clause "Simplified" License
187 stars 60 forks source link

Add Rayleigh-Taylor Instability Problem #118

Closed siuwuncheung closed 3 years ago

siuwuncheung commented 3 years ago

The Rayleigh-Taylor Instability Problem can be run as follows: ./laghos -p 7 -m data/rt2D.mesh -tf 4 -rs 1 -ok 4 -ot 3 -vis -pa The figures are taken at t = 4.0, and the formation of the ”mushroom cap” can be seen.

Velocity at t = 4.0 velocity

Energy at t = 4.0 energy

Density at t = 4.0 density

vladotomov commented 3 years ago

Could you add the sample run to README?

siuwuncheung commented 3 years ago

Could you add the sample run to README?

@vladotomov Thanks for reviewing this PR. I can add the sample run to README. However, I am not able to reproduce the results from the paper with the command line above so far. I am not sure if I misunderstood the problem setting. The amplitude of the "sliding" surface is much smaller than that in the paper, and there is a distortion of the mesh around the bottom-left corner and the top-right corner which does not appear in the paper. I will soon put a figure of the density at the final time in this PR. I am thinking that if it is because of the boundary condition. In this problem, it seems that we should enforce symmetry planes and impose v.n = 0 only for the vertical walls, unlike the other tests we do it on the whole boundary. If that is the case, it would take a while for me to figure out how to implement the boundary condition.

camierjs commented 3 years ago

There is still an issue with the CUDA run: let me look at it.

camierjs commented 3 years ago

Energy at t = 3.88 Screen Shot 2020-11-30 at 4 31 39 PM

vladotomov commented 3 years ago

@camierjs @tzanio are we merging?