Find answers, ask questions, and connect with our <br>community around the world.

Home Forums OpenFOAM Forum taking output after running simulation

  • taking output after running simulation

    Posted by Deleted User on April 12, 2023 at 6:19 am

    Hi I cant take output data from functionobject after running simulation. When I execute “sonicFoam -postProcess -func MachNo” there is no data related to Mach no.
    And same thing happened for Courant number. However I can see Courant number values of each timestep in the log file. Respond of terminal is below.
    Thaks, best regards.

    ——————————————————————————————————————
    Time = 0.0002
    Reading thermophysical properties

    Selecting thermodynamics package
    {
    type hePsiThermo;
    mixture pureMixture;
    transport const;
    thermo hConst;
    equationOfState perfectGas;
    specie specie;
    energy sensibleInternalEnergy;
    }

    Reading field U

    Reading/calculating face flux field phi

    Creating turbulence model

    Selecting turbulence model type RAS
    Selecting RAS turbulence model LaunderSharmaKE
    RAS
    {
    RASModel LaunderSharmaKE;
    turbulence on;
    printCoeffs on;
    Cmu 0.09;
    C1 1.44;
    C2 1.92;
    C3 0;
    sigmak 1;
    sigmaEps 1.3;
    }

    Creating field kinetic energy K

    No MRF models present

    No finite volume options present
    functionObjects::MachNo MachNo writing field: Ma
    ——————————————————————————————————————
    Time = 0.0002
    Reading thermophysical properties

    Selecting thermodynamics package
    {
    type hePsiThermo;
    mixture pureMixture;
    transport const;
    thermo hConst;
    equationOfState perfectGas;
    specie specie;
    energy sensibleInternalEnergy;
    }

    Reading field U

    Reading/calculating face flux field phi

    Creating turbulence model

    Selecting turbulence model type RAS
    Selecting RAS turbulence model LaunderSharmaKE
    RAS
    {
    RASModel LaunderSharmaKE;
    turbulence on;
    printCoeffs on;
    Cmu 0.09;
    C1 1.44;
    C2 1.92;
    C3 0;
    sigmak 1;
    sigmaEps 1.3;
    }

    Creating field kinetic energy K

    No MRF models present

    No finite volume options present
    functionObjects::CourantNo CourantNo writing field: Co

    Deleted User replied 1 year, 2 months ago 2 Members · 15 Replies
  • 15 Replies
  • Deleted User

    Deleted User
    April 12, 2023 at 6:44 am

    And as same when I want to see forcecoefficientts their values shows as 0, although when I saw on terminal it was not zero. I executed ‘sonicFoam -postProcess -func forceCoeffsCompressible’ command.
    Below I share 1)both running time output and 2)functionobject output after simulation.
    ——————————————————————————————————————
    1)
    Time = 0.0006262

    Courant Number mean: 0.000171948 max: 0.244198
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    PIMPLE: iteration 1
    DILUPBiCGStab: Solving for Ux, Initial residual = 8.3762e-08, Final residual = 1.54481e-13, No Iterations 1
    DILUPBiCGStab: Solving for Uy, Initial residual = 1.7342e-07, Final residual = 2.48803e-13, No Iterations 1
    DILUPBiCGStab: Solving for e, Initial residual = 1.66078e-06, Final residual = 4.14365e-13, No Iterations 1
    DILUPBiCGStab: Solving for p, Initial residual = 1.55396e-06, Final residual = 3.33243e-14, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 2.33865e-14, global = 6.28587e-16, cumulative = -9.01354e-09
    DILUPBiCGStab: Solving for p, Initial residual = 6.99048e-10, Final residual = 4.27171e-17, No Iterations 1
    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 1.85762e-16, global = 4.16408e-19, cumulative = -9.01354e-09
    DILUPBiCGStab: Solving for epsilon, Initial residual = 4.36249e-09, Final residual = 1.68404e-13, No Iterations 1
    DILUPBiCGStab: Solving for k, Initial residual = 5.33802e-07, Final residual = 9.12716e-13, No Iterations 1
    ExecutionTime = 1926.23 s ClockTime = 1929 s

    forceCoeffs forces write:
    Coefficient Total Pressure Viscous Internal
    Cd: 0.243837 0.238508 0.00532908 0
    Cd(f): -0.152576 -0.15582 0.00324387 0
    Cd(r): 0.396413 0.394327 0.00208521 0
    Cl: 0.548988 0.550147 -0.00115867 0
    Cl(f): -0.0567121 -0.0562655 -0.000446647 0
    Cl(r): 0.6057 0.606412 -0.000712019 0
    CmPitch: -0.331206 -0.331339 0.000132686 0
    CmRoll: -0.274494 -0.275073 0.000579333 0
    CmYaw: 0.121918 0.119254 0.00266454 0
    Cs: -6.42308e-22 -6.65965e-22 2.3657e-23 0
    Cs(f): 0.121918 0.119254 0.00266454 0
    Cs(r): -0.121918 -0.119254 -0.00266454 0

    ——————————————————————————————————————

    2)// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Create time

    Create mesh for time = 0

    PIMPLE: Operating solver in PISO mode

    forceCoeffs forceCoeffsCompressible:
    –> FOAM Warning :
    From Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const
    in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 910
    Cannot find any patch or group names matching patch1
    –> FOAM Warning :
    From Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const
    in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 910
    Cannot find any patch or group names matching patch2
    Not including porosity effects
    magUInf: 20
    lRef : 1
    Aref : 1
    rhoInf : 1.225
    Selecting all coefficients
    – Drag force: Cd
    – Drag force front: Cd(f)
    – Drag force rear: Cd(r)
    – Lift force: Cl
    – Lift force front: Cl(f)
    – Lift force rear: Cl(r)
    – Pitch moment: CmPitch
    – Roll moment: CmRoll
    – Yaw moment: CmYaw
    – Side force: Cs
    – Side force front: Cs(f)
    – Side force rear: Cs(r)

    Time = 0
    Reading thermophysical properties

    Selecting thermodynamics package
    {
    type hePsiThermo;
    mixture pureMixture;
    transport const;
    thermo hConst;
    equationOfState perfectGas;
    specie specie;
    energy sensibleInternalEnergy;
    }

    Reading field U

    Reading/calculating face flux field phi

    Creating turbulence model

    Selecting turbulence model type RAS
    Selecting RAS turbulence model LaunderSharmaKE
    RAS
    {
    RASModel LaunderSharmaKE;
    turbulence on;
    printCoeffs on;
    Cmu 0.09;
    C1 1.44;
    C2 1.92;
    C3 0;
    sigmak 1;
    sigmaEps 1.3;
    }

    Creating field kinetic energy K

    No MRF models present

    No finite volume options present
    forceCoeffs forceCoeffsCompressible write:
    Coefficient Total Pressure Viscous Internal
    Cd: 0 0 0 0
    Cd(f): 0 0 0 0
    Cd(r): 0 0 0 0
    Cl: 0 0 0 0
    Cl(f): 0 0 0 0
    Cl(r): 0 0 0 0
    CmPitch: 0 0 0 0
    CmRoll: 0 0 0 0
    CmYaw: 0 0 0 0
    Cs: 0 0 0 0
    Cs(f): 0 0 0 0
    Cs(r): 0 0 0 0

    writing force and moment coefficient files.

  • Barış Bicer

    Moderator
    April 12, 2023 at 9:32 pm

    Can you show me your screen shot of run folder and what kind of subfolders are available there. Please also re-run the above command and paste also as screen shot of terminal here.

    sonicFoam -postProcess -func MachNo

  • Barış Bicer

    Moderator
    April 12, 2023 at 9:32 pm

    also tell me which version of OF do you use?

  • Deleted User

    Deleted User
    April 12, 2023 at 10:10 pm

    Hi Hocam, I send the screenshots but it doesn’t show on the site.

  • Deleted User

    Deleted User
    April 12, 2023 at 10:11 pm

    I am using OpenFoam-v2212

  • Deleted User

    Deleted User
    April 12, 2023 at 10:15 pm

  • Deleted User

    Deleted User
    April 12, 2023 at 10:16 pm

  • Barış Bicer

    Moderator
    April 13, 2023 at 9:08 am

    Selam Abdulkadir,

    It looks it works. can you show the inside the output files. just type ls 0.* and put here the screen shot

  • Deleted User

    Deleted User
    April 13, 2023 at 9:54 pm

    Sure hocam, it is below.

  • Barış Bicer

    Moderator
    April 13, 2023 at 10:24 pm

    As you can see in the output folder there is Ma (Mach number) parameter already written as processing.

    • Deleted User

      Deleted User
      April 13, 2023 at 11:13 pm

      I thought the related data would be written in the terminal, thanks

  • Barış Bicer

    Moderator
    April 14, 2023 at 7:42 am

    ok you learnt 🙂

  • Deleted User

    Deleted User
    April 14, 2023 at 10:59 am

    Hi, when I execute postprocess function objects command it create a folder under the postProcessing folder for last time step (t=0.0026). But when I open the .dat file that contains the datas it shows forces and forcecoeff is zero-0. Also there is only few data of each time interval, is it normal?

    Thanks best regards.

  • Deleted User

    Deleted User
    April 14, 2023 at 11:01 am

  • Deleted User

    Deleted User
    April 14, 2023 at 11:01 am

Log in to reply.

error: Content is protected !!