[mdiazv64]

Dear nipinl and jan,

Thank you very much for your quick reply. In my case I want to estimate the permeability of a synthetic porous medium with cube shape in the z direction. The boundary conditions would be the pressure differential at z and non-flow on the other faces. But I wouldn’t know how to modify the boundary conditions in the code.
I would appreciate if you could guide me on this matter.

Best regards,

Martin

[mdiazv64]

Hello. I would appreciate support regarding the application of the Lattice Boltzmann method in porous media flow. I am using the resolvedRock3d example applied to an image following the procedure explained in the manual but it does not converge. It runs around 3% (step=4185), but from then on it gives NAN and I don’t understand why that is.

In my case I have the following parameters:

Timeframe to be simulated: 29.7014 s
MaxIT: 139908
—————– UnitConverter information —————–
–Parameters:
Resolution: N= 100
Lattice velocity: latticeU= 0.000714753
Lattice relaxation frequency: omega= 1.33333
Lattice relaxation time: tau= 0.75
Characteristic length(m): charL= 0.0509
Characteristical speed(m/s): charU= 0.00171372
Phys. kinematic viscosity(m^2/s): charNu= 0.0001017
Phys. density(kg/m^d): charRho= 963.8
Characteristical pressure(N/m^2): charPressure= 0
Mach number: machNumber= 0.00123799
Reynolds number: reynoldsNumber= 0.857704
Knudsen number: knudsenNumber= 0.00144338
— Conversion factors:
Voxel length(m): physDeltaX= 0.000509
Time step(s): physDeltaT= 0.000212292
Velocity factor(m/s): physVelocity= 2.39764
Density factor(kg/m^3): physDensity= 963.8
Mass factor(kg): physMass= 1.27098e-07
Viscosity factor(m^2/s): physViscosity= 0.0012204
Force factor(N): physForce= 0.00143546
Pressure factor(N/m^2): physPressure= 5540.59
————————————————– ———–

Please, I would appreciate any suggestions on how to proceed.

[Author]

Posts