Skip to content

Commit

Permalink
check cases
Browse files Browse the repository at this point in the history
  • Loading branch information
DongWuTUM committed Sep 3, 2023
1 parent d6630c9 commit d495f8e
Show file tree
Hide file tree
Showing 3 changed files with 6 additions and 8 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -112,8 +112,7 @@ int main(int ac, char *av[])
// Define the methods for particle relaxation for wall boundary.
//----------------------------------------------------------------------
SimpleDynamics<RandomizeParticlePosition> wall_boundary_random_particles(wall_boundary);
relax_dynamics::ShellRelaxationStepInner
relaxation_step_wall_boundary_inner(wall_boundary_inner, thickness, level_set_refinement_ratio);
relax_dynamics::ShellRelaxationStepInner relaxation_step_wall_boundary_inner(wall_boundary_inner);
relax_dynamics::ShellNormalDirectionPrediction shell_normal_prediction(wall_boundary_inner, thickness, cos(Pi / 3.75));
wall_boundary.addBodyStateForRecording<int>("UpdatedIndicator");
//----------------------------------------------------------------------
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -83,9 +83,9 @@ int main(int ac, char *av[])
//----------------------------------------------------------------------
SPHSystem sph_system(system_domain_bounds, resolution_ref);
/** Tag for running particle relaxation for the initially body-fitted distribution */
sph_system.setRunParticleRelaxation(true);
sph_system.setRunParticleRelaxation(false);
/** Tag for starting with relaxed body-fitted particles distribution */
sph_system.setReloadParticles(false);
sph_system.setReloadParticles(true);
sph_system.handleCommandlineOptions(ac, av);
IOEnvironment io_environment(sph_system);
//----------------------------------------------------------------------
Expand Down Expand Up @@ -135,8 +135,7 @@ int main(int ac, char *av[])
// Define the methods for particle relaxation for wall boundary.
//----------------------------------------------------------------------
SimpleDynamics<RandomizeParticlePosition> shell_random_particles(shell);
relax_dynamics::ShellRelaxationStepInner
relaxation_step_shell_inner(shell_inner, thickness, level_set_refinement_ratio);
relax_dynamics::ShellRelaxationStepInner relaxation_step_shell_inner(shell_inner);
relax_dynamics::ShellNormalDirectionPrediction shell_normal_prediction(shell_inner, thickness, cos(Pi / 3.75));
shell.addBodyStateForRecording<int>("UpdatedIndicator");
//----------------------------------------------------------------------
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -14,13 +14,13 @@ using namespace SPH;

Real to_rad(Real angle) { return angle * Pi / 180; }

void relax_shell(RealBody &plate_body, Real thickness, Real level_set_refinement_ratio)
void relax_shell(RealBody &plate_body, Real thickness)
{
// BUG: apparently only works if dp > thickness, otherwise ShellNormalDirectionPrediction::correctNormalDirection() throws error

InnerRelation imported_model_inner(plate_body);
SimpleDynamics<RandomizeParticlePosition> random_imported_model_particles(plate_body);
relax_dynamics::ShellRelaxationStepInner relaxation_step_inner(imported_model_inner, thickness, level_set_refinement_ratio);
relax_dynamics::ShellRelaxationStepInner relaxation_step_inner(imported_model_inner);
relax_dynamics::ShellNormalDirectionPrediction shell_normal_prediction(imported_model_inner, thickness);
//----------------------------------------------------------------------
// Particle relaxation starts here.
Expand Down

0 comments on commit d495f8e

Please sign in to comment.