Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

run isothermal Soret channel as a test #441

Draft
wants to merge 4 commits into
base: development
Choose a base branch
from
Draft

run isothermal Soret channel as a test #441

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
b1eb277
run isothermal Soret channel as a test
baperry2 Nov 20, 2024
9488262
fix warns
baperry2 Nov 20, 2024
8823cbb
remove unused input from isothermal soret
baperry2 Dec 2, 2024
db25cf2
remove more unused
baperry2 Dec 2, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions Exec/RegTests/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -21,6 +21,7 @@ else()
endif()
if(PELE_DIM EQUAL 2)
add_subdirectory(TripleFlame)
add_subdirectory(IsothermalSoretChannel)
endif()
if(PELE_DIM EQUAL 3)
add_subdirectory(TurbInflow)
Expand Down
10 changes: 5 additions & 5 deletions ...RegTests/IsothermalSoretChannel/input.inp → ...othermalSoretChannel/isothermal-soret.inp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,7 @@ prob.P_mean = 101325.0
prob.Vin = 10.0
prob.Tlow = 400.0
prob.Thigh = 800.0
prob.phi = 0.5
prob.phi = 1.0
#---------------------- Transport --------------------------------
transport.use_soret = 1

Expand All @@ -46,7 +46,7 @@ peleLM.do_species_balance = 1 # Compute species balance
amr.max_step = 2000 # Maximum number of time steps
amr.stop_time = 4.00 # final simulation physical time [s]
amr.max_wall_time = 1 # Maximum simulation run time [h]
amr.cfl = 0.5 # CFL number for hyperbolic system
amr.cfl = 0.5 # CFL number for hyperbolic system
amr.dt_shrink = 0.01 # Scale back initial timestep
amr.dt_change_max = 1.1 # Maximum dt increase btw successive steps

Expand Down Expand Up @@ -86,6 +86,6 @@ amr.gradT.adjacent_difference_greater = 200
amr.gradT.field_name = temp

#---------------------- Debug/HPC CONTROL-------------------------
amrex.fpe_trap_invalid = 1
amrex.fpe_trap_zero = 1
amrex.fpe_trap_overflow = 1
#amrex.fpe_trap_invalid = 1
#amrex.fpe_trap_zero = 1
#amrex.fpe_trap_overflow = 1
10 changes: 2 additions & 8 deletions Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -59,11 +59,8 @@ pelelmex_initdata(
const amrex::Real* prob_hi = geomdata.ProbHi();
const amrex::Real* dx = geomdata.CellSize();

const amrex::Real x = prob_lo[0] + (i + 0.5) * dx[0];
const amrex::Real y = prob_lo[1] + (j + 0.5) * dx[1];
const amrex::Real y_low = prob_lo[1];

const amrex::Real Lx = prob_hi[0] - prob_lo[0];
const amrex::Real Ly = prob_hi[1] - prob_lo[1];

auto eos = pele::physics::PhysicsType::eos();
Expand Down Expand Up @@ -98,19 +95,16 @@ AMREX_GPU_HOST_DEVICE
AMREX_FORCE_INLINE
void
bcnormal(
const amrex::Real x[AMREX_SPACEDIM],
const amrex::Real* /*x[AMREX_SPACEDIM]*/,
const int /*m_nAux*/,
amrex::Real s_ext[NVAR],
const int /*idir*/,
const int sgn,
const amrex::Real time,
const amrex::Real /*time*/,
amrex::GeometryData const& /*geomdata*/,
ProbParm const& prob_parm,
pele::physics::PMF::PmfData::DataContainer const* /*pmf_data*/)
{
auto eos = pele::physics::PhysicsType::eos();
amrex::Real massfrac[NUM_SPECIES] = {0.0};

if (sgn == 1) {
s_ext[TEMP] = prob_parm.Tlow;
} else if (sgn == -1) {
Expand Down
1 change: 1 addition & 0 deletions Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,6 +10,7 @@ PeleLM::readProbParm()
pp.query("Vin", PeleLM::prob_parm->Vin);
pp.query("Thigh", PeleLM::prob_parm->Thigh);
pp.query("Tlow", PeleLM::prob_parm->Tlow);
pp.query("phi", PeleLM::prob_parm->phi);

PeleLM::prob_parm->bathID = N2_ID;
PeleLM::prob_parm->oxidID = O2_ID;
Expand Down
1 change: 1 addition & 0 deletions Tests/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -236,6 +236,7 @@ if(NOT PELE_ENABLE_EB)
endif()
if(PELE_DIM EQUAL 2)
add_test_r(tripleflame-${PELE_DIM}d TripleFlame)
add_test_rv(isothermal-soret IsothermalSoretChannel)
endif()
if(PELE_DIM EQUAL 3)
add_test_r(hit-${PELE_DIM}d HITDecay)
Expand Down
Loading