Partial Fourier approximation of the Lamé equations in axisymmetric domains

In this paper, we study the partial Fourier method for treating the Lamé equations in three‐dimensional axisymmetric domains subjected to non‐axisymmetric loads. We consider the mixed boundary value problem of the linear theory of elasticity with the displacement û , the body force f̂ ϵ (L₂)³ and homogeneous Dirichlet and Neumann boundary conditions. The partial Fourier decomposition reduces, without any error, the three‐dimensional boundary value problem to an infinite sequence of two‐dimensional boundary value problems, whose solutions û _n (n = 0, 1, 2,…) are the Fourier coefficients of û . This process of dimension reduction is described, and appropriate function spaces are given to characterize the reduced problems in two dimensions. The trace properties of these spaces on the rotational axis and some properties of the Fourier coefficients û _n are proved, which are important for further numerical treatment, e.g. by the finite‐element method. Moreover, generalized completeness relations are described for the variational equation, the stresses and the strains. The properties of the resulting system of two‐dimensional problems are characterized. Particularly, a priori estimates of the Fourier coefficients û _n and of the error of the partial Fourier approximation are given. Copyright © 1999 John Wiley & Sons, Ltd.     

The minimax approach for a class of variable order fractional differential equation

This paper introduces an approximate solution for Liouville‐Caputo variable order fractional differential equations with order 0 < α(t) ≤ 1 . The solution is adapted using a family of fractional‐order Chebyshev functions with unknown coefficients. These coefficients have been obtained by using an optimization approach based on minimax technique and the least pth optimization function. Several linear and nonlinear fractional‐order differential equations are discussed using the proposed technique for fixed and variable order fractional‐order derivatives. Moreover, the response of RC charging circuit with variable order fractional capacitor is studied for different cases. Several comparisons with related published techniques have been added to illustrate the accuracy of the proposed approach.     

On the divergence of Fourier series of functions in several variables

The paper deals with the problems of divergence of the series from absolute values of the Fourier coefficients of functions in several variables. It is proved that as the dimension of the space increases, the absolute convergence of Fourier series with respect to any complete orthnormal system (ONS) of functions with continuous partial derivatives becomes worse. For instance, for any ? ∈ (0, 2) there exists a function in variables $k > \frac{{2(2 - \varepsilon )}} {\varepsilon }$ having all the continuous partial derivatives, however the series of absolute values of its coefficients with respect to any complete orthnormal system diverges in power 2 ? ?.     

A Padé-based algorithm for overcoming the Gibbs phenomenon

Truncated Fourier series and trigonometric interpolants converge slowly for functions with jumps in value or derivatives. The standard Fourier–Padé approximation, which is known to improve on the convergence of partial summation in the case of periodic, globally analytic functions, is here extended to functions with jumps. The resulting methods (given either expansion coefficients or function values) exhibit exponential convergence globally for piecewise analytic functions when the jump location(s) are known. Implementation requires just the solution of a linear system, as in standard Padé approximation. The new methods compare favorably in experiments with existing techniques.     

How well does the finite Fourier transform approximate the Fourier transform?

We show that the answer to the question in the title is “very well indeed.” In particular, we prove that, throughout the maximum possible range, the finite Fourier coefficients provide a good approximation to the Fourier coefficients of a piecewise continuous function. For a continuous periodic function, the size of the error is estimated in terms of the modulus of continuity of the function. The estimates improve commensurately as the functions become smoother. We also show that the partial sums of the finite Fourier transform provide essentially as good an approximation to the function and its derivatives as the partial sums of the ordinary Fourier series. Along the way we establish analogues of the Riemann‐Lebesgue lemma and the localization principle. © 2004 Wiley Periodicals, Inc.     

Pointwise Behavior of Double Fourier Series of Functions of Bounded Variation

We extend some recent results of S. A. Telyakovskii on the uniform boundedness of the partial sums of Fourier series of functions of bounded variation to periodic functions of two variables, which are of bounded variation in the sense of Hardy. As corollaries, we obtain the classical Parseval formula, the convergence theorem of the series involving the sine Fourier coefficients, and a lower estimate of the best approximation by trigonometric polynomials in the metric of L in a sharpened version.     

