# Cusomisation inertia

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

 March 15, 2010, 09:21 Cusomisation inertia #1 New Member   H Join Date: Mar 2010 Posts: 6 Rep Power: 7 Hi! I need to develop a customisation function in Flow 3D with fortran? I need to compute fluid inertia for each time step and to have the evolution as function of time ? who can send me an exemple of customisation (.f file) or a title of a book where i can find exemple ? Thanks a lot

 March 17, 2010, 01:38 #2 Senior Member   Join Date: Mar 2009 Posts: 129 Rep Power: 8 Hi Luc I also find difficulty to get any tutorial/example/reference that show user on how to conduct the customization with Flow-3D. Hope others users/expert that have the information can share with us. thanks in advance.

 March 17, 2010, 07:45 any answer ? #3 New Member   H Join Date: Mar 2010 Posts: 6 Rep Power: 7 Hi, For the moment no answer I continue to work on customisation but for the moment it isn't working.....! Who can help us ? Luc

March 22, 2010, 22:27
#4
Senior Member

michael barkhudarov
Join Date: Mar 2009
Posts: 331
Rep Power: 9
Quote:
 Originally Posted by Luc Hi! I need to develop a customisation function in Flow 3D with fortran? I need to compute fluid inertia for each time step and to have the evolution as function of time ? who can send me an exemple of customisation (.f file) or a title of a book where i can find exemple ? Thanks a lot
Hi Luc,

By inertia do you mean fluid mass or momentum? Here is a sample code to do x-direction momentum. You can put this code in the user-customizable routine QSADD:

c initialize before summing up
if (nbl.eq.1) fluid_mom_x=zero
do nn=ijklim2(nbl,1),ijklim2(nbl,2)
ijk=ijk2m1(nn)
if(vf(ijk).lt.em6 .or. icstat(ijk).le.0) cycle ! skip blocked and boundary cells
i=i_str(ijk)
j=j_str(ijk)
k=k_str(ijk)
cell_vol=vf(ijk)*delx(i)*dely(j)*delz(k)/rri(i) ! open cell volume
include '../mijk.f' ! compute neighbor indices
u_centered=(afr(imjk)*u(imjk)+afr(ijk)*u(ijk))/(afr(imjk)+afr(ijk)+em6) ! average velocity to cell center
cell_mom_x=rho(ijk)*u_centered*cell_vol ! cell momentum
if(nmat.eq.1) cell_mom_x=f(ijk)*cell_mom_x ! times fluid fraction for one-fluid problems
fluid_mom_x=fluid_mom_x+cell_mom_x! summation
enddo
c print the result
if(nbl.eq.nblcks) print*,'time=',t,'fluid x-momentum=',fluid_mom_x

 March 23, 2010, 03:37 thanks #5 New Member   H Join Date: Mar 2010 Posts: 6 Rep Power: 7 Hello, Thanks a lot I will try

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post emanuele Main CFD Forum 0 September 23, 2009 14:52 abdulnasr FLUENT 0 July 1, 2009 10:57 Miky1020 CFX 2 June 9, 2009 20:16 Sham FLUENT 0 May 1, 2005 23:57 Ard FLUENT 0 February 20, 2005 22:48

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