search

nickabattista / IB2d

An easy to use immersed boundary method in 2D, with full implementations in MATLAB and Python that contains over 75 built-in examples, including multiple options for fiber-structure models and advection-diffusion, Boussinesq approximations, and/or artificial forcing.
GNU General Public License v3.0
161 stars 88 forks source link

Naturally oscillating 2D-cylinder #2

Open Haider-BA opened 8 years ago

Haider-BA commented 8 years ago

Kindly include an example of naturally oscillating (elastically mounted / VIV) 2D cylinder.

Haider-BA commented 8 years ago

Did you try to implement the above problem. For reference, please see the code lines from 153 - 168 in the attached link. https://github.com/unicfdlab/TrainingTracks/blob/master/OpenFOAM/simpleFsi-OF3.0.0/src/basicFsiFunctionObject/basicFsi.C

nickabattista commented 8 years ago

Sorry I have not had a chance to try to implement it, yet. I have been adding more functionality to the code. I will try to get to that. Do you need it ASAP?

Haider-BA commented 8 years ago

Yes, I need it ASAP. If you could implement it.

nickabattista commented 8 years ago

I am implementing it now. I'll keep you up-to-date.

Haider-BA commented 8 years ago

Thank you. I am waiting anxiously.

nickabattista commented 7 years ago

I have an implementation ready with elastic (non-damped) springs attaching the cylinder to the channel. I am adding in functionality for damped springs now, but it may take me a day or two.

nickabattista commented 7 years ago

I have an example of the VIV cylinder with damped springs now working. Please let me know what you think. Thanks:

-> IB2d/matIB2d/Examples/Example_VIV_Cylinder/Single_VIV_Cylinder_Damped_Springs_256

nickabattista commented 7 years ago

Thoughts?

Haider-BA commented 7 years ago

I tried but cold not understand whether the circular cylinder is moving due to aerodynamic forces or still? Reporting the position (x and y) w.r.t time will be helpful to get the idea. It looks different implementation than the literature.

nickabattista commented 7 years ago

You can visualize the data that the simulation produces in either VisIt (https://wci.llnl.gov/simulation/computer-codes/visit/ ) or ParaView (http://www.paraview.org/download/ ), and look at any of the Eulerian data (velocity vectors, velocity magnitudes, vorticity, pressure, etc., along with the position of the Lagrangian Structure. You can do the explicit analysis you described using the data_analysis software within IB2d. There is an example for flow in a channel here: IB2d/data_analysis/analysis_in_matlab/Example_For_Data_Analysis/

The circular cylinder is moving solely based on the fluid dynamic forces acting on it alone, and then it is recoiling back due to the elastic deformations of the damped springs tethering the cylinder to the walls of the channel.

Haider-BA commented 7 years ago

Now I can see the movement. But the boundary looks deformed (elastic), I think it should be rigid. For data analysis same file (IB2d/data_analysis/analysis_in_matlab/Example_For_Data_Analysis/) will be used or it needs modification?

Haider-BA commented 7 years ago

How to test VIV with different Mass ratios, frequency, etc.., where to input?