CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[Publications] Darcy-Forchheimer Online Calculator

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

Like Tree4Likes
  • 2 Post By Tobi
  • 1 Post By Tobi
  • 1 Post By MKdoubleB

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 14, 2019, 19:01
Default Darcy-Forchheimer Online Calculator
  #1
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,


I just want to inform you, that I implemented an Darcy-Forchheimer Calculator on my website. It is not yet finished and I have to check the results but the next days I will finish it. Nevertheless, the first indications seems to work fine.

https://holzmann-cfd.com/community/b...cy-forchheimer
MTomić and MKdoubleB like this.
__________________
Keep foaming,
Tobias Holzmann

Last edited by Tobi; June 22, 2020 at 03:30.
Tobi is offline   Reply With Quote

Old   July 24, 2022, 07:17
Default
  #2
New Member
 
Matej Tomić
Join Date: Jul 2022
Posts: 5
Rep Power: 4
MTomić is on a distinguished road
Hello Tobias,

I checked out your calculator tool for the Darcy-Forchheimer coefficients. I am unable to obtain the Darcy contribution identical to your calculator, for a 1D case. I am trying to model a heat exchanger with a certain porosity. It would be nice to see how your tool calculates the Darcy contribution.

D_{\mathrm{Matej}} = 2.82 \cdot 10^{7}\ 1/\mathrm{m}^{2}
D_{\mathrm{Tobias}} = 2.56 \cdot 10^{-3}\ 1/\mathrm{m}^{2}


Regards,
Matej
MTomić is offline   Reply With Quote

Old   August 2, 2022, 08:56
Default
  #3
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Its still crap and buggy. I had not time to check everything in detail.

Matej = Matej Forman? So its my colleague. Do you have the calculation base from him? (community@holzmann-cfd.com).
Otherwise you can tell me the material-name and I can check it myself.

Thanks, Tobi
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   August 30, 2022, 06:18
Default
  #4
New Member
 
Matej Tomić
Join Date: Jul 2022
Posts: 5
Rep Power: 4
MTomić is on a distinguished road
Hello Tobi,

Funny coincidence, although I must state that I am in no way, shape or form affiliated with Matej Forman. Inspired by your calculator tool, I made a similar one (only with the help of MATLAB). In order to compare the results with your calculator tool, please use the following info:

v_{1} = 0.471\ \mathrm{m/s},\  \Delta p_{1} = 6.16\ \mathrm{Pa}
v_{2} = 11.26\ \mathrm{m/s},\  \  \Delta p_{2} = 1007.04\ \mathrm{Pa}



When inserting aforementioned inputs into your calculator, the Darcy contribution output differs severely.
Attached Images
File Type: jpg Screenshot_1.jpg (66.5 KB, 40 views)
__________________
Regards,
Matej
MTomić is offline   Reply With Quote

Old   September 2, 2022, 16:57
Default
  #5
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hey. Nice and thanks for your feedback.

