Nonlinear calculation of arbitrarily shaped supercavitating hydrofoils near a free surface

A nonlinear exact solution to the problem of two-dimensional gravity-free incompressible potential flow around an arbitrarily shaped supercavitating hydrofoil near a free surface is obtained. A combination of Newton’s method with a functional iterative procedure is used to solve the nonlinear integral and algebraic equations of this problem. Fast and stable convergence results by starting the iteration with a readily chosen initial solution. Some representative numerical computations are made for practical hydrofoils having both generally shaped camber and leading-edge thickness distributions. The force coefficients, pressure distribution and free-streamline shapes of the cavity are calculated for each case with an execution time on an IBM 370-158 of 200-530s depending upon the initial trial solution

Numerical approximations for population growth model by Rational Chebyshev and Hermite Functions collocation approach: A comparison

This paper aims to compare rational Chebyshev (RC) and Hermite functions (HF) collocation approach to solve the Volterra's model for population growth of a species within a closed system. This model is a nonlinear integro-differential equation where the integral term represents the effect of toxin. This approach is based on orthogonal functions which will be defined. The collocation method reduces the solution of this problem to the solution of a system of algebraic equations. We also compare these methods with some other numerical results and show that the present approach is applicable for solving nonlinear integro-differential equations.Comment: 18 pages, 5 figures; Published online in the journal of "Mathematical Methods in the Applied Sciences

R-function Theory for Bending Problem of Shallow Spherical Shells with Polygonal Boundary

The governing differential equations of the bending problem of simply supported shallow spherical shells on Winkler foundation are simplified to an independent equation of radial deflection. The independent equation of radial deflection is decomposed to two Laplace operators by intermediate variable. The R-function theory is applied to describe a shallow spherical shell on Winkler foundation with concave boundary, and then a quasi-Green’s function is established by using the fundamental solution and the normalized boundary equation. The quasi-Green’s function satisfies the homogeneous boundary condition of the problem. The Laplace operators of the problem are reduced to two simultaneous Fredholm integral equations of the second kind by the Green’s formula. The singularity of the kernel of the integral equation is eliminated by choosing a suitable form of the normalized boundary equation. The integral equations are discretized into the homogeneous linear algebraic equations to proceed numerical computing. The singular term in the discrete equation is eliminated by the integral method. Some numerical examples are given to verify the validity of the proposed method in calculating simple boundary conditions and polygonal boundary conditions. A comparison with the ANSYS finite element (FEM) solution shows a good agreement, and it demonstrates the feasibility and efficiency of the present method

High-Order Numerical Solution of Second-Order One-Dimensional Hyperbolic Telegraph Equation Using a Shifted Gegenbauer Pseudospectral Method

We present a high-order shifted Gegenbauer pseudospectral method (SGPM) to solve numerically the second-order one-dimensional hyperbolic telegraph equation provided with some initial and Dirichlet boundary conditions. The framework of the numerical scheme involves the recast of the problem into its integral formulation followed by its discretization into a system of well-conditioned linear algebraic equations. The integral operators are numerically approximated using some novel shifted Gegenbauer operational matrices of integration. We derive the error formula of the associated numerical quadratures. We also present a method to optimize the constructed operational matrix of integration by minimizing the associated quadrature error in some optimality sense. We study the error bounds and convergence of the optimal shifted Gegenbauer operational matrix of integration. Moreover, we construct the relation between the operational matrices of integration of the shifted Gegenbauer polynomials and standard Gegenbauer polynomials. We derive the global collocation matrix of the SGPM, and construct an efficient computational algorithm for the solution of the collocation equations. We present a study on the computational cost of the developed computational algorithm, and a rigorous convergence and error analysis of the introduced method. Four numerical test examples have been carried out in order to verify the effectiveness, the accuracy, and the exponential convergence of the method. The SGPM is a robust technique, which can be extended to solve a wide range of problems arising in numerous applications.Comment: 36 pages, articl

Solution for a problem of linear plane elasticity with mixed boundary conditions on an ellipse by the method of boundary integrals

AbstractA numerical boundary integral scheme is proposed for the solution of the system of field equations of plane, linear elasticity in stresses for homogeneous, isotropic media in the domain bounded by an ellipse under mixed boundary conditions. The stresses are prescribed on one half of the ellipse, while the displacements are given on the other half. The method relies on previous analytical work within the Boundary Integral Method [1,2].The considered problem with mixed boundary conditions is replaced by two subproblems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way and the problem at this stage is reduced to the solution of a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution inside the domain, and the unknown boundary values of stresses or displacements on proper parts of the boundary.On the basis of the obtained results, it is inferred that the tangential stress component on the fixed part of the boundary has a singularity at each of the two separation points, thought to be of logarithmic type. A tentative form for the singular solution is proposed to calculate the full solution in bulk directly from the given boundary conditions using the well-known Boundary Collocation Method. It is shown that this addition substantially decreases the error in satisfying the boundary conditions on some interval not containing the singular points.The obtained results are discussed and boundary curves for unknown functions are provided, as well as three-dimensional plots for quantities of practical interest. The efficiency of the used numerical schemes is discussed, in what concerns the number of boundary nodes needed to calculate the approximate solution

