[MaySea]

Hi all,

I'm using OF6.

I have been trying to use the atmBoundaryLayer package to introduce a logarithmic velocity profile at my inlets.

However, no matter how I adjust the ABLConditions file, the velocity field at the inlets is (0 0 0).

This is my ABLConditions file:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

kappa                0.40;          // 
Cmu                  0.09;          // 
flowDir              (0.332934586 -0.803775193 0); 
zDir                 (0 0 1);       
Uref                 0.87;          
Zref                 27;           
z0                   uniform 0.009525;  
d                    uniform 0.0; 
zGround				 uniform 6.56807;    
value                $internalField;

My log, clearly showing that the velocity field is 0, because the residuals are 0 too.

PIMPLE: Convergence criteria found
        U: tolerance 1e-05
        p: tolerance 0.0001
        k: tolerance 1e-05
        epsilon: tolerance 1e-05
        omega: tolerance 1e-05


PIMPLE: Operating solver in PISO mode

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model RNGkEpsilon
RAS
{
    RASModel        RNGkEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.0845;
    C1              1.42;
    C2              1.68;
    C3              0;
    sigmak          0.71942;
    sigmaEps        0.71942;
    eta0            4.38;
    beta            0.012;
}

No MRF models present

No finite volume options present
Courant Number mean: 0 max: 0

Starting time loop

wallShearStress wallShearStress1:
    processing wall patches: 
        barkley

fieldAverage fieldAverage1:
    Starting averaging at time 0

Courant Number mean: 0 max: 0
deltaT = 0.01
Time = 0.01

DILUPBiCGStab:  Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 3
GAMG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 3
time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 3
GAMG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 3
time step continuity errors : sum local = 0, global = 0, cumulative = 0
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.00049805579, Final residual = 2.4428013e-10, No Iterations 2
DILUPBiCGStab:  Solving for k, Initial residual = 1, Final residual = 7.6378181e-10, No Iterations 2
ExecutionTime = 1.4 s  ClockTime = 2 s

Now, I checked the turbineSiting tutorial a couple of times. When I run the tutorial using the allRun script, the situation is the same! Velocity residuals == 0, no velocity field.

[MaySea]

Ok, I solved the issue.

The flowDir components had too many decimal places I think.
I rounded them to 3 decimal places and it worked. Must be related to writePrecision I guess?


Sorry for littering the forum anyway. Maybe someone finds it useful one day.

[lachuktr]

If the flowDir is given as (1 0 0) flow is in x direction. Can you please explain when the digit are not whole number instead decimal digits. ( for example if flowDir is given as (0.5 0.866 0) what is the meaning?)

[MaySea]

Originally Posted by lachuktr

If the flowDir is given as (1 0 0) flow is in x direction. Can you please explain when the digit are not whole number instead decimal digits. ( for example if flowDir is given as (0.5 0.866 0) what is the meaning?)

The resultant vector's direction will be the direction of your flow. It's in Cartesian space. That vector's (0.5 0.866 0) direction will be 60 deg.

[lachuktr]

Can you pls explain ...how you got 60 degree from those numbers