Holistic discretisation of shear dispersion in a two-dimensional channel

T. MacKenzie, A. J. Roberts

Abstract


Consider the spread of a contaminant along a 2D channel or river. We directly derive the 1D discrete numerical model from the 2D advective and diffusive dynamics for the evolution of the contaminant. The holistic discretisation of the 2D advection-diffusion equation is placed within the purview of centre manifold theory by dividing the physical domain into rectangular 2D elements through introducing artificial insulating boundaries which are later removed. The resulting holistic discretisation is consistent with the 1D Taylor model for shear dispersion in the channel. This new technique accurately models the subgrid scale processes and provides a direct link between the 1D numerical discretisation and the original 2D physical dynamics. Centre manifold theory also systematically incorporates the physical inlet and outlet conditions into the 1D discretisation. This method is straightforwardly extended to nonlinear reaction-diffusion equations and more complex geometries.

Full Text:

PDF BibTeX


DOI: https://doi.org/10.21914/anziamj.v44i0.694



Remember, for most actions you have to record/upload into this online system
and then inform the editor/author via clicking on an email icon or Completion button.
ANZIAM Journal, ISSN 1446-8735, copyright Australian Mathematical Society.