Connections between discontinuous Galerkin and nonconforming finite element methods for the Stokes equations

We study a discontinuous Galerkin finite element method (DGFEM) for the Stokes equations with a weak stabilization of the viscous term. We prove that, as the stabilization parameter γ tends to infinity, the solution converges at speed γ^?1 to the solution of some stable and well‐known nonconforming finite element methods (NCFEM) for the Stokes equations. In addition, we show that an a posteriori error estimator for the DGFEM‐solution based on the reconstruction of a locally conservative H(div, Ω)‐tensor tends at the same speed to a classical a posteriori error estimator for the NCFEM‐solution. These results can be used to affirm the robustness of the DGFEM‐method and also underline the close relationship between the two approaches. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A remark on the decay of solutions to the 3‐D Navier–Stokes equations

In this paper we derive a decay rate of the L²‐norm of the solution to the 3‐D Navier–Stokes equations. Although the result which is proved by Fourier splitting method is well known, our method is new, concise and direct. Moreover, it turns out that the new method established here has a wide application on other equations. Copyright © 2007 John Wiley & Sons, Ltd.     

Mathematical and physical aspect of the solution of the Cauchy problem for a system of equations describing a nonlocal model of propagation of heat with finite speed

In our paper we present a new system of equations describing a nonlocal model of propagation of heat with finite speed in three-dimensional space. Such a system of equations is described by a system of integral – differential equations. At first using the modiffied Cagniard de Hoop method, we construct the fundamental solution of this system of equations. On the basis of the constructed fundamental solution we obtain the explicite formulate of the solution of the Cauchy problem for this system of equations and applying the method of Sobolev and Biesov spaces, we get L^p – L^q time decay estimate for the solution of the Cauchy problem. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Inviscid Limit for the Steady Incompressible Navier-Stokes Equations in the Three Dimension*

In this paper, the authors consider the zero-viscosity limit of the three dimensional incompressible steady Navier-Stokes equations in a half space R⁺×R². The result shows that the solution of three dimensional incompressible steady Navier-Stokes equations converges to the solution of three dimensional incompressible steady Euler equations in Sobolev space as the viscosity coefficient going to zero. The method is based on a new weighted energy estimates and Nash-Moser itera...     

Position vector of a time-like slant helix in Minkowski 3-space

In this paper, position vector of a time-like slant helix with respect to standard frame of Minkowski space E³₁ is studied in terms of Frenet equations. First, a vector differential equation of third order is constructed to determine position vector of an arbitrary time-like slant helix. In terms of solution, we determine the parametric representation of the slant helices from the intrinsic equations. Thereafter, we apply this method to find the representation of a time-like Salkowski and time-like anti-Salkowski curves as examples of a slant helices, by means of intrinsic equations. Moreover, we present some new characterizations of slant helices and illustrate some examples of our main results.     

Some new analysis results for a class of interface problems

Interface problems modeled by differential equations have many applications in mathematical biology, fluid mechanics, material sciences, and many other areas. Typically, interface problems are characterized by discontinuities in the coefficients and/or the Dirac delta function singularities in the source term. Because of these irregularities, solutions to the differential equations are not smooth or discontinuous. In this paper, some new results on the jump conditions of the solution across the interface are derived using the distribution theory and the theory of weak solutions. Some theoretical results on the boundary singularity in which the singular delta function is at the boundary are obtained. Finally, the proof of the convergency of the immersed boundary (IB) method is presented. The IB method is shown to be first‐order convergent in L ^∞ norm. Copyright © 2013 John Wiley & Sons, Ltd.     

A priori error estimates for interior penalty versions of the local discontinuous Galerkin method applied to transport equations

The local discontinuous Galerkin method has been developed recently by Cockburn and Shu for convection‐dominated convection‐diffusion equations. In this article, we consider versions of this method with interior penalties for the numerical solution of transport equations, and derive a priori error estimates. We consider two interior penalty methods, one that penalizes jumps in the solution across interelement boundaries, and another that also penalizes jumps in the diffusive flux across such boundaries. For the first penalty method, we demonstrate convergence of order k in the L^∞(L²) norm when polynomials of minimal degree k are used, and for the second penalty method, we demonstrate convergence of order k+1/2. Through a parabolic lift argument, we show improved convergence of order k+1/2 (k+1) in the L²(L²) norm for the first penalty method with a penalty parameter of order one (h^?1). © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 545–564, 2001     

A parameter-free stabilized finite element method for scalar advection-diffusion problems

