Nonlinear dynamical modelling of chaotic electrostatic ion cyclotron oscillations by jerk equations

Plasma being a nonlinear and complex system, is capable of sustaining a wide spectrum of waves, oscillations and instabilities. These fluctuations interact nonlinearly amongst themselves and also with particles: electrons/ions and thus lead to nonlinear wave-wave or wave-particle interaction. In the presence of coherent waves the particles are accelerated whereas irregular oscillations can give rise to particle heating which is also called stochastic heating. Particle orbits are known to be randomized by the wave fields such that their motion can also become stochastic. For fusion to be sustained one needs a very high temperature plasma for an extended duration. It quite common to deploy external waves like electron cyclotron waves or ion cyclotron waves for plasma heating and current drive. These external waves also work only in certain regimes. Conventional plasma techniques have been able to answer several of the observations of the above processes related to heating transport etc, but nonlinear dynamics as a tool has helped in comprehending the plasma oscillations better. We have for the first time obtained a Third Order nonlinear ordinary differential equation (TONLODE) also known as jerk equation to describe the electrostatic ion cyclotron plasma oscillations in a magnetic field. The interesting feature of this equation is that it does not require an external forcing term to obtain chaotic behaviour.     

On the Krook-Wu model of the Boltzmann equation

The distribution function of the Krook-Wu model of the nonlinear Boltzmann equation (elastic differential cross sections inversely proportional to the relative speed of the colliding particles) is obtained as a generalized Laguerre polynomial expansion where the only time dependence is provided by the coefficients. In a recent paper M. Barnsley and the present author have shown that these coefficients are recursively determined from the resolution of a nonlinear differential system. Here we explicitly show how to construct the solutions of the Krook-Wu model and study the properties of the corresponding Krook-Wu distribution functions.     

Numerical modeling of shock-wave instability in thermodynamically nonideal media

A numerical analysis of the nonlinear instability of shock waves is presented for solid deuterium and for a model medium described by a properly constructed equation of state. The splitting of an unstable shock wave into an absolutely stable shock and a shock that emits acoustic waves is simulated for the first time.     

Completely Integrable Equation for the Quantum Correlation Function of Nonlinear Schrödinger Equation

Correlation functions of exactly solvable models can be described by differential equations [1]. In this paper we show that for the non-free fermionic case, differential equations should be replaced by integro-differential equations. We derive an integro-differential equation, which describes a time and temperature dependent correlation function of the penetrable Bose gas. The integro-differential equation turns out to be the continuum generalization of the classical nonlinear Schr?dinger equation. Received: 15 January 1997 / Accepted: 21 March 1997     

An extended Jacobi elliptic function rational expansion method and its application to ()-dimensional dispersive long wave equation

With the aid of computerized symbolic computation, a new elliptic function rational expansion method is presented by means of a new general ansatz, in which periodic solutions of nonlinear partial differential equations that can be expressed as a finite Laurent series of some of 12 Jacobi elliptic functions, is more powerful than exiting Jacobi elliptic function methods and is very powerful to uniformly construct more new exact periodic solutions in terms of rational formal Jacobi elliptic function solution of nonlinear partial differential equations. As an application of the method, we choose a (2+1)-dimensional dispersive long wave equation to illustrate the method. As a result, we can successfully obtain the solutions found by most existing Jacobi elliptic function methods and find other new and more general solutions at the same time. Of course, more shock wave solutions or solitary wave solutions can be gotten at their limit condition.     

On the two-dimensional instability of square-wave detonations

