Accurate bidiagonal decomposition of collocation matrices of weighted ϕ‐transformed systems

Given a system of functions, we introduce the concept of weighted φ‐transformed system, which will include a very large class of useful representations in Statistics and Computer Aided Geometric Design. An accurate bidiagonal decomposition of the collocation matrices of these systems is obtained. This decomposition is used to present computational methods with high relative accuracy for solving algebraic problems with collocation matrices of weighted φ‐transformed systems such as the computation of eigenvalues, singular values, and the solution of some linear systems. Numerical examples illustrate the accuracy of the performed computations.     

Hermitian‐type generalized singular value decomposition with applications

For a pair of given Hermitian matrix A and rectangular matrix B with the same row number, we reformulate a well‐known simultaneous Hermitian‐type generalized singular value decomposition (HGSVD) with more precise structure and parameters and use it to derive some algebraic properties of the linear Hermitian matrix function A?BXB* and Hermitian solution of the matrix equation BXB* = A, and the canonical form of a partitioned Hermitian matrix and some optimization problems on its inertia and rank. Copyright © 2012 John Wiley & Sons, Ltd.     

Existence of positive solutions for nonlinear singular boundary value problem with p‐Laplacian

In this paper, we study the existence of positive solutions for the following Sturm–Liouville‐like four‐point singular boundary value problem (BVP) with p‐Laplacian where ?_p(s)=|s|^p?2 s, p>1, f is a lower semi‐continuous function. Using the fixed‐point theorem of cone expansion and compression of norm type, the existence of positive solution and infinitely many positive solutions for Sturm–Liouville‐like singular BVP with p‐Laplacian are obtained. Copyright © 2009 John Wiley & Sons, Ltd.     

Investigations of nonstandard,Mickens‐type,finite‐difference schemes for singular boundary value problems in cylindrical or spherical coordinates

It is well known that standard finite‐difference schemes for singular boundary value problems involving the Laplacian have difficulty capturing the singular (??(1/r) or ??(log r)) behavior of the solution near the origin (r = 0). New nonstandard finite‐difference schemes that can capture this behavior exactly for certain singular boundary value problems encountered in theoretical aerodynamics are presented here. These schemes are special cases of nonstandard finite differences which have been extensively researched by Professor Ronald E. Mickens of Clark Atlanta University in their most general form. Several examples of these “Mickens‐type” finite differences that illustrate both their accuracy and utility for singular boundary value problems in both cylindrical and spherical co‐ordinates are investigated. The numerical results generated by the Mickens‐type schemes are compared favorably with solutions obtained from standard finite‐difference schemes. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 380–398, 2003.     

The characteristic finite element alternating direction method with moving meshes for nonlinear convection‐dominated diffusion problems

In modern numerical simulation of prospecting and exploiting oil‐gas resources and in environmental science, it is necessary to consider numerical method of nonlinear convection‐dominated diffusion problems. This thesis, starting from actual conditions such as the three‐dimensional characteristics of large‐scale science‐engineering computation, puts forward a kind of characteristic finite element alternating direction method with moving meshes. Some techniques, such as calculus of variations, operator‐splitting, generalized L² projection, energy method, negative norm estimate, the theory of prior estimates and techniques, are adopted. Optimal order estimates in L² norm are derived to determine the errors in the approximate solution. Thus the important theoretical problem has been solved. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

An O(k2 + kh2 + h4) arithmetic average discretization for the solution of 1‐D nonlinear parabolic equations

This article develops a new two‐level three‐point implicit finite difference scheme of order 2 in time and 4 in space based on arithmetic average discretization for the solution of nonlinear parabolic equation ε u_xx = f(x, t, u, u_x, u_t), 0 < x < 1, t > 0 subject to appropriate initial and Dirichlet boundary conditions, where ε > 0 is a small positive constant. We also propose a new explicit difference scheme of order 2 in time and 4 in space for the estimates of (?u/?x). The main objective is the proposed formulas are directly applicable to both singular and nonsingular problems. We do not require any fictitious points outside the solution region and any special technique to handle the singular problems. Stability analysis of a model problem is discussed. Numerical results are provided to validate the usefulness of the proposed formulas. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Some basic boundary value problems of the plane thermoelasticity with microtemperatures

The present paper is devoted to the two‐dimensional version of statics of the linear theory of elastic materials with inner structure whose particles, in addition to the classical displacement and temperature fields, possess microtemperatures. Some results of the classical theories of elasticity and thermoelasticity are generalized. The Green's formulas in the case under consideration are obtained, basic boundary value problems are formulated, and uniqueness theorems are proved. The fundamental matrix of solutions for the governing system of the model and the corresponding single and double layer thermoelastopotentials are constructed. Properties of the potentials are studied. Applying the potential method, for the first and second boundary value problems, we construct singular integral equations of the second kind and prove the existence theorems of solutions for the bounded and unbounded domains. This paper describes the use of the LaTeX2? mmaauth.cls class file for setting papers for Mathematical Methods in the Applied Sciences. Copyright © 2012 John Wiley & Sons, Ltd.     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

Homogenizing the viscoelasticity problem with long-term memory

The system of integro-differential equations describing the small oscillations of an ?-periodic viscoelastic material with long-term memory is considered. Using the two-scale convergencemethod, we construct the systemof homogenized equations and prove the strong convergence as ? → 0 of the solutions of prelimit problems to the solution of the homogenized problem in the norm of the space L ².     

Lanczos and the Riemannian SVD in information retrieval applications

