Poisseuile Flow Test
Hello fellow Foamers,
I'm working to verify interFoam for later two phase work with a classic Poisseuile flow test case. I've followed the suggestions on other posts to use pressureInletVelocity as the inlet boundary condition. Since the only results I've gotten show a uniform zero velocity field I suspect my B.C.'s are still not correct. Please see the contents of my boundary files and my control file below. Any suggestions sill be greatly appreciated. Best Regards, Stephen FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { allWalls { type fixedValue; value uniform (0 0 0); } movingWall { type fixedValue; value uniform (0 0 0); } outlet { type pressureInletVelocity; value uniform (0 0 0); } inlet { type pressureInletVelocity; value uniform (0 0 0); } defaultFaces { type empty; } // ************************************************** *********************** // FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { allWalls { type zeroGradient; } movingWall { type zeroGradient; } inlet { type fixedValue; value uniform 3000; } outlet { type fixedValue; value uniform 0; } defaultFaces { type empty; } } // ************************************************** *********************** // FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application interFoam; startFrom startTime; startTime 0; stopAt endTime; endTime 5.0; deltaT 0.001; writeControl adjustableRunTime; writeInterval 0.05; purgeWrite 0; writeFormat ascii; writePrecision 8; writeCompression uncompressed; timeFormat general; timePrecision 8; runTimeModifiable yes; adjustTimeStep yes; maxCo 0.5; maxDeltaT 1; // ************************************************** *********************** // |
All times are GMT -4. The time now is 09:51. |