The factorization method in a rough-surface scattering problem: I

In the present paper the wave scattering problem on rough surface is considered for the Helmholtz equation with the Dirichlet boundary condition. An approximate solution is derived with using a factorization approach to the original Helmholtz equation. As a result, the system of two equations of parabolic type appears. The first system equation has an exact analytical solution whereas for the second one, an approximate solution, is considered in terms of perturbation series. It is shown that the obtained approximate solution is the modified classical small perturbation series with respect to small Rayleigh parameter. In Appendix A it is demonstrated that, when the derived perturbation series is converged, it is possible to summarize it and to represent the exact solution of original boundary problem in an analytical symbolical form.     

A new approach for modelling the damped Helmholtz oscillator:applications to plasma physics and electronic circuits

In this paper,a new approach is devoted to find novel analytical and approximate solutions to the damped quadratic nonlinear Helmholtz equation(HE)in terms of the Weiersrtrass elliptic function.The exact solution for undamped HE(integrable case)and approximate/semi-analytical solution to the damped HE(non-integrable case)are given for any arbitrary initial conditions.As a special case,the necessary and sufficient condition for the integrability of the damped HE using an elementary approach is reported.In general,a new ansatz is suggested to find a semi-analytical solution to the non-integrable case in the form of Weierstrass elliptic function.In addition,the relation between the Weierstrass and Jacobian elliptic functions solutions to the integrable case will be derived in details.Also,we will make a comparison between the semi-analytical solution and the approximate numerical solutions via using Runge-Kutta fourth-order method,finite difference method,and homotopy perturbation method for the first-two approximations.Furthermore,the maximum distance errors between the approximate/semi-analytical solution and the approximate numerical solutions will be estimated.As real applications,the obtained solutions will be devoted to describe the characteristics behavior of the oscillations in RLC series circuits and in various plasma models such as electronegative complex plasma model.     

Solution to the MHD flow over a non-linear stretching sheet by homotopy perturbation method

In this study,by means of homotopy perturbation method(HPM) an approximate solution of the magnetohydrodynamic(MHD) boundary layer flow is obtained.The main feature of the HPM is that it deforms a difficult problem into a set of problems which are easier to solve.HPM produces analytical expressions for the solution to nonlinear differential equations.The obtained analytic solution is in the form of an infinite power series.In this work,the analytical solution obtained by using only two terms from HPM soluti...     

Non-linear vibrations of a simple-simple beam with a non-ideal support in between

A simply supported Euler-Bernoulli beam with an intermediate support is considered. Non-linear terms due to immovable end conditions leading to stretching of the beam are included in the equations of motion. The concept of non-ideal boundary conditions is applied to the beam problem. In accordance, the intermediate support is assumed to allow small deflections. An approximate analytical solution of the problem is found using the method of multiple scales, a perturbation technique. Ideal and non-ideal frequencies as well as frequency-response curves are contrasted.     

Vortical and acoustical mode coupling inside a porous tube with uniform wall suction

This paper considers the oscillatory motion of gases inside a long porous tube of the closed-open type. In particular, the focus is placed on describing an analytical solution for the internal acoustico-vortical coupling that arises in the presence of appreciable wall suction. This unsteady field is driven by longitudinal oscillatory waves that are triggered by small unavoidable fluctuations in the wall suction speed. Under the assumption of small amplitude oscillations, the time-dependent governing equations are linearized through a regular perturbation of the dependent variables. Further application of the Helmholtz vector decomposition theorem enables us to discriminate between acoustical and vortical equations. After solving the wave equation for the acoustical contribution, the boundary-driven vortical field is considered. The method of matched-asymptotic expansions is then used to obtain a closed-form solution for the unsteady momentum equation developing from flow decomposition. An exact series expansion is also derived and shown to coincide with the numerical solution for the problem. The numerically verified end results suggest that the asymptotic scheme is capable of providing a sufficiently accurate solution. This is due to the error associated with the matched-asymptotic expansion being smaller than the error introduced in the Navier-Stokes linearization. A basis for comparison is established by examining the evolution of the oscillatory field in both space and time. The corresponding boundary-layer behavior is also characterized over a range of oscillation frequencies and wall suction velocities. In general, the current solution is found to exhibit features that are consistent with the laminar theory of periodic flows. By comparison to the Sexl profile in nonporous tubes, the critically damped solution obtained here exhibits a slightly smaller overshoot and depth of penetration. These features may be attributed to the suction effect that tends to attract the shear layers closer the wall.     

