[Typhain]

Hello and thanks for the answer !
Yeah I know, there is lots of step between 2.1 and 5.0 but it work for meltFoam solver post at the bigining of the conversation...
Well I wil follow your advice and climb step by step. If it work I will share my solver for 5.0 version

Bye

[cfd@kgp]

Typhain, are you saying meltfoam got compiled in OF5.0?

if so, pls send the link of the solver

Thanks,

[Typhain]

Here is the meltFoam solver for 5.0 OpenFoam.
And there is also the test case with gallium wich use meltFoam solver.
mySolvers.tar.gz

[kera]

Hallo,

If I am not mistaken, from OpenFOAM 4.0 in order to access the internal values it make uses of "primitiveField()" instead of "internalField()" (as in the earlier versions) and hence the error with gMax.

try using this:

residual = gMax(mag(alpha3New.primitiveField()-alpha3.primitiveField()));
meanResidual = gSum(mag(alpha3New.primitiveField()-alpha3.primitiveField())*mesh.V())/gSum(mesh.V().field());

hope this helps!

Regards,
Ricky

Originally Posted by Typhain

residual = gMax(mag(alpha3New.internalField()-alpha3.internalField()));
meanResidual = gSum(mag(alpha3New.internalField()-alpha3.internalField())*mesh.V())/gSum(mesh.V().field());

And the error message is the folowing

TEqn.H:37:78: error: no matching function for call to ‘gMax(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> >)’

So I am a little bit confuse about gMax function, what should be the dimension of alpha3 to use it ?

thanks

[Typhain]

Hello Ricky,
First of all I'd like to thank you for the answer, and then thank you for the help. Using primitiveField() work as expected and thanks to you the compilation of convMelt_parallel_solver work.

There is also details on simulation that I have to correct but I think soon, I will post the solver uptated for OF5.0 and a testCase.

My mission now is to change the entry of the simulation. In deed, I want to change the nature of the boundaries conditions : instead of having temperature, I would like to work with thermal flux.
Maybe you have an idea, of how to proceed ?

Regards
Typhain

[cfd@kgp]

Dear Typhain,

There are two ways to proceed for implementing heat flux boundary condition namely,

1. programming method, edit the current boundary condition file located in src directory --> (Dirichlet boundary condition) to Neumann boundary condition or Robin based on your need. This methods involves compilation and then adding the name of.lib in the controlDict

