Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn image-incorporated immersed boundary method for diffusion equations
120
0
0.0
(
0
)
Added by
Andrew Chiang
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Authors
Andrew C. Chiang
Ask ChatGPT about the research
No Arabic abstract
A novel sharp interface ghost-cell based immersed boundary method has been proposed and its parameters have been optimized against an analytical model in diffusion applications. The proposed embedded constrained moving least-squares (ECMLS) algorithm minimizes the error of the interpolated concentration at the image point of the ghost point by applying a moving least-squares method on all internal nodes, near the ghost image point, and the associated mirrored image points of these internal nodes through the corresponding boundary conditions. Using an analytical model as a reference, the ECMLS algorithm is compared to the constrained moving least-squares (CMLS) algorithm and the staircase model using various grid sizes, interpolation basis functions, weight functions, and the penalty parameter of the constraint. It is found that using ECMLS algorithm in the investigated diffusion application, the incomplete quartic basis function yields the best performance while the quadratic, cubic, and bicubic basis functions also produce results better than the staircase model. It is also found that the linear and bilinear basis functions cannot produce results better than the staircase model in diffusion applications. It is shown that the optimal radius of the region of internal nodes used for interpolation scales with the logarithm of the boundary radius of curvature. It is shown that for the diffusion application, the proposed ECMLS algorithm produces lower errors at the boundary with better numerical stability over a wider range of basis functions, weight functions, boundary radius of curvature, and the penalty parameter than the CMLS algorithm.
rate research
Read More
In this paper, we develop a robust fast method for mobile-immobile variable-order (VO) time-fractional diffusion equations (tFDEs), superiorly handling the cases of small or vanishing lower bound of the VO function. The valid fast approximation of the VO Caputo fractional derivative is obtained using integration by parts and the exponential-sum-approximation method. Compared with the general direct method, the proposed algorithm ($RF$-$L1$ formula) reduces the acting memory from $mathcal{O}(n)$ to $mathcal{O}(log^2 n)$ and computational cost from $mathcal{O}(n^2)$ to $mathcal{O}(n log^2 n)$, respectively, where $n$ is the number of time levels. Then $RF$-$L1$ formula is applied to construct the fast finite difference scheme for the VO tFDEs, which sharp decreases the memory requirement and computational complexity. The error estimate for the proposed scheme is studied only under some assumptions of the VO function, coefficients, and the source term, but without any regularity assumption of the true solutions. Numerical experiments are presented to verify the effectiveness of the proposed method.
An interface/boundary-unfitted eXtended hybridizable discontinuous Galerkin (X-HDG) method of arbitrary order is proposed for linear elasticity interface problems on unfitted meshes with respect to the interface and domain boundary. The method uses piecewise polynomials of degrees $k (>= 1)$ and $k-1$ respectively for the displacement and stress approximations in the interior of elements inside the subdomains separated by the interface, and piecewise polynomials of degree $k$ for the numerical traces of the displacement on the inter-element boundaries inside the subdomains and on the interface/boundary of the domain. Optimal error estimates in $L^2$-norm for the stress and displacement are derived. Finally, numerical experiments confirm the theoretical results and show that the method also applies to the case of crack-tip domain.
We propose an energy-stable parametric finite element method (ES-PFEM) to discretize the motion of a closed curve under surface diffusion with an anisotropic surface energy $gamma(theta)$ -- anisotropic surface diffusion -- in two dimensions, while $theta$ is the angle between the outward unit normal vector and the vertical axis. By introducing a positive definite surface energy (density) matrix $G(theta)$, we present a new and simple variational formulation for the anisotropic surface diffusion and prove that it satisfies area/mass conservation and energy dissipation. The variational problem is discretized in space by the parametric finite element method and area/mass conservation and energy dissipation are established for the semi-discretization. Then the problem is further discretized in time by a (semi-implicit) backward Euler method so that only a linear system is to be solved at each time step for the full-discretization and thus it is efficient. We establish well-posedness of the full-discretization and identify some simple conditions on $gamma(theta)$ such that the full-discretization keeps energy dissipation and thus it is unconditionally energy-stable. Finally the ES-PFEM is applied to simulate solid-state dewetting of thin films with anisotropic surface energies, i.e. the motion of an open curve under anisotropic surface diffusion with proper boundary conditions at the two triple points moving along the horizontal substrate. Numerical results are reported to demonstrate the efficiency and accuracy as well as energy dissipation of the proposed ES-PFEM.
In this work we present an adaptive boundary element method for computing the electromagnetic response of wave interactions in hyperbolic metamaterials. One unique feature of hyperbolic metamaterial is the strongly directional wave in its propagating cone, which induces sharp transition for the solution of the integral equation across the cone boundary when wave starts to decay or grow exponentially. In order to avoid a global refined mesh over the whole boundary, we employ a two-level a posteriori error estimator and an adaptive mesh refinement procedure to resolve the singularity locally for the solution of the integral equation. Such an adaptive procedure allows for the reduction of the degree of freedom significantly for the integral equation solver while achieving desired accuracy for the solution. In addition, to resolve the fast transition of the fundamental solution and its derivatives accurately across the propagation cone boundary, adaptive numerical quadrature rules are applied to evaluate the integrals for the stiff matrices. Finally, in order to formulate the integral equations over the boundary, we also derive the limits of layer potentials and their derivatives in the hyperbolic media when the target points approach the boundary.
This paper develops and analyzes a general iterative framework for solving parameter-dependent and random diffusion problems. It is inspired by the multi-modes method of [7,8] and the ensemble method of [19] and extends those methods into a more general and unified framework. The main idea of the framework is to reformulate the underlying problem into another problem with a parameter-independent diffusion coefficient and a parameter-dependent (and solution-dependent) right-hand side, a fixed-point iteration is then employed to compute the solution of the reformulated problem. The main benefit of the proposed approach is that an efficient direct solver and a block Krylov subspace iterative solver can be used at each iteration, allowing to reuse the $LU$ matrix factorization or to do an efficient matrix-matrix multiplication for all parameters, which in turn results in significant computation saving. Convergence and rates of convergence are established for the iterative method both at the variational continuous level and at the finite element discrete level under some structure conditions. Several strategies for establishing reformulations of parameter-dependent and random diffusion problems are proposed and their computational complexity is analyzed. Several 1-D and 2-D numerical experiments are also provided to demonstrate the efficiency of the proposed iterative method and to validate the theoretical convergence results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