A fast method for the solution of the Helmholtz equation

In this paper, we consider the numerical solution of the Helmholtz equation, arising from the study of the wave equation in the frequency domain. The approach proposed here differs from those recently considered in the literature, in that it is based on a decomposition that is exact when considered analytically, so the only degradation in computational performance is due to discretization and roundoff errors. In particular, we make use of a multiplicative decomposition of the solution of the Helmholtz equation into an analytical plane wave and a multiplier, which is the solution of a complex-valued advection–diffusion–reaction equation. The use of fast multigrid methods for the solution of this equation is investigated. Numerical results show that this is an efficient solution algorithm for a reasonable range of frequencies.     

Approximation of the parabolic equation and the wavefield of a point source in a layered random medium

This paper studies the wavefield of a source in a multidimensional randomly layered medium. They obtained asymptotical expressions of the wave statistical characteristics for different boundary conditions both in the framework of the parabolic equation approximation and the exact formulation of the boundary problem for the Helmholtz equation. It is shown that the presence of a small but finite absorption γ is most important for the statistics. The diffraction effects turn out to be like those of absorption, but γ cannot tend to zero in this problem. In an appendix they give the factorization formulae of the wave equation solution in a layered medium.     

Statistical characteristics of photon distributions in a semi-infinite turbid medium: Analytical expressions and their application to optical tomography

The paper focuses on the determination of statistical characteristics of photon distributions in a semi-infinite turbid medium, specifically the photon average trajectory and the root-mean-square deviation of photons from the average trajectory, with an approach based on the diffusion approximation to the radiative transfer equation. We show that the Dirichlet and Robin boundary conditions used for this purpose give close results. We derive exact analytical expressions for the case of the Dirichlet boundary condition. To demonstrate the practical value of our results we consider approximate solution of the inverse problem of time-domain diffuse optical tomography with the flat layer transmission geometry. The problem is solved with the method of photon average trajectories which are constructed with analytical expressions derived for a semi-infinite medium.     

A new approximate analytical approach for dispersion relation of the nonlinear Klein-Gordon equation

A novel approach is presented for obtaining approximate analytical expressions for the dispersion relation of periodic wavetrains in the nonlinear Klein-Gordon equation with even potential function. By coupling linearization of the governing equation with the method of harmonic balance, we establish two general analytical approximate formulas for the dispersion relation, which depends on the amplitude of the periodic wavetrain. These formulas are valid for small as well as large amplitude of the wavetrain. They are also applicable to the large amplitude regime, which the conventional perturbation method fails to provide any solution, of the nonlinear system under study. Three examples are demonstrated to illustrate the excellent approximate solutions of the proposed formulas with respect to the exact solutions of the dispersion relation. (c) 2001 American Institute of Physics.     

Approximation of the parabolic equation and the wavefield of a point source in a layered random medium

Abstract

This paper studies the wavefield of a source in a multidimensional randomly layered medium. They obtained asymptotical expressions of the wave statistical characteristics for different boundary conditions both in the framework of the parabolic equation approximation and the exact formulation of the boundary problem for the Helmholtz equation. It is shown that the presence of a small but finite absorption γ is most important for the statistics. The diffraction effects turn out to be like those of absorption, but γ cannot tend to zero in this problem. In an appendix they give the factorization formulae of the wave equation solution in a layered medium.     

Analysis of modified Van der Pol’s oscillator using He’s parameter-expanding methods

