[dancfd]

Hello all,

I receive a strange error when trying to run sonicFoam on my case. Eventually I want to run my case using sonicDyMFoam to have an oscillating airfoil. However, at this stage I only wish to test the simulation as a static, steady-state case. To do so, I ran moveDynamicMesh to get the airfoil angle where I wanted it, then copied the relevant /<time>/polyMesh/points file into the /constant/polyMesh directory, and edited the pointDisplacement file to prevent the mesh from moving further. However, when I try to run sonicFoam I get the following error:

Starting time loop

Time = 5e-08

Courant Number mean: 4.389700743e-05 max: 0.006296550206
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0

--> FOAM FATAL ERROR: 
mesh flux field does not exists, is the mesh actually moving?

    From function fvMesh::phi()
    in file fvMesh/fvMeshGeometry.C at line 392.

FOAM exiting

The simulation will run using sonicDyMFoam, though. Why is that? I do not want to solve the cellDisplacement equations if I do not need to, to speed up the simulation. sonicDyMFoam will not run with "steadyState" in the ddtSchemes subdict of system/fvSchemes, therefore I believe I need to use sonicFoam, but this error prohibits that!

Also, since this is a steady-state simulation, what is the best solver to use for compressible flow that will produce the same results as sonicDyMFoam for a steady case? As far as I can tell, sonicDyMFoam = sonicFoam + mesh motion, but which solver = sonicFoam but steady? I would guess rhoSimpleFoam. Note: these cases are all turbulent.

Thanks for the help,
Dan

[dancfd]

Hello all,

I am struggling to run rhoSimpleFoam on a case where sonicDyMFoam runs, but produces poor results. I receive the following error when running rhoSimpleFoam:

Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<eConstThermo<perfectGas>>>>>


--> FOAM FATAL ERROR: 
Not implemented

    From function basicThermo::h()
    in file basicThermo/basicThermo.C at line 259.

FOAM aborting

Can anyone suggest what the problem could be? Note that there are unfortunately no tutorials for this solver.

Thanks,

Dan

[immortality]

change epsithermo to hpsithermo