CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Slip boundary around circles

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 2, 2008, 15:49
Default Dear all, About the 1st que
  #41
Senior Member
 
ZHOU Bin
Join Date: Mar 2009
Location: Nanjing/Shanghai/Torino, Nanjing/Shanghai/Piemente, China/Italy
Posts: 164
Rep Power: 7
zhoubinwx is on a distinguished road
Send a message via ICQ to zhoubinwx Send a message via MSN to zhoubinwx Send a message via Skype™ to zhoubinwx
Dear all,

About the 1st question my post on Sep 21, that is:

The three point coordinates are (4.41226e-05,5.79779e-05), (4.71079e-05,6.24456e-05), (4.2272e-05,5.92143e-05);
The circle center coordinate is:-(4.469e-05,6.083e-05)-----this is correct, I have compared with gmsh geo file
The diamter of the circle is:5.81612e-06-----however, the diameter should be 5.93 micrometer (I am sure the equation to calculate diameter is right: 2*sqrt(....)).

In order to check why, I plot the following pictures, Pic 1 from Gmsh geometry, Pic 2 with the output information.


Pic 1


Pic 2

Finding:
I need to check if the one point coordinate of three points are the point on the circle?

I could not find them in the file constant/polyMesh/points!!! is the point the middle of the circle points? (I guess)

I have two png figures. I find point 1143 coordinate is: (4.25934e-05 5.87334e-05 0)
while our one point output (the one near point 12) is: (4.41226e-05,5.79779e-05,5e-06), while
point 12 is: (4.469e-05 5.7865e-05 0)

The way I get these three points:

const vectorField & pp = patch().Cf();
vector p0, p1, p2;
if (pp.size() >= 3)
{
p0 = pp[0];
p1 = pp[1];
p2 = pp[2];
}

I could not fully understand this numerical discrepancy. would you please give me any hint?

Thank you for your attention.

Bin
zhoubinwx is offline   Reply With Quote

Old   October 3, 2008, 04:26
Default After test, I know that my way
  #42
Senior Member
 
ZHOU Bin
Join Date: Mar 2009
Location: Nanjing/Shanghai/Torino, Nanjing/Shanghai/Piemente, China/Italy
Posts: 164
Rep Power: 7
zhoubinwx is on a distinguished road
Send a message via ICQ to zhoubinwx Send a message via MSN to zhoubinwx Send a message via Skype™ to zhoubinwx
After test, I know that my way to get the point coordinate is wrong,
const vectorField & pp = patch().Cf();

I must find a way out.
zhoubinwx is offline   Reply With Quote

Old   October 4, 2008, 12:12
Default Today I compared the Gmsh poin
  #43
Senior Member
 
ZHOU Bin
Join Date: Mar 2009
Location: Nanjing/Shanghai/Torino, Nanjing/Shanghai/Piemente, China/Italy
Posts: 164
Rep Power: 7
zhoubinwx is on a distinguished road
Send a message via ICQ to zhoubinwx Send a message via MSN to zhoubinwx Send a message via Skype™ to zhoubinwx
Today I compared the Gmsh points coordinates with the file "points", at last I am sure:

const vectorField & pp = patch().Cf();

could return the middle point of each mesh face.

Then I think since the function Cf() could get the point coordinate, I many use the way it used.

I find "makeFaceCentresAndAreas" in "primitiveMeshFaceCentresAndAreas.C",

Then I use the way it used:
-------------------
vector p0, p1, p2;
const pointField& p=points();
const faceList& fs = faces();
forAll (fs, facei)
{
const labelList& f = fs[facei];
p0=p[f[0]];
p1=p[f[1]];
p2=p[f[2]];
}
--------------------
when I compile, I get:
--------------------
error: 'points' was not declared in this scope
error: 'faces' was not declared in this scope

Since we are already in the "circleSlip" type, I do not want to use
label patchID = mesh.boundaryMesh().findPatchID("circle");
to specify the patch.Otherwise we need many findPatchID("circle1"), findPatchID("circle2")....;

This is not what I want, just like the slip already in OpenFOAM, whenever sb want circleSlip boundary, it could get the points on the circle.

