Generalized Burgers equations and Euler-Painlev transcendents. I

Initial-value problems for the generalized Burgers equation (GBE) ut+u betaux+lambdaualpha =(delta/2)uxx are discussed for the single hump type of initial data both continuous and discontinuous. The numerical solution is carried to the self-similar ``intermediate asymptotic'' regime when the solution is given analytically by the self-similar form. The nonlinear (transformed) ordinary differential equations (ODE's) describing the self-similar form are generalizations of a class discussed by Euler and Painlevé and quoted by Kamke. These ODE's are new, and it is postulated that they characterize GBE's in the same manner as the Painlev equations categorize the Kortweg-de Vries (KdV) type. A connection problem for some related ODE's satisfying proper asymptotic conditions at x=±[infinity], is solved. The range of amplitude parameter is found for which the solution of the connection problem exists. The other solutions of the above GBE, which display several interesting features such as peaking, breaking, and a long shelf on the left for negative values of the damping coefficient lambda, are also discussed. The results are compared with those holding for the modified KdV equation with damping. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Oscillation criteria for differential equations of second order

In this article, we give sufficient condition in the form of integral inequalities to establish the oscillatory nature of non linear homogeneous differential equations of the form where r, q, p, f and g are given data. We do this by separating the two cases f is monotonous and non monotonous.

Observability for a Class of Nonlinear Systems

t - N m and sufficient computable conditions are obtained for the obsemabii of systems with linear state equations and polgwmIal outputs. Based on these, initial state reconstmctors are also described.

Coupled amplitude theory of nonlinear surface acoustic waves

The nonlinear propagation characteristics of surface acoustic waves on an isotropic elastic solid have been studied in this paper. The solution of the harmonic boundary value problem for Rayleigh waves is obtained as a generalized Fourier series whose coefficients are proportional to the slowly varying amplitudes of the various harmonics. The infinite set of coupled equations for the amplitudes when solved exhibit an oscillatory slow variation signifying a continuous transfer of energy back and forth among the various harmonics. A conservation relation is derived among all the harmonic amplitudes.

Invariance group properties and exact solutions of equations describing time-dependent free surface flows under gravity

Using the method of infinitesimal transformations, a 6-parameter family of exact solutions describing nonlinear sheared flows with a free surface are found. These solutions are a hybrid between the earlier self-propagating simple wave solutions of Freeman, and decaying solutions of Sachdev. Simple wave solutions are also derived via the method of infinitesimal transformations. Incomplete beta functions seem to characterize these (nonlinear) sheared flows in the absence of critical levels.

Application of Artificial Viscosity in Establishing Supercritical Solutions to the Transonic Integra

The nonlinear singular integral equation of transonic flow is examined in the free-stream Mach number range where only solutions with shocks are known to exist. It is shown that, by the addition of an artificial viscosity term to the integral equation, even the direct iterative scheme, with the linear solution as the initial iterate, leads to convergence. Detailed tables indicating how the solution varies with changes in the parameters of the artificial viscosity term are also given. In the best cases (when the artificial viscosity is smallest), the solutions compare well with known results, their characteristic feature being the representation of the shock by steep gradients rather than by abrupt discontinuities. However, 'sharp-shock solutions' have also been obtained by the implementation of a quadratic iterative scheme with the 'artificial viscosity solution' as the initial iterate; the converged solution with a sharp shock is obtained with only a few more iterates. Finally, a review is given of various shock-capturing and shock-fitting schemes for the transonic flow equations in general, and for the transonic integral equation in particular, frequent comparisons being made with the approach of this paper.

Nonlinear mode coupling of surface acoustic waves on an isotropic solid

The nonlinear mode coupling between two co-directional quasi-harmonic Rayleigh surface waves on an isotropic solid is analysed using the method of multiple scales. This procedure yields a system of six semi-linear hyperbolic partial differential equations with the same principal part governing the slow variations in the (complex) amplitudes of the two fundamental, the two second harmonic and the two combination frequency waves at the second stage of the perturbation expansion. A numerical solution of these equations for excitation by monochromatic signals at two arbitrary frequencies, indicates that there is a continuous transfer of energy back and forth among the fundamental, second harmonic and combination frequency waves due to mode coupling. The mode coupling tends to be more pronounced as the frequencies of the interacting waves approach each other.

Nonlinear surface acoustic waves on an isotropic solid

A systematic derivation of the approximate coupled amplitude equations governing the propagation of a quasi-monochromatic Rayleigh surface wave on an isotropic solid is presented, starting from the non-linear governing differential equations and the non-linear free-surface boundary conditions, using the method of mulitple scales. An explicit solution of these equations for a signalling problem is obtained in terms of hyperbolic functions. In the case of monochromatic excitation, it is shown that the second harmonic amplitude grows initially at the expense of the fundamental and that the amplitudes of the fundamental and second harmonic remain bounded for all time.

Approximation of the perturbation equations of a quasi-linear hyperbolic system in the neighborhood of a bicharacteristic

In 1956 Whitham gave a nonlinear theory for computing the intensity of an acoustic pulse of an arbitrary shape. The theory has been used very successfully in computing the intensity of the sonic bang produced by a supersonic plane. [4.] derived an approximate quasi-linear equation for the propagation of a short wave in a compressible medium. These two methods are essentially nonlinear approximations of the perturbation equations of the system of gas-dynamic equations in the neighborhood of a bicharacteristic curve (or rays) for weak unsteady disturbances superimposed on a given steady solution. In this paper we have derived an approximate quasi-linear equation which is an approximation of perturbation equations in the neighborhood of a bicharacteristic curve for a weak pulse governed by a general system of first order quasi-linear partial differential equations in m + 1 independent variables (t, x1,…, xm) and derived Gubkin's result as a particular case when the system of equations consists of the equations of an unsteady motion of a compressible gas. We have also discussed the form of the approximate equation describing the waves propagating upsteam in an arbitrary multidimensional transonic flow.

A Finite Variable Difference Relaxation Scheme for hyperbolic-parabolic equations

Using the framework of a new relaxation system, which converts a nonlinear viscous conservation law into a system of linear convection-diffusion equations with nonlinear source terms, a finite variable difference method is developed for nonlinear hyperbolic-parabolic equations. The basic idea is to formulate a finite volume method with an optimum spatial difference, using the Locally Exact Numerical Scheme (LENS), leading to a Finite Variable Difference Method as introduced by Sakai [Katsuhiro Sakai, A new finite variable difference method with application to locally exact numerical scheme, journal of Computational Physics, 124 (1996) pp. 301-308.], for the linear convection-diffusion equations obtained by using a relaxation system. Source terms are treated with the well-balanced scheme of Jin [Shi Jin, A steady-state capturing method for hyperbolic systems with geometrical source terms, Mathematical Modeling Numerical Analysis, 35 (4) (2001) pp. 631-645]. Bench-mark test problems for scalar and vector conservation laws in one and two dimensions are solved using this new algorithm and the results demonstrate the efficiency of the scheme in capturing the flow features accurately.

A semi-analytical particle filter for identification of nonlinear oscillators

Particle filters find important applications in the problems of state and parameter estimations of dynamical systems of engineering interest. Since a typical filtering algorithm involves Monte Carlo simulations of the process equations, sample variance of the estimator is inversely proportional to the number of particles. The sample variance may be reduced if one uses a Rao-Blackwell marginalization of states and performs analytical computations as much as possible. In this work, we propose a semi-analytical particle filter, requiring no Rao-Blackwell marginalization, for state and parameter estimations of nonlinear dynamical systems with additively Gaussian process/observation noises. Through local linearizations of the nonlinear drift fields in the process/observation equations via explicit Ito-Taylor expansions, the given nonlinear system is transformed into an ensemble of locally linearized systems. Using the most recent observation, conditionally Gaussian posterior density functions of the linearized systems are analytically obtained through the Kalman filter. This information is further exploited within the particle filter algorithm for obtaining samples from the optimal posterior density of the states. The potential of the method in state/parameter estimations is demonstrated through numerical illustrations for a few nonlinear oscillators. The proposed filter is found to yield estimates with reduced sample variance and improved accuracy vis-a-vis results from a form of sequential importance sampling filter.

Testing and Modeling of MR Damper and Its Application to SDOF Systems Using Integral Backstepping Technique

Magnetorheological dampers are intrinsically nonlinear devices, which make the modeling and design of a suitable control algorithm an interesting and challenging task. To evaluate the potential of magnetorheological (MR) dampers in control applications and to take full advantages of its unique features, a mathematical model to accurately reproduce its dynamic behavior has to be developed and then a proper control strategy has to be taken that is implementable and can fully utilize their capabilities as a semi-active control device. The present paper focuses on both the aspects. First, the paper reports the testing of a magnetorheological damper with an universal testing machine, for a set of frequency, amplitude, and current. A modified Bouc-Wen model considering the amplitude and input current dependence of the damper parameters has been proposed. It has been shown that the damper response can be satisfactorily predicted with this model. Second, a backstepping based nonlinear current monitoring of magnetorheological dampers for semi-active control of structures under earthquakes has been developed. It provides a stable nonlinear magnetorheological damper current monitoring directly based on system feedback such that current change in magnetorheological damper is gradual. Unlike other MR damper control techniques available in literature, the main advantage of the proposed technique lies in its current input prediction directly based on system feedback and smooth update of input current. Furthermore, while developing the proposed semi-active algorithm, the dynamics of the supplied and commanded current to the damper has been considered. The efficiency of the proposed technique has been shown taking a base isolated three story building under a set of seismic excitation. Comparison with widely used clipped-optimal strategy has also been shown.

Numerical procedure for second order non-linear ordinary differential equations and application to heat transfer problem

In this paper, we have first given a numerical procedure for the solution of second order non-linear ordinary differential equations of the type y″ = f (x;y, y′) with given initial conditions. The method is based on geometrical interpretation of the equation, which suggests a simple geometrical construction of the integral curve. We then translate this geometrical method to the numerical procedure adaptable to desk calculators and digital computers. We have studied the efficacy of this method with the help of an illustrative example with known exact solution. We have also compared it with Runge-Kutta method. We have then applied this method to a physical problem, namely, the study of the temperature distribution in a semi-infinite solid homogeneous medium for temperature-dependent conductivity coefficient.

A nonlinear theory of wake development

A simple new series, using an expansion of the velocity profile in parabolic cylinder functions, has been developed to describe the nonlinear evolution of a steady, laminar, incompressible wake from a given arbitrary initial profile. The first term in this series is itself found to provide a very satisfactory prediction of the decay of the maximum velocity defect in the wake behind a flat plate or aft of the recirculation zone behind a symmetric blunt body. A detailed analysis, including higher order terms, has been made of the flat plate wake with a Blasius profile at the trailing edge. The same method yields, as a special case, complete results for the development of linearized wakes with arbitrary initial profile under the influence of arbitrary pressure gradients. Finally, for purposes of comparison, a simple approximate solution is obtained using momentum integral methods, and found to predict satisfactorily the decay of the maximum velocity defect. © 1970 Wolters-Noordhoff Publishing.

Application of D transforms in the analysis of nonlinear distributed parameter systems

The scope of application of Laplace transforms presently limited to the study of linear partial differential equations, is extended to the nonlinear domain by this study. This has been achieved by modifying the definition of D transforms, put forth recently for the study of classes of nonlinear lumped parameter systems. The appropriate properties of the new D transforms are presented to bring out their applicability in the analysis of nonlinear distributed parameter systems.