!-------------------------------------------------------------------------------
! Code_Saturne version 4.0.4
! --------------------------
! This file is part of Code_Saturne, a general-purpose CFD tool.
!
! Copyright (C) 1998-2015 EDF S.A.
!
! This program is free software; you can redistribute it and/or modify it under
! the terms of the GNU General Public License as published by the Free Software
! Foundation; either version 2 of the License, or (at your option) any later
! version.
!
! This program is distributed in the hope that it will be useful, but WITHOUT
! ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
! FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
! details.
!
! You should have received a copy of the GNU General Public License along with
! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
! Street, Fifth Floor, Boston, MA 02110-1301, USA.
!-------------------------------------------------------------------------------
!===============================================================================
! Purpose:
! -------
!> \file cs_user_physical_properties.f90
!> \brief Definition of physical variable laws.
!>
!> usphyv
!> \brief Definition of physical variable laws.
!>
!> \section Warning
!>
!> It is \b forbidden to modify turbulent viscosity \c visct here
!> (a specific subroutine is dedicated to that: \ref usvist)
!>
!> - icp = 1 must have been specified
!> in \ref usipsu if we wish to define a variable specific heat
!> cpro_cp (otherwise: memory overwrite).
!>
!> - the kivisl field integer key (scalar_diffusivity_id)
!> must have been specified
!> in \ref usipsu if we wish to define a variable viscosity
!> \c viscls.
!>
!>
!> \remarks
!> - This routine is called at the beginning of each time step
!> Thus, AT THE FIRST TIME STEP (non-restart case), the only
!> values initialized before this call are those defined
!> - in the GUI or \ref usipsu (cs_user_parameters.f90)
!> - density (initialized at \c ro0)
!> - viscosity (initialized at \c viscl0)
!> - in the GUI or \ref cs_user_initialization
!> - calculation variables (initialized at 0 by defaut
!> or to the value given in the GUI or in \ref cs_user_initialization)
!>
!> - We may define here variation laws for cell properties, for:
!> - density: rom kg/m3
!> - density at boundary faces: romb kg/m3)
!> - molecular viscosity: cpro_viscl kg/(m s)
!> - specific heat: cpro_cp J/(kg degrees)
!> - diffusivities associated with scalars: cpro_vscalt kg/(m s)
!>
!> \b Warning: if the scalar is the temperature, cpro_vscalt corresponds
!> to its conductivity (Lambda) in W/(m K)
!>
!>
!> The types of boundary faces at the previous time step are available
!> (except at the first time step, where arrays \c itypfb and \c itrifb have
!> not been initialized yet)
!>
!> It is recommended to keep only the minimum necessary in this file
!> (i.e. remove all unused example code)
!>
!>
!> \section usphyv_cell_id Cells identification
!>
!> Cells may be identified using the \ref getcel subroutine.
!> The syntax of this subroutine is described in the
!> \ref cs_user_boundary_conditions subroutine,
!> but a more thorough description can be found in the user guide.
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
! Arguments
!______________________________________________________________________________.
! mode name role !
!______________________________________________________________________________!
!> \param[in] nvar total number of variables
!> \param[in] nscal total number of scalars
!> \param[in] mbrom indicator of filling of romb array
!> \param[in] dt time step (per cell)
!_______________________________________________________________________________
!===============================================================================
!===============================================================================
! Purpose:
! -------
!> usvima
!> \brief User subroutine dedicated the use of ALE
!> (Arbitrary Lagrangian Eulerian Method): fills mesh viscosity arrays.
!>
!> This subroutine is called at the beginning of each time step.
!>
!> Here one can modify mesh viscosity value to prevent cells and nodes
!> from huge displacements in awkward areas, such as boundary layer for example.
!>
!> IF variable IORTVM = 0, mesh viscosity modeling is isotropic therefore VISCMX
!> array only needs to be filled.
!> IF variable IORTVM = 1, mesh viscosity modeling is orthotropic therefore
!> all arrays VISCMX, VISCMY and VISCMZ need to be filled.
!>
!> Note that VISCMX, VISCMY and VISCMZ arrays are initialized at the first time step
!> to the value of 1.
!>
!> \section usvima_cell_id Cells identification
!>
!> Cells may be identified using the \ref getcel subroutine.
!> The syntax of this subroutine is described in the
!> \ref cs_user_boundary_conditions subroutine,
!> but a more thorough description can be found in the user guide.
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
! Arguments
!______________________________________________________________________________.
! mode name role !
!______________________________________________________________________________!
!> \param[in] nvar total number of variables
!> \param[in] nscal total number of scalars
!> \param[in] dt time step (per cell)
!> \param[out] viscmx mesh viscosity in X direction
!> \param[out] viscmy mesh viscosity in Y direction
!> \param[out] viscmz mesh viscosity in Z direction
!_______________________________________________________________________________
subroutine usvima &
( nvar , nscal , &
dt , &
viscmx , viscmy , viscmz )
!===============================================================================
!===============================================================================
! Module files
!===============================================================================
use paramx
use dimens, only: ndimfb
use pointe
use numvar
use optcal
use cstphy
use entsor
use parall
use period
use albase
use field
use mesh
!===============================================================================
implicit none
! Arguments
integer nvar , nscal
double precision dt(ncelet)
double precision viscmx(ncelet), viscmy(ncelet), viscmz(ncelet)
! Local variables
integer iel, ifac
integer ilelt, nlelt
double precision rad, xr2, xcen, ycen
double precision, allocatable, dimension(:) :: xb, yb, xb0, yb0
double precision dxmax, dymax
double precision deltax, deltay
integer, allocatable, dimension(:) :: lstelt
!===============================================================================
allocate(lstelt(nfabor))
!===============================================================================
rad = (0.5)**2
dxmax = 0.0
dymax = 0.0
call getfbr('intF', nlelt, lstelt)
!==========
write(nfecra,*) 'number of bondary faces = ', nlelt
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
xb(ilelt) = cdgfbo(1,ifac)
yb(ilelt) = cdgfbo(2,ifac)
if (ntcabs.eq.0) then
xb0(ilelt) = xb(ilelt)
yb0(ilelt) = yb(ilelt)
endif
deltax = abs(xb(ilelt)-xb0(ilelt))
deltay = abs(yb(ilelt)-yb0(ilelt))
if (deltax.gt.dxmax) dxmax = deltax
if (deltay.gt.dymax) dymax = deltay
enddo
xcen = 5.025d0+dxmax
ycen = 1.000d0+dymax
do iel = 1, ncel
viscmz(iel) = 1.d10
xr2 = (xyzcen(1,iel)-xcen)**2 + (xyzcen(2,iel)-ycen)**2
if (xr2.lt.rad) then
viscmx(iel) = 1.d2
viscmy(iel) = 1.d2
else
viscmx(iel) = 0.01/xr2
viscmy(iel) = 0.01/xr2
endif
enddo
write(nfecra,1000)
!----
! Formats
!----
1000 format(/, &
' ** RESOLUTION DE LA VISCOSITY DE MAILLAGE' ,/,&
' ------------------------------------' ,/)
!----
! End
!----
return
end subroutine usvima