CFD Online Discussion Forums - melting problem: looking for appropriate solvers

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - melting problem: looking for appropriate solvers (https://www.cfd-online.com/Forums/openfoam-solving/93620-melting-problem-looking-appropriate-solvers.html)

Definition of lf

Quote:

Originally Posted by Pavithra (Post 741376)

Hello Everyone,

Finally, I have managed to modify buoyantBoussinesqPimpleFoam to simulate melting problems. I have attached the solver and Gallium melting testcase. The solver is tested in OF v6.

Any feedback or suggestions are welcome.

Also, please give me a direction or any guidance to parallelize this solver.

The solver works well when running on a single core. But, blows up while running in parallel.

Thank You.

- Pavithra.

Hi, Pavithra.. thank you for your solver. May i know what is that lf at the zero folder? my intuation says that it is alpha, the liquid fraction?

Quote:

Originally Posted by Pavithra (Post 741376)

Hi, it seems like the solution of your code is weird. I have seen a big vortex which is not appropriate. This one would be somewhat correct

https://www.cfd-online.com/Forums/op...tml#post586823

you can test it from this paper
https://www.researchgate.net/publica...th_experiments

Regards,
Mike

mesh dependent problem?

Quote:

Originally Posted by Ya_Squall2010 (Post 504403)

Thanks for reply. I am following, however, the parameters tabulated in table 1 of your <Heat Mass Transfer> paper on 2011. The only thing I have changed is that I tripled the mesh resolution used in your paper. Why different meshes are giving totally different solutions? Backed up by the experiment conducted by Gau and Viskanta, am I right to say that the converged result on the finer mesh was wrong? But how could it be?

Hi, i would like to ask if how was your study in melting of gallium? can you share your results/paper with me, if you have? your results seems strange to me too, but i would like to try it though..

Quote:

Originally Posted by fabian_roesler (Post 469557)

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

Hi Fabian,

It has been some time but did you make your melting/solidification solver with free surface interface public? I had also developed a melt pool model for arc welding application including metal transfer. However I used the error function to calculate the liquid fraction. When I compare my simulation results with the melt pool data from the experiment, it gives good results in terms of predicting the depth of the fusion penetration but has significant discrepancy when comparing the the free surface length of the melt pool. I am starting to think if it is because my liquid fraction is not converging properly.

Would it be possible to see your code specifically the part where liquid fraction and energy equation are iteratively solved.

I did infact implemented the iteratively solving of liquid fraction and energy equation, thanks to the solver uploaded in this discussion. However, it always return the following error when running in parallel mode. It has never happened before. Would you have idea what could be the reason?

Thanks !

Quote:

Fatal error in PMPI_Waitall: Message truncated, error stack:
PMPI_Waitall(405).....: MPI_Waitall(count=6, req_array=0x27479d0, status_array=0x1) failed
MPIR_Waitall_impl(181): fail failed
do_cts(649)...........: Message truncated; 67200 bytes received but buffer size is 22400
Fatal error in PMPI_Waitall: Message truncated, error stack:
PMPI_Waitall(405).................: MPI_Waitall(count=6, req_array=0x1ab5780, status_array=0x1) failed
MPIR_Waitall_impl(181)............: fail failed
MPIDI_CH3U_Receive_data_found(131): Message from rank 11 and tag 1 truncated; 48000 bytes received but buffer size is 16000
Fatal error in PMPI_Waitall: Message truncated, error stack:
PMPI_Waitall(405).................: MPI_Waitall(count=8, req_array=0x2da8670, status_array=0x1) failed
MPIR_Waitall_impl(181)............: fail failed
MPIDI_CH3U_Receive_data_found(131): Message from rank 11 and tag 1 truncated; 13440 bytes received but buffer size is 4480
mpprun info: Job terminated with error

droplet solidification

Hi，
I am doing some work on droplet solidification. I tried to add TEqn into interFoam and the solver can be compiled, but temperature is very strange in the result of my test case.
Did anyone find or build an appropriate solver?

Integrating TEqn to interFoam solver is well documented. You can find many documentation of it on the web. May be search with wolfdynamics team. In fact, the report presented in the link below also gives you a detail steps for implementaion. I suggest you try with first order scheme to run your test cases first.

http://www.tfd.chalmers.se/~hani/kur...inalReport.pdf

Quote:

Originally Posted by Madeinspace (Post 818044)

Thank you very much! I will modify my solver according the report, and then run the case with first order scheme.

DC casting - Pull velocity

Has anyone managed to apply a "pull velocity" for the solid in the meltingFoam solver? I am trying to simulate a direct chill (DC) casting process (https://www.substech.com/dokuwiki/do...ill_dc_casting) where the solidified metal is pulled out from the computational domain at a specified casting speed. Similar to what is implemented in Ansys Fluent (https://www.afs.enea.it/project/nept...th/node357.htm)

Any tips here?

Best regards

Quote:

Originally Posted by mikulo (Post 794122)

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

Can I also have a copy of the Thesis please?
All the best