setWaveForcingCoefficients
set forcing values for the wave part of the flow
Declaration
setWaveForcingCoefficients(Apbar,Ambar,options)
Parameters
Apbar
Ap ‘mean’ value to relax toAmbar
Am ‘mean’ value to relax toMAp
(optional) forcing mask, Ap. 1s at the forced modes, 0s at the unforced modes. Default is MAp = abs(Apbar) > 0MAm
(optional) forcing mask, Am. 1s at the forced modes, 0s at the unforced modes. Default is MAm = abs(Apbar) > 0tauP
(optional) relaxation time (default 0)tauM
(optional) relaxation time (default 0)
Discussion
Forcing takes the following form,
\[\frac{\partial}{\partial t} A_\pm^{klj} = \underbrace{M_{A_\pm}^{klj} \left(\bar{A}_\pm^{klj} - A_\pm^{klj} \right)/ \tau_\pm}_{F_\textrm{force}} + F_\textrm{inertial}^{klj} + F_\textrm{damp}^{klj}\]where \(M_{A_\pm}^{klj}\) are masks (1s and 0s), \(\bar{A}_\pm^{klj}\) are mean amplitudes, and \(\tau_\pm\) are time scales. If the time scale is set to 0, then the mean amplitudes remain fixed for all time. In that limit, the equations can be written as,
\[\frac{\partial}{\partial t} A_\pm^{klj} = \neg M_{A_\pm}^{klj} \left( F_\textrm{inertial}^{klj} + F_\textrm{damp}^{klj} \right)\]