Skip to content
Snippets Groups Projects
Commit 3470d0a7 authored by Philipp Suffa's avatar Philipp Suffa
Browse files

Small change, should work now

parent c106d52f
No related branches found
No related tags found
1 merge request!187Optmizations for kernel generation for the partially saturated cells method
Pipeline #78031 failed
......@@ -18,7 +18,7 @@ class PSMConfig:
Object velocity field for PSM
"""
sc: int = 1
solid_collision: int = 1
"""
Solid collision option for PSM
"""
......@@ -65,7 +65,7 @@ def get_psm_solid_collision_term(collision_rule, psm_config, particle_per_cell_c
for i, (eqFluid, eqSolid, f, offset) in enumerate(
zip(equilibrium_fluid, equilibrium_solid, pre_collision_pdf_symbols, stencil)):
inverse_direction_index = stencil.stencil_entries.index(stencil.inverse_stencil_entries[i])
if psm_config.sc == 1:
if psm_config.solid_collision == 1:
solid_collision = fraction_field.center(particle_per_cell_counter) * (
(
pre_collision_pdf_symbols[inverse_direction_index]
......@@ -73,12 +73,12 @@ def get_psm_solid_collision_term(collision_rule, psm_config, particle_per_cell_c
)
- (f - eqSolid)
)
elif psm_config.sc == 2:
elif psm_config.solid_collision == 2:
# TODO get relaxation rate vector from method and use the right relaxation rate [i]
solid_collision = fraction_field.center(particle_per_cell_counter) * (
(eqSolid - f) + (1.0 - method.relaxation_rates[0]) * (f - eqFluid)
)
elif psm_config.sc == 3:
elif psm_config.solid_collision == 3:
solid_collision = fraction_field.center(particle_per_cell_counter) * (
(
pre_collision_pdf_symbols[inverse_direction_index]
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment