CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Volumic source term

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 19, 2010, 12:11
Default Volumic source term
  #1
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 154
Rep Power: 7
Aurelien Thinat is on a distinguished road
Hello everybody,

I didn't find any topic on this subject so here I am. I would like to create a volume source term in a fluid region.

I have a fluid domain, and I need that a part of this fluid generates a scalar source term, which could be temperature or inert species.

Is there any function in OpenFoam for this purpose ? I would use reactingFoam after disabling the chemistry.

Thank you all.

Aurélien
Aurelien Thinat is offline   Reply With Quote

Old   November 29, 2010, 09:24
Default volume sources and sink in openFoam
  #2
Member
 
Antonio Liggieri
Join Date: Aug 2010
Posts: 66
Rep Power: 5
alfa_8C is an unknown quantity at this point
Hello Aurélien,

A good day to you!
I am about to face exactly the same problem. I also want to use reactingFoam without chemistry reactions, and apply volume sources and volume sinks. Did you find out how to it?

Thank you in advance,
Tony
alfa_8C is offline   Reply With Quote

Old   November 29, 2010, 10:51
Default
  #3
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
Aurelien and Tony, sources are discretized by Sp(), Su(), SuSp() operators and in direct manner (fvm::div(phi,T)==b).

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   November 29, 2010, 11:41
Default
  #4
Member
 
Antonio Liggieri
Join Date: Aug 2010
Posts: 66
Rep Power: 5
alfa_8C is an unknown quantity at this point
Dear Santiago,

thanks a lot for the input! Could you post a code example of a defined source/sink respectively? I need to implement a volume source. If you've ever done something like this, I would appreciate any advice...

Best Regards, Tony
alfa_8C is offline   Reply With Quote

Old   November 30, 2010, 13:21
Default
  #5
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
Hi, looking at epsilon equation definition from standard k-epsilon model you have an example:

Code:
00223     // Dissipation equation
00224     tmp<fvScalarMatrix> epsEqn
00225     (
00226         fvm::ddt(epsilon_)
00227       + fvm::div(phi_, epsilon_)
00228       - fvm::Sp(fvc::div(phi_), epsilon_)
00229       - fvm::laplacian(DepsilonEff(), epsilon_)
00230      ==
00231         C1_*G*epsilon_/k_
00232       - fvm::Sp(C2_*epsilon_/k_, epsilon_)
00233     );
this is the implementation of second equation in http://www.cfd-online.com/Wiki/Standard_k-epsilon_model, check the term with C2 where a quadratic source is defined.

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd 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


Similar Threads
Thread Thread Starter Forum Replies Last Post
momentum source term zwdi FLUENT 13 December 5, 2013 18:58
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51
UDF Scalar Code: HT 1 Greg Perkins FLUENT 8 October 20, 2000 12:40
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 11, 2000 03:43


All times are GMT -4. The time now is 22:58.