Skip to content

Solver parameters for conservative tracer transport - #758

Merged
jshipton merged 2 commits into
mainfrom
TBendall/ConservativeTracerSolverParams
Aug 27, 2026
Merged

Solver parameters for conservative tracer transport#758
jshipton merged 2 commits into
mainfrom
TBendall/ConservativeTracerSolverParams

Conversation

@tommbendall

Copy link
Copy Markdown
Contributor

This is an attempt to specify sensible parameters for the coupled conservative tracer transport equation:

d(rho)/dt + div(rho*u) = 0
d(m*rho)/dt + div(m*rho*u) = 0

Although this needs to use a NonlinearVariationalSolver (since the m equation has a non-linear time derivative), since we can solve for rho first, the solution should not actually require any non-linear solves.

I also do some reorganising of solver parameters for the base TimeDiscretisation and the ExplicitTimeDiscretisation

@tommbendall tommbendall added the optimisation Involves optimising performance of code label Aug 27, 2026
Comment thread gusto/time_discretisation/time_discretisation.py
@jshipton
jshipton merged commit 3330c70 into main Aug 27, 2026
5 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

optimisation Involves optimising performance of code

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants