Skip to content

Commit

Permalink
Robustify Rayleigh Bernard two phase test
Browse files Browse the repository at this point in the history
Refs failures from the combination of idaholab#26462 and idaholab#26227

To avoid negative viscosity interpolations I shrank the growth
factor from 2 to 1.5. Once I did this I got a CSVDiff which signaled
to me that the original result was solved to too loose a tolerance.
So I removed the steady state detection, which has potential issues
with it (see idaholab#26312), and extended the `end_time` until we hit
a fairly tight absolute tolerance
  • Loading branch information
lindsayad authored and schakrabortygithub committed Mar 12, 2024
1 parent 81999bb commit 0c97692
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 8 deletions.
Original file line number Diff line number Diff line change
@@ -1,2 +1,2 @@
time,average_void,max_drag_coefficient,max_x_slip_velocity,max_x_velocity,max_y_slip_velocity,max_y_velocity,min_x_velocity,min_y_velocity
2097.663,0.58378794122993,1.2036091618704,0.0089812524979983,0.12307439685091,0.0092959518687844,0.20305056360613,-0.12307439685089,-0.19525139046703
100000000,0.58378794782072,1.2036091677715,0.0089812664767032,0.1230743501087,0.0092959788633406,0.20305059272475,-0.12307435004431,-0.19525139487663
Original file line number Diff line number Diff line change
Expand Up @@ -48,7 +48,7 @@ g = -9.81
pressure = pressure
[]
[pin_pressure]
type = NSFVPressurePin
type = NSPressurePin
variable = pressure
pin_type = point-value
point = '0 0 0'
Expand Down Expand Up @@ -260,7 +260,7 @@ g = -9.81
[]
[]

[Materials]
[FunctorMaterials]
[CD]
type = NSFVDispersePhaseDragFunctorMaterial
rho = 'rho_mixture'
Expand All @@ -270,7 +270,7 @@ g = -9.81
particle_diameter = ${dp}
[]
[mixing_material]
type = NSFVMixtureMaterial
type = NSFVMixtureFunctorMaterial
phase_1_names = '${rho_d} ${mu_d}'
phase_2_names = '${rho} ${mu}'
prop_names = 'rho_mixture mu_mixture'
Expand Down Expand Up @@ -336,16 +336,15 @@ g = -9.81
type = IterationAdaptiveDT
optimal_iterations = 7
iteration_window = 2
growth_factor = 2.0
growth_factor = 1.5
cutback_factor = 0.5
dt = 1e-3
[]
nl_max_its = 10
steady_state_detection = true
steady_state_tolerance = 1e-10
nl_rel_tol = 1e-03
nl_abs_tol = 1e-11
nl_abs_tol = 1e-9
l_max_its = 5
end_time = 1e8
[]

[Outputs]
Expand Down

0 comments on commit 0c97692

Please sign in to comment.