Gas EBU example

Gas EBU example

Local variables to be added

The following local variables need to be defined for the examples in this section:

integer          iel, mode, igg, izone
double precision hinit, coefg(ngazgm)
double precision sommqf, sommqt, sommq, tentm, fmelm

double precision, dimension(:), pointer :: cvar_ygfm, cvar_fm, cvar_scalt

Allocation

Before user initialization, work arrays lstelt must be allocated, like in basic example.

Initialization

The following initialization block needs to be added for the following examples:

! Variables initialization:
!   ONLY done if there is no restart computation

if (isuite.gt.0) return

! Control output

write(nfecra,9001)

call field_get_val_s(ivarfl(isca(iygfm)), cvar_ygfm)
call field_get_val_s(ivarfl(isca(ifm)), cvar_fm)
call field_get_val_s(ivarfl(isca(iscalt)), cvar_scalt)

do igg = 1, ngazgm
  coefg(igg) = zero
enddo

! a. Preliminary calculations

sommqf = zero
sommq  = zero
sommqt = zero

! For multiple inlets
do izone = 1, nozapm
  sommqf = sommqf + qimp(izone)*fment(izone)
  sommqt = sommqt + qimp(izone)*tkent(izone)
  sommq  = sommq  + qimp(izone)
enddo

if (abs(sommq).gt.epzero) then
  fmelm = sommqf / sommq
  tentm = sommqt / sommq
else
  fmelm = zero
  tentm = t0
endif

! Calculation of the Enthalpy of the mean gas mixture
! (unburned - or fresh- gas at mean mixture fraction)
if (ippmod(icoebu).eq.1 .or. ippmod(icoebu).eq.3) then
  coefg(1) = fmelm
  coefg(2) = (1.d0-fmelm)
  coefg(3) = zero
  mode     = -1

  ! Converting the mean boundary conditions into
  ! enthalpy values
  call cothht(mode, ngazg, ngazgm, coefg, npo, npot, th, ehgazg, hinit, tentm)
endif

! b. Initialisation

do iel = 1, ncel

  ! Mass fraction of Unburned Gas

  cvar_ygfm(iel) = 5.d-1

  ! Mean Mixture Fraction

  if ( ippmod(icoebu).eq.2 .or. ippmod(icoebu).eq.3 ) then
    cvar_fm(iel) = fmelm
  endif

  ! Enthalpy

  if ( ippmod(icoebu).eq.1 .or. ippmod(icoebu).eq.3 ) then
    cvar_scalt(iel) = hinit
  endif

enddo

Finalization

At the end of the subroutine, it is recommended to deallocate the work array lstelt, like in basic example.