We formulate and study numerically a new, parameter-free stabilized finite element method for advection-diffusion problems. Using properties of compatible finite element spaces we establish connection between nodal diffusive fluxes and one-dimensional diffusion equations on the edges of the mesh. To define the stabilized method we extend this relationship to the advection-diffusion case by solving simplified one-dimensional versions of the governing equations on the edges. Then we use H(curl)-conforming edge elements to expand the resulting edge fluxes into an exponentially fitted flux field inside each element. Substitution of the nodal flux by this new flux completes the formulation of the method. Utilization of edge elements to define the numerical flux and the lack of stabilization parameters differentiate our approach from other stabilized methods. Numerical studies with representative advection-diffusion test problems confirm the excellent stability and robustness of the new method. In particular, the results show minimal overshoots and undershoots for both internal and boundary layers on uniform and non-uniform grids.     

Natural frequencies of thin rectangular plates using homotopy-perturbation method

The homotopy perturbation method (HPM) was developed to search for asymptotic solutions of nonlinear problems involving parabolic partial differential equations with variable coefficients. This paper illustrates that HPM be easily adapted to solve parabolic partial differential equations with constant coefficients. Natural frequencies of a rectangular plate of uniform thickness, simply-supported on all sides, are obtained with minimum amount of computation. The solution is shown to converge rapidly to a combination of sine and cosine functions. Truncating the series solution by using only the first three terms of the sine and cosine functions as compared to the exact solution results in an absolute error not exceeding 2 × 10⁻⁴ and 9×10⁻⁴ for the trigonometric functions respectively. HPM is then applied to solve the nonlinear problem of a rectangular plate of variable thickness. A direct expression for the eigenvalues (natural frequencies) of the rectangular plate is obtained as compared to determining its eigenvalues by solving the characteristic equation using the conventional method. Comparison of results for the frequency parameter with existing literature show that HPM is highly efficient and accurate. Natural frequencies of a simply-supported guitar soundboard were obtained using an equivalent rectangular plate with the same boundary condition.     

Generalized Alternating Direction Implicit Method for Projected Generalized Lyapunov Equations

We generalize an alternating direction implicit method for projected generalized Lyapunov equations. Low rank versions of this method is also presented that can be used to compute a low rank approximation of the solution of Lyapunov equations with symmetric, positive semidefinite right‐hand side. Numerical example is given. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A family of Newton methods for nonsmooth constrained systems with nonisolated solutions

We propose a new family of Newton-type methods for the solution of constrained systems of equations. Under suitable conditions, that do not include differentiability or local uniqueness of solutions, local, quadratic convergence to a solution of the system of equations can be established. We show that as particular instances of the method we obtain inexact versions of both a recently introduced LP-based Newton method and of a Levenberg-Marquardt algorithm for the solution of systems with nonisolated solutions, and improve on corresponding existing results.     

Lp-Lq time decay estimate to the solution of the Cauchy problem of the system of equations for nonlocal model of hyperbolic thermoelasticity theory

The aim of this paper is to present a new system of equations describing nonlocal model of hyperbolic thermoelasticity theory. We used the Papkin and Gurtin approach based on the constitutive relations for internal energy e(x), and heat flux q(x), with integral terms. Such system of equations describes the propagation of thermal perturbation with finite velocity. Using the modified Cagniard–de Hoop's method we constructed the matrix of fundamental solutions for this system of equations in three–dimensional space. Basing on the constructed matrix of fundamental solutions in the explicit formula we represent the solution of the Cauchy problem to this system of equations in the form of some kind of convolutions. Next, applying the method of Sobolev spaces, we obtain the L^p−L^q time decay estimate to the solution of the Cauchy problem. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the existence and stability of solutions for the micropolar fluids in exterior domains

We prove the existence of a global strong solution in some class of Marcinkiewicz spaces for the micropolar fluid in an exterior domain of R³, with initial conditions being a non‐smooth disturbance of a steady solution. We also analyse the large time behaviour of those solutions and apply our results in the context of the Navier–Stokes equations. Copyright © 2007 John Wiley & Sons, Ltd.     

Convergence of euler-stokes splitting of the navier-stokes equations

