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

Running dbnsFoam on a fully cyclic box

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 26, 2019, 10:21
Default Running dbnsFoam on a fully cyclic box
  #1
New Member
 
Thomas
Join Date: Jul 2018
Posts: 16
Rep Power: 7
Roxxor91 is on a distinguished road
Hi everybody,
I want to run "dbnsFoam" on the a domain that is really just a 3D rectangle. I want to apply cyclic/periodic boundary conditions for every primitive variables and I wrote the "boundary" file that way:

Code:
6                                                                                                                                                                                                                                             
(                                                                                                                                                                                                                                             
    inbound                                                                                                                                                                                                                                   
    {                                                                                                                                                                                                                                         
      type            cyclic;//patch;                                                                                                                                                                                                         
        nFaces          8;
        startFace       192;
        matchTolerance 0.001;
        neighbourPatch outbound;
    }
    top
    {
      type            cyclic;//patch;                                                                                                                                                                                                         
        nFaces          10;
        startFace       200;
        matchTolerance 0.001;
        neighbourPatch bottom;
    }
    outbound
    {
      type            cyclic;//patch;                                                                                                                                                                                                         
        nFaces          8;
        startFace       210;
        matchTolerance 0.001;
        neighbourPatch inbound;
    }
    bottom
    {
      type            cyclic;//patch;                                                                                                                                                                                                         
        nFaces          10;
        startFace       218;
        matchTolerance 0.001;
        neighbourPatch top;
    }
    left
    {
      type            cyclic;//patch;                                                                                                                                                                                                         
        nFaces          40;
        startFace       228;
        matchTolerance 0.001;
        neighbourPatch right;
    }
    right
    {
      type            cyclic;//patch;                                                                                                                                                                                                         
        nFaces          40;
        startFace       268;
        matchTolerance 0.001;
        neighbourPatch left;
    }
However, the solver stops at the very beginning of the runTime.loop and shows:

Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field rho

Reading field U

Selecting numericFlux rusanovFluxBarthJespersenLimiter
Reading field MDN

Reading analytic fields


Starting time loop

before compressibleCourantNo
Floating point exception (core dumped)
With GDB I figured out the problem has to come from "cyclicFvPatch::makeDeltaCoeffs(scalarField& dc)" (even though the execution seems to stop while reading "compressibleCourantNo.H" but I have no clue of what theses files are for at the moment.

Do you have any idea of what's causing the issue ?
Roxxor91 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
Setting up cyclic boundary condition KateEisenhower OpenFOAM Pre-Processing 6 January 11, 2017 17:17
LES in OF 3.0.1, no Output for nuSgs Mirage OpenFOAM Programming & Development 1 October 4, 2016 18:00
cyclic boundary violates mass conservation Maik OpenFOAM Running, Solving & CFD 13 August 23, 2016 09:25
Pressure instability with rhoSimpleFoam daniel_mills OpenFOAM Running, Solving & CFD 44 February 17, 2011 17:08
Immersol Simulation of a Heated Box Dong Phoenics 0 March 2, 2006 21:20


All times are GMT -4. The time now is 17:12.