Cusomisation inertia
 

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 ? :confused:

Thanks a lot

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.

any answer ?
 

Hi,

For the moment no answer I continue to work on customisation but for the moment it isn't working.....! :mad:
Who can help us ?

Luc

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

thanks
 

Hello,

Thanks a lot I will try

:)