CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[swak4Foam] a Problem with groovyBC for icoFoam

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 29, 2017, 12:42
Default a Problem with groovyBC for icoFoam
  #1
Senior Member
 
mohsen kh
Join Date: Jan 2013
Location: Iran
Posts: 125
Rep Power: 14
m5m5kh is on a distinguished road
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.
m5m5kh is offline   Reply With Quote

Old   August 29, 2017, 15:57
Default
  #2
Senior Member
 
mohsen kh
Join Date: Jan 2013
Location: Iran
Posts: 125
Rep Power: 14
m5m5kh is on a distinguished road
I am looking forward to urgent help of some professional foamers .
m5m5kh is offline   Reply With Quote

Reply


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
[swak4Foam] reactingMultiPhaseEulerFoam problems with groovyBC zanilu70 OpenFOAM Community Contributions 4 December 13, 2016 06:46
[swak4Foam] epsilon groovyBC problem Thom OpenFOAM Community Contributions 5 October 26, 2012 05:16
[swak4Foam] Problem with free surface in groovyBC Hisham OpenFOAM Community Contributions 2 August 4, 2011 08:23
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13


All times are GMT -4. The time now is 01:35.