## Abstract

We present a technique for numerically solving convection-diffusion problems in domains Ω with curved boundary. The technique consists in approximating the domain Ω by polyhedral subdomains D_{h} where a finite element method is used to solve for the approximate solution. The approximation is then suitably extended to the remaining part of the domain Ω. This approach allows for the use of only polyhedral elements; there is no need of fitting the boundary in order to obtain an accurate approximation of the solution. To achieve this, the boundary condition on the border of Ω is transferred to the border of D_{h} by using simple line integrals. We apply this technique to the hybridizable discontinuous Galerkin method and provide extensive numerical experiments showing that, whenever the distance of D_{h} to ∂ Ω is of order of the meshsize h, the convergence properties of the resulting method are the same as those for the case in which Ω = D_{h}. We also show numerical evidence indicating that the ratio of the L^{2} (Ω) norm of the error in the scalar variable computed with d > 0 to that of that computed with d = 0 remains constant (and fairly close to one), whenever the distance d is proportional to min {h, Pe^{-1}}/(k + 1)^{2}, where Pe is the so-called Péclet number.

Original language | English (US) |
---|---|

Pages (from-to) | 512-543 |

Number of pages | 32 |

Journal | Journal of Scientific Computing |

Volume | 59 |

Issue number | 2 |

DOIs | |

State | Published - May 2014 |

### Bibliographical note

Funding Information:Partially supported by the National Science Foundation (Grant DMS-1115331) and by the Minnesota Supercomputing Institute, University of Minnesota.

## Keywords

- Convection-diffusion
- Curved domains
- Discontinuous Galerkin methods
- Immersed boundary methods