From 7b4edfe41c254477775b9bfba4b5a885133750d5 Mon Sep 17 00:00:00 2001 From: danieljvickers Date: Tue, 20 Jan 2026 20:48:20 -0500 Subject: [PATCH 1/2] Now things are looking pretty --- examples/2D_kelvin_helmholtz/case.py | 1 - examples/2D_richtmyer_meshkov/case.py | 1 - examples/2D_viscous_shock_tube/case.py | 1 - examples/nD_perfect_reactor/analyze.py | 1 - examples/nD_perfect_reactor/export.py | 1 - src/simulation/m_ibm.fpp | 2 +- 6 files changed, 1 insertion(+), 6 deletions(-) diff --git a/examples/2D_kelvin_helmholtz/case.py b/examples/2D_kelvin_helmholtz/case.py index 3a12034696..e6873f8145 100644 --- a/examples/2D_kelvin_helmholtz/case.py +++ b/examples/2D_kelvin_helmholtz/case.py @@ -2,7 +2,6 @@ import json import math - eps = 1e-6 time_end = 1.0 time_save = time_end / 100.0 diff --git a/examples/2D_richtmyer_meshkov/case.py b/examples/2D_richtmyer_meshkov/case.py index efb30f17f4..fed6529c55 100644 --- a/examples/2D_richtmyer_meshkov/case.py +++ b/examples/2D_richtmyer_meshkov/case.py @@ -2,7 +2,6 @@ import json import math - mu = 1.0e-4 lambd = 1.0 time_end = 15.0 diff --git a/examples/2D_viscous_shock_tube/case.py b/examples/2D_viscous_shock_tube/case.py index e61f34cb97..7c760e696b 100644 --- a/examples/2D_viscous_shock_tube/case.py +++ b/examples/2D_viscous_shock_tube/case.py @@ -2,7 +2,6 @@ import json import math - mu = 1.0e-3 time_end = 1.0 time_save = time_end / 100.0 diff --git a/examples/nD_perfect_reactor/analyze.py b/examples/nD_perfect_reactor/analyze.py index ab0c6c0a94..b2bd4e1fa8 100644 --- a/examples/nD_perfect_reactor/analyze.py +++ b/examples/nD_perfect_reactor/analyze.py @@ -6,7 +6,6 @@ import mfc.viz from case import dt, Tend, SAVE_COUNT, sol - case = mfc.viz.Case(".", dt) sns.set_theme(style=mfc.viz.generate_cpg_style()) diff --git a/examples/nD_perfect_reactor/export.py b/examples/nD_perfect_reactor/export.py index f60f55e3fe..112622a32f 100644 --- a/examples/nD_perfect_reactor/export.py +++ b/examples/nD_perfect_reactor/export.py @@ -5,7 +5,6 @@ import mfc.viz from case import dt, NS, Tend, SAVE_COUNT, sol - case = mfc.viz.Case(".", dt) for name in tqdm(sol.species_names, desc="Loading Variables"): diff --git a/src/simulation/m_ibm.fpp b/src/simulation/m_ibm.fpp index 0b8ec1ab81..2617556ff3 100644 --- a/src/simulation/m_ibm.fpp +++ b/src/simulation/m_ibm.fpp @@ -1179,7 +1179,7 @@ contains elseif (patch_ib(ib_marker)%geometry == 8) then ! sphere patch_ib(ib_marker)%moment = 0.4*patch_ib(ib_marker)%mass*(patch_ib(ib_marker)%radius)**2 - else ! we do not have an analytic moment of inertia calculation and need to approximate it directly + else ! we do not have an analytic moment of inertia calculation and need to approximate it directly count = 0 moment = 0._wp cell_volume = (x_cc(1) - x_cc(0))*(y_cc(1) - y_cc(0)) ! computed without grid stretching. Update in the loop to perform with stretching From a89fcdb51cf185691520bc2426d9478b13477430 Mon Sep 17 00:00:00 2001 From: danieljvickers Date: Tue, 20 Jan 2026 20:54:11 -0500 Subject: [PATCH 2/2] Fix format --- src/simulation/m_ibm.fpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/simulation/m_ibm.fpp b/src/simulation/m_ibm.fpp index 2617556ff3..2c5d133060 100644 --- a/src/simulation/m_ibm.fpp +++ b/src/simulation/m_ibm.fpp @@ -1179,7 +1179,7 @@ contains elseif (patch_ib(ib_marker)%geometry == 8) then ! sphere patch_ib(ib_marker)%moment = 0.4*patch_ib(ib_marker)%mass*(patch_ib(ib_marker)%radius)**2 - else ! we do not have an analytic moment of inertia calculation and need to approximate it directly + else ! we do not have an analytic moment of inertia calculation and need to approximate it directly via a sum count = 0 moment = 0._wp cell_volume = (x_cc(1) - x_cc(0))*(y_cc(1) - y_cc(0)) ! computed without grid stretching. Update in the loop to perform with stretching