A new fast algorithm based on the augmented immersed interface method
and a fast Poisson solver is proposed to solve three dimensional elliptic interface
problems with a piecewise constant but discontinuous coefficient. In the new approach, an augmented variable along the interface, often the jump in the normal
derivative along the interface is introduced so that a fast Poisson solver can be utilized. Thus, the solution of the Poisson equation depends on the augmented variable
which should be chosen such that the original flux jump condition is satisfied. The
discretization of the flux jump condition is done by a weighted least squares interpolation using the solution at the grid points, the jump conditions, and the governing
PDEs in a neighborhood of control points on the interface. The interpolation scheme
is the key to the success of the augmented IIM particularly. In this paper, the key
new idea is to select interpolation points along the normal direction in line with the
flux jump condition. Numerical experiments show that the method maintains second order accuracy of the solution and can reduce the CPU time by 20-50%. The
number of the GMRES iterations is independent of the mesh size. 相似文献
The aim of this paper is to present a new idea to construct the nonlinear fractal interpolation function, in which we exploit the Matkowski and the Rakotch fixed point theorems. Our technique is different from the methods presented in the previous literatures. 相似文献
In this paper, we present some alternative definitions of Besov spaces of generalized smoothness, defined via Littlewood–Paley‐type decomposition, involving weak derivatives, polynomials, convolutions and generalized interpolation spaces. 相似文献
In this paper,we obtain some asymptotic behavior results for solutions to the prescribed Gaussian curvature equation.Moreover,we prove that under a con-formal metric in R2,if the total Gaussian curvature is 4π,the conformal area of R2 is finite and the Gaussian curvature is bounded,then R2 is a compact C1,α surface after completion at ∞,for any α ∈(0,1).If the Gaussian curvature has a Holder decay at in-finity,then the completed surface is C2.For radial solutions,the same regularity holds if the Gaussian curvature has a limit at infinity. 相似文献
The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix. The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities, including transcendental ones, in which the discretization process is as simple as that in solving linear problems, and only common two-term connection coefficients are needed. All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method, which does not require numerical integration in the resulting nonlinear discrete system. The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers. The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids, and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids. In addition, Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method, including the initial guess far from real solutions.
In this paper, we present two different approaches for constructing reduced‐order models (ROMs) for the two‐dimensional shallow water equation (SWE). The first one is based on the noncanonical Hamiltonian/Poisson form of the SWE. After integration in time by the fully implicit average vector field method, ROMs are constructed with proper orthogonal decomposition(POD)/discrete empirical interpolation method that preserves the Hamiltonian structure. In the second approach, the SWE as a partial differential equation with quadratic nonlinearity is integrated in time by the linearly implicit Kahan's method, and ROMs are constructed with the tensorial POD that preserves the linear‐quadratic structure of the SWE. We show that in both approaches, the invariants of the SWE such as the energy, enstrophy, mass and circulation are preserved over a long period of time, leading to stable solutions. We conclude by demonstrating the accuracy and the computational efficiency of the reduced solutions by a numerical test problem. 相似文献