In order to elucidate the basic mechanisms responsible for the cells developing on detonation fronts, the two-dimensional instability of planar Chapman-Jouguet detonations is studied by the use of a square-wave model. It has been shown that there are two main instability mechanisms: the high sensitivity of the heat release rate to temperature and the hydrodynamic effect arising from the interface (reaction front) subject to the influence of the deflection of streamlines across the perturbed shock. Since for ordinary gaseous detonations the density across the shock is very large, the instability related to this pure hydrodynamic effect is very strong. An exact dispersion relation is derived and it is shown that the instability of square-wave detonations is described by a differential-difference equation of advanced type, which has a set of an infinite number of unstable solution branches with the growth rate increasing with the number of modes. The nonlinear solution shows that the discontinuity develops from smooth initial data within a finite time. Therefore the results obtained from square-wave detonations cannot be applied directly to describe the observed instability patterns. By carrying out a Taylor development of the terms with advanced time, a dispersion relation of a polynomial form is obtained and the high-frequency instability is eliminated. On the basis of the linear results with a regular reaction model, a phenomenological nonlinear equation of the fourth order for the position of the detonation surface is obtained by carrying out a Taylor series development. One-dimensional numerical solution of the phenomenological nonlinear equation shows typical nonlinear phenomena such as periodic oscillation, period doubling and dynamic quenching.     

Acoustic shock wave propagation in a heterogeneous medium: a numerical simulation beyond the parabolic approximation

Numerical simulation of nonlinear acoustics and shock waves in a weakly heterogeneous and lossless medium is considered. The wave equation is formulated so as to separate homogeneous diffraction, heterogeneous effects, and nonlinearities. A numerical method called heterogeneous one-way approximation for resolution of diffraction (HOWARD) is developed, that solves the homogeneous part of the equation in the spectral domain (both in time and space) through a one-way approximation neglecting backscattering. A second-order parabolic approximation is performed but only on the small, heterogeneous part. So the resulting equation is more precise than the usual standard or wide-angle parabolic approximation. It has the same dispersion equation as the exact wave equation for all forward propagating waves, including evanescent waves. Finally, nonlinear terms are treated through an analytical, shock-fitting method. Several validation tests are performed through comparisons with analytical solutions in the linear case and outputs of the standard or wide-angle parabolic approximation in the nonlinear case. Numerical convergence tests and physical analysis are finally performed in the fully heterogeneous and nonlinear case of shock wave focusing through an acoustical lens.     

A Model for Coagulation with Mating

We consider in this work a model for aggregation, where the coalescing particles initially have a certain number of potential links (called arms) which are used to perform coagulations. There are two types of arms, male and female, and two particles may coagulate only if one has an available male arm, and the other has an available female arm. After a coagulation, the used arms are no longer available. We are interested in the concentrations of the different types of particles, which are governed by a modification of Smoluchowski’s coagulation equation—that is, an infinite system of nonlinear differential equations. Using generating functions and solving a nonlinear PDE, we show that, up to some critical time, there is a unique solution to this equation. The Lagrange Inversion Formula allows in some cases to obtain explicit solutions, and to relate our model to two recent models for limited aggregation. We also show that, whenever the critical time is infinite, the concentrations converge to a state where all arms have disappeared, and the distribution of the masses is related to the law of the size of some two-type Galton-Watson tree. Finally, we consider a microscopic model for coagulation: we construct a sequence of Marcus-Lushnikov processes, and show that it converges, before the critical time, to the solution of our modified Smoluchowski’s equation.     

Formation of dust-acoustic shock waves in a strongly coupled cryogenic dusty plasma

The nonlinear propagation of ultra-low-frequency dust-acoustic (DA) waves in a strongly coupled cryogenic dusty plasma has been investigated, by using the Boltzmann distributed electrons and ions, as well as modified hydrodynamic equations for strongly coupled charged dust grains. The reductive perturbation technique is used to derive the Burger equation. It is shown that strong correlations among negatively charged dust particles acts like a dissipation, which is responsible for the formation of the DA shock waves. The latter are associated with the negative potential, i.e. with the compression of negatively charged cryogenic dust particle density. It is also found that the effective dust-temperature, which arises from electrostatic interactions among negatively charged dust particles, significantly affects the height of the DA shock structures. New laboratory experiments at cryogenic temperature should be conducted to verify our theoretical prediction.     

Parametric study of cylindrical and spherical dust ion acoustic shock waves with two temperature electrons in dusty plasma relevant to Saturn's E ring

