|
[Sponsors] |
Modifications work in simpleFoam, not pisoFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 23, 2016, 12:15 |
Modifications work in simpleFoam, not pisoFoam
|
#1 |
New Member
Join Date: Feb 2016
Posts: 10
Rep Power: 10 |
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 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). |
|
Tags |
dimensions, electrohydrodynamic, incompatible, pisofoam, simplefoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |