CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Edge (http://www.cfd-online.com/Forums/edge/)
-   -   External volume forces (http://www.cfd-online.com/Forums/edge/86657-external-volume-forces.html)

merazzi March 29, 2011 09:20

External volume forces
 
Hej,

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

Best regards,

Silvio

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?

Adam

merazzi April 1, 2011 04:24

Yes, we have source code.

Silvio

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
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


All times are GMT -4. The time now is 04:05.