On surface radiation conditions for high-frequency wave scattering

A new approximation of the logarithmic derivative of the Hankel function is derived and applied to high-frequency wave scattering. We re-derive the on surface radiation condition (OSRC) approximations that are well suited for a Dirichlet boundary in acoustics. These correspond to the Engquist–Majda absorbing boundary conditions. Inverse OSRC approximations are derived and they are used for Neumann boundary conditions. We obtain an implicit OSRC condition, where we need to solve a tridiagonal system. The OSRC approximations are well suited for moderate wave numbers. The approximation of the logarithmic derivative is also used for deriving a generalized physical optics approximation, both for Dirichlet and Neumann boundary conditions. We have obtained similar approximations in electromagnetics, for a perfect electric conductor. Numerical computations are done for different objects in 2D and 3D and for different wave numbers. The improvement over the standard physical optics is verified.     

An estimate of the heat fluxes to the surface of blunt bodies moving at hypersonic velocity in the atmosphere

The problem of the chemically equilibrium three-dimensional boundary layer on a blunt body which is in motion in the atmosphere is considered. A solution of the system of equations of the boundary layer is found by the method of successive approximations, and simple analytic expressions are written in the first approximation for the surface friction and heat flux coefficients. Formulae are obtained in the final form for estimating the convective heat flux in the neighbourhood of the critical point of spherical blunting.     

Convergence of the Neumann Series in BEM for the Neumann Problem of the Stokes System

A weak solution of the Neumann problem for the Stokes system in Sobolev space is studied in a bounded Lipschitz domain with connected boundary. A solution is looked for in the form of a hydrodynamical single layer potential. It leads to an integral equation on the boundary of the domain. Necessary and sufficient conditions for the solvability of the problem are given. Moreover, it is shown that we can obtain a solution of this integral equation using the successive approximation method. Then the consequences for the direct boundary integral equation method are treated. A solution of the Neumann problem for the Stokes system is the sum of the hydrodynamical single layer potential corresponding to the boundary condition and the hydrodynamical double layer potential corresponding to the trace of the velocity part of the solution. Using boundary behavior of potentials we get an integral equation on the boundary of the domain where the trace of the velocity part of the solution is unknown. It is shown that we can obtain a solution of this integral equation using the successive approximation method.     

Galerkin approximations with embedded boundary conditions for retarded delay differential equations

Finite-dimensional approximations are developed for retarded delay differential equations (DDEs). The DDE system is equivalently posed as an initial-boundary value problem consisting of hyperbolic partial differential equations (PDEs). By exploiting the equivalence of partial derivatives in space and time, we develop a new PDE representation for the DDEs that is devoid of boundary conditions. The resulting boundary condition-free PDEs are discretized using the Galerkin method with Legendre polynomials as the basis functions, whereupon we obtain a system of ordinary differential equations (ODEs) that is a finite-dimensional approximation of the original DDE system. We present several numerical examples comparing the solution obtained using the approximate ODEs to the direct numerical simulation of the original non-linear DDEs. Stability charts developed using our method are compared to existing results for linear DDEs. The presented results clearly demonstrate that the equivalent boundary condition-free PDE formulation accurately captures the dynamic behaviour of the original DDE system and facilitates the application of control theory developed for systems governed by ODEs.     

Asymptotic determination of the formation process of nonlinear distortion of one-dimensional pulses in a layered medium : PMM vol.40, n 6, 1976, pp. 1093–1103

Nonlinear effects in the propagation, reflection, and refraction of one-dimensional pulses in a medium consisting of two layers lying on a half-space are considered and analyzed. Properties of layers and of the half-space are different, and stresses are defined by an expansion in powers of strains. The initial pulse of finite duration is specified in the form of boundary condition at the surface of the external layer either for the deformation or for the dislocation rate, and the problem of wave pattern when the initial pulse amplitude tends to zero,i.e. in the case of small nonlinear effects, is solved.Problem is solved by the method of successive integration of nonhomogeneous linear wave equations, in which the solution of the linear problem is taken as the first approximation and the subsequent approximations are derived by approximating the nonlinear terms with the use of the preceding approximation.     

