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.

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

Isothermal domain bc #684

Merged
merged 13 commits into from
Aug 30, 2023
Merged

Conversation

baperry2
Copy link
Contributor

This PR adds a robust isothermal wall capability. It addresses a need that has been brought up in various discussions/issues in the past (#263 #265 #502).

Isothermal walls are handled by separately computing the heat flux through walls after computing other diffusive fluxes (similar to how we handle EB isothermal wall fluxes). This avoids any issues with setting ghost cells that may affect the hydro calculations at the walls and also allows the values to be set directly on the face.

We see good agreement with and 2nd order convergence to the analytical solution for thermal diffusion from a wall into a semi-infinite medium Tnorm = erf(x/sqrt(4*alpha*t)) without hydro. Note the test includes diffusion from all walls in the domain but the data is taken early enough that the effects of various walls haven't started interacting with each other yet.

figure_dir0

figure_dir1

figure_convergence

Draft until documentation is added.

@baperry2 baperry2 requested review from marchdf and hsitaram August 24, 2023 23:19
Copy link
Contributor

@hsitaram hsitaram left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

looking great! thank you Bruce

@baperry2 baperry2 marked this pull request as ready for review August 26, 2023 22:59
Source/Diffusion.cpp Outdated Show resolved Hide resolved
Source/Diffusion.cpp Outdated Show resolved Hide resolved
@baperry2 baperry2 changed the title Isothermal bc Isothermal domain bc Aug 28, 2023
@baperry2 baperry2 enabled auto-merge (squash) August 28, 2023 18:25
@baperry2 baperry2 merged commit e7ccd20 into AMReX-Combustion:development Aug 30, 2023
14 checks passed
@baperry2 baperry2 deleted the isothermal-bc branch August 30, 2023 00:03
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants