[JackW]

Dear all,

In the calculation of symm(fvc::grad(U)) / T(fvc::grad(U)), <insert-anything-with-gradU>,

I am having a bit of nightmare trying to decipher this message:

#0  Foam::error::printStack(Foam::Ostream&) in "/local/software/rh53/OpenFOAM/2.3.0/gcc//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigSegv::sigHandler(int) in "/local/software/rh53/OpenFOAM/2.3.0/gcc//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  ?? in "/lib64/libc.so.6"
#3  Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in "/local/software/rh53/OpenFOAM/2.3.0/gcc//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#4  Foam::incompressible::LESModels::SSG_FSM::B() const in "/home/jw22g11/OpenFOAM/jw22g11-2.3.0/platforms/linux64GccDPOpt/lib/libSSG_FSM.so"
#5  Foam::incompressible::LESModels::SSG_FSM::correct(Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/jw22g11/OpenFOAM/jw22g11-2.3.0/platforms/linux64GccDPOpt/lib/libSSG_FSM.so"
#6  Foam::incompressible::LESModel::correct() in "/local/software/rh53/OpenFOAM/2.3.0/gcc//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so"
#7  main in "/local/software/rh53/OpenFOAM/2.3.0/gcc/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/pimpleFoam"
#8  __libc_start_main in "/lib64/libc.so.6"
#9  __gxx_personality_v0 in "/local/software/rh53/OpenFOAM/2.3.0/gcc/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/pimpleFoam"

Firstly, I am really confused as to why libIncompressibleRASModels.so is being used... I am using my own custom LESModel! (NOTE: when I use a built in LESModel, e.g. DDES, all is fine.) I am starting from a precursor URANS simulation, but the same issue persists if I start from 0/. I am using codedFixedValue at the inlet. My LESModel using k,omega,U and p. All wall resolved (none, set to 0).

The really puzzling part is when I do min and max(fvc::grad(U())),

min/max : min(grad(U)) [0 0 -1 0 0 0 0] (-636.4704112 -5824.846064 -0.04949353017 -601.214397 -4482.509972 -0.05453189162 -0.009910488156 -0.006683723019 -0.003530922467) max(grad(U)) [0 0 -1 0 0 0 0] (4453.55466 1167.020685 0.06365886192 8114.400858 635.6771783 0.04238448374 0.009990205026 0.008183356591 0.005046891309)

you can see they are both bound and issues persist when trying to transpose it, i.e. no actual calculation with the values. This has led me to believe I haven't compiled my LESModel very well. Although, I have been using it, or near enough it for 3 years. The make files look like:

SSG_FSM.C
LIB = $(FOAM_USER_LIBBIN)/libSSG_FSM
EXE_INC = \
    -I$(LIB_SRC)/turbulenceModels \
    -I$(LIB_SRC)/transportModels \
    -I$(LIB_SRC)/turbulenceModels/LES/LESdeltas/lnInclude \
    -I$(LIB_SRC)/turbulenceModels/LES/LESfilters/lnInclude \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude \
    -I$(LIB_SRC)/turbulenceModels/incompressible/LES/lnInclude

LIB_LIBS = \
    -lincompressibleTurbulenceModel \
    -lincompressibleLESModels \
    -lLESdeltas \
    -lLESfilters \
    -lfiniteVolume \
    -lmeshTools

If anyone has any idea what is going on, it would be GREATLY appreciated!

Thank you very much

Jack

[JackW]

Ok, so I have solved my problem. Issue was my fvSolution contained

cache
(
     grad(U)
);

no idea how this affects anything but there you go! Hope this helps somebody!

[fimbull]

I know this thread is old but you simply mad a syntax mistakte.
In your fvSolution you have used:

cache 
(
    grad(U) 
);

But you must use this:

cache
{
    grad(U); 
}