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

Fix radiation transport error in Y and Z directions. #633

Merged
merged 14 commits into from
May 28, 2024

Conversation

chongchonghe
Copy link
Contributor

@chongchonghe chongchonghe commented May 27, 2024

Description

Fix bug #632 . This is a dumb mistake I made on using AMREX_D_TERM in AddFluxesRK2, causing a wrong calculation of the HLL flux in 2D/3D cases. The mistake was introduced in the "Use IMEX PD-ARS ..." commit (78c5970) 6 months ago. Because of this mistake, all 2D/3D simulations since then with non-zero y/z-component radiation fluxes are wrong.

To avoid mistakes like this in the future, I added a new test of radiation streaming in the Y direction. I plan to add a better test (say, the RadShock test) in both Y and Z direction, with the following grid configuration: [4, 128, 4] and [4, 4, 128]. Or, alternatively, we should have a non-trivial 3D test with an exact solution.

Related issues

Fixed #632

Checklist

Before this pull request can be reviewed, all of these tasks should be completed. Denote completed tasks with an x inside the square brackets [ ] in the Markdown source below:

  • I have added a description (see above).
  • I have added a link to any related issues see (see above).
  • I have read the Contributing Guide.
  • I have added tests for any new physics that this PR adds to the code.
  • I have tested this PR on my local computer and all tests pass.
  • I have manually triggered the GPU tests with the magic comment /azp run.
  • I have requested a reviewer for this PR.

@chongchonghe
Copy link
Contributor Author

/azp run

Copy link

Azure Pipelines successfully started running 2 pipeline(s).

@chongchonghe
Copy link
Contributor Author

/azp run

Copy link

Azure Pipelines successfully started running 2 pipeline(s).

@chongchonghe chongchonghe requested a review from BenWibking May 27, 2024 14:09
@BenWibking
Copy link
Collaborator

/azp run

Copy link

Azure Pipelines successfully started running 2 pipeline(s).

@chongchonghe
Copy link
Contributor Author

/azp run

Copy link

Azure Pipelines successfully started running 2 pipeline(s).

@chongchonghe
Copy link
Contributor Author

/azp run

Copy link

Azure Pipelines successfully started running 2 pipeline(s).

@chongchonghe chongchonghe requested a review from BenWibking May 28, 2024 07:28
@BenWibking
Copy link
Collaborator

@chongchonghe
Copy link
Contributor Author

/azp run

Copy link

Azure Pipelines successfully started running 2 pipeline(s).

@chongchonghe
Copy link
Contributor Author

OK. All tests passed. Only have to check if the SphericalCollapse test still run on time on nvcc after setting max_time back to 0.05.

@BenWibking BenWibking added this pull request to the merge queue May 28, 2024
Merged via the queue into development with commit 7a597a5 May 28, 2024
21 checks passed
@chongchonghe chongchonghe deleted the chongchong/fix-2D-radiation-issue branch May 30, 2024 11:53
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.

Radiation transport solver gives wrong results in the Y direction after "Use IMEX ..." branch
2 participants