interFoam, setting velocity in fluid1 to zero

Martin_K_lalelu · January 19, 2015, 11:36

Hi!

To compare my VOF results of a deforming water droplet with a rigid droplet, I would like to keep using interfoam. Is there a way to set the velocity of one phase to zero, like in

Code:

forAll(mesh.cells(),celli)
{
    if(alpha1[celli]>0.0)
    {
         U[celli]==vector::zero;
    }
}

?

Unfortunately, this does not work, neither does an artificially high viscosity to simulate a rigid body.

There is a high chance that this is my fault, so any help would be really appreciated!

ThomasV · March 5, 2015, 05:17

Hi!

I'm having a similar problem: I added a temperature field to the solver and want to change the viscosity of the "water" phase depending on the temperature in each cell. I thus need access to the viscosity field...

What I figured out until now is this:
"twoPhaseProperties.nu()" returns the final viscosity field for each cell...

I didn't manage to get access to the initial viscosity fields though. All of this "unfortunately" is embedded in the incompressibleTwoPhaseMixture model. I really don't know much of this entire system that's being used for handling the phases so is there a way to get access to the respecitve fields and change them for each cell? Or with that said do functions which give access to those fields exist? I think the fields were created as protected objects inside the respective class so I won't be able to simply edit them in my solver routine...

It would be nice if someone who's a bit more familiar with this phase model could help us out here...

January 19, 2015, 11:36	interFoam, setting velocity in fluid1 to zero	#1
Martin_K_lalelu New Member Martin K Join Date: Jan 2013 Location: Germany Posts: 28 Rep Power: 13	Hi! To compare my VOF results of a deforming water droplet with a rigid droplet, I would like to keep using interfoam. Is there a way to set the velocity of one phase to zero, like in Code: forAll(mesh.cells(),celli) { if(alpha1[celli]>0.0) { U[celli]==vector::zero; } } ? Unfortunately, this does not work, neither does an artificially high viscosity to simulate a rigid body. There is a high chance that this is my fault, so any help would be really appreciated!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
interFoam vs. simpleFoam channel flow comparison	DanM	OpenFOAM Running, Solving & CFD	12	January 31, 2020 15:26
Multiple floating objects	CKH	OpenFOAM Running, Solving & CFD	14	February 20, 2019 09:08
A Question About Setting Whole Velocity Field Using "Proflie" in FLUENT	adsl17754	FLUENT	5	July 31, 2018 04:29
Setting velocity at outlet	saba_saeb	FLUENT	29	December 12, 2012 05:16
Using UDF for setting velocity for two fluids	UDF user	FLUENT	7	November 26, 2008 07:22

March 5, 2015, 05:17		#2
ThomasV Member Thomas Vossel Join Date: Aug 2013 Location: Germany Posts: 45 Rep Power: 12	Hi! I'm having a similar problem: I added a temperature field to the solver and want to change the viscosity of the "water" phase depending on the temperature in each cell. I thus need access to the viscosity field... What I figured out until now is this: "twoPhaseProperties.nu()" returns the final viscosity field for each cell... I didn't manage to get access to the initial viscosity fields though. All of this "unfortunately" is embedded in the incompressibleTwoPhaseMixture model. I really don't know much of this entire system that's being used for handling the phases so is there a way to get access to the respecitve fields and change them for each cell? Or with that said do functions which give access to those fields exist? I think the fields were created as protected objects inside the respective class so I won't be able to simply edit them in my solver routine... It would be nice if someone who's a bit more familiar with this phase model could help us out here...