[The Mad Scientist]

making appropriate required physical changes should do

[dp32]

Sir, can I know what are all the changes to be done?

[The Mad Scientist]

first check your solver for any test case.. if its working fine then all you need to do is supply physical BC & IC along with thermophyscial parameter..

[dp32]

Sir, i have checked the solver for melting problem and it's working fine. But, when it come to Solidification problem, i'm unable to capture the natural convection effects. so, could you please help me with the initial and boundary conditions.

[The Mad Scientist]

post your results and error ...

[riccardo93]

Hi everyone,
I need to simulate the solidification of steel.
At the moment I'm using the solver solidificationFoam, but there are some missing things:
-Communication between solid and liquid domain, ingot and ingot mold
-Thermal Shrinkage
-Viscosity variation as a function of temperature

What solver do you recommend?

Thanks in advance!

[Shoonya]

Originally Posted by Pavithra

Hi Everyone,

I am attaching the results of some validation cases, herewith.

I have attached the solver, cases and results. Any suggestions for improvements are welcome.

Thank You.

- Pavithra

Thanks Pavithra for sharing the solver meltingFoam. It works fine in the present version of OpenFOAM (v1912).

Regards

[wanghongjie]

Originally Posted by kuechenrole

Hey Fabian,

I'm working on a problem you might have already solved in your PhD and therefore I would like to ask you for a hint.

My Master thesis is about the simulation of casting processes including the filling, temperature transport and solidification and thus exactly the reversal of your subject.
In order to add the gas phase I'm trying to combine interFoam and convMeltFoam, as you suggested earlier.

Using the model of your PhD, everything seems to work out fine except the solution of the TEqn, which exploits whenever I include the calculation of the liquid fraction fl.

The only thing I changed so far is the update equation of the liquid fraction. I replaced the heat capacity cp of the hole cell with a cp calculated from just the liquid and solid phase.

I suggest that the new possibility of changing the fl by filling a cell instead of just melting solid is causing the trouble. By filling a cell I can produce huge differences of fl between two time steps, which the iteration algorithm for T and fl can't handle anymore.

Do you come to that problem?

My supervisor mentioned a smoothing function to restrict the time step if huge fl changes happen could be a solution. I haven't tried it yet, but will post again, if something works.

I have to excuse, if the solution is described in your PhD. We already ordered a book through inter-library loan, but it seems to be loaned. Do you know about a library online version? If I haven't got access by next week I will buy it. (Until then I'm dealing with the googleBooks prev)

Thanks,

ole

Can i have the paper？ Thanks！！！

[aerospace1365]

Originally Posted by som

I have a ball with fixed T gradient on the boundary. this ball is placed into tube full of wax. Your solver acts very well with such problem.

Now i have new problem. I want make the ball moving in the tube.
I've tried to setup the boundary of the ball to type movingWallVelocity and run the case. Everything is good but a cant see any boundary displacement in ParaFoam.
Any ideas?

Hello Som.
I want to cool a heat sink by butting it on the wax, just like your problem.
May I have your code?
would you please send it to my email?
aerospace1365@gmail.com

[aerospace1365]

Originally Posted by 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

Dear Sir
I have a case that should be solved by your solver. I have a heat sink (fins) located in wax. the wax will melt and cool the fins.
do you think this solver can solve this problem?
my Gmail is aerospace1365@gmail.com.
I am looking forward to hearing from you.
regards

[Skaiwalker]

Hello all,

I want to port the solidification solver to OF8 to utilize the tabulated thermophysicalProperties in my bachelor's thesis. "solidificaitionMeltingSource" looked like a good starting point to me, so I changed the "update" and "apply" functions according to the formulation in Fabian Röslers Solver and made some minor changes like deleting alpha1e and relax factor. The code is attached in the tar file.

---
Question 1: Is there missing a negative sign in the paper in equation 10? In the uploaded solver post 19 there is an additional negative sign in the upper numerator of TEqn, as noted in the equation snippet below. Implementation according to the paper throws a floating point exception error.


---

I ran some cases but the results are terrible as can be seen in the attached picture (as in the paper the timesteps 2, 6, 10 and 17 min are plotted). The melt front is way too fast and doesn't show curving like the results from the paper. Out of curiosity I tested the standard implemented linear fvOption with some cranked up numerics and ran the same settings again for the erf fvOption:

• linear fvOption differs from the results in the paper in melt front position and shape increasingly with time. Adjusting the model parameters helps a bit in later time steps with the shape but the initial deviation stays the same.
• erf fvOptoin with cranked up numerics moves the melt front a bit in the right direction, but is far from accurate.

---
Question 2: Does anybody have an idea where this behavior comes from? I even programmed it a second time from the scratch to avoid typos, but the results were the same. And this doesn't explain the behavior of the linear model as well.
---

For post processing of the melt front run "pvpython writeMeltFront.py" in the case directory and point plotMeltFront.py to meltFrontPosition directory in corresponding case directory. I left the results in the myCases cases for an easy first impression.

Thanks in advance and have a nice day!

[Zonk90]

Hello Kai,


Regarding your first question: the derivative of is , so I think yes there should be a minus in equation 10.

And regarding your second question. The erf-Solver is rather mediocre. Even if everything is implemented correctly you will need a very small time step, since the very strong non-linearity of the equation is not handled properly. You can take a look at: https://www.mdpi.com/1996-1073/12/5/868 for a better implementation.


Best,
Moritz

[Skaiwalker]

Hello again forum!

Moritz posted an interesting approach, I will look further into it. But in the mean time, does anyone of you have an Idea where especially the deviation of the linear model comes from?

Thanks in advance!

[fabian_roesler]

Hi,
This is my first post since years. The erf function was the first more or less running version we developed at the thermodynamics chair in Bayreuth. The more sophisticated version is in my PhD thesis. The latest version seems to be the one in Moritz's article. I had no time to reed through the article yet and I'm out of this topic since 7 years .
What deviation do you mean?
Cheers
Fabian

[Skaiwalker]

Hi Fabian,

thanks for your reply, I received your PhD thesis yesterday and it looks very helpful from what I've read until now! Especially the parts about sinking solid and literature review.

The deviation I'm talking about are the more narrow expansion in x of the linear model compared to yours and Brent and Voller's simulations and what is the most significant difference in my opinion is the shape of the phase front. The more "wave-like" character up until 360 s and the pronounced bulge from then on.

But I missed a paper Moritz sent me in an E-Mail which explains these results, shame on me Hannoun et al. published "Resolving the controversy over tin and gallium melting in a rectangular cavity heated from the side" and state, that Brent and Voller's results are not grid-converged. So the positions of vortexes are wrong and in the early stages of the simulation, when vortexes are small, they are not resolved at all.
Further the experiment of Gau and Viskanta drained the fluid portion to measure the phase position, so we have different bc's in the experiment and in the simulation.

I will look into your thesis for the theory behind phase change but most likely adopt Moritz's approach because it is the most advanced and computational cost may be reduced.

[Vincy_]

Hi:

Has anyone encountered the following problems？
I guess the boundary conditions are set incorrectly? But I have tried many, still have this problem.



--> FOAM FATAL ERROR: (openfoam-2012)
Attempt to cast type calculated to type nutWallFunction at index 2

From To& Foam::refCast(From&, Foam::label) [with To = const Foam::nutWallFunctionFvPatchScalarField; From = const Foam::fvPatchField<double>; Foam::label = int]
in file /home/user2/OpenFOAM/OpenFOAM-v2012/src/OpenFOAM/lnInclude/typeInfo.H at line 162.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::exitOrAbort(int, bool) at ??:?
#2 Foam::nutWallFunctionFvPatchScalarField::nutw(Foam ::turbulenceModel const&, int) at ??:?
#3 Foam:megaWallFunctionFvPatchScalarField::calcula te(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#4 Foam:megaWallFunctionFvPatchScalarField::calcula teTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#5 Foam:megaWallFunctionFvPatchScalarField::updateC oeffs() at ??:?
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::updateCoeffs() at ??:?
#7 Foam::kOmegaSSTBase<Foam::LESModels:ESModel<Foam ::IncompressibleTurbulenceModel<Foam::transportMod el> > >::correct() at ??:?
#8 ? at ??:?
#9 __libc_start_main in /lib64/libc.so.6
#10 ? at ??:?
Aborted (core dumped)



Thanks in advance

[mikulo]

Originally Posted by fabian_roesler

Hi steph

