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

V inner boundary solver workflow #2797

Draft
wants to merge 91 commits into
base: master
Choose a base branch
from

Conversation

Rodot-
Copy link
Contributor

@Rodot- Rodot- commented Aug 14, 2024

📝

Adds a workflow to solve for the v_inner_boundary of the simulation

Based on PR #2730 Which should be merged first

Type: 🚀 feature

This adds a complete workflow as a subclass of the SimpleSimulation workflow from PR #2730. This workflow adds convergence criteria and handling the the inner-boundary velocity value. This PR also makes small edits to the configuration schema and parser to allow accessing properties of the v_inner_boundary convergence.

📌 Resources

A notebook is provided in tardis/docs/workflows

🚦 Testing

How did you test these changes?

  • Testing pipeline
  • Other method (describe)
  • My changes can't be tested (explain why)

☑️ Checklist

  • I requested two reviewers for this pull request
  • I updated the documentation according to my changes
  • I built the documentation by applying the build_docs label

Note: If you are not allowed to perform any of these actions, ping (@) a contributor.

andrewfullard and others added 30 commits July 23, 2024 10:44
Solvers pass info to each other
Simplified some parts
Also adds docstrings to methods. Updates some methods to use new functionality of the plasma. Adds requirements for the convergence plots (still broken)
andrewfullard and others added 27 commits August 15, 2024 15:01
Also adds docstrings to methods. Updates some methods to use new functionality of the plasma. Adds requirements for the convergence plots (still broken)

class InnerVelocitySimulationSolver(SimpleSimulation):

TAU_TARGET = np.log(2.0 / 3)
Copy link
Contributor

Choose a reason for hiding this comment

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

For style reasons I would prefer 2. / 3.


TAU_TARGET = np.log(2.0 / 3)

def __init__(self, configuration, mean_optical_depth="rossland", tau=None):
Copy link
Contributor

Choose a reason for hiding this comment

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

It's spelled rosseland

Extrapolate with exponential fits

Need some way to return and inspect the optical depths for later logging"""
pass
Copy link
Contributor

Choose a reason for hiding this comment

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

This does nothing, currently?

self.plasma_solver.plasma_solver_settings.RADIATIVE_RATES_TYPE
== "dilute-blackbody"
):
j_blues = radiation_field.calculate_mean_intensity(
Copy link
Contributor

Choose a reason for hiding this comment

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

Do these not need to be masked in the same way as for detailed rates?

planckian_radiation_field = (
radiation_field.to_planckian_radiation_field()
)
j_blues = planckian_radiation_field.calculate_mean_intensity(
Copy link
Contributor

Choose a reason for hiding this comment

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

Do these not need to be masked in the same way as for detailed rates?

@andrewfullard
Copy link
Contributor

@Rodot- this needs to be rebased now

@andrewfullard andrewfullard marked this pull request as draft November 4, 2024 15:15
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants