Error analysis of corner cutting algorithms

Corner cutting algorithms are used in different fields and, in particular, play a relevant role in Computer Aided Geometric Design. Evaluation algorithms such as the de Casteljau algorithm for polynomials and the de Boor–Cox algorithm for B‐splines are examples of corner cutting algorithms. Here backward and forward error analysis of corner cutting algorithms are performed. The running error is also analyzed and as a consequence the general algorithm is modified to include the computation of an error bound. This revised version was published online in June 2006 with corrections to the Cover Date.     

用单层位势求解尖角区域上的Dirichlet外问题

采用Kress变换以及处理第一类奇异核的积分方法,运用Nystrom方法利用单层位势求解尖角区域上的Dirichlet外问题.给出具体的算法和数值例子,通过数值例子可以看出用单层位势求解尖角区域上的Dirichlet外问题与用单双层结合求解所得的结果基本上一致,说明这种方法是有效的和可行的.     

An evolution compressible Stokes system in a polygon

An evolution compressible Stokes system is studied in a polygon. A lowest order corner singularity by the Laplacian operator is sorted out and an increased regularity is established in a suitable Banach space. In the lowest level, for each time the order of the corner singularity of the system is the same as that of the Laplacian. In addition, a singular function for pressure is defined and shows that the corner singularities originated from the corners of polygon may propagate along the directed streamlines.     

Fuzzy handoff algorithms for wireless communication

In order to manage the high call density expected of future cellular systems, microcells must be used. A migration to microcells will increase the number of handoffs, and require faster handoff algorithms – in terms of decision making. In the case of line-of-sight transmission, it is important that the handoff algorithm detects the cell boundary early enough, otherwise this will lead to channel dragging into the new cell subsequently increasing the chance of co-channel interference. In the case of non-line-of-sight transmission, a mobile station on turning a street corner will experience a phenomenon known as the Manhattan corner effect that causes the received signal level to drop by 20–30 dB in 20–30 m. This corner effect problem can lead to a loss of communication if not identified early enough. This paper presents two new handoff techniques using fuzzy logic as possible solutions to microcellular handoff. The first algorithm uses an adaptive fuzzy predictor, while the second uses a fuzzy averaging technique. The results of the simulation show that fuzzy is a viable option for microcellular handoff.     

On an Algebra of a Certain Class of Operators in a Slab Domain in R2

In this paper a Dirichlet problem for the Laplacian in a domain with a corner in R² is treated and a new method of construction of a local parametrix for that Dirichlet problem at a corner point is given. In order to construct a parametrix a certain class of operators in a slab domain is introduced and it is shown that that operator class is an algebra.     

An adaptive pruning algorithm for the discrete L-curve criterion

We describe a robust and adaptive implementation of the L-curve criterion. The algorithm locates the corner of a discrete L-curve which is a log–log plot of corresponding residual norms and solution norms of regularized solutions from a method with a discrete regularization parameter (such as truncated SVD or regularizing CG iterations). Our algorithm needs no predefined parameters, and in order to capture the global features of the curve in an adaptive fashion, we use a sequence of pruned L-curves that correspond to considering the curves at different scales. We compare our new algorithm to existing algorithms and demonstrate its robustness by numerical examples.     

Some solutions of karman''s equation describing the transonic flow past a corner point on a profile with a curvilinear generatrix

Transonic flows in the neighbourhood of a corner point on a profile are investigated in a class of self-similar solutions of Karman's equation. The corner point is formed by the intersection of two smooth curves, the tangents to which make a convex angle. The generatrix, which lies in the subsonic part of the flow is assumed to be curvilinear and to vary according to a power law. Values of the self-similarity index are found for which transonic flows are possible either with a free streamline or with a rarefaction wave.     

Corner asymptotics of the magnetic potential in the eddy‐current model

In this paper, we describe the magnetic potential in the vicinity of a corner of a conducting body embedded in a dielectric medium in a bidimensional setting. We make explicit the corner asymptotic expansion for this potential as the distance to the corner goes to zero. This expansion involves singular functions and singular coefficients. We introduce a method for the calculation of the singular functions near the corner, and we provide two methods to compute the singular coefficients: the method of moments and the method of quasi‐dual singular functions. Estimates for the convergence of both approximate methods are proven. We eventually illustrate the theoretical results with finite element computations. The specific nonstandard feature of this problem lies in the structure of its singular functions: They have the form of series whose first terms are harmonic polynomials, and further terms are genuine nonsmooth functions generated by the piecewise constant zeroth order term of the operator. Copyright © 2013 John Wiley & Sons, Ltd.     