2. Second method is slightly simpler, by using the groovy boundary conditions (check the implementation details here https://openfoamwiki.net/index.php/Contrib/groovyBC).
here also u have to add "libs ( "libgroovyBC.so" );" in controlDict

3. I am not very sure whether swak4foam can also help u in doing this.

But I would go for option 2. Happy Foaming!

Regards,
cfd@kgp

[dshafer]

Fabian, would you be willing to share the solver you mention in your thesis which takes unconstrained close contact melting into account?

I realize your thesis came out some time ago - perhaps some of the other solvers mentioned on this forum have been updated to model unconstrained melting? If not, I would be very grateful for suggestions for how to implement this myself. At the moment I'm simply hoping to replicate a simple test case of Octadecane melting in a spherical enclosure.

Full disclosure: I am new to OpenFOAM and have only attempted to model this problem in COMSOL and Fluent so far, without success. I expect the solution will require me to modify the momentum equation to A*(u - u_sink), where A is the porosity function and u_sink is the sinking velocity of the solid. I am not sure how to define u_sink in this case.

I'm hoping OpenFOAM will be more suitable to this problem!

Many thanks

Originally Posted by fabian_roesler

Hi Folks

My PhD thesis on modeling and simulation of phase change processes in macro encapsulated latent heat thermal energy storage is now published. It is written in German and was published by the Logos Verlag Berlin.

Modellierung und Simulation der Phasenwechselvorgänge in makroverkapselten latenten thermischen Speichern
Thermodynamik: Energie - Umwelt - Technik, Bd. 24
Fabian Rösler
ISBN 978-3-8325-3787-6

http://www.logos-verlag.de/cgi-bin/e...87&lng=deu&id=

The convMeltFoam solver as well as the solver which takes unconstrained close contact melting and an additional gas phase into account are explained in detail. Experimental validation in an rectangular cavity is performed.
Hope you like it.

Cheers

Fabian

[dshafer]

Dear all,

Does anyone know if the solvers discussed in this thread can be used to model unconstrained close-contact melting? I am trying to replicate a simple case study of Octadecane melting in a spherical enclosure, but am not aware of any solvers which can take the sinking of the solid into account.

I should note that I am new to OpenFOAM, and have only attempted to solve this problem in COMSOL and Fluent so far (without success).

I'd love to know if anyone has had success using OpenFOAM for these purposes!

Thank you!

[ahmed1981]

Hi Typhain
is there OpenFoam version 5.0 ready for download, I search for this version in internet but I do not find it.
I now use openFoam 4.1, I ask can I use meltFoam solver for this version, If your answer is yes. How can I use this solver please .

I ask that because I am new to OpenFoam, And I want to model phase change materials
.
Thank you very much


Ahmed Almudhafar

[janghel]

Dear All,
Is there any one who has successfully modelled solidification with shrinkage. I want to do the shrinkage model in openfoam. But i dont know from where should I start.

[janghel]

Originally Posted by fabian_roesler

Hi callahance

I did exactly what you propose in my thesis. However, I am still in the PhD process and want the solver, results and thesis be public not before my defense of the doctor's thesis. You are right; I combined the compressibleInterFoam solver with my own melting solver to allow a free surface between liquid and gas.

Regards

Fabian

Hello Fabian,
I am trying to solve the three phase problem with your convFinMeltFoam solver. I have taken Fin material as air (third phase) with solid and liquid. I want to track the interface of liquid and air using VOF as solidification takes place. Due to density difference in solid and liquid, volume will shrink. I am unable to combine the solver with interFoam.
Please give some advice...

[kera]

Hallo diwaker,

I am also trying to implement shrinkage model due to density variations, were you able to find some good literature?

I am currently trying to adapt the model described in [1] to my solver:

[1] "M. Raessi and J. Mostaghimi, Three-dimensional modeling of density variation due to phase change in complex free surface flows, Numerical Heat transfer, Part B, 47: 507-531, 2005"

Thank you in advance.

Regards,
Ricky

Originally Posted by janghel

Dear All,
Is there any one who has successfully modelled solidification with shrinkage. I want to do the shrinkage model in openfoam. But i dont know from where should I start.

[janghel]

Originally Posted by kera

Hallo diwaker,

I am also trying to implement shrinkage model due to density variations, were you able to find some good literature?

I am currently trying to adapt the model described in [1] to my solver:

[1] "M. Raessi and J. Mostaghimi, Three-dimensional modeling of density variation due to phase change in complex free surface flows, Numerical Heat transfer, Part B, 47: 507-531, 2005"

Thank you in advance.

Regards,
Ricky

Dear kera,
you can look into paper by Sun and Garimella, they have done solidification shrinkage of TNT.
"Numerical and Experimental Investigation of Solidification Shrinkage" 2007
i have some more papers but they are not like Raessi's paper. Using raessi's model you can get 3D cavity.
I am also trying to implement the model by Sun and Garimella. but no progress so far..
were you able to implement shrinkage?

ThankYou

[kera]

Hallo diwaker and FOAMers,

I tried to implement the source term in "alphaEqn" as per this forum InterFoam: add a source term in alpha eq. and it seems that the source term is being calculated (checked via some print statements), but the problem is I don't see any effect in the simulation. --> This time I blindly followed Sun and Garimella [1], maybe I am doing something wrong.

Here is a snippet of the code.

volScalarField Sp 
    (
        IOobject
        (
            "Sp", 
         runTime.timeName(),
         mesh
         ),
      ((rhol - rhos)*gammas*fvc::ddt(alpha1)/rhol  - gammas*fvc::div(phiAlpha1))
     );

explicitSolve:

MULES::explicitSolve(geometricOneField(), alpha1, phi, tphiAlpha(), Sp,Su, 1, 0);

I am not sure where exactly is the problem, any help is really appreciated.

Thank you!

Regards,
Ricky

[1] " D. Sun and S.V. Garimella, Numerical and Experimental Investigation of Solidification Shrinkage, Prude University, 2007"

Originally Posted by janghel

Dear kera,
you can look into paper by Sun and Garimella, they have done solidification shrinkage of TNT.
"Numerical and Experimental Investigation of Solidification Shrinkage" 2007
i have some more papers but they are not like Raessi's paper. Using raessi's model you can get 3D cavity.
I am also trying to implement the model by Sun and Garimella. but no progress so far..
were you able to implement shrinkage?

ThankYou

[janghel]

Dear ricky,
You are One step ahead of me. I am actually working on solver of Fabian with little modification. I combined the solver with interFoam. Right now i am working on it. will update you ASAP. Meanwhile i tried to combine the solidification solver with compressibleinterFoam. But it did not work. I got an error. I will try your approach by adding source term to alphaEqn. One Question?
Did you solve continuity equation with source terms mentioned in both the papers...?

Regards
Diwakar janghel

[MSF]

Hi

Dear all,

Does anyone know if the solvers discussed in this thread can be used to model unconstrained close-contact melting? I am trying to replicate a simple case study of Octadecane melting in a spherical enclosure, but am not aware of any solvers which can take the sinking of the solid into account.

I should note that I am new to OpenFOAM, and have only attempted to solve this problem in COMSOL and Fluent so far (without success).

I'd love to know if anyone has had success using OpenFOAM for these purposes!

Thank you!

If you want to implement close-contact melting you can use multiple approaches: Variable viscosity, force balance + Darcy term or prescribe the settling velocity with forcing functions. Here are some links to papers. The first two models are in OpenFoam, the third was implemented in Matlab I think.

http://www.sciencedirect.com/science...17931017329241
Darcy-Term + force balance
https://link.springer.com/article/10.1007/s00231-016-1932-0
Variable viscosity
http://www.sciencedirect.com/science...17931017300856
Forcing functions + force balance


Best
Moritz

[Shekhar@03]

Originally Posted by AnjaMiehe

Hello David,

this I not trivial for me either as the solver was coded in 1.4.1, or at least that is the tag in the test case. Also, the test case proposes icoFoam, a PISO based solver, to solve the test case. I tried to rewrite meltFoam in 2.1.0 using the SIMPLE Algorithm but I found far more similarities in coding with the icoFoam solver in 2.1.0.

Therefore, and as my try for the meltSimpleFoam does compile but not solve the test case, here is my idea for the meltIcoFoam in 2.1.0 . For me, it compiles only with warnings and the rough test case given earlier in this thread works alright. I can't guarantee for anything.


I hope this is a help for you,
Regards, Anja

can you please provide a test case for this solver?

[Germilly]

Hello ahmmedshakil,

I am dealing with the same problem (post: #100), have you solved it?

My energy equation is:

   fvm::div(phiCpf, Tf) - fvm::laplacian(por*kf, Tf)
 - hconv*Ts + fvm::Sp(hconv, Tf)

where phiCpf is:

phiCpf = fvc::interpolate(cpf)*phi

Tf is the temperature of the fluid...cpf is function of temperature using a polynomial equation.

The solver is running, but the final result does not respect the conservation of energy.

I think I have to do some correction, because of fvm::div(phiCpf, Tf) ...

Can you help me?

Thank you

GB

[Parsa Nazmi]

Hi Everyone

I'm new to OpenFoam6 and I want to melt a two dimensional square ice with conduction. I don't know how to start. Can you help me about solver?

[hemanth17]

hey i m using this sovler and i m facing an error(heat generation at bottom left corner) intial condition upper wall at higher temperature, lower wall at lower temperature, left and right wall at zero gradient condition