We have performed numerical analysis of the one-dimensional dynamics of the cylindrical/spherical dust ion acoustic shock waves in unmagnetized dusty plasma consisting of positive ions, immobile dust particles, and nonextensive distributed cold and hot electrons. A multiple-scale expansion method is used to derive Burgers Equation (BE) and modified Burgers equation (MBE) by including higher order nonlinearity. The basic characteristics of the shock waves have been analysed numerically and graphically for different physical parameters relevant to Saturn' E ring through 2D figures. The parametric dependence of dust ion acoustic shock waves on some plasma parameters nonextensive index, density, and temperature of cold and hot electrons, concentration of dust particles, thermal effects and kinematic viscosity of ions is explored. Furthermore, it is found that the nonplanar geometry effects have an important impact on the establishment of shock waves. The amplitude of the wave decreases faster as one departs away from the axis of the cylinder or centre of the sphere. Such decaying behaviour continues as time progresses. It is also found that an increasing dust concentration decreases the amplitude of the dust ion acoustic shock waves.     

Novel solitary and shock wave solutions for the generalized variable-coefficients (2+1)-dimensional KP-Burger equation arising in dusty plasma

The 2-D generalized variable-coefficient Kadomtsev-Petviashvili-Burgers equation representing many types of acoustic waves in cosmic and/or laboratory dusty plasmas is reduced by the modified classical direct similarity reduction method to nonlinear ordinary differential equation of fourth-order. Using the extended Riccati equation mapping method for solving the reduced equation, many new shock wave, solitary wave and periodic wave solutions are obtained with some constraints between the variable coefficients. Finally, some physical interpretations for the obtained solutions as, bright and dark solitons, periodic solitary wave, and shock wave in dust plasma and quantum plasma are achieved.     

Solitary waves on a curved space-time

A nonlinear partial differential equation is derived which admits plane solitary waves on a conformally flat Riemannian space-time. The metric is determined by the amplitude of these waves. By interpreting these solitary waves as particles we arrive at the following picture: these particles are confined to regions exhibiting singular (very large) amplitudes in an otherwise continuous wavetrain. There is, thus, no distinction between the notion of a particle and that of a wave.     

Nonlinear electrostatic waves in inhomogeneous dense dusty magnetoplasmas

The nonlinear electrostatic drift waves are studied using quantum hydrodynamic model in dusty quantum magnetoplasmas. The dissipative effects due to collisions between ions and dust particles have also been taken into account. The Korteweg-de Vries Burgers (KdVB) like equation is derived and analytical solution is obtained using tanh method. The limiting cases of KdV type solitary waves, Burger type monotonic shock waves and oscillatory shock solutions are also presented. It is found that both hump and dip type solitary structures are possible in quantum dusty plasmas. However, amplitude and width of the nonlinear structure depend on the dust charge polarity and its concentration in electron-ion quantum plasmas. The monotonic shock like structure is independent of the quantum parameter. It is found that shock strength is increased in the presence of positively charged particles in comparison with negatively charged dust particles. The oscillatory shock structures are also obtained and it is found that change in dust charge polarity only shifts the phase of the oscillatory shock in plasmas. The numerical results are also presented for illustration.     

Chaotic diffusion for particles moving in a time dependent potential well

The chaotic diffusion for particles moving in a time dependent potential well is described by using two different procedures: (i) via direct evolution of the mapping describing the dynamics and; (ii) by the solution of the diffusion equation. The dynamic of the diffusing particles is made by the use of a two dimensional, nonlinear area preserving map for the variables energy and time. The phase space of the system is mixed containing both chaos, periodic regions and invariant spanning curves limiting the diffusion of the chaotic particles. The chaotic evolution for an ensemble of particles is treated as random particles motion and hence described by the diffusion equation. The boundary conditions impose that the particles can not cross the invariant spanning curves, serving as upper boundary for the diffusion, nor the lowest energy domain that is the energy the particles escape from the time moving potential well. The diffusion coefficient is determined via the equation of the mapping while the analytical solution of the diffusion equation gives the probability to find a given particle with a certain energy at a specific time. The momenta of the probability describe qualitatively the behavior of the average energy obtained by numerical simulation, which is investigated either as a function of the time as well as some of the control parameters of the problem.     

Modeling of an electrohydraulic lithotripter with the KZK equation.

