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

MotionUy and DynamicMesh with OF14

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 9, 2007, 08:04
Default Hello Community, I am tryin
  #1
bfa
Member
 
Björn Fabritius
Join Date: Mar 2009
Location: Freiberg, Germany
Posts: 31
Rep Power: 8
bfa is on a distinguished road
Hello Community,

I am trying to use automatic mesh motion with component laplacian solver and uniform diffusivity (for inner point motion).
I need to apply a non-uniform boundary motion to one of my walls, therefore I read motionUy from the dictionary and calculate a location-dependent motion velocity in y-direction. So far so good: The new motionUy values are written to the time directories as
value non-uniform List<scalar>
But the executed motion is somehow always uniform. i.e. instead of morphing into a sinus kind of shape, the whole boundary moves upwards.
How can that be? Is there anything I have to change in the laplacian motion solver?

Any hints would be welcome, I could as well send my source snippets and pointMotionU/cellMotionU files if neccessary.

Regards
Bjoern
bfa is offline   Reply With Quote

Old   July 9, 2007, 09:35
Default Hello again, I think I found
  #2
bfa
Member
 
Björn Fabritius
Join Date: Mar 2009
Location: Freiberg, Germany
Posts: 31
Rep Power: 8
bfa is on a distinguished road
Hello again,
I think I found a part of the problem:

In my application I read only cellMotionU, not pointMotionU with
volScalarField& cellMotionU =
const_cast<volscalarfield&>(objectRegistry::lookup Object<volscalarfield>("cellMo tionUy"));

then I update the cellMotion with
cellMotionU.boundaryField()[movingBottomID] ==
curMotionVel_ *
Foam::sin(5.0*M_PI*(boundaryPoints.component(vecto r::X)));

with curMotionVel_ being the amplitude of my motion dependant on simulation-time and boundaryPoints being a vectorField containing the face centres of my moving boundary. (As you can see, I want to describe a motion dependent on the x-coordinate).

When I now try to do the same with pointMotionU, namely
pointScalarField& ptMotionU =
const_cast<pointscalarfield&>(objectRegistry::look upObject<pointscalarfield>("po intMotionUy"));

I get loads of compilation errors concerning the datatype pointScalarField.

Any idea how to read the pointMotionUy values and change them similar to the cellMotionU? I thought the solver's volume-point-interpolator would do the trick for me, but I was wrong there.

Regards
Bjoern
bfa is offline   Reply With Quote

Old   July 10, 2007, 06:52
Default Hmm, looks like I'm talking to
  #3
bfa
Member
 
Björn Fabritius
Join Date: Mar 2009
Location: Freiberg, Germany
Posts: 31
Rep Power: 8
bfa is on a distinguished road
Hmm, looks like I'm talking to myself mainly. But nevertheless I made some progress. i.e. I surrounded the problem and now make it lay down its weapons and give up.

I found out that changing the cellMotionU entry in the time directories is not sufficient, since the motionSolver (in my case laplacianFvMotionSolver) only reads the value from $CASE_ROOT/0 but never updates itself.
How can I tell the motionSolver to take the new (calculated) values for cellMotionU every timestep?

Thanks for your help,
Bjoern
bfa 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
How to upgrade OF14 to OF15 ivanyao OpenFOAM Installation 7 February 15, 2009 11:19
How can i run OF14 rshssh passwordless waynezw0618 OpenFOAM Installation 2 July 11, 2008 06:41
DynamicMesh not buils using icc 100 in sgi altix 4000 arkangel OpenFOAM Installation 0 October 24, 2007 09:40
PROBLEMS WITH cfdTools in OF14 leosding OpenFOAM Running, Solving & CFD 1 August 20, 2007 13:05
Some one help!!! Write Cp location for dynamicmesh Endee FLUENT 0 July 11, 2007 05:42


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