Marangoni Convection - icoReactingMultiphaseInterFoam

Schmidtgen · July 16, 2021, 11:04

Hello,

i am trying to implement the temperature dependent surface tension (Marangoni Convection). Too validate my Implementation I want to replicate the Case published by Tan. But my Results aren't satisfying and maybe someone can have a look into my case setup.

The tangential surface tension component is implementent in UEqn:

$\vec{F}_{s}=[\sigma\kappa\vec{n}+\frac{d\sigma}{dT}(\nabla T - \vec{n}(\vec{n}\cdot \nabla T))]|\nabla \gamma|\frac{2\rho}{\rho_{1}+\rho_{2}}$

Code:

volScalarField alpha0 = fvc::average(fluid.phase(0));  
volScalarField alpha1 = fvc::average(fluid.phase(1));   
volVectorField gradAlpha = alpha1*fvc::grad(alpha0)-alpha0*fvc::grad(alpha1);
volVectorField nn = gradAlpha/(mag(gradAlpha)+deltaN);
volScalarField factor = (2*rho)/(fluid.phase(0).rho()+fluid.phase(1).rho());

fvVectorMatrix UEqn
(
    fvm::ddt(rho, U)
  + fvm::div(rhoPhi, U)
//- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
  + turbulence->divDevRhoReff(U)
  - dSigmadT*(fvc::grad(T) - nn*(nn & fvc::grad(T)))*mag(gradAlpha)*factor
  ==
    fvOptions(rho, U)
);

As normal suface tension component the constant sufaceTension Model is selected.

I attached my case setup. I look forward to any advise.

https://imgur.com/a/kD6QsQF

L. Tan, S. Leong, E. Leonardi, and T. Barber, “A numerical study of solid-liquid
phase change with marangoni effects using a multiphase approach,” Progress in
Computational Fluid Dynamics, vol. 6, no. 6, pp. 304–313, 2006.

Schmidtgen · July 17, 2021, 13:49

Here is my case.

XinXin · September 26, 2021, 03:21

hello,Schmidtgen, have you solved this problem?Could you please share your method ?

VINIT · October 20, 2022, 20:17

Quote:

Originally Posted by XinXin

hello,Schmidtgen, have you solved this problem?Could you please share your method ?

Hello XInXin

Were you able to incorporate marangoni?

Schmidtgen · October 26, 2022, 04:43

Quote:

Originally Posted by XinXin

hello,Schmidtgen, have you solved this problem?Could you please share your method ?

Sorry for the late reply. I had stopped working on the solver. And the problem was not finally solved by me.

qingqingliu · December 1, 2022, 12:42

I am also working on this solver, but different topic.
You introduced a void fraction in the U equation. Did you know how the void fraction is solved in this solver? I looked in the folder of this solver but did not find the void fraction transport equation.

July 16, 2021, 11:04	Marangoni Convection - icoReactingMultiphaseInterFoam	#1
Schmidtgen New Member Join Date: Jan 2021 Posts: 12 Rep Power: 5	Hello, i am trying to implement the temperature dependent surface tension (Marangoni Convection). Too validate my Implementation I want to replicate the Case published by Tan. But my Results aren't satisfying and maybe someone can have a look into my case setup. The tangential surface tension component is implementent in UEqn: $\vec{F}_{s}=[\sigma\kappa\vec{n}+\frac{d\sigma}{dT}(\nabla T - \vec{n}(\vec{n}\cdot \nabla T))]\|\nabla \gamma\|\frac{2\rho}{\rho_{1}+\rho_{2}}$ Code: volScalarField alpha0 = fvc::average(fluid.phase(0)); volScalarField alpha1 = fvc::average(fluid.phase(1)); volVectorField gradAlpha = alpha1fvc::grad(alpha0)-alpha0fvc::grad(alpha1); volVectorField nn = gradAlpha/(mag(gradAlpha)+deltaN); volScalarField factor = (2rho)/(fluid.phase(0).rho()+fluid.phase(1).rho()); fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) //- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U) + turbulence->divDevRhoReff(U) - dSigmadT(fvc::grad(T) - nn(nn & fvc::grad(T)))mag(gradAlpha)*factor == fvOptions(rho, U) ); As normal suface tension component the constant sufaceTension Model is selected. I attached my case setup. I look forward to any advise. https://imgur.com/a/kD6QsQF L. Tan, S. Leong, E. Leonardi, and T. Barber, “A numerical study of solid-liquid phase change with marangoni effects using a multiphase approach,” Progress in Computational Fluid Dynamics, vol. 6, no. 6, pp. 304–313, 2006.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Marangoni Convection	CneSpaulding	STAR-CCM+	2	June 6, 2014 07:38
Thermophysical properties for natural convection	Ciefdi	OpenFOAM Running, Solving & CFD	0	November 7, 2013 11:44
Modeling both radiation and convection on surfaces - Ansys Transient Thermal R13	s.mishra	ANSYS	0	March 31, 2012 04:12
Marangoni convection in CFX	shahpar73	CFX	3	January 31, 2011 06:13
Marangoni Convection	Rucy	FLUENT	0	July 6, 2005 05:17

September 26, 2021, 03:21		#3
XinXin New Member XinZhou Join Date: Apr 2018 Posts: 10 Rep Power: 8	hello,Schmidtgen, have you solved this problem?Could you please share your method ?

December 1, 2022, 12:42		#6
qingqingliu Member Ching Liu Join Date: Sep 2017 Posts: 45 Rep Power: 8	I am also working on this solver, but different topic. You introduced a void fraction in the U equation. Did you know how the void fraction is solved in this solver? I looked in the folder of this solver but did not find the void fraction transport equation.