Non UBC
DSpace
Jahandari, Hormoz
2019-03-25T05:06:08Z
2018-05-31T13:31
Discontinuous (jumping) coefficients often appear in modelling problems where the computational domain represents inhomogeneous media. An example of this is the geophysical electromagnetic (EM) modelling problem, where these jumps occur at interfaces that separate regions with different conductivities. These interfaces, along with other problem features, such as singular EM sources and pointwise solution observations, motivate mesh refinement to achieve good accuracy. The goal of this study is to investigate the combined application of h- and r-refinement to reduce numerical error in the modelling of EM data. For simplicity, aspects of this hr-adaptivity are explored in 1D. The steady-state diffusion and Helmholtz equations (which are commonly solved for the EM scalar and vector potentials, respectively) constitute the physical PDEs (PPDEs) here, while the r-refinement is based on an equidistribution principle. The PPDEs and the mesh PDE are solved alternately in an iterative manner to reduce an error estimate to a desired level. At each iteration, the old mesh is h-refined, the error estimate and its corresponding monitor function are updated and the r-refinement is performed. Various finite-element (FE) a posteriori global and local error estimates were examined: while FE residual-based estimates were cheaper to compute, hierarchical error estimates were found to be better indicators of the true errors. The solutions of adjoint problems of the PPDEs were used to construct local error estimates. These estimates were successfully used for goal-oriented refinement of the FE models.
The talk is based on a joint work with Scott MacLachlan and Ronald Haynes.
https://circle.library.ubc.ca/rest/handle/2429/69209?expand=metadata
25.0
video/mp4
Author affiliation: Memorial University of Newfoundland
Banff (Alta.)
10.14288/1.0377451
eng
Unreviewed
Vancouver : University of British Columbia Library
Banff International Research Station for Mathematical Innovation and Discovery
Attribution-NonCommercial-NoDerivatives 4.0 International
http://creativecommons.org/licenses/by-nc-nd/4.0/
Postdoctoral
BIRS Workshop Lecture Videos (Banff, Alta)
Mathematics
Numerical analysis
Partial differential equations
Combining h- and r-adaptivity for finite-element models with jumping coefficients
Moving Image
http://hdl.handle.net/2429/69209