We consider approximation by partial time steps of a smooth solution of the Navier-Stokes equations in a smooth domain in two or three space dimensions with no-slip boundary condition. For small k > 0, we alternate the solution for time k of the inviscid Euler equations, with tangential boundary condition, and the solution of the linear Stokes equations for time k, with the no-slip condition imposed. We show that this approximation remains bounded in H^2,p and is accurate to order k in L^p for p > ∞. The principal difficulty is that the initial state for each Stokes step has tangential velocity at the boundary generated during the Euler step, and thus does not satisfy the boundary condition for the Stokes step. The validity of such a fractional step method or splitting is an underlying principle for some computational methods. © 1994 John Wiley & Sons, Inc.     

An alternating direction implicit method for orthogonal spline collocation linear systems

Summary An Alternating Direction Implicit method is analyzed for the solution of linear systems arising in high-order, tensor-product orthogonal spline collocation applied to some separable, second order, linear, elliptic partial differential equations in rectangles. On anNxN partition, with Jordan's selection of the acceleration parameters, the method requiresO(N ² ln ² N) arithmetic operations to produce an approximation whose accuracy, in theH ¹-norm, is that of the collocation solution.     

Solving soliton equations with self-consistent sources by method of variation of parameters

Starting from the solutions of soliton equations and corresponding eigenfunctions obtained by Darboux transformation, we present a new method to solve soliton equations with self-consistent sources (SESCS) based on method of variation of parameters. The KdV equation with self-consistent sources (KdVSCS) is used as a model to illustrate this new method. In addition, we apply this method to construct some new solutions of the derivative nonlinear Schrödinger equation with self-consistent sources (DNLSSCS) such as phase solution, dark soliton solution, bright soliton solution and breather-type solution.     

A numerical approach for solving a class of the nonlinear Lane–Emden type equations arising in astrophysics

In this paper, a collocation method based on the Bessel polynomials is presented for the approximate solution of a class of the nonlinear Lane–Emden type equations, which have many applications in mathematical physics. The exact solution can be obtained if the exact solution is polynomial. In other cases, such as an increasing number of nodes, a good approximation can be obtained with applicable errors. In addition, the method is presented with error and stability analysis. The numerical results show the effectiveness of the method for this type of equations. Comparing the methodology with some known techniques shows that the present approach is relatively easy and highly accurate. Copyright © 2011 John Wiley & Sons, Ltd.     

A hybrid Arnoldi-Faber iterative method for nonsymmetric systems of linear equations

Summary We present here a new hybrid method for the iterative solution of large sparse nonsymmetric systems of linear equations, say of the formAx=b, whereA ^{N, N} , withA nonsingular, andb ^N are given. This hybrid method begins with a limited number of steps of the Arnoldi method to obtain some information on the location of the spectrum ofA, and then switches to a Richardson iterative method based on Faber polynomials. For a polygonal domain, the Faber polynomials can be constructed recursively from the parameters in the Schwarz-Christoffel mapping function. In four specific numerical examples of non-normal matrices, we show that this hybrid algorithm converges quite well and is approximately as fast or faster than the hybrid GMRES or restarted versions of the GMRES algorithm. It is, however, sensitive (as other hybrid methods also are) to the amount of information on the spectrum ofA acquired during the first (Arnoldi) phase of this procedure.     

Fourier–Finite-element Approximation of Elliptic Interface Problems in Axisymmetric Domains

The paper deals with the Fourier-finite-element method (FFEM), which combines the approximate Fourier method with the finite-element method, and its application to Poisson-like equations −p̂Δ₃û = f̂ in three-dimensional axisymmetric domains Ωˆ. Here, pˆ is a piecewise constant coefficient having a jump at some axisymmetric interface. Special emphasis is given to estimates of the Fourier-finite-element error in the Sobolev space H¹(Ωˆ), if the interface is smooth or if it meets the boundary of Ωˆ at some edge. In general, the solution û contains a singularity at the interface, which is described by a tensor product representation and treated numerically by appropriate mesh grading in the meridian plane of Ωˆ. The rate of convergence of the combined approximation in H¹(Ωˆ) is proved to be 𝒪(h+N⁻¹) (h, N: the parameters of the finite-element- and Fourier-approximation, with h→0, N→∞). The theoretical results are confirmed by numerical experiments.     

On the convergence of the midpoint method

The midpoint method is an iterative method for the solution of nonlinear equations in a Banach space. Convergence results for this method have been studied in [3, 4, 9, 12]. Here we show how to improve and extend these results. In particular, we use hypotheses on the second Fréchet derivative of the nonlinear operator instead of the third-derivative hypotheses employed in the previous results and we obtain Banach space versions of some results that were derived in [9, 12] only in the real or complex space. We also provide various examples that validate our results.