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

Dynamic mesh within a dynamic mesh

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

Like Tree1Likes
  • 1 Post By louisgag

Reply
 
LinkBack Thread Tools Display Modes
Old   October 8, 2013, 15:33
Default Dynamic mesh within a dynamic mesh
  #1
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 223
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hello all,

I've sucessfully used GGI and AMI in the past to simulate a rotating part within a fixed domain.

Now, I want to know if it is possible to have a rotating mesh embed within another rotating mesh. Has anyone achieved this?

Thanks


-Louis
louisgag is offline   Reply With Quote

Old   October 8, 2013, 21:24
Default
  #2
Member
 
Join Date: Apr 2009
Posts: 36
Rep Power: 9
minger is on a distinguished road
Wow, I was actually getting ready to post the same question. I have looked through the mixer2d tutorials in pimpleDyMFoam and think I have a decent grasp of what's going on there. However, I am hoping to be able to "nest" a rotating zone within another.
minger is offline   Reply With Quote

Old   October 29, 2013, 09:30
Default
  #3
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 223
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
For those wondering if it works, yes it does. I had to do a tiny C++ modification to ensure that the origin of the internal AMI zone rotates with the external one and voilà.
-Louis
louisgag is offline   Reply With Quote

Old   October 31, 2013, 21:12
Default
  #4
Member
 
Join Date: Apr 2009
Posts: 36
Rep Power: 9
minger is on a distinguished road
Would you mind attaching a case zipped up?
minger is offline   Reply With Quote

Old   November 6, 2013, 17:25
Default
  #5
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 223
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
The project is confidential so I can't put all the data online, but it should be fairly easy to replicate.

In my case I modified oscillatingRotatingMotion.C as such,
Code:
Foam::solidBodyMotionFunctions::oscillatingRotatingMotion::
transformation() const
{
    scalar t = time_.value();

    vector eulerAngles = amplitude_*sin(omega_*t);

    // Convert the rotational motion from deg to rad
    eulerAngles *= pi/180.0;

    scalar omegaBlop_ =  73.304; // outter AMI zone info
    scalar originBlop_ =  3.048;
    scalar radBlop_ =  0.6096;
    scalar omgOffsetBlop = pi/2; // offset from the right hand axis when simulation starts

    vector blop_ = vector (originBlop_+radBlop_*cos(omegaBlop_*t + omgOffsetBlop),
            originBlop_+radBlop_*sin(omegaBlop_*t + omgOffsetBlop), 0);

//    Info<< "LOUIS: " << blop_ << endl;

    quaternion R(eulerAngles.x(), eulerAngles.y(), eulerAngles.z());
//    septernion TR(septernion(origin_)*R*septernion(-origin_));
    septernion TR(septernion(blop_)*R*septernion(-blop_));

    Info<< "solidBodyMotionFunctions::oscillatingRotatingMotion::"
        << "transformation(): "
        << "Time = " << t << " transformation: " << TR << endl;
And in the dynamic mesh dictionary I ensure that the outer AMI zone angular velocity matches the modified one in the C++ file...
louisgag is offline   Reply With Quote

Old   February 26, 2014, 07:03
Default
  #6
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 223
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
here is the code for OpenFoam 2.3.x, copy in your user folder and

Code:
wmake libso
Attached Files
File Type: zip dynamicFvMesh.zip (4.9 KB, 18 views)
sylvester likes this.
louisgag is offline   Reply With Quote

Old   February 26, 2014, 07:14
Default
  #7
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 223
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Along with a test case from which I removed the confidential data... You will need to change the foil U boundary condition because it is a one I've coded myself...

http://ecbiz103.inmotionhosting.com/...leEmbedAMI.zip

any feedback is welcome !
louisgag is offline   Reply With Quote

Old   March 28, 2014, 07:34
Default
  #8
New Member
 
Diego Angeli
Join Date: Mar 2009
Posts: 26
Rep Power: 9
diego_angeli is on a distinguished road
Hi Louis,

here Diego from UNIMORE. We don't know each other but I think we will soon

I just downloaded your myOscillatingRotatingMotion and the test case, thanks for sharing and congratulations for the nice work.

If it's ok for you I will try to use them on a similar case as yours, for CROP..

I'll let you know ASAP how are we faring on that.

Regards
diego_angeli is offline   Reply With Quote

Old   April 2, 2014, 13:13
Default
  #9
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 223
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hello Diego,
It will be a pleasure to collaborate with you. Please use this thread to communicate with me so that everyone can follow. Obviously we'll have to avoid confidential project data.
Regards,
-Louis
louisgag is offline   Reply With Quote

Old   April 14, 2014, 07:47
Default
  #10
New Member
 
Diego Angeli
Join Date: Mar 2009
Posts: 26
Rep Power: 9
diego_angeli is on a distinguished road
Hi Louis,
thanks for the support.
So, Andrea and me managed to reproduce a moving mesh analogous to your case, starting from a mesh generated from us with another meshing tool. The workflow is ok on the mesh side.
Now we should work on the boundary condition for U on the foils, then maybe we will come back to you to agree on the case setup.

Best regards
diego_angeli is offline   Reply With Quote

Old   December 11, 2015, 04:23
Default
  #11
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 223
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
The link for the test case I had posted above is broken. Please use the following instead: louisgagnon.com/temporaire/exampleEmbedAMI.zip

-Louis
louisgag 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
Dynamic Mesh "Shadow Wall" thezack FLUENT 0 June 4, 2013 22:09
dynamic mesh for drop interface IndrajitW FLUENT 0 March 30, 2013 09:03
Dynamic Mesh on Pintle type injector. herntan FLUENT 15 January 4, 2012 04:31
pls help. mesh collapsed with dynamic mesh. wlt_1985 FLUENT 1 July 28, 2011 01:53
dynamic mesh on a hexa grid Manoj Kumar FLUENT 0 August 21, 2007 07:41


All times are GMT -4. The time now is 12:52.