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

Solver for superfluid helium based on Laplacian Foam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 19, 2014, 13:10
Default Solver for superfluid helium based on Laplacian Foam
  #1
New Member
 
Gennaro
Join Date: May 2014
Posts: 23
Rep Power: 8
Gennaro is on a distinguished road
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
Attached Files
File Type: gz laplacianFoamCustom.tar.gz (84.9 KB, 5 views)
Gennaro is offline   Reply With Quote

Old   October 3, 2017, 04:56
Default
  #2
New Member
 
Libor Macek
Join Date: Mar 2017
Location: Brno, Czech republic
Posts: 6
Rep Power: 5
libi.macek is on a distinguished road
Hello Gennaro,

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

Best regards. Libor
libi.macek is offline   Reply With Quote

Old   October 5, 2017, 09:01
Default
  #3
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 298
Rep Power: 14
Cyp is on a distinguished road
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
Cyp is offline   Reply With Quote

Old   October 26, 2017, 04:18
Default
  #4
New Member
 
Libor Macek
Join Date: Mar 2017
Location: Brno, Czech republic
Posts: 6
Rep Power: 5
libi.macek is on a distinguished road
Thank you for advice
libi.macek is offline   Reply With Quote

Reply

Tags
laplacianfoam, simple algorithm

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pressure based and Density based Solver Xobile Main CFD Forum 39 August 19, 2020 07:04
[OpenFOAM] Take derivative of mean velocity in paraFoam hiuluom ParaView 13 April 26, 2016 07:44
Derivative of velocity and mean velocity hiuluom OpenFOAM Post-Processing 1 May 30, 2015 00:42
[Other] mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 05:16
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 15:11


All times are GMT -4. The time now is 06:57.