Large-Eddy Simulations of cavitation in a square surface cavity

We report on the development and application of a multiphase approach to the prediction of cavitation induced by high-speed flow over and within a square surface cavity. The approach entails employing a full cavitation model in conjunction with Large-Eddy Simulations in order to capture the initiation and development of bubble formations in turbulent-flow conditions. The incipient formation of the bubble cloud, and the flow processes of vortex shedding and shear-layer oscillations are tracked using the Volume of Fluid method. The validity of the computational approach was assessed by comparisons with experiments on cavitating flow over a hydrofoil. Application to the case of flow over and within a two-dimensional square cavity with cavitation clearly reveal the presence of traveling cavitation at the corner of the cavity trailing edge, and vortex cavitation within the cavity. It is shown that the collapse of cavitation bubbles results in an impact frequency that is higher than the frequency of the shear-layer oscillations. This implies that structural damage due to cavitation is likely to be most severe at the corner formed at the intersection of the cavity’s trailing edge and the flat surface upstream of it.     

Some multilevel decoupled algorithms for a mixed navier-stokes/darcy model

In this work, several multilevel decoupled algorithms are proposed for a mixed Navier-Stokes/Darcy model. These algorithms are based on either successively or parallelly solving two linear subdomain problems after solving a coupled nonlinear coarse grid problem. Error estimates are given to demonstrate the approximation accuracy of the algorithms. Experiments based on both the first order and the second order discretizations are presented to show the effectiveness of the decoupled algorithms.     

Mathematical analysis of absorbing boundary conditions for the wave equation: the corner problem

Our goal in this work is to establish the existence and the uniqueness of a smooth solution to what we call in this paper the corner problem, that is to say, the wave equation together with absorbing conditions at two orthogonal boundaries. First we set the existence of a very smooth solution to this initial boundary value problem. Then we show the decay in time of energies of high order--higher than the order of the boundary conditions. This result shows that the corner problem is strongly well-posed in spaces smaller than in the half-plane case. Finally, specific corner conditions are derived to select the smooth solution among less regular solutions. These conditions are required to derive complete numerical schemes.

     

High Hartmann number flow through a rectangular duct with unequally and finitely conducting walls

Nonsymmetric Hartmann flow through a rectangular duct is investigated for thin duct walls with, generally, unequal but finite conductivities. A high Hartmann number is adopted. Consistent with known phenomena, both Hartmann layers transverse to the applied magnetic field are assumed to be separated from the two side boundary layers by four corner regions plus four inner corner regions. The method of singular perturbations and matched asymptotic expansions is applied to the coupled system. The equations governing the core and Hartmann layers are first partially resolved for leading terms. This is then followed by tackling equations governing one side layer and two adjacent corner regions. The latters' incorporation secures, for the former, only those boundary conditions that are compatible along the transverse walls. Both corner regions are denied access to non-required boundary conditions along the neighbouring side wall by the adjoining inner corner regions. However, the latters' boundary value problems need not be tackled for the acquirement of only dominant terms beyond all four inner corner regions. The complementary side layer and associated corners are accounted for by a non-symmetric reflection principle. Results reveal that a difference between conductivities in the transverse walls together with at least one finitely conducting side wall impart to disturbances within the core and Hartmann layers (i) a nontrivial dependence on the transverse coordinate relative to the magnetic field and flow in addition to the (usual) dependence on the field aligned coordinate, (ii) a dependence on side wall parameters in addition to the dependence on transverse wall parameters. Applications to related situations are considered. These include the case for a perfectly conducting lower wall, a finitely conducting upper wall, and equally and finitely conducting side walls.     

带尖角的障碍声波散射区域的反演

对带尖角的障碍声波散射区域进行了反演,其前提条件是整体场满足奇次Dirichlet边界条件．在用Nystrom方法解正问题的过程中,由于采用等距网格积分给尖角处带来很差的收敛性,这是因为双层位势的积分算子的核在尖角处有Mellin型奇性,不再是紧算子;为此采用梯度网格,数值例子表明该处理方法的有效可靠性．     