If you have any hints about this, i.e. how to get the point coordinates, I would like to thank you in advance.

I must conquer this even though I may not good at OpenFOAM.

Thank you for attention.

Best regards,

bin
zhoubinwx is offline   Reply With Quote

Old   October 5, 2008, 15:38
Default Could I ask in OpenFOAM, which
  #44
Senior Member
 
ZHOU Bin
Join Date: Mar 2009
Location: Nanjing/Shanghai/Torino, Nanjing/Shanghai/Piemente, China/Italy
Posts: 164
Rep Power: 7
zhoubinwx is on a distinguished road
Send a message via ICQ to zhoubinwx Send a message via MSN to zhoubinwx Send a message via Skype™ to zhoubinwx
Could I ask in OpenFOAM, which plane we are using for the 2-D models? Is it in the front(z=0)? in the back(z=1e-5)? or in the middle plane(5e-6)?

Best regards,

Bin
zhoubinwx is offline   Reply With Quote

Old   October 6, 2008, 06:12
Default It is solved about how to get
  #45
Senior Member
 
ZHOU Bin
Join Date: Mar 2009
Location: Nanjing/Shanghai/Torino, Nanjing/Shanghai/Piemente, China/Italy
Posts: 164
Rep Power: 7
zhoubinwx is on a distinguished road
Send a message via ICQ to zhoubinwx Send a message via MSN to zhoubinwx Send a message via Skype™ to zhoubinwx
It is solved about how to get the diameter for a circle patch.

If anyone needs this, just write to me: zhoubinwx@hotmail.com

let me go on with the boundary condition setting, and hope to finish soon and correctly.

Bin
zhoubinwx is offline   Reply With Quote

Old   October 9, 2008, 12:49
Default Dear foamers, After I imple
  #46
Senior Member
 
ZHOU Bin
Join Date: Mar 2009
Location: Nanjing/Shanghai/Torino, Nanjing/Shanghai/Piemente, China/Italy
Posts: 164
Rep Power: 7
zhoubinwx is on a distinguished road
Send a message via ICQ to zhoubinwx Send a message via MSN to zhoubinwx Send a message via Skype™ to zhoubinwx
Dear foamers,

After I implement my boundary, the simulation goes on for several minutes, but then I get an error, please allow me to post it here, to see if some friends have experience dealing with this.

The error is:
#0 Foam::error::printStack(Foam:stream&) in "/home/zhou/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/zhou/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/zhou/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&) in "/home/zhou/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/simpleFoam"
#5 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&) in "/home/zhou/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/simpleFoam"
#6 main in "/home/zhou/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/simpleFoam"
#7 __libc_start_main in "/lib/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/zhou/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/simpleFoam"

If you know, please let me know.

Really appreciate your attention.

Best regards,

Bin
zhoubinwx is offline   Reply With Quote

Old   June 29, 2011, 17:39
Default volVectorField_Tangential Component
  #47
New Member
 
Kelvin
Join Date: Jan 2011
Posts: 5
Rep Power: 5
Kelvin_mk is on a distinguished road
Dear Zhou Bin,

I try to find tangential component of gradient. I read all of your posts. Did you define tangential components of any vector as a volVectorField. You wrote previous post,

volVectorField U_t=U_n*(n&U)

How can you define it???

This thread is very old maybe you got it so I want to ask you?

Thank you for your helping.

Kelvin_mk is offline   Reply With Quote

Old   March 1, 2013, 06:49
Default
  #48
New Member
 
Mussie
Join Date: Feb 2011
Posts: 7
Rep Power: 5
afeworkimussie is on a distinguished road
Hi,
as we all know all of us learn from mistakes others committed while dealing with CFD. That is why we visit the site to get ideas and to learn from others mistakes and to share our mistakes and ideas. Hence, I would like you to post how you solve your problem as I encounter the same error flag message with s similar act to implement a rampedFixedValue BC.

Cheers!
afeworkimussie is offline   Reply With Quote

Old   March 1, 2013, 06:55
Default
  #49
New Member
 
