[Gennaro]

Hi all,

I'm working on a superFluid helium solver based on Laplacian Foam.

The difference between LaplacianFoam and my solver is that the thermal conductivity k is replaced with keff which is a function of the gradient of the temperature.

This model was developed a few years ago and it is used in the cryogenic industry, but no one implemented it on OpenFOAM yet. I'm trying to do it although I'm relatively new to OpenFOAM programming and C++ programming.

So I created the solver, which compiles but crashes returning the following error:

Create time
Create mesh for time = 0
Reading field T
Reading thermophysicalProperties
Selecting thermodynamics package
{
type heSolidThermo;
mixture pureMixture;
transport polynomial;
thermo hPolynomial;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}
Reading Alambda
Reading Slambda
Reading Tlambda
No finite volume options present
Radiation model not active: radiationProperties not found
Selecting radiationModel none
Calculating temperature distribution
SIMPLE: no convergence criteria found. Calculations will run for 0.01 steps.
Time = 1e-06
#0 Foam::error:rintStack(Foam::Ostream&) in "/afs/cern.ch/project/cfd/OpenFOAM/SLC6/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/afs/cern.ch/project/cfd/OpenFOAM/SLC6/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/afs/cern.ch/project/cfd/OpenFOAM/SLC6/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/afs/cern.ch/user/g/gbozza/OpenFOAM/gbozza-2.3.x/platforms/linux64Gcc46DPOpt/bin/laplacianFoamCustom"
#5
at laplacianFoamCustom.C:0
#6
in "/afs/cern.ch/user/g/gbozza/OpenFOAM/gbozza-2.3.x/platforms/linux64Gcc46DPOpt/bin/laplacianFoamCustom"
#7 __libc_start_main in "/lib64/libc.so.6"
#8
in "/afs/cern.ch/user/g/gbozza/OpenFOAM/gbozza-2.3.x/platforms/linux64Gcc46DPOpt/bin/laplacianFoamCustom"
Floating exception (core dumped)

I don't understand what SIMPLE: no convergence criteria found means.

Please find in attachment the solver and a case.

Can you please help?

Thanks

Best regards

[libi.macek]

Hello Gennaro,

Did you make any progressi in this problem?
I am interested in computation with superfluid.

Best regards. Libor

[Cyp]

Dear Libor,

please have a look at some papers we published with HellFOAM our superfluid helium simulator with OpenFOAM:

http://www.sciencedirect.com/science...10465514003415

https://journals.aps.org/prl/abstrac...ett.118.074506

http://oatao.univ-toulouse.fr/11935/...aine_11935.pdf


The first one described the algorithm we developed in OpenFOAM and some validation test cases. We uploaded the source code and you can download it here https://github.com/csoulain/HellFOAM
I haven't upgrade the code with the new OpenFOAM release, but I am sure you can do it easily.

Cheers,
Cyprien

[libi.macek]

Thank you for advice