Application of Adapted-Bubbles to the Helmholtz Equation with Large Wavenumbers in 2D

An adapted-bubbles approach which is a modification of the residual-free bubbles (RFB) method, is proposed for the Helmholtz problem in 2D. A new two-level finite element method is introduced for the approximations of the bubble functions. Unlike the other equations such as the advection-diffusion equation, RFB method when applied to the Helmholtz equation, does not depend on another stabilized method to obtain approximations to the solutions of the sub-problems. Adapted-bubbles (AB) are obtained by a simple modification of the sub-problems. This modification increases the accuracy of the numerical solution impressively. We provide numerical experiments with the AB method up to $ch = 5$ where $c$ is the wavenumber and $h$ is the mesh size. Numerical tests show that the AB method is better by far than higher order methods available in the literature.     

A new approach to the classical Stokes flow problem: Part I Methodology and first-order analytical results

The problem of determining the axisymmetric Stokes flow past an arbitrary body, the boundary shape of which can be represented by an analytic function, is examined by developing an exact method. An appropriate nonorthogonal coordinate system is introduced, and it is shown that the Hilbert space to which the stream function belongs is spanned by the set of Gegenbauer polynomials based on the physical argument that the drag on a body should be finite. The partial differential equation of the original problem is then reduced to two simultaneous vector differential equations. By the truncation of this infinite-dimensional system to the one-dimensional subspace, an explicit analytic solution to the Stokes equation valid for all bodies in question is obtained as a first approximation.     

Global existence to generalized Boussinesq equation with combined power-type nonlinearities

The Cauchy problem to the generalized Boussinesq equation with combined power-type nonlinearities is studied. Global solvability or finite time blow-up of the solutions with subcritical initial energy is proved by means of the sign preserving property of the Nehari functional. For generalized Lienard (or generalized Bernoulli) nonlinear terms the critical energy constant is explicitly evaluated. A new method, that can be considered as a modification of the potential well method, is developed. The performed numerical experiments support the theoretical results.     

Hertzian contact problem: Numerical reduction and volumetric modification

A technique for the analytical formulation and numerical implementation of an elastic contact model for rigid bodies in the framework of the Hertzian contact problem is described. The normal elastic force and the semiaxes of the contact area are computed so that the problem is sequentially reduced to a scalar transcendental equation depending on complete elliptic integrals of the first and second kinds. Based on the classical solution to the Hertzian contact problem, an invariant volumetric force function is proposed that depends on the geometric characteristics of interpenetration of two undeformed bodies. The normal forces computed using the force function agree with results obtained previously for non-Hertzian contact of elastic bodies. As an example, a ball bearing is used to compare the contact dynamics of elastic bodies simulated in the classical Hertzian model and its volumetric modification.     

A modification of the invariant imbedding method for a singular boundary value problem

We consider a dynamically-consistent analytical model of a 3D topographic vortex. The model is governed by equations derived from the classical problem of the axisymmetric Taylor–Couette flow. Using linear expansions, these equations can be reduced to a differential sixth-order equation with variable coefficients. For this differential equation, we formulate a boundary value problem, which has a number of issues for numerical solving. To avoid these issues and find the eigenvalues and eigenfunctions of the boundary value problem, we suggest a modification of the invariant imbedding method (the Riccati equation method). In this paper, we show that such a modification is necessary since the boundary conditions possess singular matrices, which sufficiently complicate the derivation of the Riccati equation. We suggest algebraic manipulations, which permit the initial problem to be reduced to a problem with regular boundary conditions. Also, we propose a method for obtaining a numerical solution of the matrix Riccati equation by means of recurrence relations, which allow us to obtain a matrizer converging to the required eigenfunction. The suggested method is tested by calculating the corresponding eigenvalues and eigenfunctions, and then, by constructing fluid particle trajectories on the basis of the eigenfunctions.     

PSE在可压缩边界层转捩问题中的应用

提出了用抛物化稳定性方程（PSE）预测层流-湍流转捩的一种新的概念．它被试用于平板可压缩边界层转捩位置的预测问题中,并将结果与直接数值模拟（DNS）所得进行比较．结果二者符合情况令人满意．二者符合的原因在于PSE方法准确地再现了层流-湍流转捩中导致breakdown过程的机理,即平均流剖面的修正导致其稳定性特性的明显改变．     

Transfer matrix method for the free and forced vibration analyses of multi-step Timoshenko beams coupled with rigid bodies on springs