I did exactly that within by PhD thesis. I defended the thesis one month ago and it will be published soon. Unfortunately it is written in German. I combined the compressibleInterFoam solver with my new melting solver based on the convMeltFoam. The density of the PCM is no longer described by means of a boussinesq approximation but by a temperature dependent density in all terms of the conservation equations. This makes it possible to simulate volume change during melting and solidification in a closed volume/capsule.
As soon as my thesis is published, I will post at least parts of the solver here.

Cheers

Fabian

Hello Fabian, is it okay to have a copy of your thesis. Would it be better if it is written in English but if not, then I will try to find a way to translate it. Thanks

[Vincy_]

Originally Posted by Vincy_

Hi:

Has anyone encountered the following problems？
I guess the boundary conditions are set incorrectly? But I have tried many, still have this problem.



--> FOAM FATAL ERROR: (openfoam-2012)
Attempt to cast type calculated to type nutWallFunction at index 2

From To& Foam::refCast(From&, Foam::label) [with To = const Foam::nutWallFunctionFvPatchScalarField; From = const Foam::fvPatchField<double>; Foam::label = int]
in file /home/user2/OpenFOAM/OpenFOAM-v2012/src/OpenFOAM/lnInclude/typeInfo.H at line 162.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::exitOrAbort(int, bool) at ??:?
#2 Foam::nutWallFunctionFvPatchScalarField::nutw(Foam ::turbulenceModel const&, int) at ??:?
#3 Foam:megaWallFunctionFvPatchScalarField::calcula te(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#4 Foam:megaWallFunctionFvPatchScalarField::calcula teTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#5 Foam:megaWallFunctionFvPatchScalarField::updateC oeffs() at ??:?
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::updateCoeffs() at ??:?
#7 Foam::kOmegaSSTBase<Foam::LESModels:ESModel<Foam ::IncompressibleTurbulenceModel<Foam::transportMod el> > >::correct() at ??:?
#8 ? at ??:?
#9 __libc_start_main in /lib64/libc.so.6
#10 ? at ??:?
Aborted (core dumped)



Thanks in advance

I have solved this problem, my boundary conditions are set incorrectly, mainly because I have to be careful about the setting of inlet, outlet and wall.

[zhaohb11_cfd]

Hi everyone
I want to cpuple interfoam and the solidification model to simulate the droplet freezing. Thus I put the T equation in the interFoam solver and compiled. The solver works fine in serial mode. But when I run in parallel. There is error:
//////////////////////////////////////////////////////////////////////////////
Fatal error in MPI_Recv: Message truncated, error stack:
MPI_Recv(224).....................: MPI_Recv(buf=0x1da60c0, count=6400, MPI_BYTE, src=8, tag=1, comm=0x84000002, status=0x7ffec83b2e40) failed
MPIDI_CH3U_Receive_data_found(131): Message from rank 8 and tag 1 truncated; 19200 bytes received but buffer size is 6400
Fatal error in MPI_Recv: Message truncated, error stack:
MPI_Recv(224).....................: MPI_Recv(buf=0x18748a0, count=6400, MPI_BYTE, src=10, tag=1, comm=0x84000002, status=0x7ffe00c99a10) failed
MPIDI_CH3U_Receive_data_found(131): Message from rank 10 and tag 1 truncated; 19200 bytes received but buffer size is 6400
Fatal error in MPI_Recv: Message truncated, error stack:
MPI_Recv(224).....................: MPI_Recv(buf=0x236adb0, count=6400, MPI_BYTE, src=11, tag=1, comm=0x84000002, status=0x7ffd32387010) failed
MPIDI_CH3U_Receive_data_found(131): Message from rank 11 and tag 1 truncated; 19200 bytes received but buffer size is 6400
Fatal error in MPI_Recv: Message truncated, error stack:
MPI_Recv(224).....................: MPI_Recv(buf=0x2a4f080, count=6400, MPI_BYTE, src=9, tag=1, comm=0x84000002, status=0x7ffd606c04e0) failed
MPIDI_CH3U_Receive_data_found(131): Message from rank 9 and tag 1 truncated; 19200 bytes received but buffer size is 6400
//////////////////////////////////////////////////////////////////////////////////
the attachmen is the T Equation file. Can someone know the reason why this error happens

[zhaohb11_cfd]

Hello, everyone, have you succeeded in running cases in parallel when using the fabian_roesler code. I have followed fabian_roesler's code and it is successful in serial running but failed in paralleling running. I do not know the reason. Have someone met and fixed this problem?