In this paper, a positivity -preserving, well-balanced finite volume scheme on a rectangular mesh is designed based on wet -dry front reconstruction to solve the shallow water equations with nonflat bottom topography....
详细信息
In this paper, a positivity -preserving, well-balanced finite volume scheme on a rectangular mesh is designed based on wet -dry front reconstruction to solve the shallow water equations with nonflat bottom topography. The crucial step is a special piecewise linear representation of the bottom. The flat bottom approximation simplifies the reconstruction of the wet -dry front, which becomes a straight line inside the partially flooded cells when the water is at rest. The continuity of the discrete bottom across the cell interface centres avoids the calculation of the source term across the cell boundaries and simplifies the numerical flux evaluation. The lake -at -rest steady state is preserved after a well-balanced source term discretization. The draining time method is used for the positivity -preserving property in partially flooded cells. Numerical experiments demonstrate the robustness of the scheme.
A macro type of biquadratic C(1) finite elements is constructed on rectangle grids. This is a rectangular version of the C(1) Powell-Sabin element, a C(1)-P(2) element on triangular grids. Here, each rectangle of the ...
详细信息
A macro type of biquadratic C(1) finite elements is constructed on rectangle grids. This is a rectangular version of the C(1) Powell-Sabin element, a C(1)-P(2) element on triangular grids. Here, each rectangle of the base grid is refined into four subrectangles. As in the case of the Powell-Sabin element, we have more constraints than the number of degrees of freedom on each macroelement. However, the extra constraints are consistent. It is shown further that the constructed finite element space is the full C(1)-Q(2) space on the grid. It is also shown that the finite element space is a tensor product space of one-dimensional C(1)-P(2) spaces, where the nodal basis is supported on four intervals. The B-spline function of P(2) is supported on three intervals. The Girault-Scott operator is extended to the element. The application and the convergence of the finite element to the biharmonic equation are presented. Numerical tests are provided.
作者:
ELLIS, JADept. of Comput. Sci.
Victoria Univ. BC Canada Abstract Authors References Cited By Keywords Metrics Similar Download Citation Email Print Request Permissions
We show that two-dimensional rectangular grids of large aspect ratio can be embedded into rectangular grids of smaller aspect ratios with small expansion and dilation. In particular, width can be reduced by a factor o...
详细信息
We show that two-dimensional rectangular grids of large aspect ratio can be embedded into rectangular grids of smaller aspect ratios with small expansion and dilation. In particular, width can be reduced by a factor of up to 2 with optimal expansion, i.e., when the host rectangle is the smallest sufficient to contain the guest, and optimal dilation, i.e., 2. A width reduction factor of 3 can be obtained with optimal expansion and dilation 3. In general, any rectangular grid can be embedded into a square grid that is no more than unity larger on the side than the minimum possible, with dilation no more than 3. These results improve on those previously obtained, in which dilation better than 18 could not be guaranteed. They might be applied to more complex grid embedding problems, such as embedding multidimensional grids into hypercubes.
We construct a family of lower-order rectangular conforming mixed finite elements, in any space dimension. In the method, the normal stress is approximated by quadratic polynomials , the shear stress by bilinear polyn...
详细信息
We construct a family of lower-order rectangular conforming mixed finite elements, in any space dimension. In the method, the normal stress is approximated by quadratic polynomials , the shear stress by bilinear polynomials , and the displacement by linear polynomials . The number of total degrees of freedom (dof) per element is 10 plus 4 in 2D, and 21 plus 6 in 3D, while the previous record of least dof for conforming element is 17 plus 4 in 2D, and 72 plus 12 in 3D. The feature of this family of elements is, besides simplicity, that shape function spaces for both stress and displacement are independent of the spatial dimension . As a result of these choices, the theoretical analysis is independent of the spatial dimension as well. The well-posedness condition and the optimal a priori error estimate are proved. Numerical tests show the stability and effectiveness of these new elements.
Let G(n,m) denote a graph with nm vertices arranged in n rows and m greater than or equal to max{n, 2} columns with an edge {u, v} between vertices u and v if they are adjacent horizontally or vertically. The bandwidt...
详细信息
Let G(n,m) denote a graph with nm vertices arranged in n rows and m greater than or equal to max{n, 2} columns with an edge {u, v} between vertices u and v if they are adjacent horizontally or vertically. The bandwidth of G(n,m) is known to equal n. We prove that the number of edges in a bandwidth-achieving linear arrangement f from the vertex set. of G(n,m) onto {1, 2,...: nm} for which \f(u) -f(v)\ = n can range from 2(n -1) + n(m -n) up to n(m -1). The lower bound is attained for a down-diagonal lexicographic linear arrangement;the upper bound is at-rained by a column-row lexicographic linear arrangement. (C) 1997 John Wiley & Sons, Inc.
Dispersion characteristics of the shallow water gravity-inertia waves are discussed for two nonstaggered rectangular grids that use the velocity potential and streamfunction, and the divergence and vorticity, respecti...
详细信息
Dispersion characteristics of the shallow water gravity-inertia waves are discussed for two nonstaggered rectangular grids that use the velocity potential and streamfunction, and the divergence and vorticity, respectively. It is shown that the simplest second-order accuracy schemes on the two grids produce identical frequencies on an infinite rotating plane, and that, therefore, the two grids are equivalent. Stability analysis of the linearized shallow water equations using leapfrog and semi-implicit time-differencing techniques indicates that viable numerical schemes can be designed on these grids.
The recently developed cross-rhombus stencil-based time-space domain finite-difference (FD) method for modeling two-dimensional acoustic equations, controls the temporal and spatial dispersions synchronously and outpe...
详细信息
The recently developed cross-rhombus stencil-based time-space domain finite-difference (FD) method for modeling two-dimensional acoustic equations, controls the temporal and spatial dispersions synchronously and outperforms cross-stencil-based FD (CS-FD). However, it is only applicable to modeling on equally spaced grids and extensive application is hindered. In this work, we extend it further and develop novel cross-rhombus stencil-based FD (CRS-FD) for modeling on arbitrarily rectangular grids. Two kinds of rhombus stencils involving grid points both on and off the axis are developed first to achieve fourth order and sixth order FD accuracies, respectively, for solving the second order temporal derivative on rectangular grids. The plane wave theory is then used to derive the time-space-domain dispersion relations, when the temporal and spatial derivatives are solved by the new rhombus-stencil-based temporal high order FD and the cross-stencil-based spatial high order FD, respectively. The extrapolation stencil is a mixture of rhombus and cross stencils, and the involved FD coefficients are determined by applying the Tayler expansion on the time-space domain dispersion relation. Our new CRS-FD is high-order accurate in both time and space, and applicable to modeling on arbitrarily rectangular grids. Dispersion analysis, stability analysis and modeling examples on rectangular grids show that the CRS-FD is more accurate and stable than the CS-FD. Meanwhile, we develop the variable-length schemes for CRS-FD to further increase efficiency and apply them to extrapolate wavefields in reverse time migration. The results validate that the CRS-FD is more efficient than the CS-FD because much larger time steps can be used while reaching a similar accuracy. The variable-length scheme further reduces the computational time in comparison with the fixed-length scheme.
In this paper, we present a general framework to construct fractal interpolation surfaces (FISs) on rectangular grids. Then we introduce bilinear FISs, which can be defined without any restriction on interpolation poi...
详细信息
In this paper, we present a general framework to construct fractal interpolation surfaces (FISs) on rectangular grids. Then we introduce bilinear FISs, which can be defined without any restriction on interpolation points and vertical scaling factors.
In this paper, a piecewise quadratic finite element method on rectangular grids for H-1 problems is presented. The proposed method can be viewed as a reduced rectangular Morley element. For the source problem, the con...
详细信息
In this paper, a piecewise quadratic finite element method on rectangular grids for H-1 problems is presented. The proposed method can be viewed as a reduced rectangular Morley element. For the source problem, the convergence rate of this scheme is proved to be O(h(2)) in the energy norm on uniform grids over a convex domain. Alower bound of the L-2-norm error is also proved, which makes the capacity of this scheme more clear. For the eigenvalue problem, the computed eigenvalues by this element are shown to be the lower bounds of the exact ones. Some numerical results are presented to verify the theoretical findings.
We construct a first-order formula phi such that all finite models of phi are non-narrow rectangular grids without using any binary relations other than the grid neighborship relations. As a corollary, we prove that a...
详细信息
We construct a first-order formula phi such that all finite models of phi are non-narrow rectangular grids without using any binary relations other than the grid neighborship relations. As a corollary, we prove that a set A subset of N is a spectrum of a formula which has only planar models if numbers n is an element of A can be recognized by a non-deterministic Turing machine (or a one-dimensional cellular automaton) in time t(n) and space s(n), where t(n)s(n) <= n and t(n), s(n) =Omega(log(n)).
暂无评论