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

Solution dependent mesh movement

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 22, 2011, 10:50
Default Solution dependent mesh movement
  #1
New Member
 
Join Date: Jun 2011
Location: Rotterdam, the Netherlands
Posts: 9
Rep Power: 6
Van de Leur is on a distinguished road
Dear users,

With dynamic meshing I want to make a boundary that moves depending on the flow velocity over this boundary. The movement of the boundary will then become solution dependent.

To achieve this I would like to use the velocityLaplacianFvMotionSolver.
However, the movement of a boundary in this solver is constant and read from the files 'pointMotionU' and 'cellMotionU'.
Instead of reading the displacement velocity from a file I would like to use the velocities in the domain.

The question is if the current input for velocityLaplacianFvMotionSolver.C:

pointMotionU_
(
IOobject
(
"pointMotionU",
fvMesh_.time().timeName(),
fvMesh_,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
pointMesh::New(fvMesh_)
),

can be replaced by something along the lines of:

pointMotionU_ = U.boundaryField()[patchi]

So far I have not been successful, does anyone have a suggestion?

Kind regards,

Kevin
Van de Leur is offline   Reply With Quote

Old   June 25, 2011, 14:41
Default
  #2
Senior Member
 
Philippose Rajan
Join Date: Mar 2009
Location: Germany
Posts: 530
Rep Power: 16
philippose will become famous soon enough
Hello Kevin,

A Good Evening to you :-)!

I had looked into this about two years ago, and adapted some of the common solvers (turbFoam and simpleFoam) to handle solution based motion of boundaries.

You need to search this forum with the following phrase "turbulent solver rigid body dynamics", and you will find a topic with three pages worth of discussions, and multiple downloads for solvers, and examples.

Since I have not looked into it for a long time now, you might have to make some small modifications to get it running on OpenFOAM-1.6-ext (Or maybe OpenFOAM-1.7.x)....

If you have any questions, I could try and help you out :-)!

Have a great weekend ahead!

Philippose
philippose is offline   Reply With Quote

Old   June 28, 2011, 15:30
Default
  #3
New Member
 
Join Date: Jun 2011
Location: Rotterdam, the Netherlands
Posts: 9
Rep Power: 6
Van de Leur is on a distinguished road
Hey Philippose,

Thanks a lot for the reference, very useful.
I'll let you know if I have some questions, or nice results.

Regards,

Kevin
Van de Leur is offline   Reply With Quote

Old   September 27, 2011, 15:26
Default
  #4
New Member
 
Join Date: Jun 2011
Location: Rotterdam, the Netherlands
Posts: 9
Rep Power: 6
Van de Leur is on a distinguished road
Philippose,

I've been following your advice and it helps me a lot.
Regarding the code you've written I do not fully understand the following line from the file solveForceBalance.H:

scalar calcVelocityX = calcVelocity * (aVector & dirVector1);

If I'm not mistaken aVector gives the velocity after multiplication with calcVelocity and dirVector1 gives the direction of this velocity. Mu question is, why is the directional vector also multiplied with calcVelocity?

Regards,

Kevin
Van de Leur is offline   Reply With Quote

Old   January 13, 2014, 08:00
Default
  #5
New Member
 
Detian Liu
Join Date: Mar 2013
Location: Beijing, P.R.China
Posts: 20
Rep Power: 4
Detian Liu is on a distinguished road
Dear Philippose and Kevin,
I am very interesting in your great work, and I want to implement the functionality in which dynamic mesh motion is based on the velocity solution, that means if U is less than U0, the mesh will become a part of boundary, if U is large than U0, the mesh will become a part of computational demain, whether I can use the Foam you mentioned before called "turbForceFoam" to implement this functionality? Or is there some inbuilt solvers in OpenFOAM2.2.2 can achieve this goal?
Thank you in advance for the help!
Best regards,

detian
2014.1.13
Detian Liu is offline   Reply With Quote

Reply

Tags
dynamic mesh

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
grid dependancy gueynard a. Main CFD Forum 19 June 27, 2014 21:22
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 06:42
Radiation and miscellaneous enhancements vtk_fan OpenFOAM Running, Solving & CFD 6 February 18, 2008 00:49
how to converge the solution for fine mesh kathiravan CD-adapco 5 August 11, 2006 01:30
Gambit problems Althea FLUENT 21 February 6, 2001 08:05


All times are GMT -4. The time now is 08:24.