|February 15, 2012, 06:01||
Join Date: Jul 2010
Posts: 4Rep Power: 6
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:
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)
|perpendicular, projection, projection tensor|
|Thread||Thread Starter||Forum||Replies||Last Post|
|projection method(3D)||ehsan_am86||FLOW-3D||1||August 28, 2010 21:45|
|BlockMeshmergePatchPairs||hjasak||OpenFOAM Native Meshers: blockMesh||11||August 15, 2008 07:36|
|projection method||sudhakar||Main CFD Forum||4||December 8, 2007 02:25|
|A Second-Order Projection Method for Incompressibl||Ivan||Main CFD Forum||5||April 29, 2001 07:27|
|Projection methods||ZY Sun||Main CFD Forum||9||September 15, 2000 18:00|