Singular boundary method for 3D time-harmonic electromagnetic scattering problems

A frequency domain singular boundary method is presented for solving 3D time-harmonic electromagnetic scattering problem from perfect electric conductors. To avoid solving the coupled partial differential equations with fundamental solutions involving hypersingular terms, we decompose the governing equation into a system of independent Helmholtz equations with mutually coupled boundary conditions. Then the singular boundary method employs the fundamental solutions of the Helmholtz equations to approximate the scattered electric field variables. To desingularize the source singularity in the fundamental solutions, the origin intensity factors are introduced. In the novel formulation, only the origin intensity factors for fundamental solutions of 3D Helmholtz equations and its derivatives need to be considered which have been derived in the paper. Several numerical examples involving various perfectly conducting obstacles are carried out to demonstrate the validity and accuracy of the present method.     

Analytical solutions of a coupled nonlinear system arising in a flow between stretching disks

In this paper, we consider the axi-symmetric flow between two infinite stretching disks. By using a similarity transformation, we reduce the governing Navier-Stokes equations to a system of nonlinear ordinary differential equations. We first obtain analytical solutions via a four-term perturbation method for small and large values of the Reynolds number R. Also, we apply the Homotopy Analysis Method (which may be used for all values of R) to obtain analytical solutions. These solutions converge over a larger range of values of the Reynolds number than the perturbation solutions. Our results agree well with the numerical results of Fang and Zhang [22]. Furthermore, we obtain the analytical solutions valid for moderate values of R by use of Homotopy Analysis.     

Lie group symmetry analysis for granular media stress equations

The Airy stress function, although frequently employed in classical linear elasticity, does not receive similar usage for granular media problems. For plane strain quasi-static deformations of a cohesionless Coulomb-Mohr granular solid, a single nonlinear partial differential equation is formulated for the Airy stress function by combining the equilibrium equations with the yield condition. This has certain advantages from the usual approach, in which two stress invariants and a stress angle are introduced, and a system of two partial differential equations is needed to describe the flow. In the present study, the symmetry analysis of differential equations is utilised for our single partial differential equation, and by computing an optimal system of one-dimensional Lie algebras, a complete set of group-invariant solutions is derived. By this it is meant that any group-invariant solution of the governing partial differential equation (provided it can be derived via the classical symmetries method) may be obtained as a member of this set by a suitable group transformation. For general values of the parameters (angle of internal friction ? and gravity g) it is found there are three distinct classes of solutions which correspond to granular flows considered previously in the literature. For the two limiting cases of high angle of internal friction and zero gravity, the governing partial differential equation admit larger families of Lie point symmetries, and from these symmetries, further solutions are derived, many of which are new. Furthermore, the majority of these solutions are exact, which is rare for granular flow, especially in the case of gravity driven flows.     

Solutions of 2nd-order linear differential equations subject to Dirichlet boundary conditions in a Bernstein polynomial basis

An algorithm for approximating solutions to 2nd-order linear differential equations with polynomial coefficients in B-polynomials (Bernstein polynomial basis) subject to Dirichlet conditions is introduced. The algorithm expands the desired solution in terms of B-polynomials over a closed interval [0, 1] and then makes use of the orthonormal relation of B-polynomials with its dual basis to determine the expansion coefficients to construct a solution. Matrix formulation is used throughout the entire procedure. However, accuracy and efficiency are dependent on the size of the set of B-polynomials, and the procedure is much simpler compared to orthogonal polynomials for solving differential equations. The current procedure is implemented to solve five linear equations and one first-order nonlinear equation, and excellent agreement is found between the exact and approximate solutions. In addition, the algorithm improves the accuracy and efficiency of the traditional methods for solving differential equations that rely on much more complicated numerical techniques. This procedure has great potential to be implemented in more complex systems where there are no exact solutions available except approximations.     

Similarity solutions for flow and heat transfer over a permeable surface with convective boundary condition

The steady laminar boundary layer flow over a permeable flat plate in a uniform free stream, with the bottom surface of the plate is heated by convection from a hot fluid is considered. Similarity solutions for the flow and thermal fields are possible if the mass transpiration rate at the surface and the convective heat transfer from the hot fluid on the lower surface of the plate vary like x^−1/2, where x is the distance from the leading edge of the solid surface. The governing partial differential equations are first transformed into ordinary differential equations, before being solved numerically. The effects of the governing parameters on the flow and thermal fields are thoroughly examined and discussed.     