Approximate solution of boundary integral equations for biharmonic problems in non-smooth domains

This paper deals with approximate solutions to integral equations arising in boundary value problems for the biharmonic equation in simply connected piecewise smooth domains. The approximation method considered demonstrates excellent convergence even in the case of boundary conditions discontinuous at corner points. In an application we obtain very accurate approximations for some characteristics of two-dimensional Stokes flow in non-smooth domains. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An advanced study on the solution of nanofluid flow problems via Adomian’s method

Nanofluid flow is one of the most important areas of research at the present time due to its wide and significant applications in industry and several scientific fields. The boundary layer flow of nanofluids is usually described by a system of nonlinear differential equations with boundary conditions at infinity. These boundary conditions at infinity cause difficulties for any of the series method, such as Adomian’s method, the variational iteration method and others.The objective of the present work is to introduce a reliable method to overcome such difficulties that arise due to an infinite domain. The proposed scheme, that we will introduce, is based on Adomian’s decomposition method, where we will solve a system of nonlinear differential equations describing the boundary layer flow of a nanofluid past a stretching sheet.     

Convergence of Finite Volume Schemes for Hamilton-Jacobi Equations with Dirichlet Boundary Conditions

We study numerical methods for time-dependent Hamilton-Jacobi equations with weak Dirichlet boundary conditions. We first propose a new class of abstract monotone approximation schemes and get a convergence rate of 1/2 . Then, according to the abstract convergence results, by newly constructing monotone finite volume approximations on interior and boundary points, we obtain convergent finite volume schemes for time-dependent Hamilton-Jacobi equations with weak Dirichlet boundary conditions. Finally give some numerical results.     

非线性积分微分方程组奇摄动边值问题

本文讨论含积分算子的非线性微分方程组Ｒｏｂｉｎ边值问题的奇摄动,在适当假设条件下通过对角化技巧,利用逐步逼近法证明了解的存在,并得到直到Ｏ（ε＾Ｎ＋１）的按范数界限的一致有效估计。     

Methods of solving spatial problems of the mechanics of a deformable solid in terms of stresses

The formulation in /1/ of a quasistatic problem of the mechanics of a deformable solid in terms of stresses is discussed, including also the variational formulation, which consists of solving six equations in six symmetric stress tensor components when six boundary conditions are satisfied. Methods of successive approximation are proposed for solving this problem and theorems on the convergence of these methods, including a “rapidly converging” method, whose rate of convergence is substantially higher than a geometric progression, are proved.     

Acoustic Properties of Continuous and Discrete Models in Fluid Dynamics

The article investigates the propagation of small perturbations in fluids whose dynamics is described by Euler and Navier–Stokes equations or by a quasi-fluid-dynamic system derived from the difference approximation of the Boltzmann equation. The problem of wave reflection from the artificial boundaries of the numerical region is solved using various boundary conditions. The analysis is repeated for difference approximations of fluid-dynamic equations. The procedure is tested for viscous subsonic flow past a plate.     

Finite difference method in the problem of diffraction of a plane acoustic wave in a half-plane with a cut

In the numerical solution of the diffraction problem for an acoustic plane wave in a half-plane with a cut, boundary conditions that are equivalent to the radiation conditions at infinity are set in a neighborhood of the points of the cut. Joining the physical boundary conditions on the cut, a closing set of equations of order 4N, where N is the number of grid points on the cut, is obtained. The so-called Green’s grid function for the half-plane is used, which makes it possible to pass from one grid layer to another one for the solution satisfying certain conditions at infinity.     

Numerical analytic method for the case of singular matrices in boundary conditions

