Skip to content

Regarding dimensions of the domain and simulation killing in Cylinder3D

OpenLB – Open Source Lattice Boltzmann Code Forums on OpenLB General Topics Regarding dimensions of the domain and simulation killing in Cylinder3D

Viewing 9 posts - 16 through 24 (of 24 total)
  • Author
    Posts
  • #9063

    And also when I open the .stl file from SolidWorks I donot get the option of Mesh Files because as per all the youtube videos I saw has this option.

    #9064

    I tried to make a new file but the simulation terminates like this:

    [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.
    [UnitConverter] —————– UnitConverter information —————–
    [UnitConverter] — Parameters:
    [UnitConverter] Resolution: N= 20
    [UnitConverter] Lattice velocity: latticeU= 0.005
    [UnitConverter] Lattice relaxation frequency: omega= 1.88679
    [UnitConverter] Lattice relaxation time: tau= 0.53
    [UnitConverter] Characteristical length(m): charL= 0.1
    [UnitConverter] Characteristical speed(m/s): charU= 0.2
    [UnitConverter] Phys. kinematic viscosity(m^2/s): charNu= 0.002
    [UnitConverter] Phys. density(kg/m^d): charRho= 1
    [UnitConverter] Characteristical pressure(N/m^2): charPressure= 0
    [UnitConverter] Mach number: machNumber= 0.00866025
    [UnitConverter] Reynolds number: reynoldsNumber= 10
    [UnitConverter] Knudsen number: knudsenNumber= 0.000866025
    [UnitConverter]
    [UnitConverter] — Conversion factors:
    [UnitConverter] Voxel length(m): physDeltaX= 0.005
    [UnitConverter] Time step(s): physDeltaT= 0.000125
    [UnitConverter] Velocity factor(m/s): physVelocity= 40
    [UnitConverter] Density factor(kg/m^3): physDensity= 1
    [UnitConverter] Mass factor(kg): physMass= 1.25e-07
    [UnitConverter] Viscosity factor(m^2/s): physViscosity= 0.2
    [UnitConverter] Force factor(N): physForce= 0.0400001
    [UnitConverter] Pressure factor(N/m^2): physPressure= 1600
    [UnitConverter] ————————————————————-
    [prepareGeometry] Prepare Geometry …
    [SuperGeometry3D] cleaned 0 outer boundary voxel(s)
    cylinder3d: ../../../src/functors/analytical/indicator/indicatorF3D.hh:643: olb::IndicatorCuboid3D<T>::IndicatorCuboid3D(olb::Vector<T, 3>, olb::Vector<T, 3>) [with S = float]: Assertion `_xLength>0 && _yLength>0 && _zLength>0′ failed.
    Aborted

    #9065
    Adrian
    Keymaster

    Likely the scaling / base unit of your new file doesn’t match the simulation setup, causing it to be voxelized as an empty volume.

    #9066

    But what is the problem with the given file. Why I am not able to edit it in Solidworks? I changed it to solid body still cannot do it.

    #9067
    Adrian
    Keymaster

    Please stop just reposting the same thing again…

    #9069

    I have edited the post Sir.

    #9070

    Do you mean I need to edit the code also also right?

    #9086
    stephan
    Moderator

    Dear atanuchaudhury,

    thank you for posting. It is likely that you have to adapt your code to the new geometry loaded.
    If this absolutely does not work, please consider using geometry primitives (indicator functors as in the 2D example) instead.

    However, we cannot offer full support on debugging your complete application.
    You might want to consider taking part in our upcoming spring school for this.
    More information is given here: https://www.openlb.net/spring-school-2025/

    BR
    Stephan

    #9157

    I used the setNewSlipboundary, its giving error while running in GPU but doesnot give any error in CPU:
    [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
    [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
    [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
    [BlockGeometryStatistics3D] WARNING: no discreteNormal is found
    [BlockGeometryStatistics3D] WARNING: no discreteNormal is found
    [BlockGeometryStatistics3D] WARNING: no discreteNormal is found
    [BlockGeometryStatistics3D] WARNING: no discreteNormal is found
    [prepareLattice] Prepare Lattice … OK
    [main] starting simulation…
    [setBoundaryValues] step=0; maxVel=0.0499999
    [Timer] step=0; percent=0; passedTime=7.513; remTime=96158.9; MLUPs=0
    [LatticeStatistics] step=0; t=0; uMax=0.0499999; avEnergy=2.49386e-06; avRho=1.0 0011
    [getResults] pressure1=0; pressure2=0; pressureDrop=0; drag=0; lift=0
    cylinder3d: ../../../src/utilities/vectorHelpers.h:204: std::vector<_Tp> olb::ut il::normalize(const std::vector<_Tp>&) [with T = float]: Assertion `scale>0′ fai led.
    Aborted (core dumped)

Viewing 9 posts - 16 through 24 (of 24 total)
  • You must be logged in to reply to this topic.