Basics of a generalized Wiman–Valiron theory for monogenic Taylor series of finite convergence radius

In this paper, we develop the basic concepts for a generalized Wiman–Valiron theory for Clifford algebra valued functions that satisfy inside an n + 1-dimensional ball the higher dimensional Cauchy-Riemann system ${\frac{\partial f}{\partial x_0} + \sum_{i=1}^n e_i\frac{\partial f}{\partial x_i}=0}In this paper, we develop the basic concepts for a generalized Wiman–Valiron theory for Clifford algebra valued functions that satisfy inside an n + 1-dimensional ball the higher dimensional Cauchy-Riemann system \frac?f?x₀ + ?_i=1ⁿ e_i\frac?f?x_i=0{\frac{\partial f}{\partial x_0} + \sum_{i=1}^n e_i\frac{\partial f}{\partial x_i}=0} . These functions are called monogenic or Clifford holomorphic inside the ball. We introduce growth orders, the maximum term and a generalization of the central index for monogenic Taylor series of finite convergence radius. Our goal is to establish explicit relations between these entities in order to estimate the asymptotic growth behavior of a monogenic function in a ball in terms of its Taylor coefficients. Furthermore, we exhibit a relation between the growth order of such a function f and the growth order of its partial derivatives.     

New approach to the Lax‐Wendroff modified differential equation for linear and nonlinear advection

