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

SetfieldsDict in the multiphaseInterFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 18, 2013, 02:54
Post SetfieldsDict in the multiphaseInterFoam
  #1
New Member
 
salame ama
Join Date: Dec 2012
Posts: 28
Rep Power: 4
salame is on a distinguished road
Hi, Formers:
It is a problem confused me a lot of days. i want to use the multiphaseInterFoam to stimulate different fluid inject in a cavity sequentially.
If one fluid injection is finished, how can the setfieldsDict be set for the next one?
If the setfieldDict not changed, i have tried this way, the calculation of the next fluid would begin from the end of previous fluid. But the result is not OK.
So, can you give me any advise?

Best regards!
salame is offline   Reply With Quote

Old   September 20, 2013, 16:50
Default
  #2
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 137
Rep Power: 7
mgdenno is on a distinguished road
If I understand what you are trying to do, I would say that maybe you need to change the fluid entering the system using a different boundary condition and not setFields. Maybe something like the uniformFixedValue with a time varying uniformValue shown here:

http://www.openfoam.org/version2.1.0...conditions.php

I have not tried this but it is where I would start.
mgdenno is offline   Reply With Quote

Old   September 23, 2013, 03:05
Unhappy
  #3
New Member
 
salame ama
Join Date: Dec 2012
Posts: 28
Rep Power: 4
salame is on a distinguished road
thank u for your reply! Matthew.
i have seen the webpage you recommended. To my issue, the change of u is not
enough, the fluid also be changed. like:First, water; second, oil.
salame is offline   Reply With Quote

Old   September 23, 2013, 12:19
Default
  #4
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 137
Rep Power: 7
mgdenno is on a distinguished road
I would think that you could use this BC to change the boundary conditions in your alphaair and alphawater files to change the fluid coming into the domain.

in alphawater:
Code:
 
inlet
{    
        type uniformFixedValue;
        uniformValue     table 
        ( 
            (  0   1.0) 
            (100  0.0) 
        ); 
}
in alphaair:
Code:
 
inlet
{    
        type uniformFixedValue;
        uniformValue     table 
        ( 
            (  0   0.0) 
            (100  1.0) 
        ); 
}
to switch from water to air at t=100.

Like I said, I have never done this but it is where I would start.

Let me know if it works.
mgdenno is offline   Reply With Quote

Old   October 22, 2013, 10:17
Post
  #5
New Member
 
salame ama
Join Date: Dec 2012
Posts: 28
Rep Power: 4
salame is on a distinguished road
sorry to reply you so later, Matthew.
your patient reply is appreciated, thank you! i have done the steps like you say. in my case, there are 3 fluids, two Non-newtonian and one Newtonian. properties file:
Code:
.phases         
           (
               oil {
                           transportModel  powerLaw;

                            nu              nu [0 2 -1 0 0 0 0] 1;
                           rho             rho [1 -3 0 0 0 0 0] 748.5;
                              powerLawCoeffs
       {
                k             k [0 2 -1 0 0 0 0] 9.065;
                n             n [0 0 0 0 0 0 0] 0.35;
                nuMin         nuMin [0 2 -1 0 0 0 0] 0.001336;
                nuMax         nuMax [0 2 -1 0 0 0 0] 8.5504342;
        }
   
} 
           mercury {
                             transportModel  powerLaw;

                         nu              nu [0 2 -1 0 0 0 0] 1;
                        rho             rho [1 -3 0 0 0 0 0] 850.5;
                              powerLawCoeffs
       {
                k             k [0 2 -1 0 0 0 0] 7.065;
                n             n [0 0 0 0 0 0 0] 0.25;
                nuMin         nuMin [0 2 -1 0 0 0 0] 0.01336;
                nuMax         nuMax [0 2 -1 0 0 0 0] 7.4504342;
        }
   
}
               air { transportModel Newtonian ; nu nu [ 0 2 -1 0 0 0 0 ] 1.48e-03 ; rho rho [ 1 -3 0 0 0 0 0 ] 1 ; } 

           );

refPhase        air;
And i set the alphaoil and alphamercy like this: oil
Code:
 inlet
    {
        type            uniformFixedValue;
        uniformValue    table
        
(
(0 1)
(0.2 0)
)
;
    }
. mer
Code:
inlet
    {
        type            uniformFixedValue;
        uniformValue    table
       (
 
         (0 0)
         (0.2 1)
)
;
the V file:
Code:
inlet
    {
        type            uniformFixedValue;
        uniformValue    table
        
       (
              (0 (1 0 0))
              (0.2 (0.5 0 0))

        );
    }
. the 3 files are used the new boundrys like you say. But, i felt uncertain about the result (0.2s) :截图00.jpg, the left is the inlet , blue is the air and just the 0.2 second result, i am puzzled about this picture like water-jetting. Does the multiphaseInterFoam solver just for newtonian fluid? By the way, do you know how is the new BC setting as we mentioned in the version 1.6-ext?


Best regards!
salame
salame is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
multiphaseInterFoam: timestep error by simulating a co-extrusion nozzle Quatschinsky OpenFOAM Running, Solving & CFD 7 March 27, 2014 06:08
SetFieldsDict: Non uniform density physics1 OpenFOAM Running, Solving & CFD 1 May 8, 2013 16:35
convergence issue with multiphaseInterFoam sachinlb OpenFOAM Running, Solving & CFD 2 October 12, 2012 11:45
alphaEqn in multiphaseInterFOAM 2.1 sina.s OpenFOAM Programming & Development 2 March 9, 2012 05:30
VOF fraction in multiphaseInterFoam mahaputra OpenFOAM Running, Solving & CFD 1 August 27, 2009 08:15


All times are GMT -4. The time now is 10:36.