A collocation method to solve higher order linear complex differential equations in rectangular domains

In this article, a collocation method is developed to find an approximate solution of higher order linear complex differential equations with variable coefficients in rectangular domains. This method is essentially based on the matrix representations of the truncated Taylor series of the expressions in equation and their derivates, which consist of collocation points defined in the given domain. Some numerical examples with initial and boundary conditions are given to show the properties of the method. All results were computed using a program written in scientific WorkPlace v5.5 and Maple v12. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Coupling of viscous and inviscid Stokes equations via a domain decomposition method for finite elements

Summary A generalized Stokes problem is addressed in the framework of a domain decomposition method, in which the physical computational domain is partitioned into two subdomains ₁ and ₂.Three different situations are covered. In the former, the viscous terms are kept in both subdomains. Then we consider the case in which viscosity is dropped out everywhere in . Finally, a hybrid situation in which viscosity is dropped out only in ₁ is addressed. The latter is motivated by physical applications.In all cases, correct transmission conditions across the interface between ₁ and ₂ are devised, and an iterative procedure involving the successive resolution of two subproblems is proposed.The numerical discretization is based upon appropriate finite elements, and stability and convergence analysis is carried out.We also prove that the iteration-by-subdomain algorithms which are associated with the various domain decomposition approaches converge with a rate independent of the finite element mesh size.This work was partially supported by CIRA S.p.A. under the contract Coupling of Euler and Navier-Stokes equations in hypersonic flowsDeceased     

A higher order difference method for singularly perturbed parabolic partial differential equations

This paper studies a higher order numerical method for the singularly perturbed parabolic convection-diffusion problems where the diffusion term is multiplied by a small perturbation parameter. In general, the solutions of these type of problems have a boundary layer. Here, we generate a spatial adaptive mesh based on the equidistribution of a positive monitor function. Implicit Euler method is used to discretize the time variable and an upwind scheme is considered in space direction. A higher order convergent solution with respect to space and time is obtained using the postprocessing based extrapolation approach. It is observed that the convergence is independent of perturbation parameter. This technique enhances the order of accuracy from first order uniform convergence to second order uniform convergence in space as well as in time. Comparative study with the existed meshes show the highly effective behavior of the present method.     

Fictitious domain method for a system of equations in elasticity theory

Two versions of the fictitious domain method are considered for the system of equilibrium equations of a nonhomogeneous anisotropic elastic solid with rigid clamping. The rate of convergence bounds are derived:.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 64, pp. 16–23, 1988.     

On a solvable system of p difference equations of higher order

Periodica Mathematica Hungarica - In this paper we present the well-defined solution of the following system of higher-order rational difference equations: $$begin{aligned}...     

An almost fourth order uniformly convergent domain decomposition method for a coupled system of singularly perturbed reaction-diffusion equations

We consider a system of M(≥2) singularly perturbed equations of reaction-diffusion type coupled through the reaction term. A high order Schwarz domain decomposition method is developed to solve the system numerically. The method splits the original domain into three overlapping subdomains. On two boundary layer subdomains we use a compact fourth order difference scheme on a uniform mesh while on the interior subdomain we use a hybrid scheme on a uniform mesh. We prove that the method is almost fourth order ε-uniformly convergent. Furthermore, we prove that when ε is small, one iteration is sufficient to get almost fourth order ε-uniform convergence. Numerical experiments are performed to support the theoretical results.     

Dynamical system solutions for homogeneous linear functional equations of higher order

We present the solution of a large class of homogeneous linear functional equations of higher order by using ideas from dynamical systems. A particularly simple example from this class is the functional equation $$f(x) = \frac{1}{2}f \left(\frac{x}{2}\right) + \frac{1}{2}f \left(\frac{x+1}{2}\right), \quad 0 < x < 1.$$ Equations such as these have found important applications in wavelet theory by Hilberdink (Aequa Math 61(1–2):179–189, 2001) where they are called dilation equations and are usually solved by Fourier methods by Daubechies (Comm Pure Appl Math 41(7):909–996, 1988) or iteration methods of Daubechies (SIAM J Math Anal 22(5):1388–1410, 1991). A recent result of Góra (Ergod Theory Dyn Syst 29(5):1549–1583, 2009) allows us to represent the solution as an infinite series that is determined by the dynamics of a map that is defined by the functional equation. In this problem the interplay between dynamical systems and solutions of functional equations is brought into sharp focus.     

Renormalization group method applied to the primitive equations