The paper presents an enhanced analysis of the Lax‐Wendroff difference scheme—up to the eighth‐order with respect to time and space derivatives—of the modified‐partial differential equation (MDE) of the constant‐wind‐speed advection equation. The modified equation has been so far derived mainly as a fourth‐order equation. The Π ‐form of the first differential approximation (differential approximation or equivalent equation) derived by expressing the time derivatives in terms of the space derivatives is used for presenting the MDE. The obtained coefficients at higher order derivatives are analyzed for indications of the character of the dissipative and dispersive errors. The authors included a part of the stencil applied for determining the modified differential equation up to the eighth‐order of the analyzed modified differential equation for the second‐order Lax‐Wendroff scheme. Neither the derived coefficients at the space derivatives of order p ∈ (7 – 8) in the modified differential equation for the Lax‐Wendroff difference scheme nor the results of analyses on the basis of these coefficients of the group velocity, phase shift errors, or dispersive and dissipative features of the scheme have been published. The MDEs for 2 two‐step variants of the Lax‐Wendroff type difference schemes and the MacCormack predictor–corrector scheme (see MacCormack's study) constructed for the scalar hyperbolic conservation laws are also presented in this paper. The analysis of the inviscid Burgers equation solution with the initial condition in a form of a shock wave has been discussed on their basis. The inviscid Burgers equation with the source is also presented. The theory of MDE started to develop after the paper of C. W. Hirt was published in 1968.     

Fractional partial differential equations and modified Riemann-Liouville derivative new methods for solution

The paper deals with the solution of some fractional partial differential equations obtained by substituting modified Riemann-Liouville derivatives for the customary derivatives. This derivative is introduced to avoid using the so-called Caputo fractional derivative which, at the extreme, says that, if you want to get the first derivative of a function you must before have at hand its second derivative. Firstly, one gives a brief background on the fractional Taylor series of nondifferentiable functions and its consequence on the derivative chain rule. Then one considers linear fractional partial differential equations with constant coefficients, and one shows how, in some instances, one can obtain their solutions on by-passing the use of Fourier transform and/or Laplace transform. Later one develops a Lagrange method via characteristics for some linear fractional differential equations with nonconstant coefficients, and involving fractional derivatives of only one order. The key is the fractional Taylor series of non differentiable functionf(x + h) =E _α (h ^α D _x ^α )f(x).     

On determination of jumps in terms of the Abel-Poisson mean of Fourier series. II

In this paper, we generalize two important results of Bagota and Móricz [1], and generalize our earlier results in [6] from one-variable to two-variable case. As special applications, we prove that the generalized jump of f(x, y) at some point (x ₀, y ₀) can be determined by the higher order mixed partial derivatives of the Abel-Poisson mean of double Fourier series and the higher order mixed partial derivatives of the Abel-Poisson means of the three conjugate double Fourier series.     

Pointwise Behavior of Double Fourier Series of Functions of Bounded Variation

We extend some recent results of S. A. Telyakovskii on the uniform boundedness of the partial sums of Fourier series of functions of bounded variation to periodic functions of two variables, which are of bounded variation in the sense of Hardy. As corollaries, we obtain the classical Parseval formula, the convergence theorem of the series involving the sine Fourier coefficients, and a lower estimate of the best approximation by trigonometric polynomials in the metric of L in a sharpened version. This research was supported by the Hungarian National Foundation for Scientific Research under Grants TS 044 782 and T 046 192.     

Analytical approximate solutions of Riesz fractional diffusion equation and Riesz fractional advection–dispersion equation involving nonlocal space fractional derivatives

In this paper, we consider the analytical solutions of fractional partial differential equations (PDEs) with Riesz space fractional derivatives on a finite domain. Here we considered two types of fractional PDEs with Riesz space fractional derivatives such as Riesz fractional diffusion equation (RFDE) and Riesz fractional advection–dispersion equation (RFADE). The RFDE is obtained from the standard diffusion equation by replacing the second‐order space derivative with the Riesz fractional derivative of order α∈(1,2]. The RFADE is obtained from the standard advection–dispersion equation by replacing the first‐order and second‐order space derivatives with the Riesz fractional derivatives of order β∈(0,1] and of order α∈(1,2] respectively. Here the analytic solutions of both the RFDE and RFADE are derived by using modified homotopy analysis method with Fourier transform. Then, we analyze the results by numerical simulations, which demonstrate the simplicity and effectiveness of the present method. Here the space fractional derivatives are defined as Riesz fractional derivatives. Copyright © 2015 John Wiley & Sons, Ltd.     

Implicit Extension of Taylor Series Method with Numerical Derivatives

The Taylor series method is one of the earliest analytic‐numeric algorithms for approximate solution of initial value problems for ordinary differential equations. The main idea of the rehabilitation of these algorithms is based on the approximate calculation of higher order derivatives using well‐known technique for the partial differential equations. The implicit extension based on a collocation term added to the explicit truncated Taylor series. This idea is different from the general collocation method construction, which led to the implicit R‐K algorithms [1].     

An algorithm for the finite difference approximation of derivatives with arbitrary degree and order of accuracy

In this paper, we introduce an algorithm and a computer code for numerical differentiation of discrete functions. The algorithm presented is suitable for calculating derivatives of any degree with any arbitrary order of accuracy over all the known function sampling points. The algorithm introduced avoids the labour of preliminary differencing and is in fact more convenient than using the tabulated finite difference formulas, in particular when the derivatives are required with high approximation accuracy. Moreover, the given Matlab computer code can be implemented to solve boundary-value ordinary and partial differential equations with high numerical accuracy. The numerical technique is based on the undetermined coefficient method in conjunction with Taylor’s expansion. To avoid the difficulty of solving a system of linear equations, an explicit closed form equation for the weighting coefficients is derived in terms of the elementary symmetric functions. This is done by using an explicit closed formula for the Vandermonde matrix inverse. Moreover, the code is designed to give a unified approximation order throughout the given domain. A numerical differentiation example is used to investigate the validity and feasibility of the algorithm and the code. It is found that the method and the code work properly for any degree of derivative and any order of accuracy.     

L p – L q Estimates To The Solution Of The Cauchy Problem For Partial Differential Equations Describing Non-Simple Thermoelastic Materials

Theory of non-simple materials is different from that of simple materials because in it the first strain gradient is taken into consideration as the constitutive variable. The consequence of this fact, from mathematical point of view, is that the equation of motion consists either of higher order derivatives of displacement (four order derivatives) and some material parameters can depend not only on the temperature and the gradient of displacement but also on the second derivative of displacement. We consider the system of partial differential equations describing non-simple thermoelastic materials. This system consists of four scalar equations, three equations of motion and one of energy balance, describing the field of displacement and the temperature in an elastic body. Using the Fourier transform, we found the L ^p – L ^q time decay estimates of the solution of the Cauchy problem for the system of equation describing the non-simple thermoelastic materials, being important for proving the global-in-time solution of this problem. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new stability number of the Bresse‐Cattaneo system

In this paper, we consider the Bresse‐Cattaneo system with a frictional damping term and prove some optimal decay results for the L²‐norm of the solution and its higher order derivatives. In fact, we show that there is a completely new stability number δ that controls the decay rate of the solution. To prove our results, we use the energy method in the Fourier space to build some very delicate Lyapunov functionals that give the desired results. We also prove the optimality of the results by using the eigenvalues expansion method. In addition, we show that for the absence of the frictional damping term, the solution of our problem does not decay at all. This result improves some early results     

On the Wl q-Regularity of Solutions to Systems of Differential Equations in the Case When the Equations Are Constructed from Discontinuous Functions

Some solution, final in a sense from the standpoint of the theory of Sobolev spaces, is obtained to the problem of regularity of solutions to a system of (generally) nonlinear partial differential equations in the case when the system is locally close to elliptic systems of linear equations with constant coefficients. The main consequences of this result are Theorems 5 and 8. According to the first of them, the higher derivatives of an elliptic C ^l -smooth solution to a system of lth-order nonlinear partial differential equations constructed from C ^l -smooth functions meet the local Hoelder condition with every exponent , 0<<1. Theorem 8 claims that if a system of linear partial differential equations of order l with measurable coefficients and right-hand sides is uniformly elliptic then, under the hypothesis of a (sufficiently) slow variation of its leading coefficients, the degree of local integrability of lth-order partial derivatives of every W ^l _q,loc-solution, q>1, to the system coincides with the degree of local integrability of lower coefficients and right-hand sides.     

A novel 2D partial unwinding adaptive Fourier decomposition method with application to frequency domain system identification

This paper proposes a two‐dimensional (2D) partial unwinding adaptive Fourier decomposition method to identify 2D system functions. Starting from Coifman in 2000, one‐dimensional (1D) unwinding adaptive Fourier decomposition and later a type called unwinding AFD have been being studied. They are based on the Nevanlinna factorization and a maximal selection. This method provides fast‐converging rational approximations to 1D system functions. However, in the 2D case, there is no genuine unwinding decomposition. This paper proposes a 2D partial unwinding adaptive Fourier decomposition algorithm that is based on algebraic transforms reducing a 2D case to the 1D case. The proposed algorithm enables rational approximations of real coefficients to 2D system functions of real coefficients. Its fast convergence offers efficient system identification. Numerical experiments are provided, and the advantages of the proposed method are demonstrated.     

Weak Galerkin finite element methods for a fourth order parabolic equation

This paper is devoted to a newly developed weak Galerkin finite element method with the stabilization term for a linear fourth order parabolic equation, where weakly defined Laplacian operator over discontinuous functions is introduced. Priori estimates are developed and analyzed in L² and an H² type norm for both semi‐discrete and fully discrete schemes. And finally, numerical examples are provided to confirm the theoretical results.     

Equations of anisotropic elastodynamics as a symmetric hyperbolic system: Deriving the time-dependent fundamental solution

The dynamic system of anisotropic elasticity from three second order partial differential equations is written in the form of the time-dependent first order symmetric hyperbolic system with respect to displacement velocity and stress components. A new method of deriving the time-dependent fundamental solution of the obtained system is suggested in this paper. This method consists of the following. The Fourier transform image of the fundamental solution with respect to a space variable is presented as a power series expansion relative to the Fourier parameters. Then explicit formulae for the coefficients of these power series are derived successively. Using these formulae the computer calculation of fundamental solution components (displacement velocity and stress components arising from pulse point forces) has been made for general anisotropic media (orthorhombic and monoclinic) and the simulation of elastic waves has been obtained. These computational examples confirm the robustness of the suggested method.