search

loganoz / horses3d

HORSES3D: A high-order discontinuous Galerkin solver for flow simulations and multi-physics applications
https://loganoz.github.io/horses3d/
MIT License
112 stars 24 forks source link

Error in IBM test case (compiled with gcc 11.2.0 with MODE=DEBUG) #84

Closed loganoz closed 11 months ago

loganoz commented 2 years ago

[grubio@flexo IBM_Cylinder]$ ./horses3d.ns IBMCylinder.control ####################################################################################################

HORSES3D High-Order (DG) Spectral Element Sequential Navier-Stokes Solver

#################################################################################################### Compiled Jun 30 2022, 15:11:41 \ Loading common physics................................................................

                          ->                  Reference length (m):  1.000E+00
                          ->                    Reference time (s):  1.467E-02

 \\ Loading Navier-Stokes physics.........................................................

           Fluid data
           ----------
                          ->                 Gas:        Air
                          ->      State constant:    287.150 I.S.
                          -> Specific heat ratio:      1.400

           Reference quantities
           --------------------
                          ->       Reference Temperature:    288.889 K.
                          ->          Reference pressure: 101325.000 Pa.
                          ->           Reference density:      1.221 kg/m^3.
                          ->          Reference velocity:     68.158 m/s.
                          ->         Reference viscosity:   2.08     Pa·s.
                          ->      Reference conductivity:  2.905E+03 W/(m·K).

           Dimensionless quantities
           ------------------------
                          ->              Mach number:      0.200
                          ->          Reynolds number:  4.000E+01
                          ->           Prandtl number:      0.720
                          -> Turbulent Prandtl number:      0.720
                          ->            Froude number: **********
                          ->        Gravity direction: [ 0.0,  0.0,  0.0]

 Reading mesh...

