CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CD-adapco

Comet - force calculation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 22, 2004, 16:12
Default Comet - force calculation
  #1
Tony Campbell
Guest
 
Posts: n/a
Using Comet

i iwanted to calculate the total force in given direction on a rectangular body floating half draft in water. without any motion, the weight should kinda balance the pressure force. i summed the forces as follows:

fx=fy=fz=0. do k=nsr(ir),ner(ir)

fx=fx+(btr(1,k)+pb(k)*sb(1,k))*sbm

fy=fy+(btr(2,k)+pb(k)*sb(2,k))*sbm

fz=fz+(btr(3,k)+pb(k)*sb(3,k))*sbm end do

example in post.f divides by /sbm to calculate force per unit area! (note i have the if loops taken care of)

the weight was calculated as density*9.81*length*geomsc * beam*geomsc * thickness*geomsc

my problem is the force could not balance the weight. i guess i have problem with calculating either of them. Can anyone point me to what i am doing wrong?

thanks you

  Reply With Quote

Old   October 25, 2004, 03:07
Default Re: Comet - force calculation
  #2
4xF
Guest
 
Posts: n/a
Here you are:

if (ir.eq.2) then Fx=0. Fy=0. Fz=0. do k=nsr(ir),ner(ir)

Fx = Fx + pb(k)*sb(1,k) + btr(1,k)

Fy = Fy + pb(k)*sb(2,k) + btr(2,k)

Fz = Fz + pb(k)*sb(3,k) + btr(3,k)

end do end if

BTW, you can check if the output is the same as the regionwise forces summary you will find at the end of the simulation in the .his2-file.

  Reply With Quote

Old   October 26, 2004, 15:07
Default Re: Comet - force calculation
  #3
Tony Campbell
Guest
 
Posts: n/a
Thanks. that was really helpful!
  Reply With Quote

Old   November 1, 2004, 12:52
Default Re: Comet - force calculation
  #4
Tony Campbell
Guest
 
Posts: n/a
4xF: i am having big swings in the pressure force values. i plotted fx or fy or fz versus time. i dont know what is going on. the forces are not sinusoidal even though we have excitation in waves. i guess it is because of the swings.. how can i get 'smooth' curves? i'd appreciate your help.

by the way i am doing fluid solid interaction with a coupled body motion.

thanks, t.c
  Reply With Quote

Old   November 2, 2004, 03:43
Default Re: Comet - force calculation
  #5
4xF
Guest
 
Posts: n/a
Because a) you do not use enough outer iterations to allow the flow to converge b) your time integration scheme may be explicit and hence prone to oscillatory behaviour in the integration of the movement when the time step used is large
  Reply With Quote

Old   November 9, 2004, 13:33
Default Re: Comet - force calculation
  #6
Tony Campbell
Guest
 
Posts: n/a
Thanks for your willingness to help! you may be right. since i use the velocity from the time integration to apply as B.C it could be possible my time integration can affect this. i am using first order Euler time int. but i average t and (t+1) force, to calculate velocity and average t and (t+1) velocities to calculate displacement but i apply current velocities as BC (kind of Crank-Nicholson). i saw this algo in a paper. i tried the av velocities but same.. do you know of any time int method that has worked well elsewhere or any other advice in this regard? thanks a bunch
  Reply With Quote

Old   November 13, 2004, 10:01
Default Re: Comet - force calculation
  #7
4xF
Guest
 
Posts: n/a
I think I know the paper (resp. the PhD thesis) you are refering to. It is stated in there that the method is a kind of Cranck-Nicholson. The fact is that the integration method is explicit in time and of first order and will be prone to give you oscillations for the forces if the fluctuations in forces due to fluid flow have a period which in about the same as the time step you use for integration. There are many solutions: - reduce the time step (which may not be desirable) - converge the fluid flow (forces should change less from one time-step to another) - smooth out any oscillations by underrelaxing the movement (which is in fact adding additional inertia to the rigid body system).

Have fun...
  Reply With Quote

Old   November 15, 2004, 02:12
Default Re: Comet - force calculation
  #8
Tony Campbell
Guest
 
Posts: n/a
Thanks. i have taken note of your advise and i will try it. thanks again!
  Reply With Quote

Old   December 2, 2004, 17:50
Default amplitude of displacement
  #9
Tony Campbell
Guest
 
Posts: n/a
the heave displacement is sinusoidal but is modulated. the peaks are not the same. i have some peaks about twice the others, and the pattern repeats.for instance with say 2.5seconds wave period, the pattern repeats like every approx 10 seconds. i tried different wave amplitudes and periods but same. also, i tried Newmark time integration scheme as well as other time integration schemes but same result. i expect the peaks to be equal or fairly so. the only thing left for me to try now is different time steps . i have been trying to figure out what is happening without much success..any clues would be appreciated!
  Reply With Quote

Old   December 3, 2004, 05:20
Default Re: amplitude of displacement
  #10
4xF
Guest
 
Posts: n/a
If you have the displacement velocity of the body in water (say V) and you calculation box has the size L, I am pretty sure that the "weird" wave period is something related to L/V, which means that you have some waves autoamplified by your computational domain (which is finite is size) and the prescribed body movement (consider it as it were yourself making waves in a bath).

Take care...
  Reply With Quote

Old   December 7, 2004, 15:53
Default Re: amplitude of displacement
  #11
Tony Campbell
Guest
 
Posts: n/a
Thanks. i think you are right. there may be reflected waves in there. i have an all wall boundary with a wavemaker on one side and no beach. i used much larger cells on the boundaries thinking there will be diffusion but apparently that did not help much. i dont know how to implement a beach or a wave absorbing boundary in Comet. from all indications it looks like that cannot be done easily using this software. is that correct?
  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
UDF for force calculation Sham FLUENT 0 May 27, 2007 22:33
Radial Force calculation sans CFX 0 July 5, 2006 09:07
Symmetry plane and force calculation Roland CFX 7 May 31, 2006 13:34
Force calculation - Urgent Please AP CD-adapco 6 January 17, 2006 10:22
Force calculation in Fluent Jason FLUENT 0 February 15, 2005 17:15


All times are GMT -4. The time now is 05:22.