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/)
-   -   How to solve this kind of continuity equation? (http://www.cfd-online.com/Forums/openfoam-programming-development/90363-how-solve-kind-continuity-equation.html)

shddx1 July 7, 2011 16:31

How to solve this kind of continuity equation?
 
Hi, I am a beginner of OpenFOAM. I try to use it to solve an equation in the mathematical form of:

div [sinh(|grad(phi)|)*div(phi)/|grad(phi)|] = 0

where phi is a volScalarField, |grad(phi)| is the magnitude of grad(phi), and sinh is the hyperbolic sine. I wrote the following codes, but it does not work.

volVectorField GradPhi=fvc::grad(phi);
surfaceScalarField MagGradPhi=mesh.Sf()*fvc::snGrad(phi);
surfaceScalarField SinhMagGradPhi=sinh(MagGradPhi)/MagGradPhi;
solve
(
fvm::div(SinhMagGradPhi, GradPhi)
);

The error is: conversion from ‘Foam::tmp<Foam::GeometricField<Foam::Vector<doubl e>, Foam::fvsPatchField, Foam::surfaceMesh> >’ to non-scalar type ‘Foam::surfaceScalarField’ requested

I don’t understand this problem. Do you know how to solve above equation? Any opinion is welcome! Thanks!


All times are GMT -4. The time now is 11:53.