At line 77 of file ZoneClass.f90 Fortran runtime warning: An array temporary was created At line 77 of file ZoneClass.f90 Fortran runtime warning: An array temporary was created At line 121 of file ZoneClass.f90 Fortran runtime warning: An array temporary was created At line 121 of file ZoneClass.f90 Fortran runtime warning: An array temporary was created

 \\ Mesh information......................................................................

           Mesh file "MESH/Cylinder.mesh".
           -------------------------------
                          ->        Number of elements:        668
                          ->           Number of faces:       2689
                          ->  Number of boundary faces:       1370
                          ->     Order of curved faces:          1
                          ->          2D extruded mesh:          T

 \\ Creating zones........................................................................

           Zone 1 for boundary: b2
           -----------------------
                          ->           Number of faces: 6
                          ->   Boundary condition type: Outflow
                          ->          Outflow pressure:  101325.00

           Zone 2 for boundary: back
           -------------------------
                          ->           Number of faces: 668
                          ->   Boundary condition type: FreeSlipWall
                          ->              Thermal type: Adiabatic

           Zone 3 for boundary: b1
           -----------------------
                          ->           Number of faces: 10
                          ->   Boundary condition type: FreeSlipWall
                          ->              Thermal type: Adiabatic

           Zone 4 for boundary: front
           --------------------------
                          ->           Number of faces: 668
                          ->   Boundary condition type: FreeSlipWall
                          ->              Thermal type: Adiabatic

           Zone 5 for boundary: b4
           -----------------------
                          ->           Number of faces: 8
                          ->   Boundary condition type: Inflow
                          ->                  Velocity:      68.16
                          ->               Mach number:       0.20
                          ->                  Pressure:  101325.00
                          ->                   Density:       1.22
                          ->                    AoaPhi:       0.00
                          ->                  AoaTheta:       0.00
                          -> Max. Vel. Fluct. in % (from      0.00

           Zone 6 for boundary: b3
           -----------------------
                          ->           Number of faces: 10
                          ->   Boundary condition type: FreeSlipWall
                          ->              Thermal type: Adiabatic

 \\ IBM parameters........................................................................

           IBM info
           --------
                          ->          Semi implicit treatment:          T
                          ->                Penalization term:         Dt
                          ->        Minimum number of objects:          3
                          ->   Number of interpolation points:         15
                          ->        Surface integration order:          3
                          ->                  Symmetry planes: [-3, 3]

 \\ Reading stl file......................................................................

           Stl file "cylinder3"
           --------------------
                          ->             Number of objects:        796
                          ->       Number of stored points:       2388

At line 560 of file OrientedBoundingBox.f90 Fortran runtime warning: An array temporary was created At line 560 of file OrientedBoundingBox.f90 Fortran runtime warning: An array temporary was created At line 560 of file OrientedBoundingBox.f90 Fortran runtime warning: An array temporary was created At line 561 of file OrientedBoundingBox.f90 Fortran runtime warning: An array temporary was created At line 561 of file OrientedBoundingBox.f90 Fortran runtime warning: An array temporary was created Error in `./horses3d.ns': double free or corruption (!prev): 0x0000000001b8d970 ======= Backtrace: ========= /lib64/libc.so.6(+0x81679)[0x7f2e1de42679] ./horses3d.ns[0xa4d8f4] ./horses3d.ns[0xa57df9] ./horses3d.ns[0xa5ae4c] ./horses3d.ns[0xa938d2] ./horses3d.ns[0x73824b] ./horses3d.ns[0x43e8e4] ./horses3d.ns[0x43f2ec] /lib64/libc.so.6(__libc_start_main+0xf5)[0x7f2e1dde3505] ./horses3d.ns[0x402fc9] ======= Memory map: ======== 00400000-00ff5000 r-xp 00000000 08:10 96636858430 /home/grubio/release_horses_v4/Solver/bin/horses3d.ns 011f4000-011f5000 r--p 00bf4000 08:10 96636858430 /home/grubio/release_horses_v4/Solver/bin/horses3d.ns 011f5000-011ff000 rw-p 00bf5000 08:10 96636858430 /home/grubio/release_horses_v4/Solver/bin/horses3d.ns 011ff000-01202000 rw-p 00000000 00:00 0 01b48000-05327000 rw-p 00000000 00:00 0 [heap] 7f2e04000000-7f2e04021000 rw-p 00000000 00:00 0 7f2e04021000-7f2e08000000 ---p 00000000 00:00 0 7f2e0c000000-7f2e0c026000 rw-p 00000000 00:00 0 7f2e0c026000-7f2e10000000 ---p 00000000 00:00 0 7f2e13ee3000-7f2e13ee4000 ---p 00000000 00:00 0 7f2e13ee4000-7f2e146e4000 rw-p 00000000 00:00 0 7f2e146e4000-7f2e146e5000 ---p 00000000 00:00 0 7f2e146e5000-7f2e14ee5000 rw-p 00000000 00:00 0 7f2e14ee5000-7f2e14ee6000 ---p 00000000 00:00 0 7f2e14ee6000-7f2e156e6000 rw-p 00000000 00:00 0 7f2e156e6000-7f2e156e7000 ---p 00000000 00:00 0 7f2e156e7000-7f2e15ee7000 rw-p 00000000 00:00 0 7f2e15ee7000-7f2e15ee8000 ---p 00000000 00:00 0 7f2e15ee8000-7f2e166e8000 rw-p 00000000 00:00 0 7f2e166e8000-7f2e166e9000 ---p 00000000 00:00 0 7f2e166e9000-7f2e16ee9000 rw-p 00000000 00:00 0 7f2e16ee9000-7f2e16eea000 ---p 00000000 00:00 0 7f2e16eea000-7f2e176ea000 rw-p 00000000 00:00 0 7f2e176ea000-7f2e176eb000 ---p 00000000 00:00 0 7f2e176eb000-7f2e17eeb000 rw-p 00000000 00:00 0 7f2e17eeb000-7f2e17eec000 ---p 00000000 00:00 0 7f2e17eec000-7f2e186ec000 rw-p 00000000 00:00 0 7f2e186ec000-7f2e186ed000 ---p 00000000 00:00 0 7f2e186ed000-7f2e18eed000 rw-p 00000000 00:00 0 7f2e18eed000-7f2e18eee000 ---p 00000000 00:00 0 7f2e18eee000-7f2e196ee000 rw-p 00000000 00:00 0 7f2e196ee000-7f2e196ef000 ---p 00000000 00:00 0 7f2e196ef000-7f2e19eef000 rw-p 00000000 00:00 0 7f2e19eef000-7f2e19ef0000 ---p 00000000 00:00 0 7f2e19ef0000-7f2e1a6f0000 rw-p 00000000 00:00 0 7f2e1a6f0000-7f2e1a6f1000 ---p 00000000 00:00 0 7f2e1a6f1000-7f2e1aef1000 rw-p 00000000 00:00 0 7f2e1aef1000-7f2e1aef2000 ---p 00000000 00:00 0 7f2e1aef2000-7f2e1b6f2000 rw-p 00000000 00:00 0 7f2e1b6f2000-7f2e1b6f3000 ---p 00000000 00:00 0 7f2e1b6f3000-7f2e1bef3000 rw-p 00000000 00:00 0 7f2e1bef3000-7f2e1bef4000 ---p 00000000 00:00 0 7f2e1bef4000-7f2e1c6f4000 rw-p 00000000 00:00 0 7f2e1c6f4000-7f2e1c6f5000 ---p 00000000 00:00 0 7f2e1c6f5000-7f2e1cef5000 rw-p 00000000 00:00 0 7f2e1cef5000-7f2e1cef6000 ---p 00000000 00:00 0 7f2e1cef6000-7f2e1dbbd000 rw-p 00000000 00:00 0 7f2e1dbbd000-7f2e1dbbf000 r-xp 00000000 fd:00 67139299 /usr/lib64/libdl-2.17.so 7f2e1dbbf000-7f2e1ddbf000 ---p 00002000 fd:00 67139299 /usr/lib64/libdl-2.17.so 7f2e1ddbf000-7f2e1ddc0000 r--p 00002000 fd:00 67139299 /usr/lib64/libdl-2.17.so 7f2e1ddc0000-7f2e1ddc1000 rw-p 00003000 fd:00 67139299 /usr/lib64/libdl-2.17.so 7f2e1ddc1000-7f2e1df84000 r-xp 00000000 fd:00 67139293 /usr/lib64/libc-2.17.so 7f2e1df84000-7f2e1e184000 ---p 001c3000 fd:00 67139293 /usr/lib64/libc-2.17.so 7f2e1e184000-7f2e1e188000 r--p 001c3000 fd:00 67139293 /usr/lib64/libc-2.17.so 7f2e1e188000-7f2e1e18a000 rw-p 001c7000 fd:00 67139293 /usr/lib64/libc-2.17.so 7f2e1e18a000-7f2e1e18f000 rw-p 00000000 00:00 0 7f2e1e18f000-7f2e1e1a6000 r-xp 00000000 fd:00 67139319 /usr/lib64/libpthread-2.17.so 7f2e1e1a6000-7f2e1e3a5000 ---p 00017000 fd:00 67139319 /usr/lib64/libpthread-2.17.so 7f2e1e3a5000-7f2e1e3a6000 r--p 00016000 fd:00 67139319 /usr/lib64/libpthread-2.17.so 7f2e1e3a6000-7f2e1e3a7000 rw-p 00017000 fd:00 67139319 /usr/lib64/libpthread-2.17.so 7f2e1e3a7000-7f2e1e3ab000 rw-p 00000000 00:00 0 7f2e1e3ab000-7f2e1e3f0000 r-xp 00000000 fd:03 413456417 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libquadmath.so.0.0.0 7f2e1e3f0000-7f2e1e5f0000 ---p 00045000 fd:03 413456417 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libquadmath.so.0.0.0 7f2e1e5f0000-7f2e1e5f1000 r--p 00045000 fd:03 413456417 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libquadmath.so.0.0.0 7f2e1e5f1000-7f2e1e5f2000 rw-p 00046000 fd:03 413456417 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libquadmath.so.0.0.0 7f2e1e5f2000-7f2e1e609000 r-xp 00000000 fd:03 413501991 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgcc_s.so.1 7f2e1e609000-7f2e1e808000 ---p 00017000 fd:03 413501991 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgcc_s.so.1 7f2e1e808000-7f2e1e809000 r--p 00016000 fd:03 413501991 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgcc_s.so.1 7f2e1e809000-7f2e1e80a000 rw-p 00017000 fd:03 413501991 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgcc_s.so.1 7f2e1e80a000-7f2e1e84a000 r-xp 00000000 fd:03 413456450 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgomp.so.1.0.0 7f2e1e84a000-7f2e1ea49000 ---p 00040000 fd:03 413456450 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgomp.so.1.0.0 7f2e1ea49000-7f2e1ea4a000 r--p 0003f000 fd:03 413456450 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgomp.so.1.0.0 7f2e1ea4a000-7f2e1ea4b000 rw-p 00040000 fd:03 413456450 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgomp.so.1.0.0 7f2e1ea4b000-7f2e1eb4c000 r-xp 00000000 fd:00 67139301 /usr/lib64/libm-2.17.so 7f2e1eb4c000-7f2e1ed4b000 ---p 00101000 fd:00 67139301 /usr/lib64/libm-2.17.so 7f2e1ed4b000-7f2e1ed4c000 r--p 00100000 fd:00 67139301 /usr/lib64/libm-2.17.so 7f2e1ed4c000-7f2e1ed4d000 rw-p 00101000 fd:00 67139301 /usr/lib64/libm-2.17.so 7f2e1ed4d000-7f2e1eff2000 r-xp 00000000 fd:03 413456425 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgfortran.so.5.0.0 7f2e1eff2000-7f2e1f1f1000 ---p 002a5000 fd:03 413456425 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgfortran.so.5.0.0 7f2e1f1f1000-7f2e1f1f2000 r--p 002a4000 fd:03 413456425 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgfortran.so.5.0.0 7f2e1f1f2000-7f2e1f1f4000 rw-p 002a5000 fd:03 413456425 /opt/ohpc/pub/compiler/gcc/11.2.0/lib64/libgfortran.so.5.0.0 7f2e1f1f4000-7f2e1f1f7000 r-xp 00000000 08:10 98785719863 /home/grubio/release_horses_v4/Solver/lib/libproblemfile_ns.so 7f2e1f1f7000-7f2e1f3f6000 ---p 00003000 08:10 98785719863 /home/grubio/release_horses_v4/Solver/lib/libproblemfile_ns.so 7f2e1f3f6000-7f2e1f3f7000 r--p 00002000 08:10 98785719863 /home/grubio/release_horses_v4/Solver/lib/libproblemfile_ns.so 7f2e1f3f7000-7f2e1f3f8000 rw-p 00003000 08:10 98785719863 /home/grubio/release_horses_v4/Solver/lib/libproblemfile_ns.so 7f2e1f3f8000-7f2e1f41a000 r-xp 00000000 fd:00 67139286 /usr/lib64/ld-2.17.so 7f2e1f5fb000-7f2e1f601000 rw-p 00000000 00:00 0 7f2e1f616000-7f2e1f619000 rw-p 00000000 00:00 0 7f2e1f619000-7f2e1f61a000 r--p 00021000 fd:00 67139286 /usr/lib64/ld-2.17.so 7f2e1f61a000-7f2e1f61b000 rw-p 00022000 fd:00 67139286 /usr/lib64/ld-2.17.so 7f2e1f61b000-7f2e1f61c000 rw-p 00000000 00:00 0 7fffd7b13000-7fffd7b35000 rw-p 00000000 00:00 0 [stack] 7fffd7b79000-7fffd7b7b000 r-xp 00000000 00:00 0 [vdso] ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:

0 0x7f2e1ddf73af in ???

1 0x7f2e1ddf7337 in ???

2 0x7f2e1ddf8a27 in ???

3 0x7f2e1de39e86 in ???

4 0x7f2e1de42678 in ???

5 0xa4d8f3 in __tessellationtypes_MOD_pointlinkedlist_destruct

at /home/grubio/release_horses_v4/Solver/src/libs/mesh/TessellationTypes.f90:470

6 0xa57df8 in __orientedboundingbox_MOD_convexhull

at /home/grubio/release_horses_v4/Solver/src/libs/mesh/OrientedBoundingBox.f90:810

7 0xa5ae4b in __orientedboundingbox_MOD_obb_construct

at /home/grubio/release_horses_v4/Solver/src/libs/mesh/OrientedBoundingBox.f90:414

8 0xa938d1 in __ibmclass_MOD_ibm_construct

at /home/grubio/release_horses_v4/Solver/src/libs/mesh/IBMClass.f90:254

9 0x73824a in __dgsemclass_MOD_constructdgsem

at /home/grubio/release_horses_v4/Solver/src/libs/discretization/DGSEMClass.f90:332

10 0x43e8e3 in horses3dmainns

at /home/grubio/release_horses_v4/Solver/src/NavierStokesSolver/main.f90:101

11 0x43f2eb in main

at /home/grubio/release_horses_v4/Solver/src/NavierStokesSolver/main.f90:3

Abortado (`core' generado)