A modification of the numerical analytic method of successive approximations is proposed for investigating the existence and constructing solutions of systems of nonlinear ordinary differential equations with two-point linear boundary conditions. The method is tailored to the case where the matrices contained in boundary conditions are singular, but their linear combination is nonsingular.Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 44, No. 5, pp. 673–681, May, 1992.     

Hydrostatic Stokes equations with non-smooth data for mixed boundary conditions

The main subject of this work is to study the concept of very weak solution for the hydrostatic Stokes system with mixed boundary conditions (non-smooth Neumann conditions on the rigid surface and homogeneous Dirichlet conditions elsewhere on the boundary). In the Stokes framework, this concept has been studied by Conca [Rev. Mat. Apl. 10 (1989)] imposing non-smooth Dirichlet boundary conditions.In this paper, we introduce the dual problem that turns out to be a hydrostatic Stokes system with non-free divergence condition. First, we obtain strong regularity for this dual problem (which can be viewed as a generalisation of the regularity results for the hydrostatic Stokes system with free divergence condition obtained by Ziane [Appl. Anal. 58 (1995)]). Afterwards, we prove existence and uniqueness of very weak solution for the (primal) problem.As a consequence of this result, the existence of strong solution for the non-stationary hydrostatic Navier-Stokes equations is proved, weakening the hypothesis over the time derivative of the wind stress tensor imposed by Guillén-González, Masmoudi and Rodríguez-Bellido [Differential Integral Equations 50 (2001)].     

A uniformly convergent hybrid scheme for singularly perturbed system of reaction-diffusion Robin type boundary-value problems

This paper deals with the study on system of reaction diffusion differential equations for Robin or mixed type boundary value problems (MBVPs). A cubic spline approximation has been used to obtain the difference scheme for the system of MBVPs, on a piecewise uniform Shishkin mesh defined in the whole domain. It has been shown that our proposed scheme, i.e., central difference approximation for outer region with cubic spline approximation for inner region of boundary layers, leads to almost second order parameter uniform convergence whereas the standard method i.e., the forward-backward approximation for mixed boundary conditions with central difference approximation inside the domain leads to almost first order convergence on Shishkin mesh. Numerical results are provided to show the efficiency and accuracy of these methods.     

The method of successive approximations for solving the problem of the decay of a discontinuity and its application to the equations of two-layer shallow water theory

The method of successive approximations for solving the problem of the decay of a small amplitude discontinuity is proposed for hyperbolic systems of conservation laws. In the linear approximation, a Cauchy problem for a linear hyperbolic system is obtained. Its solution represents lines of discontinuity separating the regions in which the solution is constant. Most attention is paid to the first and second approximations, within the limits of which the discontinuities obtained in the first approximation decay into stable shock waves and rarefaction waves. An analysis of the qualitatively different flow conditions that arise when solving the problem of the failure of a dam for a two-layer shallow water model with a free boundary is presented as an example.     

Near-circularity of the error curve in complex Chebyshev approximation

Let f(z) be analytic on the unit disk, and let p^*(z) be the best (Chebyshev) polynomial approximation to f(z) on the disk of degree at most n. It is observed that in typical problems the “error curve,” the image of the unit circle under (f − p^*)(z), often approximates to a startling degree a perfect circle with winding number n + 1. This phenomenon is approached by consideration of related problems whose error curves are exactly circular, making use of a classical theorem of Carathéodory and Fejér. This leads to a technique for calculating approximations in one step that are roughly as close to best as the best approximation error curve is close to circular, and hence to strong theorems on near-circularity as the radius of the domain shrinks to 0 or as n increases to ∞. As a computational example, very tight bounds are given for approximation of e^z on the unit disk. The generality of the near-circularity phenomenon (more general domains, rational approximation) is discussed.     

Continuum models in problems of the hypersonic flow of a rarefied gas around blunt bodiest

