Hello,
i need your feedback about the simulation of a 3 dimensional flow around a cylinder or a NACA-profile with the CompressibleFlowSolver of the NavierStokesCFE or the EulerCFE.
In every case the simulation starts, but always ends at the point of "Writing: "file_0.chk" " due to time-limit of the Cluster (limited to 1 day). I tried several times.
The parameters, Solver Info and so on were shown (see below). In 2D there is not such a problem.
So my question is, if it depends on the calculation capacity (CPU etc.) or what the reason for this? Actually iam calculating on a Cluster with 10 nodes and 120 CPUs. I think, that the calculation normally should be faster and another time-step should be done
after one day or is it normal to a 3D-simulation with the compressible solver? So please can you give me an advice about the Problem?
For example the 3D-problem of the cylinder with the CompressibleFlowSolver with the NavierStokesCFE produces this output after one day of calculation:
(The example with the cylinder is exact the same you published at this page:
http://www.nektar.info/compressible-flow/ )
Parameters:
ALPHA = 0
CINF = 294.827
CP = 1004.7
FINTIME = 2
GAMMA = 1.4
GASCONSTANT = 287.058
IO_CHECKSTEPS = 50000
IO_INFOSTEPS = 1
KAPPA = 0.1
L = 1
MACH = 0.1
MU = 0.00275172
MU0 = 0.5
NUMSTEPS = 2e+06
PINF = 22600
PR = 0.72
RE = 3900
RHOINF = 0.364
SKAPPA = -2.3
THERMALCONDUCTIVITY = 3.83981
TIMESTEP = 1e-06
TINF = 216.29
TWALL = 216.29
UINF = 29.4827
VINF = 0
WINF = 0
Solver Info:
ADVECTIONTYPE = WeakDG
DIFFUSIONTYPE = LDGNS
DRIVER = Standard
EQTYPE = NavierStokesCFE
EVOLUTIONOPERATOR = Nonlinear
GLOBALMATRIXSTORAGETYPE = SmvBSR
GLOBALSYSSOLN = IterativeStaticCond
LINEARPRECONSOLVER = Xxt
LOCALMATRIXSTORAGESTRATEGY = Sparse
PRECONDITIONER = Diagonal
PROBLEMTYPE = General
PROJECTION = DisContinuous
TIMEINTEGRATIONMETHOD = RungeKutta2_ImprovedEuler
UPWINDTYPE = ExactToro
VISCOSITYTYPE = Variable
gs_setup: 240784 unique labels shared
pairwise times (avg, min, max): 0.000148432 0.000120687 0.000157881
crystal router : 0.000947791 0.000929117 0.000967693
used all_to_all method: pairwise
handle bytes (avg, min, max): 241216 152740 313476
buffer bytes (avg, min, max): 160523 101632 208640
gs_setup: 1187072 unique labels shared
pairwise times (avg, min, max): 0.000483757 0.000410008 0.000513387
crystal router : 0.0109781 0.0108864 0.01108
used all_to_all method: pairwise
handle bytes (avg, min, max): 1.18751e+06 656164 1481220
buffer bytes (avg, min, max): 791381 437248 987136
gs_setup: 240784 unique labels shared
pairwise times (avg, min, max): 0.000113803 9.51052e-05 0.000123191
crystal router : 0.000868803 0.000849199 0.000890398
used all_to_all method: pairwise
handle bytes (avg, min, max): 241216 152740 313476
buffer bytes (avg, min, max): 160523 101632 208640
gs_setup: 1187072 unique labels shared
pairwise times (avg, min, max): 0.000463862 0.000412512 0.000496101
crystal router : 0.00953256 0.00941629 0.0096302
used all_to_all method: pairwise
handle bytes (avg, min, max): 1.18751e+06 656164 1481220
buffer bytes (avg, min, max): 791381 437248 987136
gs_setup: 240784 unique labels shared
pairwise times (avg, min, max): 0.000115333 9.5892e-05 0.000125098
crystal router : 0.00325584 0.0032357 0.0032764
used all_to_all method: pairwise
handle bytes (avg, min, max): 241216 152740 313476
buffer bytes (avg, min, max): 160523 101632 208640
gs_setup: 1187072 unique labels shared
pairwise times (avg, min, max): 0.00103487 0.00097959 0.00105419
crystal router : 0.00855757 0.00837069 0.0108672
used all_to_all method: pairwise
handle bytes (avg, min, max): 1.18751e+06 656164 1481220
buffer bytes (avg, min, max): 791381 437248 987136
gs_setup: 240784 unique labels shared
pairwise times (avg, min, max): 0.000114435 8.49009e-05 0.000136399
crystal router : 0.00204957 0.0020185 0.00207429
used all_to_all method: pairwise
handle bytes (avg, min, max): 241216 152740 313476
buffer bytes (avg, min, max): 160523 101632 208640
gs_setup: 1187072 unique labels shared
pairwise times (avg, min, max): 0.000471148 0.0004215 0.000498414
crystal router : 0.00977638 0.0085747 0.0109815
used all_to_all method: pairwise
handle bytes (avg, min, max): 1.18751e+06 656164 1481220
buffer bytes (avg, min, max): 791381 437248 987136
gs_setup: 240784 unique labels shared
pairwise times (avg, min, max): 0.000113707 9.59873e-05 0.000121403
crystal router : 0.000872626 0.000852704 0.00089221
used all_to_all method: pairwise
handle bytes (avg, min, max): 241216 152740 313476
buffer bytes (avg, min, max): 160523 101632 208640
gs_setup: 1187072 unique labels shared
pairwise times (avg, min, max): 0.000467315 0.000390005 0.000493503
crystal router : 0.010846 0.0107322 0.0109444
used all_to_all method: pairwise
handle bytes (avg, min, max): 1.18751e+06 656164 1481220
buffer bytes (avg, min, max): 791381 437248 987136
=======================================================================
EquationType: NavierStokesCFE
Session Name: Cylinder_NS_3D
Spatial Dim.: 3
Max SEM Exp. Order: 4
Expansion Dim.: 3
Riemann Solver: ExactToro
Advection Type:
Projection Type: Discontinuous Galerkin
Diffusion Type:
Advection: explicit
AdvectionType: WeakDG
Diffusion: explicit
Time Step: 1e-06
No. of Steps: 2000000
Checkpoints (steps): 50000
Integration Type: RungeKutta2_ImprovedEuler
Problem Type: General
=======================================================================
Initial Conditions:
- Field rho: rhoInf
- Field rhou: rhoInf * uInf
- Field rhov: rhoInf * vInf
- Field rhow: rhoInf * wInf
- Field E: pInf / (Gamma - 1) + 0.5 * rhoInf * (uInf * uInf + vInf * vInf + wInf * wInf)
Writing: "Cylinder_NS_3D_0.chk"
slurmstepd: *** JOB 930941 CANCELLED AT 2016-07-06T03:27:36 DUE TO TIME LIMIT ***
I run the solver by:
"mpirun -np 120 CompressibleFlowSolver -v 3D_cylinder_mesh.xml 3D_cylinder_conditions.xml"
Thanks for any hint. If you need more information, please let me know .
Best regards
Fabian Selbach
University of Siegen
_______________________________________________
Nektar-users mailing list
Nektar-users@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/nektar-users