CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   FLOW-3D (http://www.cfd-online.com/Forums/flow-3d/)
-   -   Cusomisation inertia (http://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

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

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 06:03.