Source author record

Dongwoo Sheen

Dongwoo Sheen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

$P_1$--Nonconforming Quadrilateral Finite Element Space with Periodic Boundary Conditions: Part I. Fundamental results on dimensions, bases, solvers, and error analysis

The $P_1$--nonconforming quadrilateral finite element space with periodic boundary condition is investigated. The dimension and basis for the space are characterized with the concept of minimally essential discrete boundary conditions. We show that the situation is totally different based on the parity of the number of discretization on coordinates. Based on the analysis on the space, we propose several numerical schemes for elliptic problems with periodic boundary condition. Some of these numerical schemes are related with solving a linear equation consisting of a non-invertible matrix. By courtesy of the Drazin inverse, the existence of corresponding numerical solutions is guaranteed. The theoretical relation between the numerical solutions is derived, and it is confirmed by numerical results. Finally, the extension to the three dimensional is provided.

preprint2022arXiv

$P_1$--Nonconforming Quadrilateral Finite Element Space with Periodic Boundary Conditions: Part II. Application to the Nonconforming Heterogeneous Multiscale Method

A homogenization approach is one of effective strategies to solve multiscale elliptic problems approximately. The finite element heterogeneous multiscale method (FEHMM) which is based on the finite element makes possible to simulate such process numerically. In this paper we introduce a FEHMM scheme for multiscale elliptic problems based on nonconforming spaces. In particular we use the noconforming element with the periodic boundary condition introduced in the companion paper. Theoretical analysis derives a priori error estimates in the standard Sobolev norms. Several numerical results which confirm our analysis are provided.

preprint2022arXiv

A computational macroscale model for the time fractional poroelasticity problem in fractured and heterogeneous media

In this work, we introduce a time memory formalism in poroelasticity model that couples the pressure and displacement. We assume this multiphysics process occurs in multicontinuum media. The mathematical model contains a coupled system of equations for pressures in each continuum and elasticity equations for displacements of the medium. We assume that the temporal dynamics is governed by fractional derivatives following some works in the literature. We derive an implicit finite difference approximation for time discretization based on the Caputo time fractional derivative. A Discrete Fracture Model (DFM) is used to model fluid flow through fractures and treat the complex network of fractures. We assume different fractional powers in fractures and matrix due to slow and fast dynamics. We develop a coarse grid approximation based on the Generalized Multiscale Finite Element Method (GMsFEM), where we solve local spectral problems for construction of the multiscale basis functions. We present numerical results for the two-dimensional model problems in fractured heterogeneous porous media. We investigate error analysis between reference (fine-scale) solution and multiscale solution with different numbers of multiscale basis functions. The results show that the proposed method can provide good accuracy on a coarse grid.

preprint2022arXiv

Algebraic Multiscale Method for two--dimensional elliptic problems

We introduce an algebraic multiscale method for two--dimensional problems. The method uses the generalized multiscale finite element method based on the quadrilateral nonconforming finite element spaces. Differently from the one--dimensional algebraic multiscale method, we apply the dimension reduction techniques to construct multiscale basis functions. Also moment functions are considered to impose continuity between local basis functions. Some representative numerical results are presented.

preprint2022arXiv

Simple quadrature rules for a nonparametric nonconforming quadrilateral element

We introduce simple quadrature rules for the family of nonparametric nonconforming quadrilateral element with four degrees of freedom. Our quadrature rules are motivated by the work of Meng {\it et al.} \cite{meng2018new}. First, we introduce a family of MVP (Mean Value Property)-preserving four DOFs nonconforming elements on the intermediate reference domain introduced by Meng {\it et al.}. Then we design two--points and three--points quadrature rules on the intermediate reference domain. Under the assumption on equal quadrature weights, the deviation from the quadrilateral center of the Gauss points for the two points and three points rules assumes the same quadratic polynomials with constant terms modified. Thus, the two--points rule and three--points rule are constructed at one stroke. The quadrature rules are asymptotically optimal as the mesh size is sufficiently small. Several numerical experiments are carried out, which show efficiency and convergence properties of the new quadrature rules.

preprint2020arXiv

$P_1$--nonconforming polyhedral finite elements in high dimensions

We consider the lowest--degree nonconforming finite element methods for the approximation of elliptic problems in high dimensions. The $P_1$--nonconforming polyhedral finite element is introduced for any high dimension. Our finite element is simple and cheap as it is based on the triangulation of domains into parallelotes, which are combinatorially equivalent to $d$--dimensional cube, rather than the triangulation of domains into simplices. Our nonconforming element is nonparametric, and on each polytope it contains only linear polynomials, but it is sufficient to give optimal order convergence for second--order elliptic problems.

preprint2016arXiv

Nonconforming finite element method applied to the driven cavity problem

A cheapest stable nonconforming finite element method is presented for solving the incompressible flow in a square cavity without smoothing the corner singularities. The stable cheapest nonconforming finite element pair based on $P_1\times \widetilde{{\mathcal P}_{c}^h}$ on rectangular meshes \cite{stab-cheapest} is employed with a minimal modification for the discontinuous Dirichlet data on the top boundary, where $\widetilde{{\mathcal P}_{c}^h}$ is the finite element space of piecewise constant pressures with the globally one-dimensional checker-board pattern subspace eliminated. The proposed Stokes elements have least number of degrees of freedom compared to those of known stable Stokes elements. Three accuracy indications for our elements are analyzed and numerically verified. Also, various numerous computational results obtained by using our proposed element show excellent accuracy.

preprint2015arXiv

Error Estimates for Approximations of Distributed Order Time Fractional Diffusion with Nonsmooth Data

