Question on resolvedRock3d example

[yyyding]

Hi everyone,

I’m relatively new to OpenLB and currently exploring flow simulations in image-based porous media. The resolvedRock3d example seems perfect for my work. I downloaded the image dataset and converted it to rock.vti using Paraview. However, when I ran the code, it showed “[SuperGeometry3D] error!” and “[BlockGeometryStatistics3D] WARNING: no discreteNormal is found”.

Has anyone else encountered similar issues? How to resolve them?

Thank you in advance.

The outputs are given as follows:
yan@yan:/mnt/c/Users/yan/Desktop/porousMedia/resolvedRock3d$ mpirun -np 1 ./resolvedRock3d rock.vti “Tiff Scalars” 2.5e-6 1.0 200 1.0
[ThreadPool] Sucessfully initialized, numThreads=1
[Directories] Directory ./tmp/ created.
[Directories] Directory ./tmp/imageData/ created.
[Directories] Directory ./tmp/imageData/data/ created.
[Directories] Directory ./tmp/vtkData/ created.
[Directories] Directory ./tmp/vtkData/data/ created.
[Directories] Directory ./tmp/gnuplotData/ created.
[Directories] Directory ./tmp/gnuplotData/data/ created.
[main] Using the following input:
[main] VTI file: rock.vti
[main] Array name: Tiff Scalars
[main] Pressure drop: 1 Pa
[main] Rock x-length: 0.00050125 m
[main] Rock y-length: 0.00050125 m
[main] Rock z-length: 0.00050125 m
[main] Timeframe to be simulated: 0.00251252 s
[UnitConverter] —————– UnitConverter information —————–
[UnitConverter] — Parameters:
[UnitConverter] Resolution: N= 200
[UnitConverter] Lattice velocity: latticeU= 0.0416667
[UnitConverter] Lattice relaxation frequency: omega= 1.33333
[UnitConverter] Lattice relaxation time: tau= 0.75
[UnitConverter] Characteristical length(m): charL= 0.00050125
[UnitConverter] Characteristical speed(m/s): charU= 0.199501
[UnitConverter] Phys. kinematic viscosity(m^2/s): charNu= 1e-06
[UnitConverter] Phys. density(kg/m^d): charRho= 1000
[UnitConverter] Characteristical pressure(N/m^2): charPressure= 0
[UnitConverter] Mach number: machNumber= 0.0721688
[UnitConverter] Reynolds number: reynoldsNumber= 100
[UnitConverter] Knudsen number: knudsenNumber= 0.000721688
[UnitConverter]
[UnitConverter] — Conversion factors:
[UnitConverter] Voxel length(m): physDeltaX= 2.50625e-06
[UnitConverter] Time step(s): physDeltaT= 5.23441e-07
[UnitConverter] Velocity factor(m/s): physVelocity= 4.78803
[UnitConverter] Density factor(kg/m^3): physDensity= 1000
[UnitConverter] Mass factor(kg): physMass= 1.57425e-14
[UnitConverter] Viscosity factor(m^2/s): physViscosity= 1.2e-05
[UnitConverter] Force factor(N): physForce= 1.44e-07
[UnitConverter] Pressure factor(N/m^2): physPressure= 22925.2
[UnitConverter] ————————————————————-
[prepareGeometry] Prepare Geometry …
[SuperGeometry3D] cleaned 4095788 outer boundary voxel(s)
[SuperGeometry3D] cleaned 4 outer boundary voxel(s)
[SuperGeometryStatistics3D] updated
[SuperGeometry3D] error!
[CuboidGeometry3D] —Cuboid Stucture Statistics—
[CuboidGeometry3D] Number of Cuboids: 7
[CuboidGeometry3D] Delta (min): 2.50625e-06
[CuboidGeometry3D] (max): 2.50625e-06
[CuboidGeometry3D] Ratio (min): 0.33
[CuboidGeometry3D] (max): 2.32558
[CuboidGeometry3D] Nodes (min): 1117314
[CuboidGeometry3D] (max): 1162914
[CuboidGeometry3D] Weight (min): 310505
[CuboidGeometry3D] (max): 550259
[CuboidGeometry3D] ——————————–
[SuperGeometryStatistics3D] materialNumber=0; count=4095792; minPhysR=(0,0,0); maxPhysR=(0.000498744,0.000498744,0.000498744)
[SuperGeometryStatistics3D] materialNumber=1; count=2718412; minPhysR=(2.50625e-06,0,0); maxPhysR=(0.000496238,0.000498744,0.000498744)
[SuperGeometryStatistics3D] materialNumber=2; count=1151881; minPhysR=(2.50625e-06,0,0); maxPhysR=(0.000496238,0.000498744,0.000498744)
[SuperGeometryStatistics3D] materialNumber=3; count=17616; minPhysR=(0,0,0); maxPhysR=(0,0.000498744,0.000498744)
[SuperGeometryStatistics3D] materialNumber=4; count=16299; minPhysR=(0.000498744,0,0); maxPhysR=(0.000498744,0.000498744,0.000498744)
[SuperGeometryStatistics3D] countTotal[1e6]=8
[prepareGeometry] Prepare Geometry … OK
[prepareLattice] Prepare Lattice …
[BlockGeometryStatistics3D] WARNING: no discreteNormal is found
[BlockGeometryStatistics3D] WARNING: no discreteNormal is found
[BlockGeometryStatistics3D] WARNING: no discreteNormal is found
[BlockGeometryStatistics3D] WARNING: no discreteNormal is found
[BlockGeometryStatistics3D] WARNING: no discreteNormal is found
…

Some content of converted rock.vti:
<VTKFile type=”ImageData” version=”1.0″ byte_order=”LittleEndian” header_type=”UInt64″>
<ImageData WholeExtent=”0 511 0 511 0 511″ Origin=”0 0 0″ Spacing=”1 1 1″ Direction=”1 0 0 0 1 0 0 0 1″>
<Piece Extent=”0 511 0 511 0 511″>
<PointData Scalars=”Tiff Scalars”>
<DataArray type=”UInt8″ Name=”Tiff Scalars” NumberOfComponents=”3″ format=”ascii” RangeMin=”0″ RangeMax=”441.6729559300637″>
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 255 255 255
255 255 255 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
…
<InformationKey name=”L2_NORM_RANGE” location=”vtkDataArray” length=”2″>
<Value index=”0″>
0
</Value>
<Value index=”1″>
441.67295593
</Value>
</InformationKey>
</DataArray>
</PointData>
<CellData>
</CellData>
</Piece>
</ImageData>
</VTKFile>

[FBukreev]

Hello,

if the simulation still runs after these warnings, then its ok. But if not, then it can lie on the type of boundary conditions. It is better to choose the Local type. If the problem still exists there, then you can switch to EquilibriumBoundary.

[Author]

Posts