CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Problem with running atmBoundaryLayer cases

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 13, 2021, 14:27
Default Problem with running atmBoundaryLayer cases
  #1
Member
 
Join Date: Nov 2018
Posts: 30
Rep Power: 4
MaySea is on a distinguished road
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:
Code:
/*--------------------------------*- 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.

Code:
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 is offline   Reply With Quote

Old   January 14, 2021, 08:10
Default Solved
  #2
Member
 
Join Date: Nov 2018
Posts: 30
Rep Power: 4
MaySea is on a distinguished road
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.
MaySea is offline   Reply With Quote

Reply

Tags
atmboundarylayer

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
ANSYS Licensing Problem, Processes Running but Showing as Not Running penguinman ANSYS 3 September 27, 2016 14:30
[waves2Foam] Problem in running Allrun and also setWaveFields ankitchy OpenFOAM Community Contributions 2 March 9, 2015 09:10
[Other] Problem running after decomposition paul.m.edwards OpenFOAM Meshing & Mesh Conversion 0 January 14, 2015 13:55
Problem while running in Highperformance computing environment Phanipavan STAR-CD 1 September 11, 2013 07:42
problem with running in parallel dhruv OpenFOAM 3 November 25, 2011 06:06


All times are GMT -4. The time now is 19:54.