Skip to content

Commit

Permalink
update the SDC docs
Browse files Browse the repository at this point in the history
  • Loading branch information
zingale committed Feb 5, 2024
1 parent 49fba8d commit 7b21213
Showing 1 changed file with 31 additions and 2 deletions.
33 changes: 31 additions & 2 deletions sphinx_docs/source/sdc.rst
Original file line number Diff line number Diff line change
Expand Up @@ -5,6 +5,29 @@ Spectral Deferred Corrections
Introduction
============

Spectral deferred correction (SDC) methods strongly couple reactions
and hydrodynamics, eliminating the splitting error that arises with
Strang (operator) splitting. Microphysics supports two different
SDC formulations.

.. note::

Both SDC formulations are supported by the same integrators, enabled
by the ``SDC`` preprocessor flag.

"True" SDC
----------

The true SDC implementation is described in :cite:`castro_sdc`. It divides
the timestep into temporal nodes and uses low-order approximations to update
from one temporal node to the next. Iteration is used to increase the order of accuracy.

The update from one temporal node to the next appears as:


Simplified SDC
--------------

The Simplified-SDC method provides a means to more strongly couple the
reactions to the hydrodynamics by evolving the reactions together with
an approximation of the advection over the timestep. The full details
Expand Down Expand Up @@ -128,8 +151,8 @@ system we are integrating, including the advective terms.
:math:`\dot{Y}_k` and the nuclear energy release, :math:`\dot{S}`.

#. Convert back to the integrator’s internal representation via ``rhs_to_int``
This converts the ``ydot``s to mass fractions and adds the advective terms
to all ``ydots``.
This converts the ``ydot`` to mass fractions and adds the advective terms
to ``ydot``.

Jacobian
--------
Expand All @@ -147,3 +170,9 @@ the Jacobian as:
where :math:`{\bf w} = (X_k, T)^\intercal` are the more natural variables
for a reaction network.

.. note::

In the original "true SDC" paper (:cite:`castro_sdc`), the matrix
system as more complicated, and we included density in ${\bf w}$.
This is not needed, and we use the Jacobian defined in
:cite:`castro_simple_sdc` instead.

0 comments on commit 7b21213

Please sign in to comment.