Mussie
Join Date: Feb 2011
Posts: 7
Rep Power: 5
afeworkimussie is on a distinguished road
Quote:
Originally Posted by zhoubinwx View Post
Hi, in order to know the step to implement simple boundary condition on models, I use ramped boundary condition, on the base of the boundary "oscillatingFixedValue".

However, when I wmake, I get the error:
Making dependency list for source file icoFoamRamped.C
could not open file rampedFixedValueFVPatchFiedls.H for source file icoFoamRamped.C
make: *** No rule to make target `rampedFixedValueFvPatchFields.dep', needed by `Make/linuxGccDPOpt/dependencies'. Stop.

I do not know how could I solve this problem. Is here anyone who could get me out?

Thank you
Hello, Please , inform what was the problem and how did you solve it!

Thanks!
afeworkimussie is offline   Reply With Quote

Old   March 1, 2013, 06:59
Default
  #50
New Member
 
Mussie
Join Date: Feb 2011
Posts: 7
Rep Power: 5
afeworkimussie is on a distinguished road
Quote:
Originally Posted by zhoubinwx View Post
Hi, the above problem is solved.

Now when I do as follows:
--------------------------
cp -r $FOAM_TUTORIALS/icoFoam/cavity ~/OpenFOAM/zhou-1.5/applications/icoFoamRamped/cavityRamped
cd cavityRamped/
cd 0
gedit U &
--------------------------
movingWall
{
type rampedFixValue
refValueLow (0 0 0);
refValueHigh (10 0 0);
startRamp 0.1;
endRamp 0.4;
}
--------------------------
I have changed the movingwall boundary with our modified ramped boundary, then
icoFoamRamped> log &
I get error:
----------------------------
Cannot find 'value' entry on patch movingWall of field U in file "/home/zhou/OpenFOAM/zhou-1.5/applications/icoFoamRamped/cavityRamped/0/U"
which is required to set the values of the generic patch field.
(Actual type rampedFixValue)

Please add the 'value' entry to the write function of the user-defined boundary-condition
or link the boundary-condition into libfoamUtil.so

file: /home/zhou/OpenFOAM/zhou-1.5/applications/icoFoamRamped/cavityRamped/0/U::moving Wall from line 37 to line 41.

From function genericFvPatchField::genericFvPatchField(const fvPatch&, const Field&, const dictionary&)
in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 72.

FOAM exiting
--------------------------

I go deep into the file "rampedFixedValueFvPatchField.C", and could not find anything wrong.

Could I ask for any suggestion about this?

Thank you.

Best regards,

Bin

I think you make a mistake in writing the boundary type you just created as rampedFixValue rather than rampedFixedValue , hopefully that may be one cause from what can be sen from the displayed code.

Cheers!
afeworkimussie is offline   Reply With Quote

Old   March 2, 2013, 06:59
Default
  #51
New Member
 
Mussie
Join Date: Feb 2011
Posts: 7
Rep Power: 5
afeworkimussie is on a distinguished road
Quote:
Originally Posted by zhoubinwx View Post
Hi, in order to know the step to implement simple boundary condition on models, I use ramped boundary condition, on the base of the boundary "oscillatingFixedValue".

However, when I wmake, I get the error:
Making dependency list for source file icoFoamRamped.C
could not open file rampedFixedValueFVPatchFiedls.H for source file icoFoamRamped.C
make: *** No rule to make target `rampedFixedValueFvPatchFields.dep', needed by `Make/linuxGccDPOpt/dependencies'. Stop.

I do not know how could I solve this problem. Is here anyone who could get me out?

Thank you
Hi there,
what was the problem and how did you solve it?
afeworkimussie is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Boundary details, free slip or no slip? Tiago CFX 5 August 20, 2013 22:54
Circle of circles Andore Main CFD Forum 1 February 8, 2008 14:45
Modelling with many circles in a rectangle gdbaldw OpenFOAM Pre-Processing 37 November 27, 2007 04:30
Slip Boundary Condition for Moving Boundary Shukla Main CFD Forum 3 November 11, 2005 15:02
Slip boundary for N.S. Rasmus Main CFD Forum 0 June 6, 2002 08:49


All times are GMT -4. The time now is 17:58.