CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

lowerPtr_ or upperPtr_ unallocated

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

Like Tree1Likes
  • 1 Post By hjasak

Reply
 
LinkBack Thread Tools Display Modes
Old   May 31, 2011, 02:42
Default lowerPtr_ or upperPtr_ unallocated
  #1
New Member
 
Gregory
Join Date: Nov 2010
Location: Dresden
Posts: 14
Rep Power: 7
catapult is on a distinguished road
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.
catapult is offline   Reply With Quote

Old   May 31, 2011, 11:38
Default
  #2
Senior Member
 
David Gaden
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 436
Rep Power: 14
marupio is on a distinguished road
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();
Then use the reference for all operations on the matrix:

Code:
epsEqnRef.relax();
solve(epsEqnRef);
// etc...
[/CODE]
marupio is offline   Reply With Quote

Old   May 21, 2014, 13:01
Default
  #3
New Member
 
Stefano Puggelli
Join Date: Sep 2013
Posts: 4
Rep Power: 5
StefanoP is on a distinguished road
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.
StefanoP is offline   Reply With Quote

Old   May 30, 2014, 04:03
Default
  #4
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 10
Tushar@cfd is on a distinguished road
I don't know much about turbulence.

By looking at the error which you have mentioned in your post.

Quote:
From function lduMatrix::lower() const
in file matrices/lduMatrix/lduMatrix/lduMatrix.C at line 205.
When a user tries to use a faulty solver file, this error message is displayed.

Re-check your code somewhere something you might have missed out.

Best Luck!
Tushar@cfd is offline   Reply With Quote

Old   May 31, 2014, 06:03
Default
  #5
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,783
Rep Power: 22
hjasak will become famous soon enough
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
ngj likes this.
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   August 16, 2014, 07:27
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,648
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings to all!

@Stefano:
Quote:
Originally Posted by StefanoP View Post
Code:
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 ?
Unfortunately it's not enough information for me to try to reproduce the same error. If possible, provide more details so that I/we can confirm in which OpenFOAM versions this problem is present in.


@Hrvoje:
Quote:
Originally Posted by hjasak View Post
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.
I'm glad that OpenFOAM isn't as security critical as OpenSSL Otherwise this should have been reported as a CVE: http://en.wikipedia.org/wiki/Common_..._and_Exposures

Best regards,
Bruno
__________________
wyldckat 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



All times are GMT -4. The time now is 07:19.