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

Relative Fluxes in icoDyMFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 22, 2008, 15:59
Default I'm wondering why fluxes are m
  #1
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 368
Rep Power: 14
deepsterblue will become famous soon enough
I'm wondering why fluxes are made relative to the mesh-motion at the end of the PISO loop in icoDyMFoam:


for (int corr=0; corr<nCorr; corr++)
{
rUA = 1.0/UEqn.A();

U = rUA*UEqn.H();

...
...
...

// Make the fluxes relative to the mesh motion - why?
fvc::makeRelative(phi, U);

U -= rUA*fvc::grad(p);
U.correctBoundaryConditions();
}

I'm guessing that this is required for the correctBoundaryConditions call for velocity. Can anyone confirm/explain why this is so?
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   September 23, 2008, 01:27
Default Hy Sandeep, If I remember r
  #2
Senior Member
 
Jens Klostermann
Join Date: Mar 2009
Posts: 117
Rep Power: 7
jens_klostermann is on a distinguished road
Hy Sandeep,

If I remember right is this because of the space conservation law and the additional flux which is associated with a moving boundary (mesh).

Jens
jens_klostermann is offline   Reply With Quote

Old   September 23, 2008, 01:57
Default Hi Sandeep, Read any paper
  #3
Member
 
lillberg's Avatar
 
Eric Lillberg
Join Date: Mar 2009
Location: Stockholm
Posts: 78
Rep Power: 7
lillberg is on a distinguished road
Send a message via Skype™ to lillberg
Hi Sandeep,

Read any paper on Arbitrary Lagrangian Eulerian (ALE) methods and it will all be clear. For example Charbel Farhat, "The Discrete Geometric Conservation Law and the Nonlinear Stability of ALE Schemes for the Solution of Flow Problems on Moving Grids", Journal of Computational Physics 174, 669&ndash;694 (2001) or more recent papers on mesh motion by ˇZeljko Tukovi'c.

/Eric
lillberg is offline   Reply With Quote

Old   September 23, 2008, 10:27
Default Hello Jens and Eric, Thanks
  #4
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 368
Rep Power: 14
deepsterblue will become famous soon enough
Hello Jens and Eric,

Thanks for the help, but it still isn't clear. The makeRelative function merely subtracts meshPhi from phi. However, for the next iteration of the nCorr loop, phi is over-written by:
phi = (fvc::interpolate(U) & mesh.Sf());
The line:
U -= rUA*fvc::grad(p);
doesn't require phi in any way, and the only reason phi should be relative to mesh-motion (in my opinion) is to satisfy a particular boundary-condition (at U.correctBoundaryConditions()) which requires it to be so. If none of the boundary conditions require it, is the line really necessary?

Comments?
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   September 25, 2008, 05:49
Default Hi Sandeep, Yes you are rig
  #5
New Member
 
Zeljko Tukovic
Join Date: Mar 2009
Posts: 22
Rep Power: 7
ztukovic is on a distinguished road
Hi Sandeep,

Yes you are right. You will need relative flux at that stage only if velocity boundary condition requires it. But you will also need relative flux in the turbDyMFoam solver because the transport equations of a turbulence model will require relative flux.

Regards,
Zeljko
ztukovic is offline   Reply With Quote

Old   October 6, 2009, 10:38
Default
  #6
New Member
 
Mark
Join Date: Oct 2009
Posts: 1
Rep Power: 0
mgeiger is on a distinguished road
Hello!
This thread has been quiet for a while. But what I don't understand about makeRelative and makeAbsolute is, how makeRelative can be called twice without makeAbsolute in between?
The interDyMFoam solver calls:
Code:
Abs
Rel
     if ...{
       Rel
       Abs
     }
   Rel
Obviously, the flux is made relative either inside or after the if{} although it is already relative. How should I understand this?

Here are the relevant calls with all for{} loops and if{} conditions:
Code:
while (runTime.run())
    {
...
        // Make the fluxes absolute
        fvc::makeAbsolute(phi, U);
...
        mesh.update();
...
        // Make the fluxes relative to the mesh motion
        fvc::makeRelative(phi, U);
...
        // --- PIMPLE loop
        for (int ocorr=0; ocorr<nOuterCorr; ocorr++)
        {
...
            // --- PISO loop
            for (int corr=0; corr<nCorr; corr++)
            {
...
                if (p.needReference())
                {
                    fvc::makeRelative(phi, U);
                    adjustPhi(phi, U, p);
                    fvc::makeAbsolute(phi, U);
                }
...
                // Make the fluxes relative to the mesh motion
                fvc::makeRelative(phi, U);
...
            }
        }
    }
regards
Mark
mgeiger is offline   Reply With Quote

Old   April 2, 2012, 04:46
Default
  #7
Member
 
fisch
Join Date: Feb 2010
Posts: 97
Rep Power: 6
fisch is on a distinguished road
Hello Mark,

did you find an answer to your question?

I'm interested.

thanks
fisch is offline   Reply With Quote

Old   July 20, 2012, 02:39
Default
  #8
Member
 
Yogesh Bapat
Join Date: Oct 2010
Posts: 33
Rep Power: 6
ybapat is on a distinguished road
Hello Mark,

Did you find answer to your question? Please let me know .
ybapat 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
Update Fluxes luca OpenFOAM Running, Solving & CFD 2 March 29, 2006 09:47
Fluxes selen CD-adapco 1 August 24, 2003 16:23
Fluxes validation Mavinakere FLUENT 0 May 22, 2002 23:24
Help need on fluxes Zhu Main CFD Forum 7 August 12, 2001 10:14
help needed on fluxes zhu Phoenics 0 August 7, 2001 21:35


All times are GMT -4. The time now is 05:50.