车路耦合非线性振动高阶Galerkin截断研究

将移动车辆模型化为运动的两自由度质量-弹簧-阻尼系统,道路模型化为立方非线性黏弹性地基上的弹性梁,并将路面不平度设定为简谐函数．通过受力分析,建立车路非线性耦合振动高阶偏微分方程．采用高阶Galerkin截断结合数值方法求解耦合系统的动态响应．首次研究不同截断阶数对车路耦合非线性振动动态响应的影响,确定Galerkin截断研究车路耦合振动的收敛性．研究结果表明,对于软土地基的沥青路面,耦合振动的动态响应,需要150阶以上的截断才能达到收敛效果．并通过高阶收敛的Galerkin截断研究了系统参数对车路耦合非线性振动动态响应的影响.     

A general solution for three dimensional steady-state transversely isotropic hygro-thermo-magneto-piezoelectric media

In this article, we extract the general solution of three dimensional (3D) equations using potential theory method (PTM) for steady-state, transversely isotropic, hygro-thermo-magneto-piezoelectric media (HTMPM). The governing equations are simplified by introducing the displacement functions. A general solution is completely determined by advantage of the superposition principle and operator theory, which is connected in terms of two functions, fulfilling a second-order and twelfth-order homogeneous partial differential equation (PDE), separately. With the help of Almansi’s theorem, the general solution can be further shortened, which is stated by seven harmonic functions only. The acquired general solutions are straightforward structure and helpful in boundary value problems of HTMPM. Further, we apply the 3D fundamental solutions inside an infinite and on the surface of semi-infinite of a steady point heat source united with a steady point moisture source transversely isotropic HTMPM. Comprehensive and exact solutions are given in the form of elementary functions, which appear as a standard for various types of approximate solutions and numerical codes. Some numerical simulation is conducted based on the obtained general solutions.     

Similarity solutions to nonlinear heat conduction and Burgers/Korteweg–deVries fractional equations

We analyze self-similar solutions to a nonlinear fractional diffusion equation and fractional Burgers/Korteweg–deVries equation in one spatial variable. By using Lie-group scaling transformation, we determined the similarity solutions. After the introduction of the similarity variables, both problems are reduced to ordinary nonlinear fractional differential equations. In two special cases exact solutions to the ordinary fractional differential equation, which is derived from the diffusion equation, are presented. In several other cases the ordinary fractional differential equations are solved numerically, for several values of governing parameters. In formulating the numerical procedure, we use special representation of a fractional derivative that is recently obtained.     

A model of multiphase flow and transport in porous media applied to gas migration in underground nuclear waste repository

We prove existence of solutions for a new model of two compressible and partially miscible phase flow in porous media, applied to gas migration in an underground nuclear waste repository. This model, modeling fully and partially water saturated situations, consist of a coupled system of quasilinear parabolic partial differential equations. We seek a new set of variables in order to obtain a system which belongs to the class of equations considered by Alt and Luckhaus such that it would be possible to use their existence theorem. A simulation of a numerical test case is performed in order to numerically demonstrate the ability of this model to take in account the appearance of one phase. To cite this article: F. Smaï, C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

Coupled vortex equations and moduli: deformation theoretic approach and Kähler geometry

We investigate differential geometric aspects of moduli spaces parametrizing solutions of coupled vortex equations over a compact Kähler manifold X. These solutions are known to be related to polystable triples via a Kobayashi–Hitchin type correspondence. Using a characterization of infinitesimal deformations in terms of the cohomology of a certain elliptic double complex, we construct a Hermitian structure on these moduli spaces. This Hermitian structure is proved to be Kähler. The proof involves establishing a fiber integral formula for the Hermitian form. We compute the curvature tensor of this Kähler form. When X is a Riemann surface, the holomorphic bisectional curvature turns out to be semi-positive. It is shown that in the case where X is a smooth complex projective variety, the Kähler form is the Chern form of a Quillen metric on a certain determinant line bundle.     

