[cfdonline2mohsen]

Dear Foamers
Hi!

I want to add a constant volumetric source term to the scalar transport equation which acts only in a particular volumetric region of the domain.

I modified the scalar transport equation as follows:

solve
	(
	  fvm::ddt(T)
	  + fvm::div(phi, T)
	  - fvm::laplacian(DT, T)
	  ==
	  sourceValue
	);

Then I defined the sourceValue in the createFields.H file as follows:

dimensionedScalar sourceValue(twoPhaseProperties.lookup("sourceValue"));

It will read the sourceValue from the transportProperties file of my case. In the transportProperties file,I gave the sourceValue as follows:

sourceValue      sourceValue [0 0 -1 0 0 0 0] 0.42;

By using this approach, the sourceValue=0.42 is used for the entire domain but I want to use the sourceValue=0.42 for only a specific volumetric region inside my domain.

I have two questions:

1) How can I define this volumetric region inside my domain? is there any utility?

2) How can I change the transportProperties file to consider sourceValue=0.42 only for this specific volumetric region and sourceValue=0 for the rest of the domain?


Any help in this regard would be much appreciated.
Thank you so much in advance.