For the stationary compressible viscous Navier-Stokes equations with no-slip condition on a convex polygon

Our concern is with existence and regularity of the stationary compressible viscous Navier-Stokes equations with no-slip condition on convex polygonal domains. Note that [u,p]=[0,c], c a constant, is the eigenpair for the singular value λ=1 of the Stokes problem on the convex sector. It is shown that, except the pair [0,c], the leading order of the corner singularities for the nonlinear equations is the same as that of the Stokes problem. We split the leading corner singularity from the solution and show an increased regularity for the remainder. As a consequence the pressure solution changes the sign at the convex corner and its derivatives blow up.     

Isotropic moments over integer lattices

Many modern edge and corner detection algorithms use moment transforms, which convolve images with tensor-valued filters, namely the product of a window function with a monomial. Over continuous domains, one may easily show that such transforms are isotropic. We generalize these continuous results to digital images, that is, to functions over the canonical integer lattice in a finite-dimensional real space. In particular, we first introduce a mathematically well-behaved method for the dilation and rotation of digital images, and then show these operations commute with discrete moment transforms in a manner consistent with the continuous results.     

Reflected brownian motion in a wedge: Semimartingale property

In this paper, the object of study is reflected Brownian motion in a two-dimensional wedge with constant direction of reflection on each side of the wedge. The basic question considered here is “When is this process a semimartingale?”. It is first shown that a related process, defined by specifying the corner of the wedge to be an absorbing state, rather than an instantaneous one, is a semimartingale. Conditions for the existence and uniqueness of the process for which the corner is an instantaneous state were given by Vardhan and Williams (“Brownian motion in a wedge with oblique reflection”, Comm. Pure Appl. Math., to appear). Under these conditions, it is shown that starting away from the corner, the process is a semimartingale if and only if there is a convex combination of the directions of reflection that points into the wedge. This equivalence is also shown to hold starting from the corner, except in one unresolved case for which the wedge angle exceeds π and the directions of reflection are exactly opposed.     

For the inf‐sup condition on mesh‐dependent norms on a nonconvex polygon

It is shown that the inf‐sup condition, called the Babuska–Brezzi condition, is valid for certain mesh‐dependent norms on a nonconvex polygonal domain. A bilinear form that is derived by inserting the corner singularity expansion into the Laplace equation is considered. A mesh‐dependent fractional norm related to the least order of the corner singularity at a concave vertex is considered. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

自适应有限元方法在凹角域线性椭圆方程的应用(续)

作为序列文章自适应有限元方法在凹角域线性椭圆方程的应用的第三篇,在本文我们将给出并详细论证一个重要结论即 |?(u(x)-U(x))|≤Ch(x)|x|^β-2,|x|≥C′h且进一步分析说明在本序列文章的第一部分和地二部分得出方法都是以此为基础作出的。     

Remark on a trace theorem for transmission problems

In this paper we discuss elliptic transmission problems of second order in plane non-smooth domains with non-homogeneous boundary data. It is known that if the boundary data are homogeneous, then the variational solution of the transmission problem admits a representation as a sum of a regular function and certain singular functions, analogous to that encountered in boundary value problems on corner domains. We determine for which domains arbitrary non-homogeneous boundary data can be reduced to the homogeneous ones preserving availability of the representation formula. In the remaining cases we find the compatibility conditions for the data which also yield such a reduction.     

Constrained multi-degree reduction of triangular Bézier surfaces

This paper proposes and applies a method to sort two-dimensional control points of triangular Bezier surfaces in a row vector. Using the property of bivariate Jacobi basis functions, it further presents two algorithms for multi-degree reduction of triangular Bezier surfaces with constraints, providing explicit degree-reduced surfaces. The first algorithm can obtain the explicit representation of the optimal degree-reduced surfaces and the approximating error in both boundary curve constraints and corner constraints. But it has to solve the inversion of a matrix whose degree is related with the original surface. The second algorithm entails no matrix inversion to bring about computational instability, gives stable degree-reduced surfaces quickly, and presents the error bound. In the end, the paper proves the efficiency of the two algorithms through examples and error analysis.