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

Segmentation Fault, calculating symm(fvc::grad(U)) -custom LES Model library problem?

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 2, 2015, 09:21
Default Segmentation Fault, calculating symm(fvc::grad(U)) -custom LES Model library problem?
  #1
Member
 
Jack
Join Date: Aug 2012
Posts: 47
Rep Power: 7
JackW is on a distinguished road
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:

Code:
#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())),
Code:
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:

Code:
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 is offline   Reply With Quote

Old   March 2, 2015, 17:23
Default
  #2
Member
 
Jack
Join Date: Aug 2012
Posts: 47
Rep Power: 7
JackW is on a distinguished road
Ok, so I have solved my problem. Issue was my fvSolution contained

Code:
cache
(
     grad(U)
);
no idea how this affects anything but there you go! Hope this helps somebody!
JackW is offline   Reply With Quote

Old   August 1, 2017, 08:04
Default
  #3
New Member
 
Felix
Join Date: Sep 2015
Location: Hannover, Berlin
Posts: 23
Rep Power: 4
fimbull is on a distinguished road
I know this thread is old but you simply mad a syntax mistakte.
In your fvSolution you have used:
Code:
cache 
(
    grad(U) 
);
But you must use this:
Code:
cache
{
    grad(U); 
}
fimbull is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
problem loading UDF library in parallel cluster Veera Gutti FLUENT 8 July 26, 2016 07:24
LES and combustion model Margherita Cadorin CFX 0 October 29, 2008 06:24
LES convergence problem: Ning FLUENT 2 March 14, 2007 22:08
problem with MFR model for multiphase mixing tanks Srinivas Main CFD Forum 1 November 7, 2005 15:16
multiphase mixing Problem with MRF model in MixSim Srinivas FLUENT 0 October 17, 2005 06:35


All times are GMT -4. The time now is 09:50.