Commit b1349a8d authored by Jonas Kusch's avatar Jonas Kusch
Browse files

waterphantom first results

parent a8ae768b
Pipeline #103361 passed with stages
in 24 minutes and 33 seconds
......@@ -2,17 +2,17 @@ OUTPUT_DIR = ../result
OUTPUT_FILE = example_csd
LOG_DIR = ../result/logs
MESH_FILE = linesource_pseudo_1D.su2
PROBLEM = LINESOURCE_PSEUDO_1D_PHYSICS
PROBLEM = WATERPHANTOM
SOLVER = CSD_SN_SOLVER
CONTINUOUS_SLOWING_DOWN = NO
HYDROGEN_FILE = ENDL_H.txt
OXYGEN_FILE = ENDL_O.txt
KERNEL = ISOTROPIC_1D
CFL_NUMBER = 0.1
CFL_NUMBER = 1.5
TIME_FINAL = 1.0
CLEAN_FLUX_MATRICES = NO
BC_DIRICHLET = ( void )
BC_NEUMANN = ( adiabatic )
QUAD_TYPE = GAUSS_LEGENDRE_1D
QUAD_ORDER = 10
QUAD_ORDER = 21
......@@ -17,8 +17,8 @@ VectorVector WaterPhantom::SetupIC() {
auto cellMids = _mesh->GetCellMidPoints();
double s = 0.1;
for( unsigned j = 0; j < cellMids.size(); ++j ) {
double x = cellMids[j][0];
psi[j] = 1.0 / ( s * sqrt( 2 * M_PI ) ) * std::exp( -x * x / ( 2 * s * s ) );
double x = cellMids[j][0];
psi[j][_settings->GetNQuadPoints() - 1] = 1.0 / ( s * sqrt( 2 * M_PI ) ) * std::exp( -x * x / ( 2 * s * s ) );
}
return psi;
}
......
......@@ -16,7 +16,7 @@ CSDSNSolver::CSDSNSolver( Config* settings ) : SNSolver( settings ) {
_angle = Vector( _settings->GetNQuadPoints(), 0.0 ); // my
_energies = Vector( _nEnergies, 0.0 ); // equidistant
double energyMin = 1e-1;
double energyMax = 1e0;
double energyMax = 5e0;
// write equidistant energy grid
_dE = ComputeTimeStep( settings->GetCFL() );
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment