Saved in:
Bibliographic Details
Main Authors: Carney, Sean P., Dussinger, Milica, Engquist, Bjorn
Format: Preprint
Published: 2023
Subjects:
Online Access:https://arxiv.org/abs/2308.07563
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866911822283538432
author Carney, Sean P.
Dussinger, Milica
Engquist, Bjorn
author_facet Carney, Sean P.
Dussinger, Milica
Engquist, Bjorn
contents Numerical homogenization of multiscale equations typically requires taking an average of the solution to a microscale problem. Both the boundary conditions and domain size of the microscale problem play an important role in the accuracy of the homogenization procedure. In particular, imposing naive boundary conditions leads to a $\mathcal{O}(ε/η)$ error in the computation, where $ε$ is the characteristic size of the microscopic fluctuations in the heterogeneous media, and $η$ is the size of the microscopic domain. This so-called boundary, or ``cell resonance" error can dominate discretization error and pollute the entire homogenization scheme. There exist several techniques in the literature to reduce the error. Most strategies involve modifying the form of the microscale cell problem. Below we present an alternative procedure based on the observation that the resonance error itself is an oscillatory function of domain size $η$. After rigorously characterizing the oscillatory behavior for one dimensional and quasi-one dimensional microscale domains, we present a novel strategy to reduce the resonance error. Rather than modifying the form of the cell problem, the original problem is solved for a sequence of domain sizes, and the results are averaged against kernels satisfying certain moment conditions and regularity properties. Numerical examples in one and two dimensions illustrate the utility of the approach.
format Preprint
id arxiv_https___arxiv_org_abs_2308_07563
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle On the nature of the boundary resonance error in numerical homogenization and its reduction
Carney, Sean P.
Dussinger, Milica
Engquist, Bjorn
Numerical Analysis
Numerical homogenization of multiscale equations typically requires taking an average of the solution to a microscale problem. Both the boundary conditions and domain size of the microscale problem play an important role in the accuracy of the homogenization procedure. In particular, imposing naive boundary conditions leads to a $\mathcal{O}(ε/η)$ error in the computation, where $ε$ is the characteristic size of the microscopic fluctuations in the heterogeneous media, and $η$ is the size of the microscopic domain. This so-called boundary, or ``cell resonance" error can dominate discretization error and pollute the entire homogenization scheme. There exist several techniques in the literature to reduce the error. Most strategies involve modifying the form of the microscale cell problem. Below we present an alternative procedure based on the observation that the resonance error itself is an oscillatory function of domain size $η$. After rigorously characterizing the oscillatory behavior for one dimensional and quasi-one dimensional microscale domains, we present a novel strategy to reduce the resonance error. Rather than modifying the form of the cell problem, the original problem is solved for a sequence of domain sizes, and the results are averaged against kernels satisfying certain moment conditions and regularity properties. Numerical examples in one and two dimensions illustrate the utility of the approach.
title On the nature of the boundary resonance error in numerical homogenization and its reduction
topic Numerical Analysis
url https://arxiv.org/abs/2308.07563