From: michael Date: Mon, 23 Nov 2020 11:57:05 +0000 (+0100) Subject: Saved the numerical results if the spherical explosion in an hexagonal cavity X-Git-Tag: V9_7_0~140 X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=61e3d8e6f2e9efd17d35f622d0ff34c6f9664645;p=tools%2Fsolverlab.git Saved the numerical results if the spherical explosion in an hexagonal cavity --- diff --git a/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/FinalState_HEXAGON300.png b/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/FinalState_HEXAGON300.png deleted file mode 100644 index 22355b1..0000000 Binary files a/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/FinalState_HEXAGON300.png and /dev/null differ diff --git a/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/SinglePhase_2DSphericalExplosion_HEXAGON.py b/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/SinglePhase_2DSphericalExplosion_HEXAGON.py deleted file mode 100755 index 526bcf8..0000000 --- a/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/SinglePhase_2DSphericalExplosion_HEXAGON.py +++ /dev/null @@ -1,85 +0,0 @@ -#!/usr/bin/env python -# -*-coding:utf-8 -* - -import CoreFlows as cf -import cdmath - - -def SinglePhase_2DSphericalExplosion_HEXAGON(): - - inputfile="./meshHexagonWithTriangles300.med"; - my_mesh=cdmath.Mesh(inputfile); - spaceDim=2 - - # Initial field data - nVar=2+spaceDim; - radius=0.5; - Center=cdmath.Vector(spaceDim);#default value is (0,0,0) - Vout=cdmath.Vector(nVar) - Vin =cdmath.Vector(nVar) - Vin[0]=155e5; - Vin[1]=0; - Vin[2]=0; - Vin[3]=563; - Vout[0]=154e5; - Vout[1]=0; - Vout[2]=0; - Vout[3]=563; - - myProblem = cf.SinglePhase(cf.Liquid,cf.around155bars600K,spaceDim); - - # Initial field creation - print ("Setting mesh and initial data" ) ; - myProblem.setInitialFieldSphericalStepFunction( my_mesh, Vout, Vin, radius, Center); - - # set the boundary conditions - wallVelocityX=0; - wallVelocityY=0; - wallTemperature=563; - - myProblem.setWallBoundaryCondition("boundaries", wallTemperature, wallVelocityX, wallVelocityY); - - # set the numerical method - myProblem.setNumericalScheme(cf.upwind, cf.Explicit); - - # name file save - fileName = "2DSphericalExplosion_HEXAGON"; - - # parameters calculation - MaxNbOfTimeStep = 3000 ; - freqSave = 20; - cfl = 0.5; - maxTime = 5; - precision = 1e-6; - - myProblem.setCFL(cfl); - myProblem.setPrecision(precision); - myProblem.setMaxNbOfTimeStep(MaxNbOfTimeStep); - myProblem.setTimeMax(maxTime); - myProblem.setFreqSave(freqSave); - myProblem.setFileName(fileName); - myProblem.setNewtonSolver(precision,20); - myProblem.saveConservativeField(False); - if(spaceDim>1): - myProblem.saveVelocity(); - pass - - - # evolution - myProblem.initialize(); - - ok = myProblem.run(); - if (ok): - print( "Simulation python " + fileName + " is successful !" ); - pass - else: - print( "Simulation python " + fileName + " failed ! " ); - pass - - print( "------------ End of calculation !!! -----------" ); - - myProblem.terminate(); - return ok - -if __name__ == """__main__""": - SinglePhase_2DSphericalExplosion_HEXAGON() diff --git a/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/pressureWaveHEXAGON300.ogv b/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/pressureWaveHEXAGON300.ogv deleted file mode 100644 index 84ca9ac..0000000 Binary files a/CoreFlows/Documentation/Figures/Simulations/SphericalExploxion_HEXAGON/pressureWaveHEXAGON300.ogv and /dev/null differ diff --git a/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/FinalState_HEXAGON300.png b/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/FinalState_HEXAGON300.png new file mode 100644 index 0000000..22355b1 Binary files /dev/null and b/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/FinalState_HEXAGON300.png differ diff --git a/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/SinglePhase_2DSphericalExplosion_HEXAGON.py b/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/SinglePhase_2DSphericalExplosion_HEXAGON.py new file mode 100755 index 0000000..526bcf8 --- /dev/null +++ b/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/SinglePhase_2DSphericalExplosion_HEXAGON.py @@ -0,0 +1,85 @@ +#!/usr/bin/env python +# -*-coding:utf-8 -* + +import CoreFlows as cf +import cdmath + + +def SinglePhase_2DSphericalExplosion_HEXAGON(): + + inputfile="./meshHexagonWithTriangles300.med"; + my_mesh=cdmath.Mesh(inputfile); + spaceDim=2 + + # Initial field data + nVar=2+spaceDim; + radius=0.5; + Center=cdmath.Vector(spaceDim);#default value is (0,0,0) + Vout=cdmath.Vector(nVar) + Vin =cdmath.Vector(nVar) + Vin[0]=155e5; + Vin[1]=0; + Vin[2]=0; + Vin[3]=563; + Vout[0]=154e5; + Vout[1]=0; + Vout[2]=0; + Vout[3]=563; + + myProblem = cf.SinglePhase(cf.Liquid,cf.around155bars600K,spaceDim); + + # Initial field creation + print ("Setting mesh and initial data" ) ; + myProblem.setInitialFieldSphericalStepFunction( my_mesh, Vout, Vin, radius, Center); + + # set the boundary conditions + wallVelocityX=0; + wallVelocityY=0; + wallTemperature=563; + + myProblem.setWallBoundaryCondition("boundaries", wallTemperature, wallVelocityX, wallVelocityY); + + # set the numerical method + myProblem.setNumericalScheme(cf.upwind, cf.Explicit); + + # name file save + fileName = "2DSphericalExplosion_HEXAGON"; + + # parameters calculation + MaxNbOfTimeStep = 3000 ; + freqSave = 20; + cfl = 0.5; + maxTime = 5; + precision = 1e-6; + + myProblem.setCFL(cfl); + myProblem.setPrecision(precision); + myProblem.setMaxNbOfTimeStep(MaxNbOfTimeStep); + myProblem.setTimeMax(maxTime); + myProblem.setFreqSave(freqSave); + myProblem.setFileName(fileName); + myProblem.setNewtonSolver(precision,20); + myProblem.saveConservativeField(False); + if(spaceDim>1): + myProblem.saveVelocity(); + pass + + + # evolution + myProblem.initialize(); + + ok = myProblem.run(); + if (ok): + print( "Simulation python " + fileName + " is successful !" ); + pass + else: + print( "Simulation python " + fileName + " failed ! " ); + pass + + print( "------------ End of calculation !!! -----------" ); + + myProblem.terminate(); + return ok + +if __name__ == """__main__""": + SinglePhase_2DSphericalExplosion_HEXAGON() diff --git a/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/pressureWaveHEXAGON300.ogv b/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/pressureWaveHEXAGON300.ogv new file mode 100644 index 0000000..84ca9ac Binary files /dev/null and b/CoreFlows/Documentation/NumericalResults/SphericalExploxion_HEXAGON/pressureWaveHEXAGON300.ogv differ