recon1d_plm_cw module reference
Piecewise Linear Method 1D reconstruction.
Data Types
PLM reconstruction following Colella and Woodward, 1984. |
Functions/Subroutines
Initialize a 1D PLM reconstruction for n cells. |
|
Calculate a 1D PLM reconstructions based on h(:) and u(:) |
|
Value of PLM_CW reconstruction at a point in cell k [A]. |
|
Derivative of PLM_CW reconstruction at a point in cell k [A]. |
|
Solver for x such that f(x)=t. |
|
Average between xa and xb for cell k of a 1D PLM reconstruction [A]. |
|
Deallocate the PLM reconstruction. |
|
Checks the PLM_CW reconstruction for consistency. |
|
Runs PLM reconstruction unit tests and returns True for any fails, False otherwise. |
Detailed Description
This implementation of PLM follows Colella and Woodward, 1984 [17], with cells resorting to PCM for extrema including the first and last cells in column. The cell-wise reconstructions are limited so that the edge values (which are also the extrema in a cell) are bounded by the neighboring cell means. This does not yield monotonic profiles for the general remapping problem.
Type Documentation
- type recon1d_plm_cw/plm_cw
PLM reconstruction following Colella and Woodward, 1984.
- Type fields:
% ul ::
real, dimension(:), allocatable, privateLeft edge value [A].% ur ::
real, dimension(:), allocatable, privateRight edge value [A].% init ::
procedure, privateImplementation of the PLM_CW initialization.% reconstruct ::
procedure, privateImplementation of the PLM_CW reconstruction.% average ::
procedure, privateImplementation of the PLM_CW average over an interval [A].% f ::
procedure, privateImplementation of evaluating the PLM_CW reconstruction at a point [A].% dfdx ::
procedure, privateImplementation of the derivative of the PLM_CW reconstruction at a point [A].% x ::
procedure, privateImplementation of solver for x: f(x)=t.% destroy ::
procedure, privateImplementation of deallocation for PLM_CW.% check_reconstruction ::
procedure, privateImplementation of check reconstruction for the PLM_CW reconstruction.% unit_tests ::
procedure, privateImplementation of unit tests for the PLM_CW reconstruction.% init_parent ::
procedure, privateDuplicate interface to% reconstruct_parent ::
procedure, privateDuplicate interface to
Function/Subroutine Documentation
- subroutine recon1d_plm_cw/init(this, n, h_neglect, check)
Initialize a 1D PLM reconstruction for n cells.
- Parameters:
this ::
this[out] This reconstructionn ::
n[in] Number of cells in this columnh_neglect ::
h_neglect[in] A negligibly small width used in cell reconstructions [H]check ::
check[in] If true, enable some consistency checking
- subroutine recon1d_plm_cw/reconstruct(this, h, u)
Calculate a 1D PLM reconstructions based on h(:) and u(:)
- Parameters:
this ::
this[inout] This reconstructionh ::
h[in] Grid spacing (thickness) [typically H]u ::
u[in] Cell mean values [A]
- function recon1d_plm_cw/f(this, k, x)
Value of PLM_CW reconstruction at a point in cell k [A].
- Parameters:
this ::
this[in] This reconstructionk ::
k[in] Cell numberx ::
x[in] Non-dimensional position within element [nondim]
- Call to:
- function recon1d_plm_cw/dfdx(this, k, x)
Derivative of PLM_CW reconstruction at a point in cell k [A].
- Parameters:
this ::
this[in] This reconstructionk ::
k[in] Cell numberx ::
x[in] Non-dimensional position within element [nondim]
- function recon1d_plm_cw/x(this, k, t)
Solver for x such that f(x)=t.
- Parameters:
this ::
this[in] This reconstructionk ::
k[in] Cell numbert ::
t[in] Value to solve for [A]
- Called from:
- function recon1d_plm_cw/average(this, k, xa, xb)
Average between xa and xb for cell k of a 1D PLM reconstruction [A].
- Parameters:
this ::
this[in] This reconstructionk ::
k[in] Cell numberxa ::
xa[in] Start of averaging interval on element (0 to 1)xb ::
xb[in] End of averaging interval on element (0 to 1)
- subroutine recon1d_plm_cw/destroy(this)
Deallocate the PLM reconstruction.
- Parameters:
this ::
this[inout] This reconstruction
- function recon1d_plm_cw/check_reconstruction(this, h, u)
Checks the PLM_CW reconstruction for consistency.
- Parameters:
this ::
this[in] This reconstructionh ::
h[in] Grid spacing (thickness) [typically H]u ::
u[in] Cell mean values [A]
- function recon1d_plm_cw/unit_tests(this, verbose, stdout, stderr)
Runs PLM reconstruction unit tests and returns True for any fails, False otherwise.
- Parameters:
this ::
this[inout] This reconstructionverbose ::
verbose[in] True, if verbosestdout ::
stdout[in] I/O channel for stdoutstderr ::
stderr[in] I/O channel for stderr