In this work, a powerful analytical method, called He’s parameter-expanding methods (HPEM) is used to obtain the exact solutions of non-linear modified Van der Pol’s oscillator. The classical Van der Pol equation with delayed feedback and a modified equation where a delayed term provides the damping are considered. It is shown that one term in series expansions is sufficient to obtain a highly accurate solution, which is valid for the whole solution domain. Comparison of the obtained solution with those obtained using perturbation method shows that this method is effective and convenient to solve this problem. This method introduces a capable tool to solve this kind of non-linear problems.     

Self-similar factor approximants for evolution equations and boundary-value problems

The method of self-similar factor approximants is shown to be very convenient for solving different evolution equations and boundary-value problems typical of physical applications. The method is general and simple, being a straightforward two-step procedure. First, the solution to an equation is represented as an asymptotic series in powers of a variable. Second, the series are summed by means of the self-similar factor approximants. The obtained expressions provide highly accurate approximate solutions to the considered equations. In some cases, it is even possible to reconstruct exact solutions for the whole region of variables, starting from asymptotic series for small variables. This can become possible even when the solution is a transcendental function. The method is shown to be more simple and accurate than different variants of perturbation theory with respect to small parameters, being applicable even when these parameters are large. The generality and accuracy of the method are illustrated by a number of evolution equations as well as boundary value problems.     

Studies of scattering theory using numerical methods

Abstract

Numerical simulations, using both exact and approximate methods, are used to study rough surface scattering in both the smd and large roughness regimes. This study is limited lo scattcring lrom rough one-dimensional surfaces that obey the Dirichlet boundary condition and have a Gaussian roughness spectrum. For surfdces with small roughness (kh?1, where k is the radiation wavenumber and h is the root-mean-square (RMS) Surface height), perturbation theory is known to be valid. However, it is shown numerically that when kh?1 and kl?6 (where I is the surface correlation length) the Kirchhoffapprorimation is valid except at low grazing angles, and one must sum the first three orders of perturbation theory obtain the correct result. For kh?1 and kl?1, first-order perturbation theory is accurate. In this region, the accuracy of the first two terms of the iterative series solution of the exact integral equation is examined; the first term a1 this series is the Kirchhoff approximation, It is shown numerically that lor very small kh these first two terms reduce to first-order perturbation theory. However, lor this reduction to occur, kh must be made smaller than necessdry lor first-order perturbation theory to be accurate. In the regime of large roughness (kh?1) backscattering enhancement occurs when the RMS slope is on the order of unity. Several investigators have recently shown that the second term of the iterative series solution (the double-scattering term) replicates the properties of backscattering enhancement reasonably well. However, the double-scattering term has a lundamental flaw: predictions lor the scattering cross section per unit length based on the double-scattering term increase as the surfdce length is increased. This is shown here with numerical simulations and with an approximate analytical result based on the high frequency limit. The physical significance of this finding is also discussed. The final topic is the use of the double-scattering approximation to study the mechanism for backscattering enhancement with the Dirichlet boundary condition. This mechanism is usually assumed to be interference between reciprocal scattering paths. When the interlerence between reciprocal scattering paths is removed, the enhancement is eliminated. This shows that interference between reciprocal paths is almost certainly the dominant mechanism for backscattering enhancement in the scattering regime studied.     

New solutions for conformable fractional Nizhnik–Novikov–Veselov system via $$G'/G$$ expansion method and homotopy analysis methods

The main purpose of this paper is to find the exact and approximate analytical solution of Nizhnik–Novikov–Veselov system which may be considered as a model for an incompressible fluid with newly defined conformable derivative by using \(G'/G\) expansion method and homotopy analysis method (HAM) respectively. Authors used conformable derivative because of its applicability and lucidity. It is known that, the NNV system of equations is an isotropic Lax integrable extension of the well-known KdV equation and has physical significance. Also, NNV system of equations can be derived from the inner parameter-dependent symmetry constraint of the KP equation. Then the exact solutions obtained by using \(G'/G\) expansion method are compared with the approximate analytical solutions attained by employing HAM.     

