External volume forces

merazzi · March 29, 2011, 10:20

Hej,

How does one add external volume forces in Edge? Is there a documentation on this topic?

Best regards,

Silvio

jka · March 31, 2011, 03:32

I think you have to go to source code to do it. Do you have access to the source code?

Adam

merazzi · April 1, 2011, 05:24

Yes, we have source code.

Silvio

jka · April 1, 2011, 15:10

It depends where you define your external force - in mesh points or mesh edges.

For edges:
Check how is it done for other cases. Have a look at the routine called gvortexgen_m.f90. It is called from FSOURCE only at the highest level of the multigrid cycle which I think is the best way. Then in gvortexgen_m.f90 the source term is added to the viscous residuals (line 222 - 235)
Because I define source term on edge I redistributed it to point on the edge with indeces IP1 and IP2. If you define your force directly in points, you do not have to loop over edges and instead make a loop directly over
points.

For points:

Check path fsource_m.f90 -> gsoukz_m.f90 which calls
source terms for turbulence models. They do loop over points.

Hope it helps

Adam

March 29, 2011, 10:20	External volume forces	#1
merazzi New Member Silvio Merazzi Join Date: Mar 2011 Posts: 2 Rep Power: 0	Hej, How does one add external volume forces in Edge? Is there a documentation on this topic? Best regards, Silvio

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
CGNS Compiling	Diego	Main CFD Forum	17	December 21, 2014 02:40
[blockMesh] BlockMesh FOAM warning	gaottino	OpenFOAM Meshing & Mesh Conversion	7	July 19, 2010 15:11
On the damBreak4phaseFine cases	paean	OpenFOAM Running, Solving & CFD	0	November 14, 2008 22:14
fluent add additional zones for the mesh file	SSL	FLUENT	2	January 26, 2008 12:55
[blockMesh] Axisymmetrical mesh	Rasmus Gjesing (Gjesing)	OpenFOAM Meshing & Mesh Conversion	10	April 2, 2007 15:00

March 31, 2011, 03:32		#2
jka New Member Adam Jirasek Join Date: Mar 2011 Posts: 18 Rep Power: 15	I think you have to go to source code to do it. Do you have access to the source code? Adam

April 1, 2011, 05:24		#3
merazzi New Member Silvio Merazzi Join Date: Mar 2011 Posts: 2 Rep Power: 0	Yes, we have source code. Silvio

April 1, 2011, 15:10		#4
jka New Member Adam Jirasek Join Date: Mar 2011 Posts: 18 Rep Power: 15	It depends where you define your external force - in mesh points or mesh edges. For edges: Check how is it done for other cases. Have a look at the routine called gvortexgen_m.f90. It is called from FSOURCE only at the highest level of the multigrid cycle which I think is the best way. Then in gvortexgen_m.f90 the source term is added to the viscous residuals (line 222 - 235) Because I define source term on edge I redistributed it to point on the edge with indeces IP1 and IP2. If you define your force directly in points, you do not have to loop over edges and instead make a loop directly over points. For points: Check path fsource_m.f90 -> gsoukz_m.f90 which calls source terms for turbulence models. They do loop over points. Hope it helps Adam