Implementation of a porous domain

megacrout · December 19, 2011, 11:19

Hi Foamers,

I am simulating reactions through a porous domain using porousExplicitSourceReactingParcelFoam. The header of the .C file of this solver reads:

Note: ddtPhiCorr not used here when porous zones are active
- not well defined for porous calculations

I guess this means that the solver is not well defined for porous domains because of the given reason, not the opposite (i.e. ddtPhiCorr is not well defined for porous domains, reason why it is not used here), right?

If so, my question is: why does ddtPhiCorr make the solver ill-defined and how can it be improved? I guess the programmers would have done it if simply using ddtPhiCorr was the answer...

Thx.

Tibo

megacrout · January 12, 2012, 07:02

Ok, I figured out what it means by opening the corresponding pEqn.H file.

However, I still do not know what the REASON for ddtPhiCorr not to be well defined is.

Furthermore, I´d be glad to know if the fact of not using this function leads to the "right" solution or simply avoids the numerical concern but returns a solution that is not as good as it could be.

Thx.

Tibo

December 19, 2011, 11:19	Implementation of a porous domain	#1
megacrout Member Tibo Join Date: Jun 2011 Posts: 68 Rep Power: 14	Hi Foamers, I am simulating reactions through a porous domain using porousExplicitSourceReactingParcelFoam. The header of the .C file of this solver reads: Note: ddtPhiCorr not used here when porous zones are active - not well defined for porous calculations I guess this means that the solver is not well defined for porous domains because of the given reason, not the opposite (i.e. ddtPhiCorr is not well defined for porous domains, reason why it is not used here), right? If so, my question is: why does ddtPhiCorr make the solver ill-defined and how can it be improved? I guess the programmers would have done it if simply using ddtPhiCorr was the answer... Thx. Tibo

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Radiation at interface between fluid and porous domain	Hitch8	CFX	19	April 20, 2015 06:24
Porous domain set-up from single pressure loss value	siw	CFX	1	December 8, 2011 16:36
Should I be using True or Superficial Velocity for Porous Domain?	siw	CFX	2	November 17, 2010 07:52
Porous Domain	ijk	CFX	1	August 30, 2010 19:03
Source within a porous domain	Hitch8	CFX	4	July 2, 2010 03:01

January 12, 2012, 07:02		#2
megacrout Member Tibo Join Date: Jun 2011 Posts: 68 Rep Power: 14	Ok, I figured out what it means by opening the corresponding pEqn.H file. However, I still do not know what the REASON for ddtPhiCorr not to be well defined is. Furthermore, I´d be glad to know if the fact of not using this function leads to the "right" solution or simply avoids the numerical concern but returns a solution that is not as good as it could be. Thx. Tibo