Variations of the latent semantic indexing (LSI) method in information retrieval (IR) require the computation of singular subspaces associated with the k dominant singular values of a large m × n sparse matrix A, where k?min(m,n). The Riemannian SVD was recently generalized to low‐rank matrices arising in IR and shown to be an effective approach for formulating an enhanced semantic model that captures the latent term‐document structure of the data. However, in terms of storage and computation requirements, its implementation can be much improved for large‐scale applications. We discuss an efficient and reliable algorithm, called SPK‐RSVD‐LSI, as an alternative approach for deriving the enhanced semantic model. The algorithm combines the generalized Riemannian SVD and the Lanczos method with full reorthogonalization and explicit restart strategies. We demonstrate that our approach performs as well as the original low‐rank Riemannian SVD method by comparing their retrieval performance on a well‐known benchmark document collection. Copyright 2004 John Wiley & Sons, Ltd.     

Existence of solution for the n-dimension second order semilinear hyperbolic equations

In the paper we establish the local and global existence of solution for the n-dimensional second order semilinear hyperbolic equation with a strongly singular coefficient which appears in the boundary-value problems of fluid dynamics. Based on the analysis about the loss of regularity on the line t=0 for the solution of the corresponding linear equation and the decay at infinity which caused by the singular coefficient, we obtain the existence of a small solution for the semilinear equation by use of fixed point theorem.     

Well‐posedness for the incompressible magneto‐hydrodynamic system

This paper is concerned with well‐posedness of the incompressible magneto‐hydrodynamics (MHD) system. In particular, we prove the existence of a global mild solution in BMO^?1 for small data which is also unique in the space C([0, ∞); BMO^?1). We also establish the existence of a local mild solution in bmo^?1 for small data and its uniqueness in C([0, T); bmo^?1). In establishing our results an important role is played by the continuity of the bilinear form which was proved previously by Kock and Tataru. In this paper, we give a new proof of this result by using the weighted L^p‐boundedness of the maximal function. Copyright © 2006 John Wiley & Sons, Ltd.     

A layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011     

A parameter‐uniform scheme for the parabolic singularly perturbed problem with a delay in time

In this paper, a parameter‐uniform numerical scheme for the solution of singularly perturbed parabolic convection–diffusion problems with a delay in time defined on a rectangular domain is suggested. The presence of the small diffusion parameter ? leads to a parabolic right boundary layer. A collocation method consisting of cubic B ‐spline basis functions on an appropriate piecewise‐uniform mesh is used to discretize the system of ordinary differential equations obtained by using Rothe's method on an equidistant mesh in the temporal direction. The parameter‐uniform convergence of the method is shown by establishing the theoretical error bounds. The numerical results of the test problems validate the theoretical error bounds.     

Determination of residual stresses in glass shells by a computational-experimental method

A method is proposed for determining the technological residual stresses in glass shells of revolution based on the solution of inverse problems of shell theory with proper stresses and the use of experimental data obtained by the photoelastic method. Translated fromMatematicheskie Metody i Fiziko-Mekhanicheskie Polya, Issue 35, 1992, pp. 100–105.     

Correctors for some asymptotic problems

In the theory of anisotropic singular perturbation boundary value problems, the solution u _ɛ does not converge, in the H ¹-norm on the whole domain, towards some u ₀. In this paper we construct correctors to have good approximations of u _ɛ in the H ¹-norm on the whole domain. Since the anisotropic singular perturbation problems can be connected to the study of the asymptotic behaviour of problems defined in cylindrical domains becoming unbounded in some directions, we transpose our results for such problems.     

A new algorithm for generalized Sylvester‐observer equation and its application to state and velocity estimations in vibrating systems

We propose a new algorithm for block‐wise solution of the generalized Sylvester‐observer equation XA?FXE = GC, where the matrices A, E, and C are given, the matrices X, F, and G need to be computed, and matrix E may be singular. The algorithm is based on an orthogonal decomposition of the triplet (A, E, C) into the observer‐Hessenberg‐triangular form. It is a natural generalization of the widely known observer‐Hessenberg algorithm for the Sylvester‐observer equation: XA?FX = GC, which arises in state estimation of a standard first‐order state‐space control system. An application of the proposed algorithm is made to state and velocity estimations of second‐order control systems modeling a wide variety of vibrating structures. For dense un‐structured data, the proposed algorithm is more efficient than the recently proposed SVD‐based algorithm of the authors. Copyright © 2010 John Wiley & Sons, Ltd.     

Singular perturbations of some nonlinear problems

We deal with singular perturbations of nonlinear problems depending on a small parameter ε > 0. First we consider the abstract theory of singular perturbations of variational inequalities involving some nonlinear operators, defined in Banach spaces, and describe the asymptotic behavior of these solutions as ε → 0. Then these abstract results are applied to some boundary value problems. Bibliography: 15 titles.     

Quantitative unique continuation for Neumann problems of elliptic equations with weight

Let P(ω, ?) be an elliptic operator with weight ω, and let u be a solution in some Lipschitz domains to ?P(ω, ?u)+W?u+Vu=0 with sharp singular potentials W and V. The weighted doubling estimates, the weighted three‐ball inequalities and the unique continuation at the boundary for solution u of the Neumann problem are established in this paper. Copyright © 2009 John Wiley & Sons, Ltd.     

Finite element superconvergence approximation for one‐dimensional singularly perturbed problems

Superconvergence approximations of singularly perturbed two‐point boundary value problems of reaction‐diffusion type and convection‐diffusion type are studied. By applying the standard finite element method of any fixed order p on a modified Shishkin mesh, superconvergence error bounds of (N^?1 ln (N + 1))^p+1 in a discrete energy norm in approximating problems with the exponential type boundary layers are established. The error bounds are uniformly valid with respect to the singular perturbation parameter. Numerical tests indicate that the error estimates are sharp; in particular, the logarithmic factor is not removable. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 374–395, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10001