From e85d09653b2884dc89bc7bb16292da0a2743e426 Mon Sep 17 00:00:00 2001
From: dy94rovu <ravi.k.ayyala@fau.de>
Date: Tue, 15 Apr 2025 16:38:31 +0200
Subject: [PATCH] set continuous_equilibrium to false and force in z-direction
---
.../thermal_PSM_particles/thermalPSMCodegenParticles.py | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/apps/benchmarks/MaterialTransport/thermal_PSM/thermal_PSM_particles/thermalPSMCodegenParticles.py b/apps/benchmarks/MaterialTransport/thermal_PSM/thermal_PSM_particles/thermalPSMCodegenParticles.py
index 0add80ff7..8e903f83d 100644
--- a/apps/benchmarks/MaterialTransport/thermal_PSM/thermal_PSM_particles/thermalPSMCodegenParticles.py
+++ b/apps/benchmarks/MaterialTransport/thermal_PSM/thermal_PSM_particles/thermalPSMCodegenParticles.py
@@ -107,8 +107,8 @@ with CodeGeneration() as ctx:
B = ps.fields(f"b({1}): {data_type}[3D]", layout=layout)
- force_concentration_on_fluid = sp.Matrix([0, (rho_0)*alpha*(concentration_field.center - T0)*gravity_LBM,0])
- #force_concentration_on_fluid = sp.Matrix([0, 0,(rho_0)*alpha*(concentration_field.center - T0)*gravity_LBM])
+ #force_concentration_on_fluid = sp.Matrix([0, (rho_0)*alpha*(concentration_field.center - T0)*gravity_LBM,0])
+ force_concentration_on_fluid = sp.Matrix([0, 0,(rho_0)*alpha*(concentration_field.center - T0)*gravity_LBM])
# Fluid LBM optimisation
lbm_fluid_opt = LBMOptimisation(
@@ -190,6 +190,8 @@ with CodeGeneration() as ctx:
output={"density": concentration_field},
compressible=True,
psm_config=psm_config_C,
+ continuous_equilibrium=False,
+ zero_centered=False,
)
if config_tokens[0] == "srt-smagorinsky" or config_tokens[0] == "trt-smagorinsky":
--
GitLab