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

FluxRequired entry

Register Blogs Community New Posts Updated Threads Search

Like Tree13Likes
  • 10 Post By hjasak
  • 3 Post By Topcresix

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 26, 2007, 13:39
Default Can someone explain briefly wh
  #1
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 21
msrinath80 is on a distinguished road
Can someone explain briefly what this entry is used for?
msrinath80 is offline   Reply With Quote

Old   April 26, 2007, 13:48
Default Yup, easy. Thin of the calcul
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33
hjasak will become famous soon enough
Yup, easy. Thin of the calculation of the fluxes in the pressure equation. Once everything is done properly, the pressure part of the flux is equal to the off-diagonal matrix coefficient related to the face, multiplied by the pressure field (difference, or for compressible flow, a slightly more complex).

Now, think of discretisation of a Laplacian: this one also got the non-orthogonal correction, which is explicit. The correction therefore also appears in the flux and must be fully consistent: if not, you get a mass continuity error.

The best way to accumulate all corrections of this kind is to simply collect the stuff during matrix assembly. Because this costs money (storage and operations), you don't want to do this for every matrix assembly. Thus, a fluxRequired flag means "collect the explicit contributions", after which you can ask the matrix for the flux (or flux correction). If you don't need the flux, you don't do fluxRequired and all is well!

Clear?

Hrv
kaifu, fumiya, salehi144 and 7 others like this.
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   May 5, 2009, 07:27
Default
  #3
New Member
 
Join Date: Apr 2009
Posts: 29
Rep Power: 17
chapman is on a distinguished road
Currently I am trying to use turbDyMFoam, but I get the following error:

Code:
BiCGStab:  Solving for Ux, Initial residual = 1, Final residual = 8.72088e-14, No Iterations 1
BiCGStab:  Solving for Uy, Initial residual = 1, Final residual = 1.22064e-13, No Iterations 1
BiCGStab:  Solving for Uz, Initial residual = 1, Final residual = 2.69711e-12, No Iterations 1
BiCGStab:  Solving for p, Initial residual = 1, Final residual = 8.52076e-08, No Iterations 404
BiCGStab:  Solving for p, Initial residual = 0.0442768, Final residual = 6.64096e-08, No Iterations 302


flux requested but p not specified in the fluxRequired sub-dictionary of fvSchemes.

    From function fvMatrix<Type>::flux()
    in file /home/chapman/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude/fvMatrix.C at line 777.
The corresponding entry in fvSchemes looks like this:

Code:
fluxRequired
{
        default no;
        pcorr;
        p;
}
Does anyone have a clue what is going wrong here?

Regards, Andy

Edit:
I found this thread: http://www.cfd-online.com/Forums/ope...ictionary.html
and checked my fvSchemes. However, the problem was not fixed for me.

Last edited by chapman; May 5, 2009 at 08:48. Reason: Added link to another thread.
chapman is offline   Reply With Quote

Old   January 23, 2015, 05:57
Default
  #4
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 839
Rep Power: 17
sharonyue is on a distinguished road
Quote:
Originally Posted by hjasak View Post
Yup, easy. Thin of the calculation of the fluxes in the pressure equation. Once everything is done properly, the pressure part of the flux is equal to the off-diagonal matrix coefficient related to the face, multiplied by the pressure field (difference, or for compressible flow, a slightly more complex).

Now, think of discretisation of a Laplacian: this one also got the non-orthogonal correction, which is explicit. The correction therefore also appears in the flux and must be fully consistent: if not, you get a mass continuity error.

The best way to accumulate all corrections of this kind is to simply collect the stuff during matrix assembly. Because this costs money (storage and operations), you don't want to do this for every matrix assembly. Thus, a fluxRequired flag means "collect the explicit contributions", after which you can ask the matrix for the flux (or flux correction). If you don't need the flux, you don't do fluxRequired and all is well!

Clear?

Hrv
Anyone can explain this more? Thanks.
sharonyue is offline   Reply With Quote

Old   May 1, 2015, 11:30
Default
  #5
New Member
 
Diego
Join Date: Mar 2015
Location: Turin,Italy
Posts: 2
Rep Power: 0
Topcresix is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
Anyone can explain this more? Thanks.
Yes, if I'm not wrong:

1) go to Santiago Marquez PhD Thesis here:
https://docs.google.com/file/d/0B2lp...VLb3lGekk/edit

2) go to page 95 (MULES + PISO Algorithm explanation)

You will see that in order to solve the flux at point 1.C, you need the to create a flux from 1.B. (fluxRequired: alpha.water).

In the same manner, in order to solve the flux at point 3.D, you need first the flux of the gradient of the pressure calculated from 3.C (fluxRequired: p_rgh).

So in this case we have to set:

fluxRequired
{
default no;
p_rgh;
alpha.water;
}

in order to storage the flux created from the discretization of alpha equation in the MULES algorithm and the flux created from the pressure equation in PISO loop.

This is what I understood from openFoam user guide
minh khang, m.b. and Teresa.Z like this.
Topcresix is offline   Reply With Quote

Reply

Tags
fluxrequired


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
space shuttle re-entry CFD Nimesh Meghpara FLUENT 0 July 5, 2008 18:34
Entry point lucian OpenFOAM Pre-Processing 1 March 13, 2007 04:04
Tangential Entry in 2D probelms Kirru FLUENT 1 October 16, 2006 08:26
Tangential Entry in 2D structure? Kiran FLUENT 0 December 6, 2005 04:42
table entry carno Siemens 1 May 28, 2005 11:29


All times are GMT -4. The time now is 01:02.