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

Marangoni bc

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

Like Tree6Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 19, 2017, 03:21
Default
  #21
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,379
Rep Power: 26
akidess will become famous soon enough
Quote:
Originally Posted by gtarang View Post
Hello Thomas,
I am implementing the same boundary condition. However I want to ask, do we really need to update the solver application? There must be another way of adding this IOobject externally I mean through codeStream or function object. Pl. suggest.
-
Tarang
That's right, you can replicate the functionality using a groovyBC boundary condition. Probably also with codeStream, but I haven't used that yet.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   March 1, 2019, 04:31
Default
  #22
New Member
 
Navid
Join Date: Nov 2015
Posts: 17
Rep Power: 6
navidamin is on a distinguished road
Hey guys!
the links to the files that tehache posted are apparently corrupt now. Can any one please upload the files and post the link here?

Thanks



Quote:
Originally Posted by tehache View Post
I updated the boundary condition a little (its now writing also its value, was missing before..), and its working for me now with 1.6.x.
As the testcase was too large to upload here, please find it (and the B.C.) at the links given below.
Under the link named MarangoniTest I add a small testcase. Gravity is set to 0, Marangoni-Number is small in order to obtain a stationary and axisymmetric solution (although calculation is 3D), which I have compared to results from our own 2D-Code Crysmas, and with Fluent, results are (within some percent, did not check grid dependence..) the same.
Therefore, I personally trust it, but dont give any guarantee

I have run the testcase using a slightly modified buoyantBoussinesqSimpleFoam:

add this to createFields to store temperature gradients:

Info<< "Generating/reading field gradT\n" << endl;
volVectorField gradT
(
IOobject
(
"gradT",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
mesh,
dimensionSet(0,-1,0,1,0,0,0)
);
gradT=fvc::grad(T);


and add this in the solver code:

turbulence->correct();
-> gradT=fvc::grad(T);
runTime.write();


specify the B.C. in U as follows:
outerwall
{
type marangoni;
marangonicoeff -9.9e-06;
}

the marangonicoeff to specifiy is the usual marangonicoefficient (d sigma/dT) divided by DYNAMIC viscosity!




http://servww6.ww.uni-erlangen.de/~jung/marangoniFvPatchVectorField-1.6.x.tar.gz
http://servww6.ww.uni-erlangen.de/~jung/MarangoniTest.tar.gz


Thomas
navidamin is offline   Reply With Quote

Old   March 1, 2019, 05:13
Default Marangoni Condition
  #23
Member
 
Tarang
Join Date: Feb 2011
Location: Delhi, India
Posts: 42
Rep Power: 11
gtarang is on a distinguished road
Here, use this. I am also right now seeing the feasibility if gradT can be calculated inside boundary patch only. Let's see when I am able to do it. For the time being you calculate gradT as mentioned in the post. Further you can use
Code:
IOobject::NO_READ,
IOobject::NO_WRITE
in place of
Code:
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
There is no need to read or write the gradT.


PS: The file will work for Openfoam-4.1 onwards

-
Tarang
Attached Files
File Type: gz marangoniGradientCondition.tar.gz (3.0 KB, 20 views)
navidamin and fjgwd002 like this.
gtarang is offline   Reply With Quote

Old   March 1, 2019, 05:42
Default
  #24
New Member
 
Navid
Join Date: Nov 2015
Posts: 17
Rep Power: 6
navidamin is on a distinguished road
Thanks Tarang!
could you please also upload your sample case?
navidamin is offline   Reply With Quote

Old   November 1, 2019, 00:49
Default
  #25
New Member
 
Nikhil Chitnavis
Join Date: Sep 2019
Location: India
Posts: 7
Rep Power: 2
Nikhil2511 is on a distinguished road
hello,

I am using the marangoni Bc for the 2D square cavity with temperature and concentration given to side wall but the results i am getting is not correct, the flow is not developing. I am using Marangoni BC with simpleFoam.
please suggest something on this.

thanks
Attached Images
File Type: png streamlines1.png (35.0 KB, 12 views)
Nikhil2511 is offline   Reply With Quote

Old   July 12, 2020, 14:59
Default Marangoni Boundary condition
  #26
Member
 
Tarang
Join Date: Feb 2011
Location: Delhi, India
Posts: 42
Rep Power: 11
gtarang is on a distinguished road
I have updated the marangoni condition to work with both compressible as well as incompressible solvers. The boundary condition calculates gradient of temperature internally. I have tested it in my problems and it is working well. Let me know if you face any problem.


You can also check out my github page
https://github.com/gargtarang/OpenFOAM_Extensions_TG


This I maintain mainly as an hobby or interest.


--

Thanks
gtarang
Attached Files
File Type: zip marangoniGradientCondition.zip (6.1 KB, 0 views)
gtarang is offline   Reply With Quote

Reply

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
Marangoni flow Tao CFX 0 September 4, 2007 17:14
Marangoni stress, help!! dong FLUENT 1 January 25, 2007 11:54
Marangoni Convection Rucy FLUENT 0 July 6, 2005 05:17
Marangoni Stress At the Wall Peter Nagy FLUENT 4 January 26, 2004 15:06
Marangoni VS Rayleigh-Bernard Peter Nagy FLUENT 2 January 12, 2004 09:53


All times are GMT -4. The time now is 15:37.