Cracked elastic layer under compressive mechanical loads

AbstractWe consider boundary value problem in which an elastic layer containing a finite length crack is under compressive loading. The crack is parallel to the layer surfaces and the contact between crack surfaces are either frictionless or with adhesive friction or Coulomb friction.Based on fourier integral transformation techniques the solution of the formulated problems is reduced to the solution of a singular integral equation, then, using Chebyshev’s orthogonal polynomials, to an infinite system of linear algebraic equations. The regularity of these equations is established. The expressions for stress and displacement components in the elastic layer are presented. Based on the developed analytical algorithm, extensive numerical investigations have been conducted.The results of these investigations are illustrated graphically, exposing some novel qualitative and quantitative knowledge about the stress field in the cracked layer and their dependence on geometric and applied loading parameters. It can be seen from this study that the crack tip stress field has a mode II type singularity

Combining Nyström methods for a fast solution of Fredholm integral equations of the second kind

In this paper, we propose a suitable combination of two different Nyström methods, both using the zeros of the same sequence of Jacobi polynomials, in order to approximate the solution of Fredholm integral equations on [−1, 1]. The proposed procedure is cheaper than the Nyström scheme based on using only one of the described methods . Moreover, we can successfully manage functions with possible algebraic singularities at the endpoints and kernels with different pathologies. The error of the method is comparable with that of the best polynomial approximation in suitable spaces of functions, equipped with the weighted uniform norm. The convergence and the stability of the method are proved, and some numerical tests that confirm the theoretical estimates are given

Rigorous Approach to Analysis of Two-Dimensional Potential Problems, Wave Propagation and Scattering for Multi-conductor Systems

The research described in this chapter analyses two-dimensional potential problems for the multi-body systems, transverse electromagnetic wave propagation along multi-conductor transmission lines and two-dimensional plane wave scattering by various arrays. All conductors may be of arbitrary cross-sections; the only restriction on the system geometry is a smooth parameterization. These problems are mathematically modelled by Dirichlet boundary value problems for either the Laplace or the Helmholtz equation, with the classical integral representation of the solutions in the form of single-layer potential. The analytical-numerical algorithm presented here is based on the method of analytical regularization. The key idea behind this technique is an analytical transformation of the initial ill-posed integral equations to a well-conditioned Fredholm second kind matrix equation. The resulting system of infinite linear algebraic equations is effectively solved using the truncation method: the solution of the truncated system converges to the solution of the infinite system with the guaranteed accuracy that only depends on the truncation number and thus may be pre-specified. The solution obtained is applied to the accurate analysis of 2-D electrostatic- and electrodynamic-field problems for multi-conductor systems with arbitrary profiled conductors. Examples of some conceptual shielded transmission lines incorporating various configurations of conductors and scattering problems for the arrays of thick strips establish the utility of our method and its reliability in various situation

Homotopy Analysis And Legendre Multi-Wavelets Methods For Solving Integral Equations

Due to the ability of function representation, hybrid functions and wavelets have a special position in research. In this thesis, we state elementary definitions, then we introduce hybrid functions and some wavelets such as Haar, Daubechies, Cheby- shev, sine-cosine and linear Legendre multi wavelets. The construction of most wavelets are based on stepwise functions and the comparison between two categories of wavelets will become easier if we have a common construction of them. The properties of the Floor function are used to and a function which is one on the interval [0; 1) and zero elsewhere. The suitable dilation and translation parameters lead us to get similar function corresponding to the interval [a; b). These functions and their combinations enable us to represent the stepwise functions as a function of floor function. We have applied this method on Haar wavelet, Sine-Cosine wavelet, Block - Pulse functions and Hybrid Fourier Block-Pulse functions to get the new representations of these functions. The main advantage of the wavelet technique for solving a problem is its ability to transform complex problems into a system of algebraic equations. We use the Legendre multi-wavelets on the interval [0; 1) to solve the linear integro-differential and Fredholm integral equations of the second kind. We also use collocation points and linear legendre multi wavelets to solve an integro-differential equation which describes the charged particle motion for certain configurations of oscillating magnetic fields. Illustrative examples are included to reveal the sufficiency of the technique. In linear integro-differential equations and Fredholm integral equations of the second kind cases, comparisons are done with CAS wavelets and differential transformation methods and it shows that the accuracy of these results are higher than them. Homotopy Analysis Method (HAM) is an analytic technique to solve the linear and nonlinear equations which can be used to obtain the numerical solution too. We extend the application of homotopy analysis method for solving Linear integro- differential equations and Fredholm and Volterra integral equations. We provide some numerical examples to demonstrate the validity and applicability of the technique. Numerical results showed the advantage of the HAM over the HPM, SCW, LLMW and CAS wavelets methods. For future studies, some problems are proposed at the end of this thesis