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

Problems with volScalarField

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 28, 2005, 04:30
Default Why not? volScalarField pf =
  #1
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
Why not?
volScalarField pf = pow(10,25.0)*pow(max(0.0, alpha - alphaMax), 10.0);

different results?
does the pf field differ or what do you mean?

N
niklas is offline   Reply With Quote

Old   November 28, 2005, 14:42
Default I'm still working on twoPhaseE
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
I'm still working on twoPhaseEulerFoam to make granular pressure, radial distribution functions, and frictional stress models runtime selectable.

The problem happens with frictional stress models, where pf is the frictional pressure.

If I use the first piece of code, I get solution singularity at the second time step due to an abnormal value in the granular pressure, which is related to pf.
If I use the second piece of code, the calculation goes on regularly.
However, it seems the initialization of the pf field has some effect. I mean the forAll cycle doesn't seem to overwrite all initial values also in the second piece of code.

Alberto
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   November 28, 2005, 17:05
Default You were right. The problem wa
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
You were right. The problem was that I forgot the Foam::max(alpha-alphaMax, 0.0).

Thank you again

P.S. How are you?

Alberto :-)
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto 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
Reading a volScalarField from File schnitzlein OpenFOAM Running, Solving & CFD 5 November 23, 2022 13:40
Is there a way to change the name a volScalarField liu OpenFOAM Running, Solving & CFD 2 October 18, 2007 18:49
What is the difference between tmp%2360volScalarField and volScalarField haibinli OpenFOAM Running, Solving & CFD 0 November 3, 2006 05:06
VolScalarField dimensionSet luca OpenFOAM Running, Solving & CFD 0 July 10, 2006 10:55
Incorporate known VolScalarfield in PDE stephan OpenFOAM Running, Solving & CFD 0 November 29, 2005 08:13


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