Hyper viscosity model
Hi all,
I would like to implement an hyper viscosity model which mean that I need to add the following term in the momentum equation : -nuHV*d⁴U/dx⁴ or nuHV*laplacian(laplacian(U)) where nuHV is the hyper viscosity coefficient. I try this: Code:
volVectorField lap = fvc::laplacian(nuHV/nuHV, U); Code:
[1] --> FOAM FATAL ERROR: Then I tried this: Code:
fvVectorMatrix UEqn Thank you, Pascal |
Hi,
about the second source code you tried: You want to solve the equation for U (what you do in the first two terms), but in the "new" term, you say that you want to solve the equation for "U-fvm::laplacian(...)". About the first source code you tried: the same error as above, you want to solve the equation for U, not lap. Try to exchange the fvm and fvc for the lap equation and the laplacian term of the first source code Regards, Christian |
Thank you Christian,
It works when I tried this: Code:
volVectorField lap = fvc::laplacian(nuHV/nuHV, U); Code:
volVectorField lap = fvm::laplacian(nuHV/nuHV, U); Regards, Pascal |
All times are GMT -4. The time now is 22:11. |