Skip to content

Commit 2b34c70

Browse files
committed
Robustify Rayleigh Bernard two phase test
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
1 parent 41d8b30 commit 2b34c70

File tree

2 files changed

+7
-8
lines changed

2 files changed

+7
-8
lines changed
Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,2 +1,2 @@
11
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
2-
2097.663,0.58378794122993,1.2036091618704,0.0089812524979983,0.12307439685091,0.0092959518687844,0.20305056360613,-0.12307439685089,-0.19525139046703
2+
100000000,0.58378794782072,1.2036091677715,0.0089812664767032,0.1230743501087,0.0092959788633406,0.20305059272475,-0.12307435004431,-0.19525139487663

modules/navier_stokes/test/tests/finite_volume/two_phase/mixture_model/rayleigh-bernard-two-phase.i

Lines changed: 6 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -48,7 +48,7 @@ g = -9.81
4848
pressure = pressure
4949
[]
5050
[pin_pressure]
51-
type = NSFVPressurePin
51+
type = NSPressurePin
5252
variable = pressure
5353
pin_type = point-value
5454
point = '0 0 0'
@@ -260,7 +260,7 @@ g = -9.81
260260
[]
261261
[]
262262

263-
[Materials]
263+
[FunctorMaterials]
264264
[CD]
265265
type = NSFVDispersePhaseDragFunctorMaterial
266266
rho = 'rho_mixture'
@@ -270,7 +270,7 @@ g = -9.81
270270
particle_diameter = ${dp}
271271
[]
272272
[mixing_material]
273-
type = NSFVMixtureMaterial
273+
type = NSFVMixtureFunctorMaterial
274274
phase_1_names = '${rho_d} ${mu_d}'
275275
phase_2_names = '${rho} ${mu}'
276276
prop_names = 'rho_mixture mu_mixture'
@@ -336,16 +336,15 @@ g = -9.81
336336
type = IterationAdaptiveDT
337337
optimal_iterations = 7
338338
iteration_window = 2
339-
growth_factor = 2.0
339+
growth_factor = 1.5
340340
cutback_factor = 0.5
341341
dt = 1e-3
342342
[]
343343
nl_max_its = 10
344-
steady_state_detection = true
345-
steady_state_tolerance = 1e-10
346344
nl_rel_tol = 1e-03
347-
nl_abs_tol = 1e-11
345+
nl_abs_tol = 1e-9
348346
l_max_its = 5
347+
end_time = 1e8
349348
[]
350349

351350
[Outputs]

0 commit comments

Comments
 (0)