CFD Online Discussion Forums

CFD Online Discussion Forums (
-   Edge (
-   -   External volume forces (

merazzi March 29, 2011 09:20

External volume forces

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

Best regards,


jka March 31, 2011 02:32

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


merazzi April 1, 2011 04:24

Yes, we have source code.


jka April 1, 2011 14: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

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


All times are GMT -4. The time now is 02:53.