CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   perpendicular projection (http://www.cfd-online.com/Forums/openfoam-programming-development/97339-perpendicular-projection.html)

Janez Turk February 15, 2012 06:01

perpendicular projection
 
Hello,

I'm developing solver for vector field (denoted n) that has only direction. It's length is always a unit. Therefore I must perpendiculary project any forces that act on the vector field.

The source term is projected by multiplication with projection tensor (I-n*n):
(I-n*n)&(source term of force).

Laplacian term can maintain implicit nature because it can take volTensorField as argument:
fvm:laplacian((I-n*n),n)

I would like to have the advection term also projected, as this is crucial to the solution, and maintain implicit nature. What shoud I do with fvm:div(phi,n) term ?

The equation looks like this:

fvm::ddt(n) + fvm::div(phi,n) - fvm::laplacian((I-n*n),n) == (I-n*n)&(many source terms)

Thanks,
regards,
Janez


All times are GMT -4. The time now is 09:34.