CFD Online Logo CFD Online URL
Home > Forums > Phoenics

Source terms linearisation

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

LinkBack Thread Tools Display Modes
Old   February 10, 2006, 11:01
Default Source terms linearisation
Posts: n/a
Hi every body,

Is it a paper or a document for the linearisation (patch & coval) of complex source terms in PHOENICS equations.

  Reply With Quote

Old   February 10, 2006, 13:52
Default Re: Source terms linearisation
Posts: n/a
Document, papers books. Since Spalding inventing this formulation years ago everyone uses it.

Book: Patankar and Spalding "Numerical Heat Transfer"

Lecture : E:\phoenics\d_polis\d_lecs\GENERAL\BCOND.HTM

  Reply With Quote

Old   February 13, 2006, 12:13
Default Re: Source terms linearisation
Posts: n/a
Yes, but I need some special examples. I want to linearise a source term for a water wave generation sourec which contains third order of stokes expression.

Can I send expressions (JPG Format) on PHOENICS User Forum or to PHOENICS experts.

Seriously yours.
  Reply With Quote

Old   February 16, 2006, 05:14
Default Re: Source terms linearisation
Posts: n/a
Yes, the text book, "Numerical Heat Transfer" by Patankar and Spalding is a great book for understanding this concept, its easy to read and was the text that made the whole thing understandable to everyone.

In Polis, (the PHOENICS encyclopedia), there is a lecture regarding the mathematics behind CFD finite Volume, this has a boundary condition section: Polis.htm \phoenics\d_polis\d_lecs\GENERAL\BCOND.HTM

Now quickly, 'boundary conditions' and 'source terms' are handled the same way in CFD codes.

They need to be set as BC = T.Co(Val-Phi)

T is a geometric multiplier usually area or volume Co stands for coefficient and Val for Value and PHI is the value of the solved variable at the cell in question.

So as an example : ---------------- Take a radiation source sigma.emmisivity(Tv^4-Tp^4) sigma = stephan boltzman emmisivity(e) = emmissivity Tv^4 is 'receiving temperature' i.e. that responsible for heating/cooling Tp^4 is temperature in the cell, that being solved for.

step 1/ work out the units of this source = Watts/m^2, so the source needs to have a geometric multiplier of AREA type to make it into WATTS for the energy equation.

step2/ Rearrange source in linear form if possible i.e. T.Co(Val-Phi)form --> (Tv^4-Tp^4)= (Tv+Tp)(Tv^2+Tp^2)[Tv-Tp]

step3/ Source = Area.sigma.e.(Tv+Tp)(Tv^2+Tp^2)[Tv-Tp] So T=Area ; Co = sigma.e.(Tv+Tp)(Tv^2+Tp^2); Val=Tv

step4/ apply in PHOENICS as PATCH and COVAL, through sources panel in menu, or inform or by hand, many types of sources are automatically handled by PHOENICS.

so PATCH -> defines the region over which this source will exist

COVAL -> defines the source

step5/ IN this case the source will be , (assuming north facing wall) :

Patch(RADIAT1,NORTH,1,nx,ny,ny,1,nz,1,lstep) Coval(RADIAT,TEM1,sigma.e.(Tv+Tp)(Tv^2+Tp^2),Tv)

This should give a nett source in WATTS at the bottom of the RESULT file for you to check.

Write down your source and see if you can linearise it, but remember some sources can not be linearised, and certainly if it does not contain the variable for which you are solving.

In these cases we just lump the whole source in using Co=FIXFLU and VAL='the SOURCE', with an appropriate geometric multiplier to give the correct units for the equation.

For momentum = Newtons For Temperature = Watts

Look at the lectures and read the Patankar/Spalding book.

  Reply With Quote


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

Similar Threads
Thread Thread Starter Forum Replies Last Post
Source Terms Babar FLUENT 0 September 11, 2007 01:11
Source Terms Rajan CFX 1 March 14, 2005 10:26
Source Term Linearisation N.P CFX 0 October 24, 2004 23:35
UDF source terms Andrew Garrard FLUENT 1 October 30, 2003 11:48
UDF for source terms Tom FLUENT 1 April 26, 2002 01:13

All times are GMT -4. The time now is 14:21.