ABL case files

group abl

Contains: user specified routines for the ABL case.

  • uservp : variable properties

  • userf : local acceleration term for fluid

  • userq : local source term for scalars

  • userbc : boundary conditions

  • useric : initial conditions

  • userchk : general purpose routine for checking errors etc.

  • userqtl : thermal divergence for lowMach number flows

  • usrdat : modify element vertices

  • usrdat2 : modify mesh coordinates

  • usrdat3 : general purpose routine for initialization

Functions

subroutine uservp(ix, iy, iz, ieg)

Set variable properties, does not call any subroutines / functions.

Parameters

  • ix: x-index

  • iy: y-index

  • iz: z-index

  • ieg: element index

subroutine userf(ix, iy, iz, eg)
subroutine userq(ix, iy, iz, eg)
subroutine userbc(ix, iy, iz, iside, eg)
subroutine useric(ix, iy, iz, eg)

Set up initial conditions.

subroutine userchk()

Compute the turbulent stress tensors and write statistics.

subroutine userqtl
subroutine usrdat()
subroutine usrdat2()
subroutine usrdat3()
subroutine eddy_visc(e)

Compute eddy viscosity using dynamic smagorinsky model.

subroutine set_grid_spacing(dg2)

Compute D^2, the grid spacing used in the DS sgs model.

subroutine set_obj

Define objects for surface integrals.

subroutine set_forcing(f_new, u, idir)

Compute inflow / outflow conditions a.k.a. driving force. Ubar=1.

subroutine coriolis(ix, iy, iz, ieg)

Compute Coriolis acceleration

\[ F_x, F_z = - f\hat{e}_y \times \vec{u} \]
See https://lists.mcs.anl.gov/pipermail/nek5000-users/2014-May/002798.html.

Parameters

  • [in] ix[in] iy[in] iz: GLL point index

  • [in] ieg: global element number