Dear all.
I'm modeling a dr
Dear all.
I'm modeling a drop of the water drop. First calculations were done using ACE program. Now I have recalculate this problem using OpenFoam/interFoam. Results between ACE and interFoam are different, despite fact that both configurations are similar. All physical units are in mm, ms, g. Configuration of my project: ## /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: 1.4.1 | | \ / A nd | Web: http://www.openfoam.org | | \/ M anipulation | | \*---------------------------------------------------------------------------*/ // FoamX Case Dictionary. FoamFile { version 2.0; format ascii; root "/opt/OpenFOAM/caelinux-1.4.1/run"; case "3d"; instance "constant/polyMesh"; local ""; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // arguments "/opt/OpenFOAM/caelinux-1.4.1/run/3d" off; convertToMeters 1; vertices ( (-7e-2 0 0) ( 0 0 0) ( 0 4e-2 0) (-7e-2 4e-2 0) (-7e-2 0 4e-2) ( 0 0 4e-2) ( 0 4e-2 4e-2) (-7e-2 4e-2 4e-2) ); blocks ( hex (0 1 2 3 4 5 6 7) (35 20 20) simpleGrading (1 1 1) ); edges ( ); patches ( symmetryPlane front ( (0 1 2 3) ) symmetryPlane bottom ( (0 1 5 4) ) wall sciana ( (1 2 6 5) ) patch atm ( (4 5 6 7) (3 7 6 2) (0 3 7 4) ) ); mergePatchPairs ( ); // ************************************************** *********************** // ## ## /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: 1.4.1 | | \ / A nd | Web: http://www.openfoam.org | | \/ M anipulation | | \*---------------------------------------------------------------------------*/ // FoamX Case Dictionary. FoamFile { version 2.0; format ascii; root "/opt/OpenFOAM/caelinux-1.4.1/run"; case "3d"; instance "constant"; local ""; class dictionary; object environmentalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // g g [0 1 -2 0 0 0 0] (9.81e-3 0 0); // ************************************************** *********************** // ## ## /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: 1.4.1 | | \ / A nd | Web: http://www.openfoam.org | | \/ M anipulation | | \*---------------------------------------------------------------------------*/ // FoamX Case Dictionary. FoamFile { version 2.0; format ascii; root "/opt/OpenFOAM/caelinux-1.4.1/run"; case "3d"; instance "constant"; local ""; class dictionary; object transportProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // phase1 { transportModel Newtonian; rho rho [1 -3 0 0 0 0 0] 0.000997; nu nu [0 2 -1 0 0 0 0] 1e-03; } phase2 { transportModel Newtonian; rho rho [1 -3 0 0 0 0 0] 1.1614e-6; nu nu [0 2 -1 0 0 0 0] 1.589e-02; } sigma sigma [1 0 -2 0 0 0 0] 3e-4; // ************************************************** *********************** // ## ## Initial gamma configuration: ## x0=-3e-2 y0=0.0 z0=0.0 r=2e-2 xmin=$x0-$r xmax=$x0+$r ymin=$y0-$r ymax=$y0+$r zmin=$z0-$r zmax=$z0+$r funkySetFields -time 0 -field gamma -keepPatches -expression "(sqr(pos().x-$x0)+sqr(pos().y-$y0) +sqr(pos().z-$z0)) <sqr($r)>= $y0 && pos().z >= $z0 ? 1 : 0" funkySetFields -time 0 -field U -keepPatches -expression "(sqr(pos().x-$x0)+sqr(pos().y-$y0) +sqr(pos().z-$z0)) <sqr($r)>= $y0 && pos().z >= $z0 ? vector(6,0,0) : vector(0,0,0)" ## Boundary conditions: gamma ## front { type symmetryPlane; } bottom { type symmetryPlane; } sciana { type constantGammaContactAngle; theta0 90; value uniform 0; } atm { type zeroGradient; } ## U ## front { type symmetryPlane; } bottom { type symmetryPlane; } sciana { type fixedValue; value uniform (0 0 0); } atm { type zeroGradient; } ## pd ## front { type symmetryPlane; } bottom { type symmetryPlane; } sciana { type zeroGradient; } atm { type fixedValue; value uniform 0; } ## When I introduce ConvertToMeters 0.1, and rescale the radius of sphere by factor 10 I obtain correct physical behaviour. Maybe some physical parameters is missing ? Thank You. Jarek |
All times are GMT -4. The time now is 20:10. |