|
[Sponsors] |
May 31, 2011, 02:42 |
lowerPtr_ or upperPtr_ unallocated
|
#1 |
New Member
Gregory
Join Date: Nov 2010
Location: Dresden
Posts: 14
Rep Power: 15 |
Hi there,
I am trying to test solve a modified epsilon equation within the kEpsilon model. For the sake of simplicity I first consider the following: --------------- volScalarField y_ = wallDist(mesh_).y(); volScalarField Ry_ = y_ * sqrt(k_) / nu(); scalar kappa_= 0.42; scalar Cl_ = kappa_ / pow(Cmu_.value(),0.75); scalar Aeps_ = 2*Cl_; volScalarField Leps_ = Cl_ * y_ * (1 - exp( -Ry_ / Aeps_ ) ); dimensionedScalar dummyValue_ ( "dummyValue", dimensionSet(0,0,-1,0,0,0,0), 1 ); tmp<fvScalarMatrix> epsEqn ( fvm::ddt(epsilon_) == dummyValue_ * (pow(k_, 1.5)/Leps_ - epsilon_) ); epsEqn().relax(); solve(epsEqn); bound(epsilon_, epsilon0_); ----------------- but when I try to solve it i get the message: --> FOAM FATAL ERROR: lowerPtr_ or upperPtr_ unallocated From function lduMatrix::lower() const in file matrices/lduMatrix/lduMatrix/lduMatrix.C at line 205. FOAM aborting Does anyone know where the problem lie ? Thanks a lot in advance. |
|
May 31, 2011, 11:38 |
|
#2 |
Senior Member
David Gaden
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 437
Rep Power: 21 |
I believe the problem is your use of tmp. Do you really need it? (Are you using it to reduce peak memory?) See here for a more detailed explanation:
http://openfoamwiki.net/index.php/OpenFOAM_guide/tmp If you do need to use tmp, assign a reference to its enclosed object: Code:
fvScalarMatrix& epsEqnRef = epsEqn(); Code:
epsEqnRef.relax(); solve(epsEqnRef); // etc... |
|
May 21, 2014, 13:01 |
|
#3 |
New Member
Stefano Puggelli
Join Date: Sep 2013
Posts: 4
Rep Power: 12 |
Hi,
I am having the same problem. I would like to write down an equation for liquid volume fraction with a source term due to evaporation. fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) == - Rate_alpha1_ees ); When I try to run I have this error: -> FOAM FATAL ERROR: lowerPtr_ or upperPtr_ unallocated From function lduMatrix::lower() const in file matrices/lduMatrix/lduMatrix/lduMatrix.C at line 202. Does anyone know where the problem lie ? Thanks a lot in advance. |
|
May 30, 2014, 04:03 |
|
#4 | |
Senior Member
|
I don't know much about turbulence.
By looking at the error which you have mentioned in your post. Quote:
Re-check your code somewhere something you might have missed out. Best Luck! |
||
May 31, 2014, 06:03 |
|
#5 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33 |
This is an old bug in matrix classes: you only have the diag and matrix algebra croaks. I fixed it ~5 years ago. foam-extend, of course, works.
Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
August 16, 2014, 07:27 |
|
#6 | ||
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Greetings to all!
@Stefano: Quote:
@Hrvoje: Quote:
Best regards, Bruno
__________________
|
|||
February 15, 2020, 09:18 |
|
#7 | |
Member
Join Date: Dec 2009
Location: Kanpur, India
Posts: 54
Rep Power: 16 |
Quote:
Code:
solve ( fvm::ddt(alpha1) == - Rate_alpha1_ees ); Code:
fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) == - Rate_alpha1_ees ); alpha1Eqn.solve(); Last edited by Shoonya; February 15, 2020 at 09:20. Reason: correction |
||
|
|