Skip to content

Arbitrary-Order SGH 3D solver (AO-SGH) #333

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
steven-t-walton wants to merge 1 commit into
base: main
Choose a base branch
from
Draft

Arbitrary-Order SGH 3D solver (AO-SGH) #333

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2,476 changes: 2,476 additions & 0 deletions apps/multiphysics/example_inputs/box_4x4x4_p3.vtu
View file
Open in desktop

Large diffs are not rendered by default.

2,779 changes: 2,779 additions & 0 deletions apps/multiphysics/example_inputs/box_8x8x1_p3.vtu
View file
Open in desktop

Large diffs are not rendered by default.

109 changes: 109 additions & 0 deletions apps/multiphysics/example_inputs/input-ao-sgh-noh.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,109 @@
# AO_SGH3D Noh implosion (octant). rho = 1, v = -r_hat, sie ~ 0, gamma = 5/3.
# Exact at t = 0.6: shock at r = 0.2, post-shock rho = 64, upstream
# rho = (1 + t/r)^2. Slip only on the origin planes; outer faces free.
# dt_max capped because the CFL is sound-speed based and c ~ 0 at t = 0.
# Regenerate the mesh from build/ via:
# ./app/ao_sgh_box_vtu_gen --nx 12 --ny 12 --nz 12 --p 2 \
# --origin 0,0,0 --length 1,1,1 \
# --out box_12x12x12_p2_unit.vtu

dynamic_options:
time_final: 0.6
dt_min: 1.0e-9
dt_max: 5.0e-4
dt_start: 1.0e-5
dt_cfl: 0.25
cycle_stop: 100000

mesh_options:
source: file
file_path: box_12x12x12_p2_unit.vtu
num_dims: 3
polynomial_order: 2

output_options:
timer_output_level: thorough
output_file_format: viz
graphics_time_step: 0.1
node_field_outputs:
- vel

solver_options:
- solver:
method: ao_dynx_FE
id: 0
time_integration: rk3hc

boundary_conditions:
- boundary_condition:
solver_id: 0
surface:
type: x_plane
plane_position: 0.0
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 1
- 0
- 0
- boundary_condition:
solver_id: 0
surface:
type: y_plane
plane_position: 0.0
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 0
- 1
- 0
- boundary_condition:
solver_id: 0
surface:
type: z_plane
plane_position: 0.0
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 0
- 0
- 1

materials:
- material:
id: 0
eos_model_type: decoupled
eos_model: gamma_law_gas
dissipation_model: MARS
dissipation_global_vars:
- 1.0
- 1.0
- 1.333
- 1.333
- 0
- 1.0
eos_global_vars:
- 1.666666666666667
- 1.0e-14
- 1.0

regions:
- region:
id: 0
volume:
type: global
solver_id: 0

initial_conditions:
- initial_condition:
region_id: 0
material_id: 0
density:
type: uniform
value: 1.0
specific_internal_energy:
type: uniform
value: 1.0e-10
velocity:
type: spherical
speed: -1.0
166 changes: 166 additions & 0 deletions apps/multiphysics/example_inputs/input-ao-sgh-sedov.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,166 @@
# AO_SGH3D Sedov blast (octant). Kamm-calibrated energy for gamma = 5/3:
# E_total = 0.493390 -> E_octant = 0.0616738 deposited in the corner element
# (V = 0.1^3 = 1e-3) -> sie = 61.67375. Shock reaches r = 1 at t = 1.
# Slip on all six faces. Regenerate the mesh from build/ via:
# ./app/ao_sgh_box_vtu_gen --nx 12 --ny 12 --nz 12 --p 2 \
# --origin 0,0,0 --length 1.2,1.2,1.2 \
# --out box_12x12x12_p2.vtu

dynamic_options:
time_final: 1.0
dt_min: 1.0e-8
dt_max: 1.0e-2
dt_start: 1.0e-5
dt_cfl: 0.25
cycle_stop: 100000

mesh_options:
source: file
file_path: box_12x12x12_p2.vtu
num_dims: 3
polynomial_order: 2

output_options:
timer_output_level: thorough
output_file_format: viz
graphics_time_step: 0.1
node_field_outputs:
- vel

solver_options:
- solver:
method: ao_dynx_FE
id: 0
time_integration: rk3hc

boundary_conditions:
- boundary_condition:
solver_id: 0
surface:
type: x_plane
plane_position: 0.0
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 1
- 0
- 0
- boundary_condition:
solver_id: 0
surface:
type: x_plane
plane_position: 1.2
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 1
- 0
- 0
- boundary_condition:
solver_id: 0
surface:
type: y_plane
plane_position: 0.0
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 0
- 1
- 0
- boundary_condition:
solver_id: 0
surface:
type: y_plane
plane_position: 1.2
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 0
- 1
- 0
- boundary_condition:
solver_id: 0
surface:
type: z_plane
plane_position: 0.0
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 0
- 0
- 1
- boundary_condition:
solver_id: 0
surface:
type: z_plane
plane_position: 1.2
tolerance: 1.0e-5
velocity_model: reflected
velocity_bc_global_vars:
- 0
- 0
- 1

materials:
- material:
id: 0
eos_model_type: decoupled
eos_model: gamma_law_gas
dissipation_model: MARS
dissipation_global_vars:
- 1.0
- 1.0
- 1.333
- 1.333
- 0
- 1.0
eos_global_vars:
- 1.666666666666667
- 1.0e-14
- 1.0

regions:
- region:
id: 0
volume:
type: global
solver_id: 0

- region:
id: 1
volume:
type: sphere
radius1: 0.0
radius2: 0.09
solver_id: 0

initial_conditions:
- initial_condition:
region_id: 0
material_id: 0
density:
type: uniform
value: 1.0
specific_internal_energy:
type: uniform
value: 1.0e-10
velocity:
type: cartesian
u: 0.0
v: 0.0
w: 0.0

- initial_condition:
region_id: 1
material_id: 0
density:
type: uniform
value: 1.0
specific_internal_energy:
type: uniform
value: 61.67375002
velocity:
type: cartesian
u: 0.0
v: 0.0
w: 0.0
Loading