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

Boundary Conditions for a Wind-Tunnel Calculation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 8, 2021, 06:54
Default Boundary Conditions for a Wind-Tunnel Calculation
  #1
New Member
 
Alex
Join Date: Feb 2021
Posts: 5
Rep Power: 5
sepbub is on a distinguished road
Hallo together!
I'm just trying to calculate my first lager project, so maybe I've made some basic mistakes.
So my goal is to calculate an opened Windtunnel (see photo). All walls should have a boundary layer and the model is symmetric. So in my oppinion it should be ok, to calculate a quarter of the tunnel.
When using simpleFoam I'll get the following Error:
Code:
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 15.5344, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 1.26609, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.44793, No Iterations 1000
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#8  Foam::fvMatrix<double>::solve() in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#9  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
Gleitkomma-Ausnahme (Speicherabzug geschrieben)
I assume the reason for that unconvergence are wrong start conditions?!
Here are they:
U:

Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (19 0 0);

boundaryField
{
    symmetry
    {
        type            symmetry;
    }

    wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            fixedValue;
        value           uniform (-8 0 0);
        
    }
    inlet
    {
        type            fixedValue;
        value           uniform (-1.2666 0 0);        
//type            zeroGradient;
    }
}
p:
Code:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    symmetry
    {
        type            symmetry;
    }
    wall
    {
        type            zeroGradient;
    }

    inlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    outlet
    {
        type            zeroGradient;
    }
}
nut:
Code:
dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    symmetry
    {
        type            symmetry;
    }
    wall
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    inlet
    {
        type            calculated;
        value           uniform 0;
    }
    outlet
    {
        type            calculated;
        value           uniform 0;
    }
}
and nutTilda:
Code:
dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0.14;

boundaryField
{
    symmetry
    {
        type            symmetry;
    }
    wall
    {
        type            zeroGradient; 
    }

    inlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    outlet
    {
        type            zeroGradient; 
    }
}
For the completeness my controlDict:
Code:
application     simpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         100000;

deltaT          1;

writeControl    timeStep;

writeInterval   5000;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;
For tubulence modeling I use SpalartAlmaras.

So are there obvious mistakes in my Files, or shall I have a look on some other stuff I didn't considere right now?
Thanks a lot!


(All Mesh Checks went well)
Attached Images
File Type: png Bildschirmfoto vom 2021-03-08 12-42-48.png (16.0 KB, 7 views)
sepbub is offline   Reply With Quote

Reply


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
Problem in setting Boundary Condition Madhatter92 CFX 12 January 12, 2016 04:39
Waterwheel shaped turbine inside a pipe simulation problem mshahed91 CFX 3 January 10, 2015 11:19
Difficulty In Setting Boundary Conditions Moinul Haque CFX 4 November 25, 2014 17:30
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 15:55


All times are GMT -4. The time now is 14:06.