CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Modifications work in simpleFoam, not pisoFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 23, 2016, 12:15
Default Modifications work in simpleFoam, not pisoFoam
  #1
New Member
 
Join Date: Feb 2016
Posts: 10
Rep Power: 10
Aicharem is on a distinguished road
Hi Foamers,

I am trying to create a solver to model EHD thrusters in OpenFOAM. To do this I need to solve electrostatic equations (Gauss' Law and conservation of charge), and then the fluid problem. The equations in OpenFOAM format are as follows:

Code:
        solve
        (
            fvm::laplacian(phiE) + rhoE/epsilon0
        );

        rhoEFlux = -kE*mesh.magSf()*fvc::snGrad(phiE);

        solve
        (
            fvm::ddt(rhoE) + fvm::div(rhoEFlux, rhoE) + fvm::div(phi, rhoE) - fvm::laplacian(DE, rhoE)
        );

when I implement this into simpleFoam I can run cases fine, but when trying to implement them in pisoFOAM I get the following error:

Code:
--> FOAM FATAL ERROR: 
incompatible dimensions for operation 
    [rhoE[0 -3 0 0 0 1 0] ] - [rhoE[1 -8 1 0 0 1 0] ]

    From function checkMethod(const fvMatrix<Type>&, const fvMatrix<Type>&)
    in file /opt/openfoam30/src/finiteVolume/lnInclude/fvMatrix.C at line 1295.

FOAM aborting
The pisoFoam based solver will run succesfully with every term except "- fvm::laplacian(DE, rhoE)".

I cannot see where I am going wrong, is it something to do with the difference in iterative schemes? I know these equations have been successfully implemented in pisoFoam in the past. Also the fact that the same equations will run when implemented in simpleFoam makes me even more confused!

I understand this is quite a difficult case to jump into, but help would really be appreciated! I have attached both solvers, the working one (thrusterFoam) and the one that does not (pisoThrust).
Attached Files
File Type: gz thrusterFoam11.tar.gz (183.7 KB, 6 views)
File Type: gz pisoThrust4.tar.gz (178.1 KB, 3 views)
Aicharem is offline   Reply With Quote

Reply

Tags
dimensions, electrohydrodynamic, incompatible, pisofoam, simplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How can I get stedaystate solution from pisoFoam? mykkujinu2201 OpenFOAM Running, Solving & CFD 1 December 8, 2015 08:37
MPI error with simpleFoam blaise OpenFOAM Running, Solving & CFD 0 November 7, 2015 15:01
First order in fvSchemes does not seem to work gerritgroot OpenFOAM Running, Solving & CFD 0 September 30, 2015 21:06
simpleFoam parallel solver & Fluent polyhedral mesh Zlatko OpenFOAM Running, Solving & CFD 3 September 26, 2014 07:53
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45


All times are GMT -4. The time now is 17:37.