Using your darcy and forchheimer values by using the equation given on my website (first one right hand side: https://holzmann-cfd.com/community/b...cy-forchheimer), we end up with the graph attached:

  • Darcy Holzmann: 0.0013535060215094123
  • Forchheimer Holzmann: 171.6287292707513


  • Darcy Matej: 2.927e7
  • Forchheimer Matej: 258.2

dp = (D*U+F*U^2)*dx using dx = 0.073 m





So with your extreme linear contribution, its clear that you will get such an offset. However, my curve is not fitting too good either. Shame on me that I still have that buggy guy available.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   September 3, 2022, 10:12
Default
  #6
New Member
 
Matej Tomić
Join Date: Jul 2022
Posts: 5
Rep Power: 4
MTomić is on a distinguished road
Hello Tobi,

Thanks for taking time and looking into this topic.

The deviation of the Forchheimer contribution is self-explanatory, as your tool only takes into account two data points. Therefore, the curve fit is of lower fidelity.

However, I think that there is some underlying mystery about the Darcy contribution.

If we assume that the \Delta p = aU + bU^{2}, then the Darcy and Forchheimer contributions are determined in the following manner:

D = \frac{a}{\mu \Delta x}
F = \frac{2b}{\rho \Delta x}

Just by looking at the expression for the Darcy contribution, we can note that in the case of air (dynamic viscosity in the vicinity of 1.73e-5 Pa s) that the Darcy contribution should be a rather large number.

So to return to the actual point of my question, how does your tool determine the Darcy contribution?
__________________
Regards,
Matej
MTomić is offline   Reply With Quote

Old   September 3, 2022, 15:22
Default
  #7
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
The main problem in the whole calculation is that I need (or we need) to avoid negative numbers as they are treated differently by OpenFOAM (a negative number means blockage).

I second your statement and if you only take the Darcy-Calculation into account (in my tool), you will get D to be around 119246891. Hence, a very large number as you have.

My tool is doing some magic
No, its simply some derivation - need to recheck but its in the javascript that you can checkout manually if you go through my source-code of the website (view-source:https://holzmann-cfd.com/modules/mod...Forchheimer.js)

Code:
                /* OLD 
                var f =
                    parseFloat(dp1)/parseFloat(dx)
                  /(
                      (
                          (
                              parseFloat(dp1)*parseFloat(u2)*parseFloat(u2)
                             -parseFloat(dp2)*parseFloat(u1)*parseFloat(u1)
                          )
                          /
                          (
                            parseFloat(dp2)*parseFloat(u1)
                           -parseFloat(dp1)*parseFloat(u2)
                          )
                      )*u1+parseFloat(u1)
                  );

                var r =
                    (
                        parseFloat(dp2)/parseFloat(dx)
                       -parseFloat(f)*parseFloat(u2)*parseFloat(u2)
                    )/parseFloat(u2);
                */

                var f =
                    (
                        parseFloat(dp2)*parseFloat(u1)*parseFloat(dx)
                      - parseFloat(dp1)*parseFloat(u2)*parseFloat(dx)
                    )
                    /
                    (
                        (-1 * parseFloat(u1)*parseFloat(dx)*parseFloat(u2)*parseFloat(dx))
                      + (parseFloat(u2)*parseFloat(u2)*parseFloat(dx)*parseFloat(dx)*parseFloat(u1))
                    );

                var r =
                    (
                        parseFloat(dp1) - parseFloat(f)*parseFloat(u1)*parseFloat(dx)
                    )
                    /
                    (
                        parseFloat(u1)*parseFloat(dx)
                    );

                // 10 % more than the given datapoint
                var largerU = 0;
                var largerdp = 0;

                if (parseFloat(dp1) > parseFloat(dp2)) { largerdp = dp1 } else { largerdp = dp2; }
                if (parseFloat(u1) > parseFloat(u2)) { largerU = u1 } else { largerU = u2; }

                xGiven = largerU;
                yGiven = largerdp;

                console.log("xGiven = " + xGiven);
                console.log("yGiven = " + yGiven);

                xTotal = xGiven * 1.1;
                deltax = xTotal / 20;

                var i;
                for (i = 0; i < 50; i++)
                {
                    xValue.push(i*deltax);
                    yValue.push(r*i*deltax*dx+f*(i*deltax)*(i*deltax)*dx);
                }
            }
I wanted to introduce some square-root calculation (while giving the opportunity to have a lot of measurement data available and built the best fit here). However, doing so lead to negative values for the linear contribution (D) which is not supported by FOAM at the moment. I had some discussions about that behavior with my colleague Kutalmis and he would also agree. Nevertheless, having only the Darcy contribution, a negative value does not make sense at all.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   May 11, 2023, 06:58
Default
  #8
New Member
 
Markus Bähner
Join Date: May 2023
Posts: 4
Rep Power: 7
MKdoubleB will become famous soon enough
Quote:
Originally Posted by Tobi View Post

My tool is doing some magic
No, its simply some derivation - need to recheck but its in the javascript that you can checkout manually if you go through my source-code of the website (view-source:https://holzmann-cfd.com/modules/mod...Forchheimer.js)

Hey Tobi,


first of all, thanks for making your contributions and the tool available to us all!
I just checked out the source code, since i was confused comparing the online calculators small result for the darcy contribution with the very large value i came up with in my own calculations. In Line 397 of the code the coefficients are calculated as:

D.value = r*mu.value;
F.value = 2*f/rho.value;


Isn't there a little mistake here? r gets multiplied by mu, yet per definition the Darcy coefficient is:
D = r/mu
as Matej stated above. F seems correct to me. Or is this the "magic" you speak of? Im not very experienced with Porous zones in OpenFOAM so please let me know if i am missing anything!


Best regards,


Markus
MKdoubleB is offline   Reply With Quote

Old   May 11, 2023, 12:59
Default
  #9
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hey Markus,


you are excellent. I love people who are investigating into the code and know how to get to the simple javascript. (EDIT: I saw that I posted the link - but anyway - nice that you checked it).


I agree - as it is already postulated on my website - D = r / mu.

In the code I have D = r * mu for the Darcy-Forchheimer calculation. We also see that the pure Darcy calculation is correct hence, for any reason, I made a mistake here.

Thanks for your hint and investigation. Code is already corrected.
Tobi
MKdoubleB likes this.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   May 11, 2023, 17:38
Default
  #10
New Member
 
Markus Bähner
Join Date: May 2023
Posts: 4
Rep Power: 7
MKdoubleB will become famous soon enough
Stuff like this happens to me all the time - youre a legend for already having the corrected one up and running.



Good yard,


Markus
MKdoubleB is offline   Reply With Quote

Old   May 12, 2023, 16:55
Default
  #11
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Do your results look better now?
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   May 24, 2023, 12:24
Default
  #12
New Member
 
Markus Bähner
Join Date: May 2023
Posts: 4
Rep Power: 7
MKdoubleB will become famous soon enough
Hey,


scusi for the late reply. Up to that point i hadnt even set up my porous model. However, initial tests with the coefficients seem to produce results which are in the expected ball park. At a later time, i will need to check against some velocity measurements of the prototype to see how far we are off.


Cheers
Tobi likes this.
MKdoubleB is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Most useful online calculator Coris Main CFD Forum 0 July 15, 2019 10:38
Darcy Forchheimer: trying to understand the code and how to run with it? student666 OpenFOAM Running, Solving & CFD 13 February 21, 2018 08:54
New Servers Launched for CFD Online Website pete Site News & Announcements 0 June 12, 2015 07:09
CFD Online Celebrates 20 Years Online jola Site News & Announcements 22 January 31, 2015 01:30
CFD Online Down due to Hardware Failure CFD Online Team Main CFD Forum 2 November 19, 2007 06:06


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