From 28834f029b2c2527c3a18cdfe7fe1fac80fe0f61 Mon Sep 17 00:00:00 2001 From: Philip Mueller <philip.paul.mueller@bluemain.ch> Date: Sun, 2 Feb 2020 09:58:52 +0100 Subject: [PATCH] Added configuration files for the const ball setting. --- .../corello/demoIni/constBallDrag/README.md | 20 ++ .../constBallDrag/test_dragConstBall_jm.ini | 197 ++++++++++++++++++ 2 files changed, 217 insertions(+) create mode 100644 code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/README.md create mode 100644 code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/test_dragConstBall_jm.ini diff --git a/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/README.md b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/README.md new file mode 100644 index 00000000..7f50d19b --- /dev/null +++ b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/README.md @@ -0,0 +1,20 @@ +# Info +This is a similar scenario than the linear momentum test was. +The main difference however is, that this time the ball is +powered. Meaning that a special rehology is used to pull the +ball with constant velocity. Such a setting could be used for +measuring the drag. + +The viscosity of the surrounding media is not zero, but very low. +This is important for measuring the drag properly, since we can +not measure the viscus drag. + +## Integrators +We will only use the jm integrator. + + + + + + + diff --git a/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/test_dragConstBall_jm.ini b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/test_dragConstBall_jm.ini new file mode 100644 index 00000000..d728e13c --- /dev/null +++ b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/test_dragConstBall_jm.ini @@ -0,0 +1,197 @@ +; Test for conservation of the constant ball. +; +; A ball is pulled with constant velocity through a +; low viscosity media. We can then measure +; the drag. The szenario builds on the linearMomentum +; test. +; + + +[Misc] +hdf5_file=corello_constBallDrag_jm.hdf5 + +typeOnlyFirst=yes +NoStrainRateDump=yes +NoStressDump=yes +DumpOnlyRep=yes + +; Probability is larger than one, so dump all. +MarkerSelectionProb=1.2 + +; No density stabilization. +RhoStab=false + +; Dump in intervals of one mili seconds +StepBasedDumping=no +DumpingStepInterval=5 +TimeBasedDumping=yes +MinDumpTimeInterval=1ms + + + +[TimeStepping] + +FinalTime=600ms + +minDeltaTime=40us +maxDeltaTime=0.5ms + +; Setup restriction +nSuperSmallTimeSteps=10 +nRampSteps=40 + +maxDisplacementX=0.1 +maxDisplacementY=0.1 + + +[Integrator] +; Use the min mod lin integrator. +type=jmMod + + +[Domain] +; This will lead of an equal grid spacing in each direction. +; A cell is 1cm in wide. +Nx=301 +Ny=101 + +nMarkerX=5 +nMarkerY=5 +MarkerRandDisplace=yes + +; Only use the full heatingterm +WriteOnylFullHeatTerm=yes + + +[SetUp] +; Use the normal ball scenario +type=constball + +; No gravity is acting +gravY=0.0 + +; The velocitoies that is read in +; negative x means moving to the left. +vx=-4.0 +vy=0.0 + + +[StokesSolver] +; Use the extended Navier-Stokes solver +type=nsext + + +[TemperatureSolver] +; No not use one. +type=nullext +dumpSol=no + + +[StokesBoundary] +; In this section the condition for the stokes boundary conditions +; this are velocities (x & y) and pressure are configured. +; Note that not all options are aviable. + +; Here the conditions on the top can be configured. +; This sets the x velocity condition at the top to +; free slip. +vel_x_top_type=freeSlip +;vel_x_top_value= + +; Here we can controll the bottom condition for the x velocity.. +; We also set it to free slip. +vel_x_bot_type=freeSlip +;vel_x_bot_value= + +; Here we can controll the x velocity at the right boundary, +; these are lower value of x. It is set to dirichlet. +; We force it to zero. +vel_x_rig_type=dirichlet +vel_x_rig_value=0 + +; Here we controll the left side of the domain, these are larger x values. +vel_x_lef_type=dirichlet +vel_x_lef_value=0 + +; Here the conditions on the top can be configured. +; This sets the y velocity condition at the top to +; dirichlet and force it to zero. +vel_y_top_type=dirichlet +vel_y_top_value=0 + +; Here we can controll the bottom condition for the y velocity.. +; We also set it to dirichlet with zero velocity. +vel_y_bot_type=dirichlet +vel_y_bot_value=0 + +; Here we can controll the y velocity at the right boundary, +; these are lower value of x. It is set to free slip. +; We force it to zero. +vel_y_rig_type=freeSlip +;vel_y_rig_value= + +; Here we controll the y velocity at the left side of the domain, +; these are larger x values. +vel_y_lef_type=freeSlip +;vel_y_lef_value= + +; Here we controll the pressure value. This is the linear offset. +; This value is needed since we can only solve for the derivative +; and we need an integration constant. +pressure_value=1.0 + + +[TemperatureBoundary] +; These section controles the temperature boundary of +; the problem. + +; Here we control the temperature condition that is applied +; at the top. We use a dirichlet condition andforce the value to 273K. +; The value here is interpreted in Kelvin. +temp_top_type=dirichlet +temp_top_value=273 + +; Here we controll the temperature that is applied to the +; bottom of the domain. We use a dirichlet condition, the +; value that is used is interpreted as Kelvin. +temp_bot_type=dirichlet +temp_bot_value=1573 + +; Here we controll the temperature condition that is applied +; to the left boundary, these is smaller x values. +; It is set to insulation, this means that the derivation in +; x direction is set to zero, and no temperature flux happens. +temp_lef_type=insulation +;temp_lef_value= + +; Here we controll the temperature condition that is applied +; to the roight boundary. It is set to insulation meaning that +; no temperature flux is allowed. +temp_rig_type=insulation +;temp_rig_value= + + +[Rheology] +; No not use a rheology, also constant marker properties. +type=constball + + +[MarkerToGrid] +; Do not use navive interpolation, interpolate markers directly +; to all grids that are needed. +type=full + + +[GridToMarker] +; Only map velocity back to the marker. +type=velo + +; Use change interpolation. +JobType=c + + +[ApplyBC] +; The standard implementation uses the original conditions. +type=pipe + + -- GitLab