The acoustic pressure field of an electrohydraulic extracorporeal shock wave lithotripter is modeled with a nonlinear parabolic wave equation (the KZK equation). The model accounts for diffraction, nonlinearity, and thermoviscous absorption. A numerical algorithm for solving the KZK equation in the time domain is used to model sound propagation from the mouth of the ellipsoidal reflector of the lithotripter. Propagation within the reflector is modeled with geometrical acoustics. It is shown that nonlinear distortion within the ellipsoidal reflector can play an important role for certain parameters. Calculated waveforms are compared with waveforms measured in a clinical lithotripter and good agreement is found. It is shown that the spatial location of the maximum negative pressure occurs pre-focally which suggests that the strongest cavitation activity will also be in front of the focus. Propagation of shock waves from a lithotripter with a pressure release reflector is considered and because of nonlinear propagation the focal waveform is not the inverse of the rigid reflector. Results from propagation through tissue are presented; waveforms are similar to those predicted in water except that the higher absorption in the tissue decreases the peak amplitude and lengthens the rise time of the shock.     

Numerical Investigation Into the Highly Nonlinear Heat Transfer Equation with Bremsstrahlung Emission in the Inertial Confinement Fusion Plasmas

A highly nonlinear parabolic partial differential equation that models the electron heat transfer process in laser inertial fusion has been solved numerically. The strong temperature dependence of the electron thermal conductivity and heat loss term (Bremsstrahlung emission) makes this a highly nonlinear process. In this case, an efficient numerical method is developed for the energy transport mechanism from the region of energy deposition into the ablation surface by a combination of the Crank‐Nicolson scheme and the Newton‐Raphson method. The quantitative behavior of the electron temperature and the comparison between analytic and numerical solutions are also investigated. For more clarification, the accuracy and conservation of energy in the computations are tested. The numerical results can be used to evaluate the nonlinear electron heat conduction, considering the released energy of the laser pulse at the Deuterium‐Tritium (DT) targets and preheating by heat conduction ahead of a compression shock in the inertial confinement fusion (ICF) approach. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Heat Transfe for Power Law Non-Newtonian Fluids

We present a theoretical analysis for heat transfer in power law non-Newtonian fluid by assuming that the thermal diffusivity is a function of temperature gradient. The laminar boundary layer energy equation is considered as an example to illustrate the application. It is shown that the boundary layer energy equation subject to the corresponding boundary conditions can be transformed to a boundary value problem of a nonlinear ordinary differential equation when similarity variables are introduced. Numerical solutions of the similarity energy equation are presented.     

Ion-acoustic shock waves in multi-electron temperature collisional plasma

The nonlinear propagation of ion-acoustic waves in a collision-dominated double electron temperature plasma is considered. Accounting for the ion viscosity and the ion heat conductivity, it is shown by means of two-warm fluid equations that the nonlinear evolution of the ion-acoustic waves is governed by the Korteweg—de Vries—Burgers equation. Stationary shock solution of the KdV—Burgers equation is presented.     

非线性LC电路方程的显式精确行波解

理论上考察了具有耗散的非线性LC电路中的行波. 借助于作者最近发展的精确求解非线性偏微分方程的扩展的双曲函数方法解析地研究了模拟非线性电路中冲击波的四阶耗散非线性波动方程. 一致地获得了丰富的显式精确解析行波解, 包括精确冲击波解和奇异的行波解, 和三角函数有理形式的周期波解. 关键词： LC电路')" href="#">非线性LC电路 非线性耗散波动方程 冲击波 周期波     

Hurst exponents for interacting random walkers obeying nonlinear Fokker-Planck equations

Anomalous diffusion of random walks has been extensively studied for the case of non-interacting particles. Here we study the evolution of nonlinear partial differential equations by interpreting them as Fokker-Planck equations arising from interactions among random walkers. We extend the formalism of generalized Hurst exponents to the study of nonlinear evolution equations and apply it to several illustrative examples. They include an analytically solvable case of a nonlinear diffusion constant and three nonlinear equations which are not analytically solvable: the usual Fisher equation which contains a quadratic nonlinearity, a generalization of the Fisher equation with density-dependent diffusion constant, and the Nagumo equation which incorporates a cubic rather than a quadratic nonlinearity. We estimate the generalized Hurst exponents.