AMReX-Astro / AMReX-Astro/Microphysics
Add new components for radiation energy density
Nobody has claimed this yet.
- Dominant language
- C++
- Stars
- 43
- Forks
- 46
- Avg merge
- 2d 18h
- Merged PRs (30d)
- 15
Description
To implement a photoionization network, we will need radiation energy density as a state components, call it state.rad. It has a RHS (ydot) that depends on T, state.xn, and state.rad, and it should evolve along with state.xn and T. However, it differs from state.xn in three different ways:
- it does not have aion and zion;
- it won't affect the EOS
- it's not limited to 0 -- 1.
Is there an existing component that can work this way? Any suggestion on how I should go about adding this new component?
It seems to me that the auxiliary composition won't work, because it does not have a RHS and Jacobian, but I might be missing something.
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by reviewing the existing state.xn component and the auxiliary composition implementation, especially how RHS and Jacobian support are defined. Determine which component interface can represent state.rad without aion, zion, EOS coupling, or 0–1 bounds. Done means a documented implementation path is agreed for evolving radiation energy density with state.xn and T.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- cpp
- Domain
- backend
- Issue type
- Feature
- Difficulty
- 5/5
- Estimated time
- Over a week
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100