[Thomas Baumann]

Hi everybody,

I try to implement a RSTM in OpenFOAM.

I have problems to implement correctly a laplacian-term of a volSymmTensorField:

my Reynolds stresses R_ are defined as volSymmTensorField

and I need the term:
d/dx_n (constants*R_nm dR_ij/dx_m)

I have implemented it in OF by:
volSymmTensorField laplacianR = fvc::laplacian(constants*R_,R_);

The code can be compiled without any errors,
but the simulation stops always at this position with errors...

I have seen that the RSTM: LaunderGibsonRSTM includes the line
//- fvm::laplacian(Cs*(k_/epsilon_)*R_, R_)
but it is commanded out,
Blending it in results in the same errors:


#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::fv::laplacianScheme<Foam::SymmTensor<double> , Foam::SymmTensor<double> >::New(Foam::fvMesh const&, Foam::Istream&) in "/user/hi204/OpenFOAM/hi204-1.7.0/lib/linux64GccDPOpt/libRSMEBM.so"
#4 Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::laplacian<Foam::SymmTensor<double>, Foam::SymmTensor<double> >(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/user/hi204/OpenFOAM/hi204-1.7.0/lib/linux64GccDPOpt/libRSMEBM.so"
#5 Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::laplacian<Foam::SymmTensor<double>, Foam::SymmTensor<double> >(Foam::tmp<Foam::GeometricField<Foam::SymmTensor< double>, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/user/hi204/OpenFOAM/hi204-1.7.0/lib/linux64GccDPOpt/libRSMEBM.so"
#6 Foam::incompressible::RASModels::RSMEBM::correct() in "/user/hi204/OpenFOAM/hi204-1.7.0/lib/linux64GccDPOpt/libRSMEBM.so"
#7
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/simpleFoam"
#8 __libc_start_main in "/lib/libc.so.6"
#9
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/simpleFoam"
Speicherzugriffsfehler



Has anybody an idea? Do I maybe need a special laplacianSchemes?

Thanks, Thomas