- **OpenFOAM Programming & Development**
(*https://www.cfd-online.com/Forums/openfoam-programming-development/*)

- - **perpendicular projection**
(*https://www.cfd-online.com/Forums/openfoam-programming-development/97339-perpendicular-projection.html*)

perpendicular projectionHello,
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 08:59. |