All possible continuum (hydrodynamic) models in the case of two-dimensional problems of supersonic and hypersonic flows around blunt bodies in the two-layer model (a viscous shock layer and shock-wave structure) over the whole range of Reynolds numbers, Re, from low values (free molecular and transitional flow conditions) up to high values (flow conditions with a thin leading shock wave, a boundary layer and an external inviscid flow in the shock layer) are obtained from the Navier-Stokes equations using an asymptotic analysis. In the case of low Reynolds numbers, the shock layer is considered but the structure of the shock wave is ignored. Together with the well-known models (a boundary layer, a viscous shock layer, a thin viscous shock layer, parabolized Navier-Stokes equations (the single-layer model) for high, moderate and low Re numbers, respectively), a new hydrodynamic model, which follows from the Navier-Stokes equations and reduces to the solution of the simplified (“local”) Stokes equations in a shock layer with vanishing inertial and pressure forces and boundary conditions on the unspecified free boundary (the shock wave) is found at Reynolds numbers, and a density ratio, k, up to and immediately after the leading shock wave, which tend to zero subject to the condition that (k/Re)^1/2 → 0. Unlike in all the models which have been mentioned above, the solution of the problem of the flow around a body in this model gives the free molecular limit for the coefficients of friction, heat transfer and pressure. In particular, the Newtonian limit for the drag is thereby rigorously obtained from the Navier-Stokes equations. At the same time, the Knudsen number, which is governed by the thickness of the shock layer, which vanishes in this model, tends to zero, that is, the conditions for a continuum treatment are satisfied. The structure of the shock wave can be determined both using continuum as well as kinetic models after obtaining the solution in the viscous shock layer for the weak physicochemical processes in the shock wave structure itself. Otherwise, the problem of the shock wave structure and the equations of the viscous shock layer must be jointly solved. The equations for all the continuum models are written in Dorodnitsyn--Lees boundary layer variables, which enables one, prior to solving the problem, to obtain an approximate estimate of second-order effects in boundary-layer theory as a function of Re and the parameter k and to represent all the aerodynamic and thermal characteristic; in the form of a single dependence on Re over the whole range of its variation from zero to infinity.

An efficient numerical method of global iterations, previously developed for solving viscous shock-layer equations, can be used to solve problems of supersonic and hypersonic flows around the windward side of blunt bodies using a single hydrodynamic model of a viscous shock layer for all Re numbers, subject to the condition that the limit (k/Re)^1/2 → 0 is satisfied in the case of small Re numbers. An aerodynamic and thermal calculation using different hydrodynamic models, corresponding to different ranges of variation Re (different types of flow) can thereby, in fact, be replaced by a single calculation using one model for the whole of the trajectory for the descent (entry) of space vehicles and natural cosmic bodies (meteoroids) into the atmosphere.     

Generalized rational approximation in interpolating subspaces of Lp(μ)

We extend the geometric approach of Cheney and Loeb in [2] to the problem of approximation in L_p(μ) by “admissable” generalized rational functions. We obtain a characterization for locally best approximations and find the interpolating condition sufficient for their local unicity. Our results are comparable to those for the linear approximation problem as investigated by Singer and Ault, Deutsch, Morris, and Olson.     

Application of homotopy analysis method for natural convection of Darcian fluid about a vertical full cone embedded in pours media prescribed surface heat flux

In this paper, a powerful analytical method, called homotopy analysis method (HAM) is used to obtain the analytical solution for a nonlinear ordinary deferential equation that often appear in boundary layers problems arising in heat and mass transfer which these kinds of the equations contain infinity boundary condition. The boundary layer approximations of fluid flow and heat transfer of vertical full cone embedded in porous media give us the similarity solution for full cone subjected to surface heat flux boundary conditions. Nonlinear ODE which is obtained by similarity solution has been solved through homotopy analysis method (HAM). The main objective is to propose alternative methods of solution, which do not require small parameters and avoid linearization and physically unrealistic assumptions. The obtained analytical solution in comparison with the numerical ones represents a remarkable accuracy. The results also indicate that HAM can provide us with a convenient way to control and adjust the convergence region.