VERB4D
|
Calculates the convection in 1D given a 1D matrix of Phase Space Densities, boundary conditions, diffusion, sources and losses. More...
Go to the source code of this file.
Functions | |
bool | Convection_1D_ULTIMATE_QUICKEST6 (Matrix1D< double > &PSD, Matrix1D< double > x, int x_size, double x_LBC, double x_UBC, string x_LBC_type, string x_UBC_type, Matrix1D< double > Ux, Matrix1D< double > Sources, Matrix1D< double > Losses, double dt) |
Calculates the convection in 1D given a 1D matrix of Phase Space Densities, boundary conditions, diffusion, sources and losses.
All equations and formulas for these calculations can be found at http://www.hadian.ir/teaching/CompHydr/3.pdf . Mostly coming from 3.5 Simplified Ultimate Quickest strategy from B.P Leonard the Ultimate conservative difference scheme. In the source code mathematical equations are numbered corresponding to the numbering of equations found in the paper
bool Convection_1D_ULTIMATE_QUICKEST6 | ( | Matrix1D< double > & | PSD, |
Matrix1D< double > | x, | ||
int | x_size, | ||
double | x_LBC, | ||
double | x_UBC, | ||
string | x_LBC_type, | ||
string | x_UBC_type, | ||
Matrix1D< double > | Ux, | ||
Matrix1D< double > | Sources, | ||
Matrix1D< double > | Losses, | ||
double | dt | ||
) |
Function that implements 1D Convection Quickest and returns a bool upon completion as referenced in the following publication
All equations and formulas for this function can be found at http://www.hadian.ir/teaching/CompHydr/3.pdf Mostly coming from 3.5 Simplified Ultimate Quickest strategy from B.P Leonard the Ultimate conservative difference scheme. In the source code mathematical equations are numbered corresponding to the numbering of equations found in the paper
where c = Courant number
where
and
refer to the right and left faces (between (i and i-1), and (i and i+1)) for PSD
PSD | - 1D matrix of phase space densities |
x | - matrix for determining step size dx = x[1] - x[0] |
x_size | - size of matrix x |
x_LBC,x_UBC,_type | - The upper and lower boundary conditions including the type of boundary condition |
Ux | - diffusion matrix |
Sources | - Sources matrix |
Losses | - Losses (loss cone) |