|
[Sponsors] |
[swak4Foam] a Problem with groovyBC for icoFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 29, 2017, 13:42 |
a Problem with groovyBC for icoFoam
|
#1 |
Senior Member
mohsen kh
Join Date: Jan 2013
Location: Iran
Posts: 125
Rep Power: 14 |
Dear Foamers
I use OF v4.1. I installed the swak4foam. my geometry is a 2D channel with a pulsatile inlet velocity which behaves as below: u (y,t)= u(y,t)=2U(1-(y/D)**2)(1+0.2*sin(t)) where U=2 m/s and D = 1 m. my U file is as follows: --------------------------------------------------------------------------------------------------------------------------- dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { outlet { type zeroGradient; } inlet { type groovyBC; // variables "Uinlet=1 ;D=1;" valueExpression "vector(2*Uinlet*(1-pow((pos().y)/D,2))*(1+0.2*sin(time())))"; } upperWall { type noSlip; } lowerWall { type noSlip; } frontAndBackPlanes { type empty; } } --------------------------------------------------------------------------------------------------------------------------- my P file is: --------------------------------------------------------------------------------------------------------------------------- dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { outlet { type fixedValue; value uniform 0; } inlet { type zeroGradient; } upperWall { type zeroGradient; } lowerWall { type zeroGradient; } frontAndBackPlanes { type empty; } } --------------------------------------------------------------------------------------------------------------------------- my controlDict is: --------------------------------------------------------------------------------------------------------------------------- application icoFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 10; deltaT 0.0001; writeControl timeStep; writeInterval 500; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable true; libs ( "libOpenFOAM.so" "libsimpleSwakFunctionObjects.so" "libswakFunctionObjects.so" "libgroovyBC.so" ); ---------------------------------------------------------------------------------------------------------------------------- following error occurred after runing icoFoam in terminal mohsen@mohsen-SVE14122CVW:~/OpenFOAM/mohsen-4.1/run/dilate100p$ icoFoam /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 4.1 Exec : icoFoam Date : Aug 29 2017 Time : 21:08:56 Host : "mohsen-SVE14122CVW" PID : 3268 Case : /home/mohsen/OpenFOAM/mohsen-4.1/run/dilate100p nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 PISO: Operating solver in PISO mode Reading transportProperties Reading field p Reading field U --> FOAM Warning : From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict) in file groovyBCFvPatchField.C at line 131 No value defined for U on inlet therefore using the internal field next to the patch Reading/calculating face flux field phi Starting time loop Time = 0.0001 Courant Number mean: 0 max: 0 smoothSolver: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 DICPCG: Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0, global = 0, cumulative = 0 DICPCG: Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0, global = 0, cumulative = 0 ExecutionTime = 0.18 s ClockTime = 0 s Time = 0.0002 Courant Number mean: 0 max: 0 smoothSolver: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 DICPCG: Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0, global = 0, cumulative = 0 DICPCG: Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0, global = 0, cumulative = 0 ExecutionTime = 0.19 s ClockTime = 0 s ------------------------------------------------------------------------------------------------------------------ And this happens for all time steps. I do not know where the mistake is. If anyone experienced a similar problem I would be happy to have advices. Kind regards. |
|
August 29, 2017, 16:57 |
|
#2 |
Senior Member
mohsen kh
Join Date: Jan 2013
Location: Iran
Posts: 125
Rep Power: 14 |
I am looking forward to urgent help of some professional foamers .
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] reactingMultiPhaseEulerFoam problems with groovyBC | zanilu70 | OpenFOAM Community Contributions | 4 | December 13, 2016 07:46 |
[swak4Foam] epsilon groovyBC problem | Thom | OpenFOAM Community Contributions | 5 | October 26, 2012 06:16 |
[swak4Foam] Problem with free surface in groovyBC | Hisham | OpenFOAM Community Contributions | 2 | August 4, 2011 09:23 |
natural convection problem for a CHT problem | Se-Hee | CFX | 2 | June 10, 2007 07:29 |
Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 20:13 |