Scattering from a rough layer of a random medium

This paper deals with scattering from a random-medium layer with rough boundaries. The fluctuations of the surface heights and medium permittivity are assumed to be small and smooth. All random quantities are assumed to be stationary and independent of each other. After the introduction of approximate boundary conditions, the system of partial differential equations is transformed into an integral equation where the fluctuations of the problem are represented as a zero-mean random operator. Employing smoothing, integral equations for the coherent fields are obtained. Use of the Helmholtz operator leads to solution for the coherent propagation constant while the boundary operators lead to coherent Fresnel coefficients. The characteristics of the results are illustrated by considering several examples.     

Perturbed Lie Symmetry and Systems of Non-Linear Diffusion Equations

Abstract

The method of one parameter, point symmetric, approximate Lie group invariants is applied to the problem of determining solutions of systems of pure one-dimensional, diffusion equations. The equations are taken to be non-linear in the dependent variables but otherwise homogeneous. Moreover, the matrix of diffusion coefficients is taken to differ from a constant matrix by a linear perturbation with respect to an infinitesimal parameter. The conditions for approximate Lie invariance are developed and are applied to the coupled system. The corresponding prolongation operator is derived and it is shown that this places a power law and logarithmic constraints on the nature of the perturbed diffusion matrix. The method is used to derive an approximate solution of the perturbed diffusion equation corresponding to impulsive boundary conditions.     

激光输出功率随耦合率及放大率的变化规律

建立了涵盖稳定腔和非稳腔的激光有源理论模型和激光开式谐振腔的边界条件。根据定态薛定谔方程与亥姆霍兹方程的等价性,采用量子力学方法求解谐振腔的亥姆霍兹波动方程,得到满足方程、边界条件、矩阵光学、稳定条件的本征解和本征方程。根据本征方程推导出横模数目随耦合率变化的规律,进而推导出激光输出功率随耦合率变化公式,以及激光输出功率随放大率的变化规律。该理论模型能够从稳定腔自然过渡到非稳腔,在小耦合率情况下退化到与传统公式基本一致的形式,又能在大耦合率情形下与实验结果符合得很好。     

Analytical approximations to the dynamics of an array of coupled DC SQUIDs

Coupled dynamical systems that operate near the onset of a bifurcation can lead, under certain conditions, to strong signal amplification effects. Over the past years we have studied this generic feature on a wide range of systems, including: magnetic and electric fields sensors, gyroscopic devices, and arrays of loops of superconducting quantum interference devices, also known as SQUIDs. In this work, we consider an array of SQUID loops connected in series as a case study to derive asymptotic analytical approximations to the exact solutions through perturbation analysis. Two approaches are considered. First, a straightforward expansion in which the non-linear parameter related to the inductance of the DC SQUID is treated as the small perturbation parameter. Second, a more accurate procedure that considers the SQUID phase dynamics as non-uniform motion on a circle. This second procedure is readily extended to the series array and it could serve as a mathematical framework to find approximate solutions to related complex systems with high-dimensionality. To the best of our knowledge, an approximate analytical solutions to an array of SQUIDs has not been reported yet in the literature.     

On the theory of electromagnetic radiation absorption during electron intersubband transitions in superlattices

An exact solution to the Schrödinger equation for electrons in superlattices with rectangular potential barriers and a periodic potential was obtained. An exact analytical expression for the electromagnetic radiation absorbance during electron intersubband transitions in such superlattices was derived in the first-order perturbation theory with arbitrary parameters. A number of extreme cases were considered.     

三维Helmholtz方程外问题的自然积分方程及其数值解

用文[2,3]提出的自然边界归化方法来处理三维Helmholtz方程的外边值问题。在简要介绍如何用球谐展开的方法得到Helmholtz问题在外球域上的自然积分方程后,给出求解该自然积分方程的一种数值方法及相应的数值算例。