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

how to express a heat source term which depends on coordinate

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

Like Tree4Likes
  • 1 Post By bryant_k
  • 2 Post By chegdan
  • 1 Post By a19910112a

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 22, 2012, 09:34
Default how to express a heat source term which depends on coordinate
  #1
New Member
 
kob
Join Date: Nov 2011
Posts: 28
Rep Power: 14
bryant_k is on a distinguished road
Hello everyone
I am doing a heat transfer problem with source term. Source term itself depends on the x,y coordinates. The expression of the volumetric heat source is:
Code:
when -0.2<x<0.2,-0.3<y<0.3
q=20*exp(-10*y)
otherwise 
q=0
should I use the functionSetFields? If I don't use the function,can you tell me how to write the code when declaring the heat soure q?

There is another question which is related to the x,y coordinate.
Is it correct to express it like:
Code:
volScalarField x = mesh.C().component(vector::X);
volScalarField y = mesh.C().component(vector::Y);
or
Code:
pos().x
pos().y
what is the difference between them?

regards!

bryant
raj kumar saini likes this.
bryant_k is offline   Reply With Quote

Old   August 25, 2012, 08:37
Default
  #2
New Member
 
kob
Join Date: Nov 2011
Posts: 28
Rep Power: 14
bryant_k is on a distinguished road
Can you tell me about the above things if you know about that?
Thank you very much!

regards!

bryant
bryant_k is offline   Reply With Quote

Old   April 8, 2013, 00:54
Default
  #3
Member
 
Neilson Whit
Join Date: Aug 2011
Posts: 74
Rep Power: 14
wolfindark is on a distinguished road
Dear bryant_k

did you find the solution? I am also working on it...
pleased if you can inform me.
wolfindark is offline   Reply With Quote

Old   April 8, 2013, 03:22
Default
  #4
Senior Member
 
Mohammad Shakil Ahmmed
Join Date: Oct 2012
Location: AUS
Posts: 137
Rep Power: 14
ahmmedshakil is on a distinguished road
@ wolfindK
you can easily use funkySetFields utility, it's easy to apply, http://openfoamwiki.net/index.php/Co...funkySetFields
ahmmedshakil is offline   Reply With Quote

Old   April 9, 2013, 16:10
Default
  #5
Senior Member
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 621
Rep Power: 0
chegdan will become famous soon enoughchegdan will become famous soon enough
@wolfindark

Have you looked at using the new fvOptions capabilities? Here I have defined a source names "energySource1" with a time duration at a single point. The injection rate coefficients are specified as pairs of Su-Sp coefficients e.g. h (10 0) where the 10 is the explicit component and 0 is the implicit component. This is a linearized source term (read Patankar's book for more details and a source of this method)

Code:
	energySource1
	{
		type            scalarSemiImplicitSource;
		active          true;
		timeStart       0.2;
		duration        2.0;
		selectionMode   points;
		points
		(
		    (2.75 0.5 0)
		);

		scalarSemiImplicitSourceCoeffs
		{
		    volumeMode      absolute;
		    injectionRateSuSp
		    {
		        h           (10 0);
		    }
		}
	}
for selectionMode you can use cellZone and then choose a cellZone where you want a heat source to be. For more information, look at the source files located at

Code:
$FOAM_SRC/fvOptions/sources/general/semiImplicitSource
good luck
raj kumar saini and adiraman9 like this.

Last edited by chegdan; April 9, 2013 at 18:29.
chegdan is offline   Reply With Quote

Old   April 9, 2013, 23:57
Default
  #6
Member
 
Neilson Whit
Join Date: Aug 2011
Posts: 74
Rep Power: 14
wolfindark is on a distinguished road
Thank you chegdan,

For my case I needed a constant heat source, therefore, I used another method which was mentioned in this forum sometime ago. I just added explicit source term to temperature equation (in TEqn.H). Then I set the field values with setFields, a spherical heat source covering a particular region in fluid region.
So it works fine until now.

I am not sure but if time dependency is needed, it can be implemented in TEqn.H, What do you think..?
wolfindark is offline   Reply With Quote

Old   April 10, 2013, 10:14
Default
  #7
Senior Member
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 621
Rep Power: 0
chegdan will become famous soon enoughchegdan will become famous soon enough
Quote:
I am not sure but if time dependency is needed, it can be implemented in TEqn.H, What do you think..?
Yes, anything is possible but it will require programming and if are comfortable with that then that is a good option.
chegdan is offline   Reply With Quote

Old   November 29, 2013, 07:19
Default
  #8
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 18
samiam1000 is on a distinguished road
Dear All,

pardon the interruption. I have a question about fvOption and I think you can help me.

I would like to add a body force on certain cells.

I have tried to write something like

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

momentumSource
{
    type            vectorExplicitSetValue;
    active          on;            //on/off switch
    selectionMode   all;       //cellSet // points //cellZone

    vectorExplicitSetValueCoeffs
    {
        injectionRate
        {
            F    ( 0.1335 0 0 );
        }
    }
}
But I cannot und how to tell OpenFOAM that F is a force. Is that possible?

Thanks a lot,
Samuele
samiam1000 is offline   Reply With Quote

Old   December 1, 2014, 16:32
Default
  #9
New Member
 
zech
Join Date: Oct 2014
Location: Cambridge,England
Posts: 22
Rep Power: 11
a19910112a is on a distinguished road
Does scalarSemiImplicitSourceCoeffs{ h (10,0)} mean that the source is S=10+0*x, where S is the magnitude of the source and x is the coordinate. If that is true, then:
1.this method is applicable when the source is linearly related to the coordinate. What if it is not a linear function. Such as a cosine or experiential function?
2.Even worse, it is not even an function: it's just a matrix of random numbers.
3. what if the value is changing in both x and y direction?
raj kumar saini likes this.
a19910112a is offline   Reply With Quote

Reply

Tags
heat source

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
GPU Linear Solvers for OpenFOAM gocarts OpenFOAM Announcements from Other Sources 37 August 17, 2022 15:22
Source Term UDF VS Porous Media Model pchoopanya Fluent UDF and Scheme Programming 1 August 28, 2013 07:12
UDF source term jerome_ FLUENT 2 July 11, 2011 12:55
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 02:24
Source Term in Cylindrical Coordinate. Apurv Main CFD Forum 0 January 26, 2000 08:15


All times are GMT -4. The time now is 00:55.