Sinusoidal initialisation of initial holdup field.

balkrishna · December 13, 2011, 09:29

Hi,

I am using VOF to simulate a 2d case of bubble rising in a liquid. I wish to initialise the gas holdup field in a sinusoidal manner. Precisely ,

Code:

If y < Y 
alpha =1 

Else 
alpha =0

Is there any way to do this using setFields ?

Thanks,
Balkrishna.

Bernhard · December 13, 2011, 11:54

Using setFields you have the boxToCell method for initialization. Did you try this one?

balkrishna · December 13, 2011, 14:20

Yes , but i think that allows one to give only initialisations in the form of quadrilaterals ... thats why the name BOXtoCell . Pls correct me if I am wrong ...

Bernhard · December 13, 2011, 16:16

Ok, I misread your pseudocode, but I am quite sure you can do this with funkySetFields of the swak4Foam utility: http://openfoamwiki.net/index.php/Contrib/swak4Foam

ngj · December 13, 2011, 16:51

Hi Balkrishna

You would also be able to do it with the waves2Foam utility found in

http://openfoamwiki.net/index.php/Contrib/waves2Foam

If you need a sinusoidal interface with a quiescent fluid, then you could initialise the computational domain with a 'stokesFirstStanding' type wave. However, it might be more easily merely to use swak4Foam.

Kind regards,

Niels

December 13, 2011, 09:29	Sinusoidal initialisation of initial holdup field.	#1
balkrishna Senior Member Balkrishna Patankar Join Date: Mar 2009 Location: Pune Posts: 123 Rep Power: 17	Hi, I am using VOF to simulate a 2d case of bubble rising in a liquid. I wish to initialise the gas holdup field in a sinusoidal manner. Precisely , Code: If y < Y alpha =1 Else alpha =0 Is there any way to do this using setFields ? Thanks, Balkrishna.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Extrusion with OpenFoam problem No. Iterations 0	Lord Kelvin	OpenFOAM Running, Solving & CFD	8	March 28, 2016 12:08
Moving mesh	Niklas Wikstrom (Wikstrom)	OpenFOAM Running, Solving & CFD	122	June 15, 2014 07:20
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3	bookie56	OpenFOAM Installation	8	August 13, 2011 05:03
Orifice Plate with a fully developed flow - Problems with convergence	jonmec	OpenFOAM Running, Solving & CFD	3	July 28, 2011 06:24
Differences between serial and parallel runs	carsten	OpenFOAM Bugs	11	September 12, 2008 12:16

December 13, 2011, 11:54		#2
Bernhard Senior Member Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 22	Using setFields you have the boxToCell method for initialization. Did you try this one?

December 13, 2011, 14:20		#3
balkrishna Senior Member Balkrishna Patankar Join Date: Mar 2009 Location: Pune Posts: 123 Rep Power: 17	Yes , but i think that allows one to give only initialisations in the form of quadrilaterals ... thats why the name BOXtoCell . Pls correct me if I am wrong ...

December 13, 2011, 16:16		#4
Bernhard Senior Member Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 22	Ok, I misread your pseudocode, but I am quite sure you can do this with funkySetFields of the swak4Foam utility: http://openfoamwiki.net/index.php/Contrib/swak4Foam

December 13, 2011, 16:51		#5
ngj Senior Member Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Copenhagen, Denmark Posts: 1,903 Rep Power: 37	Hi Balkrishna You would also be able to do it with the waves2Foam utility found in http://openfoamwiki.net/index.php/Contrib/waves2Foam If you need a sinusoidal interface with a quiescent fluid, then you could initialise the computational domain with a 'stokesFirstStanding' type wave. However, it might be more easily merely to use swak4Foam. Kind regards, Niels