CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Using version 15 for oscillating cylinder

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

Like Tree1Likes
  • 1 Post By mr.sattari

Reply
 
LinkBack Thread Tools Display Modes
Old   November 24, 2008, 22:54
Default Hi former: I used version
  #1
New Member
 
liu yang
Join Date: Mar 2009
Posts: 7
Rep Power: 7
liuyang is on a distinguished road
Hi former:
I used version 1.5 for flow around oscillating cylinder,but it seems that two files:
cellMotionU and pointMotionU.
I used a file pointMotionU like this:
dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
cylinder
{
type oscillatingVelocity;
omega 0.1;
amplitude (0 -0.628318531 0);
value uniform (0 0 0);
}
up
{
type fixedValue;
value uniform (0 0 0);
}
down
{
type fixedValue;
value uniform (0 0 0);
}
inlet
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBackPlanes
{
type empty;
}
}
and the same for cellMotionU. There seems to be something wrong, the error message is:
gradientInternalCoeffs cannot be called for a genericFvPatchField (actual type oscillatingVelocity)
on patch cylinder of field cellMotionU in file "/home/ly/OpenFOAM/ly-1.5/run/os_cyclinder/0/cellMotionU"
You are probably trying to solve for a field with a generic boundary condition.

Then what should the file cellMotionU should look like?
Thanks a lot.
liuyang is offline   Reply With Quote

Old   November 25, 2008, 01:21
Default Hi, I think removing the ce
  #2
Member
 
Mathieu Olivier
Join Date: Mar 2009
Location: Quebec City, Canada
Posts: 76
Rep Power: 7
mathieu is on a distinguished road
Hi,

I think removing the cellMotionU file from the "0" directory will solve the problem.

Mathieu
mathieu is offline   Reply With Quote

Old   November 25, 2008, 03:22
Default Hi: Following your instruct
  #3
New Member
 
liu yang
Join Date: Mar 2009
Posts: 7
Rep Power: 7
liuyang is on a distinguished road
Hi:
Following your instruction,it works now.Thank you very much.

LiuYang
liuyang is offline   Reply With Quote

Old   February 24, 2009, 06:16
Default Hello, for an oscillating s
  #4
New Member
 
patrick riemer
Join Date: Mar 2009
Location: Uni Stuttgart, Germany
Posts: 19
Rep Power: 7
patriem is on a distinguished road
Hello,

for an oscillating sine motion with vel_amp_x=0.226m/s and frequency of 128Hz i used oscillatingVelocity and a pointMotionU with


type oscillatingVelocity;
omega 804.096; //=2*pi*f
amplitude (0.226 0 0);
value uniform (0 0 0);

But the oscillating amplitude seems to be much too high. what do i do wrong? The oscill....C file under fvMotionSolver/pointPatchFields/derived/oscillatingVelocity
is standard with
(p0_ + amplitude_*sin(omega_*t.value()) - p.localPoints())
/t.deltaT().value()
);

Thank you for a hint
patriem is offline   Reply With Quote

Old   February 24, 2009, 12:16
Default Patrick, You need oscillati
  #5
ssp
New Member
 
Sopheak Seng
Join Date: Mar 2009
Posts: 5
Rep Power: 7
ssp is on a distinguished road
Patrick,

You need oscillatingDisplacement and not oscillatingVelocity. Are you using OF 1.5? Have you had any problem with oscillatingFixedValue?

Sopheak
ssp is offline   Reply With Quote

Old   February 24, 2009, 15:16
Default Hello, I donīt think that i
  #6
New Member
 
patrick riemer
Join Date: Mar 2009
Location: Uni Stuttgart, Germany
Posts: 19
Rep Power: 7
patriem is on a distinguished road
Hello,

I donīt think that i MUST use oscillatingDisplacemt. I can also describe the motion by oscillatiingVelocity, with the maximum value of velocity as the amplitude, and omega as 2*pi*f.

I want to use it for pointMotionU, the inputfile for velocitylaplacian (dynamic mesh). The oscillatingFixedValue b.c. is not supported here, so i need oscillatingVelocity.

COuld you tell me the correct entries and meanings for omega and value uniform? By the way, what is p0?
patriem is offline   Reply With Quote

Old   February 25, 2009, 09:36
Default canīt anybody help me? tried i
  #7
New Member
 
patrick riemer
Join Date: Mar 2009
Location: Uni Stuttgart, Germany
Posts: 19
Rep Power: 7
patriem is on a distinguished road
canīt anybody help me? tried it again today, but the meshmotions isnīt the expected. the amplitude is much too high!
patriem is offline   Reply With Quote

Old   February 25, 2009, 10:01
Default Patrick, Try to use oscilla
  #8
ssp
New Member
 
Sopheak Seng
Join Date: Mar 2009
Posts: 5
Rep Power: 7
ssp is on a distinguished road
Patrick,

Try to use oscillatingDisplacemt on the velocity field and you will be suprise. Just change oscillatingVelocity to oscillatingDisplacemt.

Sopheak
ssp is offline   Reply With Quote

Old   February 25, 2009, 10:12
Default like this: type oscillating
  #9
ssp
New Member
 
Sopheak Seng
Join Date: Mar 2009
Posts: 5
Rep Power: 7
ssp is on a distinguished road
like this:

type oscillatingDisplacement;
omega 804.096; //=2*pi*f
amplitude (0.226 0 0);
value uniform (0 0 0);
ssp is offline   Reply With Quote

Old   February 25, 2009, 11:21
Default Thank you very much, i am one
  #10
New Member
 
patrick riemer
Join Date: Mar 2009
Location: Uni Stuttgart, Germany
Posts: 19
Rep Power: 7
patriem is on a distinguished road
Thank you very much, i am one step further right now!

I could change the sin to a cos-motion by editing the oscillatingDisplacement...C file. And it still works.

Patrick
patriem is offline   Reply With Quote

Old   February 25, 2009, 11:43
Default I have to mention that there w
  #11
New Member
 
patrick riemer
Join Date: Mar 2009
Location: Uni Stuttgart, Germany
Posts: 19
Rep Power: 7
patriem is on a distinguished road
I have to mention that there were of course changes in the dynamicmeshdict (to displacementlaplacian)
fvsolution and fvschemes necessary. and the file in 0 should be called pointDisplacement.
patriem is offline   Reply With Quote

Old   July 7, 2013, 03:49
Default which changes should we do for (fvsolution and fvschemes)?!?
  #12
New Member
 
Pedram Pakad
Join Date: Jul 2013
Posts: 5
Rep Power: 3
pedrampakzad2006 is on a distinguished road
which changes should have for
(fvsolution and fvschemes)?!?
please help me!!
pedrampakzad2006 is offline   Reply With Quote

Old   July 25, 2013, 18:02
Default
  #13
New Member
 
Gecamp
Join Date: Oct 2010
Posts: 13
Rep Power: 6
desmoge is on a distinguished road
Prolly he refers to the fact that you have to specify a solver for 'cellDisplacement' in fvSolution and a scheme for the diffusivity under 'laplacianSchemes' in fvSchemes.
I personally use 'PCG' for the first and 'Gauss linear corrected' for the second.
Bye
desmoge is offline   Reply With Quote

Old   July 26, 2013, 07:49
Default
  #14
New Member
 
Pedram Pakad
Join Date: Jul 2013
Posts: 5
Rep Power: 3
pedrampakzad2006 is on a distinguished road
Quote:
Originally Posted by desmoge View Post
Prolly he refers to the fact that you have to specify a solver for 'cellDisplacement' in fvSolution and a scheme for the diffusivity under 'laplacianSchemes' in fvSchemes.
I personally use 'PCG' for the first and 'Gauss linear corrected' for the second.
Bye
hi Gecamp,
thanx for your hint.
i have another question??!!
what is your idea about dynamicMeshDict for angular Oscillating displacement.
which method, solver and diffusivity should i use?!?
pedrampakzad2006 is offline   Reply With Quote

Old   July 26, 2013, 10:18
Default
  #15
New Member
 
Gecamp
Join Date: Oct 2010
Posts: 13
Rep Power: 6
desmoge is on a distinguished road
Hi Pedram!
Well, I'm quite new regarding dynamic mesh movement. Here below you can find my setting which has worked so far but more runs need to go before I can develop a critical point of on it

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solver displacementLaplacian;

diffusivity inverseDistance (solidGeometryName);


Now I've got a question for you, do you know how to run with OF the only mesh movement without solving the NS equations? I just need to check the motion quickly and I dont want to solve the whole flow fields.
desmoge is offline   Reply With Quote

Old   July 26, 2013, 11:05
Default
  #16
New Member
 
Pedram Pakad
Join Date: Jul 2013
Posts: 5
Rep Power: 3
pedrampakzad2006 is on a distinguished road
Hi desmoge.
I have done the this and unfortunately in large amplitude it got an error!!!

And about your question...
You can solve your problem and after that, you can see the mesh changes in paraView.
In paraView you can set solid color and change surface to surface with edges.
Now you can run paraview and see the mesh motion.
pedrampakzad2006 is offline   Reply With Quote

Old   August 20, 2013, 03:22
Default oscillate problem
  #17
New Member
 
Pedram Pakad
Join Date: Jul 2013
Posts: 5
Rep Power: 3
pedrampakzad2006 is on a distinguished road
Hi everyone. I want to analyze heat transfer of oscillating cylinder.
When I use inversedistance as diffusivity I didn't get a correct answer for nusselt number. It is same as the fixed cylinder.
What should I do?!?
pedrampakzad2006 is offline   Reply With Quote

Old   November 3, 2013, 15:55
Default
  #18
New Member
 
mohammadreza sattari
Join Date: Oct 2013
Posts: 3
Rep Power: 3
mr.sattari is on a distinguished road
Quote:
Originally Posted by desmoge View Post
Hi Pedram!
Well, I'm quite new regarding dynamic mesh movement. Here below you can find my setting which has worked so far but more runs need to go before I can develop a critical point of on it

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solver displacementLaplacian;

diffusivity inverseDistance (solidGeometryName);


Now I've got a question for you, do you know how to run with OF the only mesh movement without solving the NS equations? I just need to check the motion quickly and I dont want to solve the whole flow fields.
hi Gecamp
In order to only deform the mesh without doing any flow calculations the "moveMesh" utility can be used.
alirezahd likes this.
mr.sattari 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
problem with oscillating cylinder ammi FLUENT 0 December 18, 2006 05:06
Oscillating cylinder MC Main CFD Forum 5 September 6, 2006 22:41
UDF for oscillating cylinder fredD FLUENT 0 July 16, 2004 10:32
flow around oscillating cylinder gs Main CFD Forum 0 May 31, 2004 07:34
need help,rotatary oscillating cylinder Filly FLUENT 3 October 4, 2000 10:11


All times are GMT -4. The time now is 11:57.