In this article we study the limit, as the Rossby number ε goes to zero, of the primitive equations of the atmosphere and the ocean. From the mathematical viewpoint we study the averaging of a penalization problem displaying oscillations generated by an antisymmetric operator and by the presence of two time scales.     

A fourth-order numerical method for the planetary geostrophic equations with inviscid geostrophic balance

The planetary geostrophic equations with inviscid balance equation are reformulated in an alternate form, and a fourth-order finite difference numerical method of solution is proposed and analyzed in this article. In the reformulation, there is only one prognostic equation for the temperature field and the velocity field is statically determined by the planetary geostrophic balance combined with the incompressibility condition. The key observation is that all the velocity profiles can be explicitly determined by the temperature gradient, by utilizing the special form of the Coriolis parameter. This brings convenience and efficiency in the numerical study. In the fourth-order scheme, the temperature is dynamically updated at the regular numerical grid by long-stencil approximation, along with a one-sided extrapolation near the boundary. The velocity variables are recovered by special solvers on the 3-D staggered grid. Furthermore, it is shown that the numerical velocity field is divergence-free at the discrete level in a suitable sense. Fourth order convergence is proven under mild regularity requirements. R. Samelson was supported by NSF grant OCE04-24516 and Navy ONR grant N00014-05-1-0891. R. Temam was supported by NSF grant DMS-0604235 and the research fund of Indiana University. S. Wang was supported by NSF grant DMS-0605067 and Navy ONR grant N00014-05-1-0218.     

On limiting for higher order discontinuous Galerkin method for 2D Euler equations

We present an implementation of discontinuous Galerkin method for 2-D Euler equations on Cartesian meshes using tensor product Lagrange polynomials based on Gauss nodes. The scheme is stabilized by a version of the slope limiter which is adapted for tensor product basis functions together with a positivity preserving limiter. We also incorporate and test shock indicators to determine which cells need limiting. Several numerical results are presented to demonstrate that the proposed approach is capable of computing complex discontinuous flows in a stable and accurate fashion.     

Bounds for analytic solutions to integral equations in the complex domain

We establish the existence of analytic solutions to integral equations in the complex domain when the nonlinearity satisfies either a growth condition or a monotonic type condition.     

A higher order method for input-affine uncertain systems

Uncertainty is unavoidable in modeling dynamical systems and it may be represented mathematically by differential inclusions. In the past, we proposed an algorithm to compute validated solutions of differential inclusions; here we provide several theoretical improvements to the algorithm, including its extension to piecewise constant and sinusoidal approximations of uncertain inputs, updates on the affine approximation bounds and a generalized formula for the analytical error. The approach proposed is able to achieve higher order convergence with respect to the current state-of-the-art. We implemented the methodology in Ariadne, a library for the verification of continuous and hybrid systems. For evaluation purposes, we introduce ten systems from the literature, with varying degrees of nonlinearity, number of variables and uncertain inputs. The results are hereby compared with two state-of-the-art approaches to time-varying uncertainties in nonlinear systems.     

Analytic solutions to integral equations in the complex domain

Schauder's fixed point theorem and the Leray–Schauder alternative are used to establish the existence of analytic solutions to integral equations in the complex domain.     

Laplace transform for the solution of higher order deformation equations arising in the homotopy analysis method

A new analytic approach for solving nonlinear ordinary differential equations with initial conditions is proposed. First, the homotopy analysis method is used to transform a nonlinear differential equation into a system of linear differential equations; then, the Laplace transform method is applied to solve the resulting linear initial value problems; finally, the solutions to the linear initial value problems are employed to form a convergent series solution to the given problem. The main advantage of the new approach is that it provides an effective way to solve the higher order deformation equations arising in the homotopy analysis method.     

First order decoupled method of the primitive equations of the ocean I: Time discretization

In this article, we study the time discretization of the first order decoupled semi-implicit scheme of the 3D primitive equations of the ocean in the case of the non-slip and non-heat flux boundary conditions on the side. The almost unconditional stability of the scheme and the optimal error estimates of the time discrete velocity, pressure and density are provided.     

A HIGH ORDER METHOD FOR NON-SMOOTH FREDHOLM EQUATIONS

1.IntroductionConsidertheequationwherek(s,t)=k(f)tandf(s)aregiven,uistheunknownsolution.SinceitisrelatedcloselytoWiener-Hopfequationsandisveryimportantinpractice,therearemanynumericalresultsaboutit(e.g.[1--11]).Itiswellknownthattheaccuracyoftheapproximati…