search

CEED / Laghos

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

A new ROM implementation of Rayleigh-Taylor instability problem #142

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.

The following results are generated using the main directory which mirrors the most recent master.

Density at t = 4.0

master_density

Velocity at t = 4.0

master_velocity

Energy at t = 4.0

master_energy

siuwuncheung commented 3 years ago

All regression tests passed up to this point.

Density at t = 4.0 density

Velocity at t = 4.0 velocity

Energy at t = 4.0 energy

siuwuncheung commented 3 years ago

After a change in boundary condition, regression tests are broken as expected, but with negligible error.

Density at t = 4.0 density

Velocity at t = 4.0 velocity

Energy at t = 4.0 energy

dylan-copeland commented 3 years ago

The test ./laghos -p 7 -m data/rt2D.mesh -tf 4 -rs 1 -ok 4 -ot 3 -vis -pa works in parallel (e.g. mpirun -n 4) for the main version, but for the rom/ version the timestep size goes to zero. @siuwuncheung have you looked into the parallel case?

siuwuncheung commented 3 years ago

The test ./laghos -p 7 -m data/rt2D.mesh -tf 4 -rs 1 -ok 4 -ot 3 -vis -pa works in parallel (e.g. mpirun -n 4) for the main version, but for the rom/ version the timestep size goes to zero. @siuwuncheung have you looked into the parallel case?

I have been testing the serial runs only. I will look into the parallel runs.

siuwuncheung commented 3 years ago

The test ./laghos -p 7 -m data/rt2D.mesh -tf 4 -rs 1 -ok 4 -ot 3 -vis -pa works in parallel (e.g. mpirun -n 4) for the main version, but for the rom/ version the timestep size goes to zero. @siuwuncheung have you looked into the parallel case?

I have been testing the serial runs only. I will look into the parallel runs.

The parallel runs have been fixed.