Multi-step Timoshenko beams coupled with rigid bodies on springs can be regarded as a generalized model to investigate the dynamic characteristics of many structures and mechanical systems in engineering. This paper presents a novel transfer matrix method for the free and forced vibration analyses of the hybrid system. It is modeled as a chain system, where each beam and each rigid body with its supporting spring are dealt with one element, respectively. The transfer equation of each element is deduced based on separation of variables method. The system overall transfer equation is obtained by substituting an element transfer equation into another. Then, the free vibration characteristics are acquired by solving exact homogeneous linear equations. To compute the forced vibration response with modal superposition method, the body dynamic equations and augmented eigenvectors are established, and the orthogonality of augmented eigenvectors is mathematically proved. Without high-order global dynamic equation or approximate spatial discretization, the free and forced vibration analyses of the hybrid system are achieved efficiently and accurately in this study. As an analytical approach, the present method is easy, highly stylized, robust, powerful and general for the complex hybrid systems containing any number of Timoshenko beams and rigid bodies. Four numerical examples are implemented, and the results show that this method is computationally efficient with high precision.     

An analytical approximation for solving nonlinear Blasius equation by NHPM

In this article, an analytic approximation to the solution of Blasius equation is obtained by using a new modification of homotopy perturbation method. The Blasius equation is a nonlinear ordinary differential equation which arises in the boundary layer flow. The comparison with Howart's numerical solution shows that the new homotopy perturbation method is an effective mathematical method with high accuracy. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Propagation of a plane wave to a materially uniform but inhomogeneous body

We produce the equations of small deformations superimposed upon large for materially uniform but inhomogeneous bodies and specialize to an isotropic material and to a homogeneous finite elastic deformation. By assuming the small deformation to be a plane wave, a set of equations for the amplitude of the wave is produced which is accompanied by an additional set of conditions. By requiring a non-trivial solution for the amplitude, we obtain the secular equation and from it a set of necessary and sufficient conditions for having a real wave speed. The second set of conditions that have to be satisfied is due to the materials inhomogeneity. Essentially, the present analysis enhances the approach of Hayes and Rivlin for materially uniform but inhomogeneous bodies. The outcome is that for such bodies the restrictions on the constitutive law for having real wave speeds for an isotropic material subjected to a pure homogeneous deformation involves the field of the inhomogeneity as well.     

A method of solving the integral equation of plane contact problems for semi-bounded bodies

A typical integral equation, which arises when solving linear plane contact problems for semi-bounded bodies, is considered. By using a special representation of the kernel of this equation, an approximate method is developed for solving it that is effective over a wide range of variation of the dimensionless geometrical parameter occurring in the kernel. The method is tested on the problem of the symmetrical compression of an elastic strip along its boundaries by two similar punches.     

Modification of Block Pulse Functions and their application to solve numerically Volterra integral equation of the first kind

In this paper a modification of Block Pulse Functions is introduced and used to solve Volterra integral equation of the first kind. Some theorems are included to show convergence and advantage of the method. Some examples show accuracy of the method.     

Three-dimensional generalization for <Emphasis Type="Italic">W</Emphasis> modification of a Godunov method

A high-accuracy modification of Godunov’s method for three-dimensional unsteady ideal gas flows is proposed. For the linear advection equation, a fully three-dimensional second-order accurate monotone scheme is designed with corrections computed on a variable stencil whose orientation depends on the signs of the equation coefficients. For the linear scalar advection equation, the scheme is proved to possess the positive approximation property. The method is tested by computing the flow in a three-dimensional Ludwieg tube with a square cross section.     

通过可渗透近球体的轴对称流动

提出了可渗透近球体轴对称流动的分析方法.用修正边界条件的办法反映可渗透性.用正规摄动法求解了Stokes方程,达到ε的2阶修正.ε是描述不变形球体半径偏差的小参数.计算了阻力和流量,并从几何方面和表面渗透性方面考查了计算结果.还尝试将此理论应用于过滤供水问题.小型的生态学上重要的水生生物体的过滤器,被模型化为轴对称可渗透物体,用扁球体或近球体建立了该问题的初级模型.     

Economic difference scheme for a parabolic equation with a mixed spatial derivative

A modification of a well-known locally one-dimensional method for a parabolic equation is proposed. The method remains economic even if the equation involves a mixed derivative with respect to spatial variables. A model case study of the method is presented. Numerical results are given that demonstrate the efficiency of the method.     

Green's function of the Brinkman equation in a 2D anisotropic case