In this work, we consider the numerical solution of an initial boundary value problem for the distributed order time fractional diffusion equation. The model arises in the mathematical modeling of ultra-slow diffusion processes observed in some physical problems, whose solution decays only logarithmically as the time $t$ tends to infinity. We develop a space semidiscrete scheme based on the standard Galerkin finite element method, and establish error estimates optimal with respect to data regularity in $L^2(D)$ and $H^1(D)$ norms for both smooth and nonsmooth initial data. Further, we propose two fully discrete schemes, based on the Laplace transform and convolution quadrature generated by the backward Euler method, respectively, and provide optimal convergence rates in the $L^2(D)$ norm, which exhibits exponential convergence and first-order convergence in time, respectively. Extensive numerical experiments are provided to verify the error estimates for both smooth and nonsmooth initial data, and to examine the asymptotic behavior of the solution.

preprint2015arXiv

Stable cheapest nonconforming finite elements for the Stokes equations

We introduce two pairs of stable cheapest nonconforming finite element space pairs to approximate the Stokes equations. One pair has each component of its velocity field to be approximated by the $P_1$ nonconforming quadrilateral element while the pressure field is approximated by the piecewise constant function with globally two-dimensional subspaces removed: one removed space is due to the integral mean--zero property and the other space consists of global checker--board patterns. The other pair consists of the velocity space as the $P_1$ nonconforming quadrilateral element enriched by a globally one--dimensional macro bubble function space based on $DSSY$ (Douglas-Santos-Sheen-Ye) nonconforming finite element space; the pressure field is approximated by the piecewise constant function with mean--zero space eliminated. We show that two element pairs satisfy the discrete inf-sup condition uniformly. And we investigate the relationship between them. Several numerical examples are shown to confirm the efficiency and reliability of the proposed methods.

preprint2013arXiv

A class of nonparametric DSSY nonconforming quadrilateral elements

A new class of nonparametric nonconforming quadrilateral finite elements is introduced which has the midpoint continuity and the mean value continuity at the interfaces of elements simultaneously as the rectangular DSSY element [J.Douglas, Jr., J. E. Santos, D. Sheen, and X. Ye. Nonconforming {G}alerkin methods based on quadrilateral elements for second order elliptic problems. ESAIM--Math. Model. Numer. Anal., 33(4):747--770, 1999]. The parametric DSSY element for general quadrilaterals requires five degrees of freedom to have an optimal order of convergence [Z. Cai, J. Douglas, Jr., J. E. Santos, D. Sheen, and X. Ye. Nonconforming quadrilateral finite elements: A correction. Calcolo, 37(4):253--254, 2000], while the new nonparametric DSSY elements require only four degrees of freedom. The design of new elements is based on the decomposition of a bilinear transform into a simple bilinear map followed by a suitable affine map. Numerical results are presented to compare the new elements with the parametric DSSY element.

preprint2009arXiv

F. John's stability conditions vs. A. Carasso's SECB constraint for backward parabolic problems

In order to solve backward parabolic problems F. John [{\it Comm. Pure. Appl. Math.} (1960)] introduced the two constraints "$\|u(T)\|\le M$" and $\|u(0) - g \| \le δ$ where $u(t)$ satisfies the backward heat equation for $t\in(0,T)$ with the initial data $u(0).$ The {\it slow-evolution-from-the-continuation-boundary} (SECB) constraint has been introduced by A. Carasso in [{\it SIAM J. Numer. Anal.} (1994)] to attain continuous dependence on data for backward parabolic problems even at the continuation boundary $t=T$. The additional "SECB constraint" guarantees a significant improvement in stability up to $t=T.$ In this paper we prove that the same type of stability can be obtained by using only two constraints among the three. More precisely, we show that the a priori boundedness condition $\|u(T)\|\le M$ is redundant. This implies that the Carasso's SECB condition can be used to replace the a priori boundedness condition of F. John with an improved stability estimate. Also a new class of regularized solutions is introduced for backward parabolic problems with an SECB constraint. The new regularized solutions are optimally stable and we also provide a constructive scheme to compute. Finally numerical examples are provided.

preprint2009arXiv

Laplace transformation method for the Black-Scholes equation

In this paper we apply the innovative Laplace transformation method introduced by Sheen, Sloan, and Thomée (IMA J. Numer. Anal., 2003) to solve the Black-Scholes equation. The algorithm is of arbitrary high convergence rate and naturally parallelizable. It is shown that the method is very efficient for calculating various options. Existence and uniqueness properties of the Laplace transformed Black-Scholes equation are analyzed. Also a transparent boundary condition associated with the Laplace transformation method is proposed. Several numerical results for various options under various situations confirm the efficiency, convergence and parallelization property of the proposed scheme.

preprint2009arXiv

Turing Instability for a Ratio-Dependent Predator-Prey Model with Diffusion

Ratio-dependent predator-prey models have been increasingly favored by field ecologists where predator-prey interactions have to be taken into account the process of predation search. In this paper we study the conditions of the existence and stability properties of the equilibrium solutions in a reaction-diffusion model in which predator mortality is neither a constant nor an unbounded function, but it is increasing with the predator abundance. We show that analytically at a certain critical value a diffusion driven (Turing type) instability occurs, i.e. the stationary solution stays stable with respect to the kinetic system (the system without diffusion). We also show that the stationary solution becomes unstable with respect to the system with diffusion and that Turing bifurcation takes place: a spatially non-homogenous (non-constant) solution (structure or pattern) arises. A numerical scheme that preserve the positivity of the numerical solutions and the boundedness of prey solution will be presented. Numerical examples are also included.