External volume forces
 

Hej,

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

Best regards,

Silvio

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

Adam

Yes, we have source code.

Silvio

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