The purpose of this paper is to obtain the Green's functionof the Brinkman equation in a 2D case of hydrodynamic anisotropywith respect to the permeability. The anisotropic nature ofthe permeability is assumed to be not space or time dependent.We use the method of Fourier transform which reduces the computationof the Green's function to the computation of the fundamentalsolution of a fourth-order partial differential equation. Thisresearch work has several applications in engineering and medicineto the motion of bodies in anisotropic porous media.     

The successive approximations method for the airfoil equation

The successive approximations (or Neumann iterations) method for the solution of Fredholm integral equations of the second kind is applied here for the first time, after an appropriate modification, to a Cauchy-type singular integral equation of the first kind, the airfoil equation. The convergence of the method is investigated and three simple applications are made. The numerical implementation of the method (by using Gaussian quadrature rules) is also described in detail and numerical results verifying the accuracy and convergence of the method are displayed.     

Application of new advanced CNN structure with adaptive thresholds to color edge detection

Color edge detection is much more efficient than gray scale detection when edges exist at the boundary between regions of different colors with no change in intensity. This paper presents adaptive templates, which are capable of detecting various color and intensity changes in color image. To avoid conception of multilayer proposed in literatures, modification has been done to the CNN structure. This modified structure allows a matrix C, which carries the change information of pixels, to replace the control parts in the basic CNN equation. This modification is necessary because in multilayer structure, it faces the challenge of how to represent the intrinsic relationship among each primary layer. Additionally, in order to enhance the accuracy of edge detection, adaptive detection threshold is employed. The adaptive thresholds are considered to be alterable criteria in designing matrix C. The proposed synthetic system not only avoids the problem which is engendered by multi-layers but also exploits full information of pixels themselves. Experimental results prove that the proposed method is efficient.     

Numerical solution of a secular equation

Summary. A method is proposed for the solution of a secular equation, arising in modified symmetric eigenvalue problems and in several other areas. This equation has singularities which make the application of standard root-finding methods difficult. In order to solve the equation, a class of transformations of variables is considered, which transform the equation into one for which Newton's method converges from any point in a certain given interval. In addition, the form of the transformed equation suggests a convergence accelerating modification of Newton's method. The same ideas are applied to the secant method and numerical results are presented. Received July 1, 1994     

A new spectral-homotopy analysis method for solving a nonlinear second order BVP

A modification of the homotopy analysis method (HAM) for solving nonlinear second-order boundary value problems (BVPs) is proposed. The implementation of the new approach is demonstrated by solving the Darcy–Brinkman–Forchheimer equation for steady fully developed fluid flow in a horizontal channel filled with a porous medium. The model equation is solved concurrently using the standard HAM approach and numerically using a shooting method based on the fourth order Runge–Kutta scheme. The results demonstrate that the new spectral homotopy analysis method is more efficient and converges faster than the standard homotopy analysis method.     

On the linearized theory op flow around bodies. Method of force sources

The method of force sources is proposed for solving linear problems related to the interaction between rigid bodies, and fluids, or gases. Method is based on the introduction of perturbation force sources into equation of motion of fluid media. Boundary conditions at the rigid body surface make it possible to reduce the problem of hydrodynamic reactions to an integral equation defining the function of force sources. Method is illustrated by the solution of three simple problems in the field of acoustics, and of viscous, and compressible media flow around bodies.

In the linearized theory of flow around rigid bodies, as well as in acoustics, an important part of the sound wave generation analysis concerns the determination of hydrodynamic reactions of the medium on moving, pulsating, or oscillating bodies. Such reactions make themselves felt as constant, or variable mechanical forces, such as drag and lift, or in the case of sound wave emitters, as the wave resistance. Various methods had been proposed for the computation of such forces, as for example, in the monographs [1 to 6].

Here, a different approach to the problem of determination of surface forces exerted by liquids and gases on the rigid body is proposed. By resorting to the formalism of the generalized functions it is possible to introduce into the equations of motion of fluid media a perturbation source in the form volume density of forces exercised by the body on the gas. The distribution of surface tension entering into the expression of this force is selected in such a manner as to satisfy boundary conditions at the body surface. It becomes possible with the use of this device to reduce the problem of determination of forces acting on the body surface to the solution of certain Integral equations. The proposed method is in all respects completely analogous to the well-known method of sources and sinks [1 to 1]. Both methods reduce the problem of interaction between body and gas to the solution of Integral equations. The method of sources and sinks, however, leads to an integral equation which describes the distribution of fictitious sources and sinks in the volume of the body having the density of the medium, while the method of force sources yields an integral equation which directly defines the distribution of mechanical forces over the surface of the body (*).

We may note that the method of force sources had to a certain extent been already used in papers [6 and 7] for the determination of sound radiation by means of point-force sources.