Nonlinear analysis of thin rectangular plates on Winkler–Pasternak elastic foundations by DSC–HDQ methods

This article introduces a coupled methodology for the numerical solution of geometrically nonlinear static and dynamic problem of thin rectangular plates resting on elastic foundation. Winkler–Pasternak two-parameter foundation model is considered. Dynamic analogues Von Karman equations are used. The governing nonlinear partial differential equations of the plate are discretized in space and time domains using the discrete singular convolution (DSC) and harmonic differential quadrature (HDQ) methods, respectively. Two different realizations of singular kernels such as the regularized Shannon’s kernel (RSK) and Lagrange delta (LD) kernel are selected as singular convolution to illustrate the present DSC algorithm. The analysis provides for both clamped and simply supported plates with immovable inplane boundary conditions at the edges. Various types of dynamic loading, namely a step function, a sinusoidal pulse, an N-wave pulse, and a triangular load are investigated and the results are presented graphically. The effects of Winkler and Pasternak foundation parameters, influence of mass of foundation on the response have been investigated. In addition, the influence of damping on the dynamic analysis has been studied. The accuracy of the proposed DSC–HDQ coupled methodology is demonstrated by the numerical examples.     

Coupled nonlinear stochastic differential equations

Systems of n coupled linear or nonlinear differential equations which may be deterministic or stochastic are solved by methods of the first author and his co-workers. Examples include two coupled Riccati equations, coupled linear equations, stochastic coupled equations with product terms, and n coupled stochastic differential equations.     

Stagnation-point flow and heat transfer over an exponentially shrinking sheet

An analysis is carried out to investigate the stagnation-point flow and heat transfer over an exponentially shrinking sheet. Using the boundary layer approximation and a similarity transformation in exponential form, the governing mathematical equations are transformed into coupled, nonlinear ordinary differential equations which are then solved numerically by a shooting method with fourth order Runge-Kutta integration scheme. The analysis reveals that a solution exists only when the velocity ratio parameter satisfies the inequality −1.487068 ? c/a. Also, the numerical calculations exhibit the existence of dual solutions for the velocity and the temperature fields; and it is observed that their boundary layers are thinner for the first solution (in comparison with the second). Moreover, the heat transfer from the sheet increases with an increase in c/a for the first solution, while the heat transfer decreases with increasing c/a for the second solution, and ultimately heat absorption occurs.     

Asymptotic solutions and stability analysis for generalized non-homogeneous Mathieu equation

The asymptotic solutions and transition curves for the generalized form of the non-homogeneous Mathieu differential equation are investigated in this paper. This type of governing differential equation of motion arises from the dynamic behavior of a pendulum undergoing a butterfly-type end support motion. The strained parameter technique is used to obtain periodic asymptotic solutions. The transition curves for some special cases are presented and their corresponding periodic solutions with the periods of 2π and 4π are evaluated. The stability analyses of those transition curves in the ε–δ plane are carried out, analytically, using the multiple scales method. The numerical simulations for some typical points in the ε–δ plane are performed and the dynamic characteristics of the resulting phase plane trajectories are discussed.     

Analyzing of the effect of the dusty fluid,anisotropy and layered of the wall on the wave propagation

This study considers the propagation of time harmonic waves in, prestressed, anisotropic elastic tubes filled with viscous fluid containing dusty particles. The fluid is assumed to be incompressible and Newtonian. The tube material is considered to be incompressible, anisotropic, and elastic. The tube is subjected to a static inner pressure P_i and an axial stretch λ. Utilizing the theory of “Superposing small deformations on large initial static deformations”, differential equations governing wave propagation inside the tube are obtained in terms of cylindrical coordinates. Analytical solutions for the equations of motion for the dust and the fluid are obtained, and expressed numerically. The dispersion relation is obtained as a function of the stretch, the thickness ratio and the parameters for dusty particles.     

Propagation of harmonic waves in prestressed fiber viscoelastic thick tubes filled with a viscous dusty fluid

