CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   FLOW-3D (https://www.cfd-online.com/Forums/flow-3d/)
-   -   Cusomisation inertia (https://www.cfd-online.com/Forums/flow-3d/73713-cusomisation-inertia.html)

 Luc March 15, 2010 09:21

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

 ck5285 March 17, 2010 01:38

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.

 Luc March 17, 2010 07:45

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

 MuxaB March 22, 2010 22:27

Quote:
 Originally Posted by Luc (Post 250036) 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,

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

 Luc March 23, 2010 03:37

thanks

Hello,

Thanks a lot I will try

:)

 All times are GMT -4. The time now is 18:23.