|
[Sponsors] |
November 18, 2013, 13:41 |
6dof model - floating body
|
#1 |
New Member
Florian
Join Date: Jun 2013
Posts: 8
Rep Power: 13 |
Hi,
I would like to use Fluent to simulate the movement of a floating body. I am currently experimenting with dynamic mesh and a "simple" model (2D): I use a rigid rectangular body of 1.54*0.08m in a box of 7.7*0.34m. For the boundary conditions I use the same parameters as Fluent tutorial "Solving a 2D Box Falling into Water". The only difference are the dimensions and the fact that my body is already immerse (by 0.04m) in the water (water height is 0.18m). I meshed the whole domain with trias and used the following UDF to parameter the 6dof model: #include "udf.h" DEFINE_SDOF_PROPERTIES(sdof_properties, prop, dt, time, dtime) { prop[SDOF_MASS] = 61.49; prop[SDOF_IXX] = 129.6296; prop[SDOF_IYY] = 111.1111; prop[SDOF_IZZ] = 129.6296; prop[SDOF_ZERO_TRANS_X] = TRUE; prop[SDOF_ZERO_ROT_Z] = TRUE; printf ("\n updated 6DOF properties"); } Before using dynamic mesh I started the calculation with the body fixed in order to get hydrostatic pressure. Then I allowed the body to move but after roughly 5 time steps, I get the error: "Warning: non-positive volumes exist. Error: Update-Dynamic-Mesh failed. Negative cell volume detected. Error Object: #f" And when i display the mesh, I can see that the body has moved out of the box... If I stop the calculation after 2 or 3 time steps and display the pressure, I can see that I get huge values. If I calculate the lift force on the body, the force is also bigger than what it should be. If i look at the lift force history, when the body is fixed, I get 610.13984 N but after I allow the body to move it is oscillating with amplitude such as -2468606 N, 1.2394356e+08 N, ... until I get the error message. From what could this instability come from? I really don't understand why it does something like that... Could someone help me solve this problem? Thanks Last edited by flinde; November 20, 2013 at 09:42. |
|
December 9, 2013, 04:40 |
|
#2 |
New Member
Florian
Join Date: Jun 2013
Posts: 8
Rep Power: 13 |
Could someone help me?
Thanks. |
|
December 9, 2013, 04:54 |
|
#3 |
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 12 |
Hi,
I'm kind of new to fluent and deal with the same issue myself at times. Here are a few things you can try (have given me some success): 1. Lower your timestep. Try 0.001 or 0.0001 if you haven't already. 2. Enable remeshing under dynamic mesh. Then press on "use defaults" (make sure local cell is enabled). This is under the remeshing tab. 3. Use relaxation factors. In the tutorial you mentioned, they used relaxation for the body forces (if i'm not misstaking). Try to lower the values for, say, pressure etc and see if that helps. Dont go overboard with lowering these as the calculations will/can take alot longer. Again, I'm new to this and hope someone more experienced can help you. |
|
December 11, 2013, 14:33 |
|
#4 |
New Member
Florian
Join Date: Jun 2013
Posts: 8
Rep Power: 13 |
Hi,
First, thanks for your answer. I actually already tried all the things you suggested me and it did not work. I also tried to use the udf allowing to take water compressibility into account (given in the "Solving a 2D Box Falling into Water" tutorial) and I obtained the stability. However, this is stable only for small time steps (1e-5 s) and if I try to increase the time step, it is unstable again. I also tried to understand from what this instability could come from. I took the fluent tutorial case but instead of letting the box fall, I put water up to the middle height of the box and then I run the 6dof model without water compressibility UDF for a 1e-3s time step. In this case the calculations worked. However, I modified the case by taking a box twice larger (did not change any other geometrical parameters) and kept the same numerical parameters, but this time the calculation was not stable. It would seem that decreasing the distance between the moving body and the box wall causes this instability. What I don't understand is how increasing the moving body size can create an instability on the pressure calculation. Another thing that I also tried with this case is using a CGmotion udf and giving the moving body a constant speed. And with a constant speed, there no pressure instability. Well actually you can see that after the first time step the pressure field is wrong but after a few more time steps, the pressure field is corrected. However, if I use the same udf with a varying speed, each time the velocity changes, the pressure field becomes wrong... What could possibly cause that? Thanks again for your help. |
|
March 26, 2015, 21:49 |
Hi,I got exactly the same problem of instability of pressure
|
#5 |
New Member
Wei Liu
Join Date: Mar 2015
Posts: 5
Rep Power: 11 |
In my simulation case, I set a solid sphere (lighter) immersed in liquid, which is supposed to float. The pressure is true when the mesh motion is turned down while the opposite is a mess. I think I got the same question as you. Have you got an solution to this?
|
|
May 4, 2016, 12:19 |
2d-falling-box udf
|
#6 |
New Member
hamid
Join Date: Apr 2014
Posts: 3
Rep Power: 12 |
Hi everyone
I confused about falling-box tutorial udf. DEFINE_SDOF_PROPERTIES(test_box, prop, dt, time, dtime) { prop[SDOF_MASS] = 666.66; prop[SDOF_IXX] = 129.6296; prop[SDOF_IYY] = 111.1111; prop[SDOF_IZZ] = 129.6296; printf ("\n2d_test_box: Updated 6DOF properties"); } as i know moments of inertia of cubic is ((m.d^2)/6) about all axis. since d=1 in this case then moments of inertia is equal 111.111. then why IZZ= IXX = 129.6296. is there anyone can help me.......... thanks a lot |
|
January 11, 2017, 10:56 |
|
#7 | |
Member
Join Date: Jun 2015
Posts: 46
Rep Power: 11 |
Quote:
|
||
June 29, 2017, 02:55 |
|
#8 |
Senior Member
Join Date: Jun 2014
Location: Taiwan
Posts: 100
Rep Power: 12 |
I wonder that since SDOF_IXX, SDOF_IYY and SDOF_IZZ are constant all the time, will the object continue to rotate all the time? @@
That is to say: will the value of SDOF_IXX,, SDOF_IYY and SDOF_IZZ not change? In this way, will the object not stop rotating? -- Sorry... in 2D case, Izz should be given value. Last edited by SJSW; July 6, 2017 at 23:47. |
|
July 1, 2017, 21:29 |
|
#9 | |
New Member
Join Date: Apr 2017
Posts: 1
Rep Power: 0 |
Quote:
|
||
January 11, 2018, 17:06 |
|
#10 | |
New Member
Nidhi
Join Date: Oct 2017
Posts: 10
Rep Power: 9 |
Quote:
|
||
January 12, 2018, 07:26 |
|
#11 | |
Member
ram
Join Date: Apr 2013
Posts: 37
Rep Power: 13 |
Quote:
|
||
June 11, 2018, 05:08 |
|
#12 |
New Member
Enrique Pérez Heredia
Join Date: May 2018
Location: Madrid, Spain.
Posts: 16
Rep Power: 8 |
Hello, i'm trying to impose a fixed velocity to the body but i am not able.
My UDF is the following DEFINE_CG_MOTION(OW_fall, dt, vel, omega, time, dtime) { vel[1] = -5; } but when the simulation begins the body does not move at any time. Cny suggestions? Thank you.
__________________
Loading signature... |
|
June 11, 2018, 20:09 |
|
#13 |
Member
Join Date: Jun 2015
Posts: 46
Rep Power: 11 |
I used the 'Dynamic Mesh Zones' window in Ansys GUI for inserting the velocity and location of the particle. Please note that it would be the initial velocity and location and both of these parameters will start changing according to the problem definition as the simulation proceeds. Please see the attachment.
|
|
March 6, 2021, 14:32 |
half geometry
|
#14 |
Member
Deutschland
Join Date: Sep 2020
Posts: 69
Rep Power: 6 |
Hey,
I am simulating a ship. Due to symmetry I am only working with half of the geometry. for Dynamic mesh I have following doubts 1. centre of mass - the centre of half or full geometry 2. moment of inertia -moment of inertia of half geometry or (moment of inertia of full geometry)/2 Thanks in advance vava10 |
|
Tags |
6dof, rigid body, udf |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Beginner questions - modelling lift and drag on a towed body | ABF | Main CFD Forum | 15 | October 19, 2014 23:27 |
6DOF with floating body bounce back. | paka | FLUENT | 4 | June 6, 2013 11:32 |
calculating buoyancy of a floating body | raghav777 | CFX | 0 | June 27, 2011 02:57 |
6DOF model with sliding mesh | Laurent | FLUENT | 0 | January 27, 2009 05:00 |
floating point exception in turbulent model | marco | Siemens | 0 | December 19, 2008 11:02 |