In this work, propagation of harmonic waves in initially stressed cylindrical viscoelastic thick tubes filled with a Newtonian fluid is studied. The tube, subjected to a static inner pressure P_i and a positive axial stretch λ, will be considered as an incompressible viscoelastic and fibrous material. The fluid is assumed as an incompressible, viscous and dusty fluid. The field equations for the fluid are obtained in the cylindrical coordinates. The governing differential equations of the tube’s viscoelastic material are obtained also in the cylindrical coordinates utilizing the theory of small deformations superimposed on large initial static deformations. For the axially symmetric motion the field equations are solved by assuming harmonic wave solutions. A closed form solution can be obtained for equations governing the fluid body, but due to the variability of the coefficients of resulting differential equations of the solid body, such a closed form solution is not possible to obtain. For that reason, equations for the solid body and the boundary conditions are treated numerically by the finite-difference method to obtain the effects of the thickness of the tube on the wave characteristics. Dispersion relation is obtained using the long wave approximation and, the wave velocities and the transmission coefficients are computed.     

Convective heat transfer in a conducting fluid over a permeable stretching surface with suction and internal heat generation/absorption

We consider a convective flow in a porous medium of an incompressible viscous conducting fluid impinging on a permeable stretching surface with suction, and internal heat generation/absorption. Using a similarity transformation the governing equations of the problem are reduced to a coupled third-order nonlinear ordinary differential equations. We first examine a number of special cases for which we may obtain exact solutions. We then obtain analytical solutions (by the Homotopy Analysis Method) and numerical solutions (by a boundary value problem solver), in order to further study the behavior of the nonlinear differential equations, for various values of the physical parameters. Our numerical solutions are shown to agree with the available results in the literature. We then employ the numerical results to bring out the effects of the suction parameter, heat source/sink parameter, stretching parameter, porosity parameter, the Prandtl number and the free convection parameter on the flow and heat transfer characteristics. In the absence of suction and free convection, our findings are in agreement with the corresponding numerical results of Attia [H.A. Attia, On the effectiveness of porosity on stagnation point flow towards a stretching surface with heat generation, Comput. Mater. Sci. 38 (2007) 741-745].     

Poroelastodynamic Boundary Element Method in Time Domain: Numerical Aspects

Based on Biot's theory the governing equations for a poroelastic continuum are given as a coupled set of partial differential equations (PDEs) for the unknowns solid displacements and pore pressure. Using the Convolution Quadrature Method (CQM) proposed by Lubich a boundary time stepping procedure is established based only on the fundamental solutions in Laplace domain. To improve the numerical behavior of the CQM-based Boundary Element Method (BEM) dimensionless variables are introduced and different choices studied. This will be performed as a numerical study at the example of a poroelastic column. Summarizing the results, the normalization to time and spatial variable as well as on Young's modulus yields the best numerical behavior. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mixed convection in the stagnation-point flow of a Maxwell fluid towards a vertical stretching surface

In the present analysis, we study the steady mixed convection boundary layer flow of an incompressible Maxwell fluid near the two-dimensional stagnation-point flow over a vertical stretching surface. It is assumed that the stretching velocity and the surface temperature vary linearly with the distance from the stagnation-point. The governing nonlinear partial differential equations have been reduced to the coupled nonlinear ordinary differential equations by the similarity transformations. Analytical and numerical solutions of the derived system of equations are developed. The homotopy analysis method (HAM) and finite difference scheme are employed in constructing the analytical and numerical solutions, respectively. Comparison between the analytical and numerical solutions is given and found to be in excellent agreement. Both cases of assisting and opposing flows are considered. The influence of the various interesting parameters on the flow and heat transfer is analyzed and discussed through graphs in detail. The values of the local Nusselt number for different physical parameters are also tabulated. Comparison of the present results with known numerical results of viscous fluid is shown and a good agreement is observed.     

On the oscillation of some linear difference equations with periodic coefficients

Some linear difference equations with periodic coefficients (not necessarily nonnegative) are considered. Necessary conditions and sufficient conditions for the oscillation of the solutions are established. Conditions under which all nonoscillatory solutions tend to zero at ∞ are also presented. The results obtained are the discrete analogues of the oscillation results for some linear delay differential equations with periodic coefficients, which were given earlier by the second author [Oscillations of some delay differential equations with periodic coefficients, J. Math. Anal. Appl. 162 (1991) 452–475].