Large Time Behavior of the Vlasov-Poisson-Boltzmann System

Directory of Open Access Journals (Sweden)

Li Li

2013-01-01

Full Text Available The motion of dilute charged particles can be modeled by Vlasov-Poisson-Boltzmann system. We study the large time stability of the VPB system. To be precise, we prove that when time goes to infinity, the solution of VPB system tends to global Maxwellian state in a rate Ot−∞, by using a method developed for Boltzmann equation without force in the work of Desvillettes and Villani (2005. The improvement of the present paper is the removal of condition on parameter λ as in the work of Li (2008.

From the Hartree dynamics to the Vlasov equation

DEFF Research Database (Denmark)

Benedikter, Niels Patriz; Porta, Marcello; Saffirio, Chiara

2016-01-01

We consider the evolution of quasi-free states describing N fermions in the mean field limit, as governed by the nonlinear Hartree equation. In the limit of large N, we study the convergence towards the classical Vlasov equation. For a class of regular interaction potentials, we establish precise...

Solving the Vlasov equation in two spatial dimensions with the Schrödinger method

Science.gov (United States)

Kopp, Michael; Vattis, Kyriakos; Skordis, Constantinos

2017-12-01

We demonstrate that the Vlasov equation describing collisionless self-gravitating matter may be solved with the so-called Schrödinger method (ScM). With the ScM, one solves the Schrödinger-Poisson system of equations for a complex wave function in d dimensions, rather than the Vlasov equation for a 2 d -dimensional phase space density. The ScM also allows calculating the d -dimensional cumulants directly through quasilocal manipulations of the wave function, avoiding the complexity of 2 d -dimensional phase space. We perform for the first time a quantitative comparison of the ScM and a conventional Vlasov solver in d =2 dimensions. Our numerical tests were carried out using two types of cold cosmological initial conditions: the classic collapse of a sine wave and those of a Gaussian random field as commonly used in cosmological cold dark matter N-body simulations. We compare the first three cumulants, that is, the density, velocity and velocity dispersion, to those obtained by solving the Vlasov equation using the publicly available code ColDICE. We find excellent qualitative and quantitative agreement between these codes, demonstrating the feasibility and advantages of the ScM as an alternative to N-body simulations. We discuss, the emergence of effective vorticity in the ScM through the winding number around the points where the wave function vanishes. As an application we evaluate the background pressure induced by the non-linearity of large scale structure formation, thereby estimating the magnitude of cosmological backreaction. We find that it is negligibly small and has time dependence and magnitude compatible with expectations from the effective field theory of large scale structure.

Multi-scale approximation of Vlasov equation

International Nuclear Information System (INIS)

Mouton, A.

2009-09-01

One of the most important difficulties of numerical simulation of magnetized plasmas is the existence of multiple time and space scales, which can be very different. In order to produce good simulations of these multi-scale phenomena, it is recommended to develop some models and numerical methods which are adapted to these problems. Nowadays, the two-scale convergence theory introduced by G. Nguetseng and G. Allaire is one of the tools which can be used to rigorously derive multi-scale limits and to obtain new limit models which can be discretized with a usual numerical method: this procedure is so-called a two-scale numerical method. The purpose of this thesis is to develop a two-scale semi-Lagrangian method and to apply it on a gyrokinetic Vlasov-like model in order to simulate a plasma submitted to a large external magnetic field. However, the physical phenomena we have to simulate are quite complex and there are many questions without answers about the behaviour of a two-scale numerical method, especially when such a method is applied on a nonlinear model. In a first part, we develop a two-scale finite volume method and we apply it on the weakly compressible 1D isentropic Euler equations. Even if this mathematical context is far from a Vlasov-like model, it is a relatively simple framework in order to study the behaviour of a two-scale numerical method in front of a nonlinear model. In a second part, we develop a two-scale semi-Lagrangian method for the two-scale model developed by E. Frenod, F. Salvarani et E. Sonnendrucker in order to simulate axisymmetric charged particle beams. Even if the studied physical phenomena are quite different from magnetic fusion experiments, the mathematical context of the one-dimensional paraxial Vlasov-Poisson model is very simple for establishing the basis of a two-scale semi-Lagrangian method. In a third part, we use the two-scale convergence theory in order to improve M. Bostan's weak-* convergence results about the finite

Response to Comment on 'On Higher-Order Corrections to Gyrokinetic Vlasov-Poisson Equations in the Long Wavelength Limit [Phys. Plasmas 16,044506 (2009)]'

International Nuclear Information System (INIS)

Lee, W.W.; Kolesnikov, R.A.

2009-01-01

We show in this Response that the nonlinear Poisson's equation in our original paper derived from the drift kinetic approach can be verified by using the nonlinear gyrokinetic Poisson's equation of Dubin et al. (Phys. Fluids 26, 3524 (1983)). This nonlinear contribution in φ 2 is indeed of the order of k # perpendicular# 4 in the long wavelength limit and remains finite for zero ion temperature, in contrast to the nonlinear term by Parra and Catto (Plasma Phys. Control. Fusion 50, 065014 (2008)), which is of the order of k # perpendicular# 2 and diverges for T i → 0. For comparison, the leading term for the gyrokinetic Poisson's equation in this limit is of the order of k # perpendicular# 2 φ.

The Arrow of Time in the Collapse of Collisionless Self-gravitating Systems: Non-validity of the Vlasov-Poisson Equation during Violent Relaxation

Science.gov (United States)

Beraldo e Silva, Leandro; de Siqueira Pedra, Walter; Sodré, Laerte; Perico, Eder L. D.; Lima, Marcos

2017-09-01

The collapse of a collisionless self-gravitating system, with the fast achievement of a quasi-stationary state, is driven by violent relaxation, with a typical particle interacting with the time-changing collective potential. It is traditionally assumed that this evolution is governed by the Vlasov-Poisson equation, in which case entropy must be conserved. We run N-body simulations of isolated self-gravitating systems, using three simulation codes, NBODY-6 (direct summation without softening), NBODY-2 (direct summation with softening), and GADGET-2 (tree code with softening), for different numbers of particles and initial conditions. At each snapshot, we estimate the Shannon entropy of the distribution function with three different techniques: Kernel, Nearest Neighbor, and EnBiD. For all simulation codes and estimators, the entropy evolution converges to the same limit as N increases. During violent relaxation, the entropy has a fast increase followed by damping oscillations, indicating that violent relaxation must be described by a kinetic equation other than the Vlasov-Poisson equation, even for N as large as that of astronomical structures. This indicates that violent relaxation cannot be described by a time-reversible equation, shedding some light on the so-called “fundamental paradox of stellar dynamics.” The long-term evolution is well-described by the orbit-averaged Fokker-Planck model, with Coulomb logarithm values in the expected range 10{--}12. By means of NBODY-2, we also study the dependence of the two-body relaxation timescale on the softening length. The approach presented in the current work can potentially provide a general method for testing any kinetic equation intended to describe the macroscopic evolution of N-body systems.

Integration of the three-dimensional Vlasov equation for a magnetized plasma

International Nuclear Information System (INIS)

Cheng, C.Z.

1976-04-01

A second order splitting scheme is developed to integrate the three dimensional Vlasov equation for a plasma in a magnetic field. The integration of the Vlasov equation is divided into a series of intermediate steps and Fourier interpolation and the ASD method with a third order Taylor expansion are used to integrate the fractional equations. Numerical experiments related to cyclotron waves in 2 and 2 1 / 2 D are demonstrated with high accuracy and efficiency. The computer storage requirements are modest; for example, a typical 2D nonlinear electron plasma simulation requires only 4000 ''particles.''

Nonlinear stationary solutions of the Wigner and Wigner-Poisson equations

OpenAIRE

Haas, F.; Shukla, P. K.

2008-01-01

Exact nonlinear stationary solutions of the one-dimensional Wigner and Wigner-Poisson equations in the terms of the Wigner functions that depend not only on the energy but also on position are presented. In this way, the Bernstein-Greene-Kruskal modes of the classical plasma are adapted for the quantum formalism in the phase space. The solutions are constructed for the case of a quartic oscillator potential, as well as for the self-consistent Wigner-Poisson case. Conditions for well-behaved p...

Kinetic description of intense nonneutral beam propagation through a periodic solenoidal focusing field based on the nonlinear Vlasov-Maxwell equations

International Nuclear Information System (INIS)

Davidson, R.C.; Chen, C.

1997-08-01

A kinetic description of intense nonneutral beam propagation through a periodic solenoidal focusing field B sol (rvec x) is developed. The analysis is carried out for a thin beam with characteristic beam radius r b much-lt S, and directed axial momentum γ b mβ b c (in the z-direction) large compared with the transverse momentum and axial momentum spread of the beam particles. Making use of the nonlinear Vlasov-Maxwell equations for general distribution function f b (rvec x,rvec p,t) and self-consistent electrostatic field consistent with the thin-beam approximation, the kinetic model is used to investigate detailed beam equilibrium properties for a variety of distribution functions. Examples are presented both for the case of a uniform solenoidal focusing field B z (z) = B 0 = const. and for the case of a periodic solenoidal focusing field B z (z + S) = B z (z). The nonlinear Vlasov-Maxwell equations are simplified in the thin-beam approximation, and an alternative Hamiltonian formulation is developed that is particularly well-suited to intense beam propagation in periodic focusing systems. Based on the present analysis, the Vlasov-Maxwell description of intense nonneutral beam propagation through a periodic solenoidal focusing field rvec B sol (rvec x) is found to be remarkably tractable and rich in physics content. The Vlasov-Maxwell formalism developed here can be extended in a straightforward manner to investigate detailed stability behavior for perturbations about specific choices of beam equilibria

Nonlinear stationary solutions of the Wigner and Wigner-Poisson equations

International Nuclear Information System (INIS)

Haas, F.; Shukla, P. K.

2008-01-01

Exact nonlinear stationary solutions of the one-dimensional Wigner and Wigner-Poisson equations in the terms of the Wigner functions that depend not only on the energy but also on position are presented. In this way, the Bernstein-Greene-Kruskal modes of the classical plasma are adapted for the quantum formalism in the phase space. The solutions are constructed for the case of a quartic oscillator potential, as well as for the self-consistent Wigner-Poisson case. Conditions for well-behaved physically meaningful equilibrium Wigner functions are discussed.

Stability analysis of cylindrical Vlasov equilibria

International Nuclear Information System (INIS)

Short, R.W.

1979-01-01

A general method of stability analysis is described which may be applied to a large class of such problems, namely those which are described dynamically by the Vlasov equation, and geometrically by cylindrical symmetry. The method is presented for the simple case of the Vlasov-Poisson (electrostatic) equations, and the results are applied to a calculation of the lower-hybrid-drift instability in a plasma with a rigid rotor distribution function. The method is extended to the full Vlasov-Maxwell (electromagnetic) equations. These results are applied to a calculation of the instability of the extraordinary electromagnetic mode in a relativistic E-layer interacting with a background plasma

Stability analysis of cylindrical Vlasov equilibria

International Nuclear Information System (INIS)

Short, R.W.

1979-01-01

A general method of stability analysis is described which may be applied to a large class of such problems, namely those which are described dynamically by the Vlasov equation, and geometrically by clindrical symmetry. The method is presented for the simple case of the Vlasov-Poisson (electrostatic) equations, and the results are applied to a calculation of the lower-hybrid-drift instability in a plasma with a rigid rotor distribution function. The method is extended to the full Vlasov-Maxwell (electromagnetic) equations. These results are applied to a calculation of the instability of the extraordinary electromagnetic mode in a relativistic E-layer interacting with a background plasma

Kinetic Boltzmann, Vlasov and Related Equations

CERN Document Server

Sinitsyn, Alexander; Vedenyapin, Victor

2011-01-01

Boltzmann and Vlasov equations played a great role in the past and still play an important role in modern natural sciences, technique and even philosophy of science. Classical Boltzmann equation derived in 1872 became a cornerstone for the molecular-kinetic theory, the second law of thermodynamics (increasing entropy) and derivation of the basic hydrodynamic equations. After modifications, the fields and numbers of its applications have increased to include diluted gas, radiation, neutral particles transportation, atmosphere optics and nuclear reactor modelling. Vlasov equation was obtained in

Comment on 'On higher order corrections to gyrokinetic Vlasov-Poisson equations in the long wavelength limit' [Phys. Plasmas 16, 044506 (2009)

International Nuclear Information System (INIS)

Parra, Felix I.; Catto, Peter J.

2009-01-01

A recent publication [F. I. Parra and P. J. Catto, Plasma Phys. Controlled Fusion 50, 065014 (2008)] warned against the use of the lower order gyrokinetic Poisson equation at long wavelengths because the long wavelength, radial electric field must remain undetermined to the order the equation is obtained. Another reference [W. W. Lee and R. A. Kolesnikov, Phys. Plasmas 16, 044506 (2009)] criticizes these results by arguing that the higher order terms neglected in the most common gyrokinetic Poisson equation are formally smaller than the terms that are retained. This argument is flawed and ignores that the lower order terms, although formally larger, must cancel without determining the long wavelength, radial electric field. The reason for this cancellation is discussed. In addition, the origin of a nonlinear term present in the gyrokinetic Poisson equation [F. I. Parra and P. J. Catto, Plasma Phys. Controlled Fusion 50, 065014 (2008)] is explained.

Expeditious 3D poisson vlasov algorithm applied to ion extraction from a plasma

International Nuclear Information System (INIS)

Whealton, J.H.; McGaffey, R.W.; Meszaros, P.S.

1983-01-01

A new 3D Poisson Vlasov algorithm is under development which differs from a previous algorithm, referenced in this paper, in two respects: the mesh lines are cartesian, and the Poisson equation is solved iteratively. The resulting algorithm has been used to examine the same boundary value problem as considered in the earlier algorithm except that the number of nodes is 2 times greater. The same physical results were obtained except the computational time was reduced by a factor of 60 and the memory requirement was reduced by a factor of 10. This algorithm at present restricts Neumann boundary conditions to orthogonal planes lying along mesh lines. No such restriction applies to Dirichlet boundaries. An emittance diagram is shown below where those points lying on the y = 0 line start on the axis of symmetry and those near the y = 1 line start near the slot end

Theoretical and numerical study of the equations of Vlasov-Maxwell in the covariant formalism

International Nuclear Information System (INIS)

Back, A.

2011-11-01

A new point of view is proposed for the simulation of plasmas using the kinetic model which links the equations of Vlasov for the distribution of particles and the equations of Maxwell for the electromagnetic contribution of fields. We use the following principle: the equations of Physics are mathematical objects which put in relation geometrical objects. To preserve the geometrical properties of the various objects in an equation, we use, for the theoretical and numerical study, the differential geometry. All the equations of Physics can be written with differential forms and this point of view is not dependent on the choice of coordinates. We propose then a discretization of the differential forms by using B-Splines. To be coherent with the theory, we also propose a discretization of the various operations of the differential geometry. We test our scheme, first on the equations of Maxwell with several boundary conditions and since it does not depend on the system of coordinates, we also test it when we change coordinates. Finally, we apply the same method to the equations of Vlasov-Poisson in one-dimension and we propose several numerical schemes. (author)

POISSON SUPERFISH, Poisson Equation Solver for Radio Frequency Cavity

International Nuclear Information System (INIS)

Colman, J.

2001-01-01

1 - Description of program or function: POISSON, SUPERFISH is a group of (1) codes that solve Poisson's equation and are used to compute field quality for both magnets and fixed electric potentials and (2) RF cavity codes that calculate resonant frequencies and field distributions of the fundamental and higher modes. The group includes: POISSON, PANDIRA, SUPERFISH, AUTOMESH, LATTICE, FORCE, MIRT, PAN-T, TEKPLOT, SF01, and SHY. POISSON solves Poisson's (or Laplace's) equation for the vector (scalar) potential with nonlinear isotropic iron (dielectric) and electric current (charge) distributions for two-dimensional Cartesian or three-dimensional cylindrical symmetry. It calculates the derivatives of the potential, the stored energy, and performs harmonic (multipole) analysis of the potential. PANDIRA is similar to POISSON except it allows anisotropic and permanent magnet materials and uses a different numerical method to obtain the potential. SUPERFISH solves for the accelerating (TM) and deflecting (TE) resonant frequencies and field distributions in an RF cavity with two-dimensional Cartesian or three-dimensional cylindrical symmetry. Only the azimuthally symmetric modes are found for cylindrically symmetric cavities. AUTOMESH prepares input for LATTICE from geometrical data describing the problem, (i.e., it constructs the 'logical' mesh and generates (x,y) coordinate data for straight lines, arcs of circles, and segments of hyperbolas). LATTICE generates an irregular triangular (physical) mesh from the input data, calculates the 'point current' terms at each mesh point in regions with distributed current density, and sets up the mesh point relaxation order needed to write the binary problem file for the equation-solving POISSON, PANDIRA, or SUPERFISH. FORCE calculates forces and torques on coils and iron regions from POISSON or PANDIRA solutions for the potential. MIRT optimizes magnet profiles, coil shapes, and current densities from POISSON output based on a

KEEN Wave Simulations: Comparing various PIC to various fixed grid Vlasov to Phase-Space Adaptive Sparse Tiling & Effective Lagrangian (PASTEL) Techniques

Science.gov (United States)

Afeyan, Bedros; Larson, David; Shadwick, Bradley; Sydora, Richard

2017-10-01

We compare various ways of solving the Vlasov-Poisson and Vlasov-Maxwell equations on rather demanding nonlinear kinetic phenomena associated with KEEN and KEEPN waves. KEEN stands for Kinetic, Electrostatic, Electron Nonlinear, and KEEPN, for electron-positron or pair plasmas analogs. Because these self-organized phase space structures are not steady-state, or single mode, or fluid or low order moment equation limited, typical techniques with low resolution or too much noise will distort the answer too much, too soon, and fail. This will be shown via Penrose criteria triggers for instability at the formation stage as well as particle orbit statistics in fully formed KEEN waves and KEEN-KEEN and KEEN-EPW interacting states. We will argue that PASTEL is a viable alternative to traditional methods with reasonable chances of success in higher dimensions. Work supported by a Grant from AFOSR PEEP.

Hamiltonian formalism of two-dimensional Vlasov kinetic equation.

Science.gov (United States)

Pavlov, Maxim V

2014-12-08

In this paper, the two-dimensional Benney system describing long wave propagation of a finite depth fluid motion and the multi-dimensional Russo-Smereka kinetic equation describing a bubbly flow are considered. The Hamiltonian approach established by J. Gibbons for the one-dimensional Vlasov kinetic equation is extended to a multi-dimensional case. A local Hamiltonian structure associated with the hydrodynamic lattice of moments derived by D. J. Benney is constructed. A relationship between this hydrodynamic lattice of moments and the two-dimensional Vlasov kinetic equation is found. In the two-dimensional case, a Hamiltonian hydrodynamic lattice for the Russo-Smereka kinetic model is constructed. Simple hydrodynamic reductions are presented.

Discrete nonlinear Schrodinger equations with arbitrarily high-order nonlinearities

DEFF Research Database (Denmark)

Khare, A.; Rasmussen, Kim Ø; Salerno, M.

2006-01-01

-Ladik equation. As a common property, these equations possess three kinds of exact analytical stationary solutions for which the Peierls-Nabarro barrier is zero. Several properties of these solutions, including stability, discrete breathers, and moving solutions, are investigated.......A class of discrete nonlinear Schrodinger equations with arbitrarily high-order nonlinearities is introduced. These equations are derived from the same Hamiltonian using different Poisson brackets and include as particular cases the saturable discrete nonlinear Schrodinger equation and the Ablowitz...

Second order oscillations of a Vlasov-Poisson plasma in the Fourier transformed space

International Nuclear Information System (INIS)

Sedlacek, Z.; Nocera, L.

1991-05-01

The Vlasov-Poisson system of equations in the Fourier-transformed velocity space is studied. At first some results of the linear theory are reformulated: in the new representation the Van Kampen eigenmodes and their adjoint are found to be ordinary functions with convenient piece-wise continuity properties. A transparent derivation is given of the free-streaming temporal echo in terms of the kinematics of wave packets in the Fourier-transformed velocity space. This analysis is further extended to include Coulomb interactions which allows to establish a connection between the echo theory, the second order oscillations of Best and the phenomenon of linear sidebands. The calculation of the time evolution of the global second order electric field is performed in detail in the case of a Maxwellian equilibrium distribution function. It is concluded that the phenomenon of linear sidebands may be properly explained in terms of the intrinsic features of the equilibrium distribution function. (author) 5 figs., 32 refs

Global Vlasov simulation on magnetospheres of astronomical objects

International Nuclear Information System (INIS)

Umeda, Takayuki; Ito, Yosuke; Fukazawa, Keiichiro

2013-01-01

Space plasma is a collisionless, multi-scale, and highly nonlinear medium. There are various types of self-consistent computer simulations that treat space plasma according to various approximations. We develop numerical schemes for solving the Vlasov (collisionless Boltzmann) equation, which is the first-principle kinetic equation for collisionless plasma. The weak-scaling benchmark test shows that our parallel Vlasov code achieves a high performance and a high scalability. Currently, we use more than 1000 cores for parallel computations and apply the present parallel Vlasov code to various cross-scale processes in space plasma, such as a global simulation on the interaction between solar/stellar wind and magnetospheres of astronomical objects

Dielectric energy versus plasma energy, and Hamiltonian action-angle variables for the Vlasov equation

International Nuclear Information System (INIS)

Morrison, P.J.

1992-04-01

Expressions for the energy content of one-dimensional electrostatic perturbations about homogeneous equilibria are revisited. The well-known dielectric energy, var-epsilon D , is compared with the exact plasma free energy expression, δ 2 F, that is conserved by the Vlasov-Poisson system. The former is an expression in terms of the perturbed electric field amplitude, while the latter is determined by a generating function, which describes perturbations of the distribution function that respect the important constraint of dynamical accessibility of the system. Thus the comparison requires solving the Vlasov equation for such a perturbations of the distribution function in terms of the electric field. This is done for neutral modes of oscillation that occur for equilibria with stationary inflection points, and it is seen that for these special modes δ 2 F = var-epsilon D . In the case of unstable and corresponding damped modes it is seen that δ 2 F ≠ var-epsilon D ; in fact δ 2 F ≡ 0. This failure of the dielectric energy expression persists even for arbitrarily small growth and damping rates since var-epsilon D is nonzero in this limit, whereas δ 2 F remains zero. The connection between the new exact energy expression and the at-best approximate var-epsilon D is described. The new expression motivates natural definitions of Hamiltonian action variables and signature. A general linear integral transform is introduced that maps the linear version of the noncanonical Hamiltonian structure, which describes the Vlasov equation, to action-angle (diagonal) form

Deterministic methods for the relativistic Vlasov-Maxwell equations and the Van Allen belts dynamics

International Nuclear Information System (INIS)

Le Bourdiec, S.

2007-03-01

Artificial satellites operate in an hostile radiation environment, the Van Allen radiation belts, which partly condition their reliability and their lifespan. In order to protect them, it is necessary to characterize the dynamics of the energetic electrons trapped in these radiation belts. This dynamics is essentially determined by the interactions between the energetic electrons and the existing electromagnetic waves. This work consisted in designing a numerical scheme to solve the equations modelling these interactions: the relativistic Vlasov-Maxwell system of equations. Our choice was directed towards methods of direct integration. We propose three new spectral methods for the momentum discretization: a Galerkin method and two collocation methods. All of them are based on scaled Hermite functions. The scaling factor is chosen in order to obtain the proper velocity resolution. We present in this thesis the discretization of the one-dimensional Vlasov-Poisson system and the numerical results obtained. Then we study the possible extensions of the methods to the complete relativistic problem. In order to reduce the computing time, parallelization and optimization of the algorithms were carried out. Finally, we present 1Dx-3Dv (mono-dimensional for x and three-dimensional for velocity) computations of Weibel and whistler instabilities with one or two electrons species. (author)

Fourier–Hermite spectral representation for the Vlasov–Poisson system in the weakly collisional limit

KAUST Repository

Parker, Joseph T.

2015-02-03

Copyright © Cambridge University Press 2015. We study Landau damping in the 1+1D Vlasov-Poisson system using a Fourier-Hermite spectral representation. We describe the propagation of free energy in Fourier-Hermite phase space using forwards and backwards propagating Hermite modes recently developed for gyrokinetic theory. We derive a free energy equation that relates the change in the electric field to the net Hermite flux out of the zeroth Hermite mode. In linear Landau damping, decay in the electric field corresponds to forward propagating Hermite modes; in nonlinear damping, the initial decay is followed by a growth phase characterized by the generation of backwards propagating Hermite modes by the nonlinear term. The free energy content of the backwards propagating modes increases exponentially until balancing that of the forward propagating modes. Thereafter there is no systematic net Hermite flux, so the electric field cannot decay and the nonlinearity effectively suppresses Landau damping. These simulations are performed using the fully-spectral 5D gyrokinetics code SpectroGK, modified to solve the 1+1D Vlasov-Poisson system. This captures Landau damping via Hou-Li filtering in velocity space. Therefore the code is applicable even in regimes where phase mixing and filamentation are dominant.

Irreversible energy flow in forced Vlasov dynamics

KAUST Repository

Plunk, Gabriel G.; Parker, Joseph T.

2014-01-01

© EDP Sciences, Società Italiana di Fisica, Springer-Verlag. The recent paper of Plunk [G.G. Plunk, Phys. Plasmas 20, 032304 (2013)] considered the forced linear Vlasov equation as a model for the quasi-steady state of a single stable plasma wavenumber interacting with a bath of turbulent fluctuations. This approach gives some insight into possible energy flows without solving for nonlinear dynamics. The central result of the present work is that the forced linear Vlasov equation exhibits asymptotically zero (irreversible) dissipation to all orders under a detuning of the forcing frequency and the characteristic frequency associated with particle streaming. We first prove this by direct calculation, tracking energy flow in terms of certain exact conservation laws of the linear (collisionless) Vlasov equation. Then we analyze the steady-state solutions in detail using a weakly collisional Hermite-moment formulation, and compare with numerical solution. This leads to a detailed description of the Hermite energy spectrum, and a proof of no dissipation at all orders, complementing the collisionless Vlasov result.

Irreversible energy flow in forced Vlasov dynamics

KAUST Repository

Plunk, Gabriel G.

2014-10-01

© EDP Sciences, Società Italiana di Fisica, Springer-Verlag. The recent paper of Plunk [G.G. Plunk, Phys. Plasmas 20, 032304 (2013)] considered the forced linear Vlasov equation as a model for the quasi-steady state of a single stable plasma wavenumber interacting with a bath of turbulent fluctuations. This approach gives some insight into possible energy flows without solving for nonlinear dynamics. The central result of the present work is that the forced linear Vlasov equation exhibits asymptotically zero (irreversible) dissipation to all orders under a detuning of the forcing frequency and the characteristic frequency associated with particle streaming. We first prove this by direct calculation, tracking energy flow in terms of certain exact conservation laws of the linear (collisionless) Vlasov equation. Then we analyze the steady-state solutions in detail using a weakly collisional Hermite-moment formulation, and compare with numerical solution. This leads to a detailed description of the Hermite energy spectrum, and a proof of no dissipation at all orders, complementing the collisionless Vlasov result.

The energy of perturbations for Vlasov plasmas

International Nuclear Information System (INIS)

Morrison, P.J.

1994-02-01

The energy content of electrostatic perturbations about homogeneous equilibria is discussed. The calculation leading to the well-known dielectric (or as it is sometimes called the wave) energy is revisited and interpreted in light of Vlasov theory. It is argued that this quantity is deficient because resonant particles are not correctly handled. A linear integral transform is presented that solves the linear Vlasov-Poisson equation. This solution together with the Kruskal-Oberman energy [Phys. Fluids 1, 275 (1958)] is used to obtain an energy expression in terms of the electric field [Phys. Fluids B 4, 3038 (1992)]. It is described how the integral transform amounts to a change to normal coordinates in an infinite dimensional Hamiltonian system

Deterministic methods for the relativistic Vlasov-Maxwell equations and the Van Allen belts dynamics; Methodes deterministes de resolution des equations de Vlasov-Maxwell relativistes en vue du calcul de la dynamique des ceintures de Van Allen

Energy Technology Data Exchange (ETDEWEB)

Le Bourdiec, S

2007-03-15

Artificial satellites operate in an hostile radiation environment, the Van Allen radiation belts, which partly condition their reliability and their lifespan. In order to protect them, it is necessary to characterize the dynamics of the energetic electrons trapped in these radiation belts. This dynamics is essentially determined by the interactions between the energetic electrons and the existing electromagnetic waves. This work consisted in designing a numerical scheme to solve the equations modelling these interactions: the relativistic Vlasov-Maxwell system of equations. Our choice was directed towards methods of direct integration. We propose three new spectral methods for the momentum discretization: a Galerkin method and two collocation methods. All of them are based on scaled Hermite functions. The scaling factor is chosen in order to obtain the proper velocity resolution. We present in this thesis the discretization of the one-dimensional Vlasov-Poisson system and the numerical results obtained. Then we study the possible extensions of the methods to the complete relativistic problem. In order to reduce the computing time, parallelization and optimization of the algorithms were carried out. Finally, we present 1Dx-3Dv (mono-dimensional for x and three-dimensional for velocity) computations of Weibel and whistler instabilities with one or two electrons species. (author)

Nonadiabatic quantum Vlasov equation for Schwinger pair production

International Nuclear Information System (INIS)

Kim, Sang Pyo; Schubert, Christian

2011-01-01

Using Lewis-Riesenfeld theory, we derive an exact nonadiabatic master equation describing the time evolution of the QED Schwinger pair-production rate for a general time-varying electric field. This equation can be written equivalently as a first-order matrix equation, as a Vlasov-type integral equation, or as a third-order differential equation. In the last version it relates to the Korteweg-de Vries equation, which allows us to construct an exact solution using the well-known one-soliton solution to that equation. The case of timelike delta function pulse fields is also briefly considered.

Non-linear free streaming in Vlasov plasma

Czech Academy of Sciences Publication Activity Database

Sedláček, Zdeněk

2004-01-01

Roč. 54, suppl.C (2004), C82-C88 ISSN 0011-4626. [Symposium on Plasma Physics and Technology/21th./. Prague, 14.06.2004-17.06.2004] Institutional research plan: CEZ:AV0Z2043910 Keywords : plasma oscillations * Vlasov equation Subject RIV: BL - Plasma and Gas Discharge Physics Impact factor: 0.292, year: 2004

Mathematical and numerical methods for Vlasov-Maxwell equations: the contributions of data mining

International Nuclear Information System (INIS)

Assous, F.; Chaskalovic, J.

2014-01-01

There exist a lot of formulations that can model plasma physics or particle accelerators problems as the Vlasov- Maxwell equations. This paper deals with the applications of data mining techniques in the evaluation of numerical solutions of Vlasov-Maxwell models. This is part of the topic of characterizing the model and approximation errors via learning techniques. We give two examples of application. The first one aims at comparing two Vlasov-Maxwell approximate models. In the second one, a scheme based on data mining techniques is proposed to characterize the errors between a P1 and a P2 finite element Particle-In-Cell approach. Beyond these examples, this original approach should operate in all cases where intricate numerical simulations like for the Vlasov-Maxwell equations take a central part. (authors)

On the Magnetic Shield for a Vlasov-Poisson Plasma

Science.gov (United States)

Caprino, Silvia; Cavallaro, Guido; Marchioro, Carlo

2017-12-01

We study the screening of a bounded body Γ against the effect of a wind of charged particles, by means of a shield produced by a magnetic field which becomes infinite on the border of Γ . The charged wind is modeled by a Vlasov-Poisson plasma, the bounded body by a torus, and the external magnetic field is taken close to the border of Γ . We study two models: a plasma composed by different species with positive or negative charges, and finite total mass of each species, and another made of many species of the same sign, each having infinite mass. We investigate the time evolution of both systems, showing in particular that the plasma particles cannot reach the body. Finally we discuss possible extensions to more general initial data. We show also that when the magnetic lines are straight lines, (that imposes an unbounded body), the previous results can be improved.

On Poisson Nonlinear Transformations

Directory of Open Access Journals (Sweden)

Nasir Ganikhodjaev

2014-01-01

Full Text Available We construct the family of Poisson nonlinear transformations defined on the countable sample space of nonnegative integers and investigate their trajectory behavior. We have proved that these nonlinear transformations are regular.

Relativistic simulation of the Vlasov equation for plasma expansion into vacuum

Directory of Open Access Journals (Sweden)

H Abbasi

2012-12-01

Full Text Available   In this study, relativistic Vlasov simulation of plasma for expansion of collisionless plasma for into vacuum is presented. The model is based on 1+1 dimensional phase space and electrostatic approximation. For this purpose, the electron dynamics is studied by the relativistic Vlasov equation. Regardless of the ions temperature, fluid equations are used for their dynamics. The initial electrons distribution function is the relativistic Maxwellian. The results show that due to the electrons relativistic temperature, the process of the plasma expansion takes place faster, the resulting electric field is stronger and the ions are accelerated to higher velocities, in comparison to the non-relativistic case.

Superexponentially damped Vlasov plasma oscillations in the Fourier transformed velocity space

International Nuclear Information System (INIS)

Sedlacek, Z.; Nocera, L.

2002-01-01

The Landau (exponentially) damped solutions of the Vlasov-Poisson equation Fourier transformed with respect to velocity are genuine eigenmodes corresponding to complex eigenvalues. In addition there exist solutions decaying faster than exponentially which exhibit no oscillatory behaviour. A new characterization is given of the initial conditions that give rise to these solutions together with a numerical demonstration

Nonlinear conformally invariant generalization of the Poisson equation to D>2 dimensions

International Nuclear Information System (INIS)

Milgrom, M.

1997-01-01

I propound a nonlinear generalization of the scalar-field Poisson equation [(var-phi , i var-phi ,i ) D/2-1 var-phi ; k ] ;k ∝ρ, in curved D-dimensional space. It is derivable from the Lagrangian density L D =L f D -Aρ var-phi, with L f D ∝-(var-phi , i var-phi ,i ) D/2 , and ρ the distribution of sources. Specializing to Euclidean spaces, where the field equation is ∇·(|∇ var-phi | D-2 ∇ var-phi)∝ρ, I find that L f D is the only conformally invariant (CI) Lagrangian in D dimensions, containing only first derivatives of var-phi, beside the free Lagrangian (∇ var-phi) 2 , which underlies the Laplace equation. When var-phi is coupled to the sources in the above manner, L D is left as the only CI Lagrangian. The symmetry is one's only recourse in solving this nonlinear theory for some nontrivial configurations. Systems comprising N point charges are special and afford further application of the symmetry. In spite of the CI, the energy function for such a system is not invariant under conformal transformations of the charges' positions. The anomalous transformation properties of the energy stem from effects of the self-energies of the charges. It follows from these that the forces F i on the charges q i at positions r i must satisfy certain constraints beside the vanishing of the net force and net moment: e.g., summation i r i ·F i must equal some given function of the charges. The constraints total (D+1)(D+2)/2, which tallies with the dimension of the conformal group in D dimensions. Among other things I use all these to derive exact expressions for the following quantities: (1) The general two-point-charge force. (Abstract Truncated)

Global existence of classical solutions to the Vlasov-Poisson system in a three dimensional, cosmological setting

OpenAIRE

Rein, Gerhard; Rendall, Alan D.

1993-01-01

The initial value problem for the Vlasov-Poisson system is by now well understood in the case of an isolated system where, by definition, the distribution function of the particles as well as the gravitational potential vanish at spatial infinity. Here we start with homogeneous solutions, which have a spatially constant, non-zero mass density and which describe the mass distribution in a Newtonian model of the universe. These homogeneous states can be constructed explicitly, and we consider d...

Reduced Vlasov-Maxwell simulations

International Nuclear Information System (INIS)

Helluy, P.; Navoret, L.; Pham, N.; Crestetto, A.

2014-01-01

The Maxwell-Vlasov system is a fundamental model in physics. It can be applied to plasma simulations, charged particles beam, astrophysics, etc. The unknowns are the electromagnetic field, solution to the Maxwell equations and the distribution function, solution to the Vlasov equation. In this paper we review two different numerical methods for Vlasov-Maxwell simulations. The first method is based on a coupling between a Discontinuous Galerkin (DG) Maxwell solver and a Particle-In-Cell (PIC) Vlasov solver. The second method only uses a DG approach for the Vlasov and Maxwell equations. The Vlasov equation is first reduced to a space-only hyperbolic system thanks to the finite-element method. The two numerical methods are implemented using OpenCL in order to achieve high performance on recent Graphic Processing Units (GPU). We obtained interesting speedups, but we also observe that the PIC method is the most expensive part of the computation. Therefore we propose another fully Eulerian approach. Thanks to a decomposition of the distribution function on velocity basis functions, we obtain a reduced Vlasov model, which appears to be a hyperbolic system of conservation laws written only in the (x,t) space. We can thus adapt very easily our DG solver to the reduced model

An implicit meshless scheme for the solution of transient non-linear Poisson-type equations

KAUST Repository

Bourantas, Georgios

2013-07-01

A meshfree point collocation method is used for the numerical simulation of both transient and steady state non-linear Poisson-type partial differential equations. Particular emphasis is placed on the application of the linearization method with special attention to the lagging of coefficients method and the Newton linearization method. The localized form of the Moving Least Squares (MLS) approximation is employed for the construction of the shape functions, in conjunction with the general framework of the point collocation method. Computations are performed for regular nodal distributions, stressing the positivity conditions that make the resulting system stable and convergent. The accuracy and the stability of the proposed scheme are demonstrated through representative and well-established benchmark problems. © 2013 Elsevier Ltd.

An implicit meshless scheme for the solution of transient non-linear Poisson-type equations

KAUST Repository

Bourantas, Georgios; Burganos, Vasilis N.

2013-01-01

A meshfree point collocation method is used for the numerical simulation of both transient and steady state non-linear Poisson-type partial differential equations. Particular emphasis is placed on the application of the linearization method with special attention to the lagging of coefficients method and the Newton linearization method. The localized form of the Moving Least Squares (MLS) approximation is employed for the construction of the shape functions, in conjunction with the general framework of the point collocation method. Computations are performed for regular nodal distributions, stressing the positivity conditions that make the resulting system stable and convergent. The accuracy and the stability of the proposed scheme are demonstrated through representative and well-established benchmark problems. © 2013 Elsevier Ltd.

Instability of the filtering method for Vlasov's equation

International Nuclear Information System (INIS)

Figua, H.; Bouchut, F.; Fijalkow, E.

1999-01-01

Klimas has introduced a smoothed Fourier-Fourier method. This method consists in convolving the original distribution function with a Gaussian distribution function, and, next, in solving the new system with a transformed splitting algorithm. Unfortunately, a second-order term appears in the new equation. In this work, it is studied how this term affects the numerical equation. In particular it is proven that instability occurs in the linear version of the Vlasov equation obtained by considering only free non-interacting particles. It is proved that the use of Fourier-Fourier transform is a fundamental requirement to solve this new equation. An important property is pointed out concerning the filtered distribution function in the transformed space. (K.A.)

On the relativistic Vlasov equation in guiding-center coordinates

International Nuclear Information System (INIS)

Salimullah, M.; Chaudhry, M.B.; Hassan, M.H.A.

1989-11-01

The relativistic Vlasov equation has been expressed in terms of the guiding-center coordinates in a hot magnetized plasma. It is noted that the relativistic effect reduces the cyclotron resonance frequency for electrostatic and electromagnetic waves propagating transverse to the direction of the static magnetic field in the plasma. (author). 4 refs

On invariant measures for the Vlasov equation with a regular potential

International Nuclear Information System (INIS)

Zhidkov, P.E.

2003-01-01

We consider a Vlasov equation with a smooth bounded potential of interaction between particles in a class of measure-valued solutions and construct a measure which is invariant for this problem in a sense

Nonlinear wave evolution in VLASOV plasma: a lie-transform analysis

International Nuclear Information System (INIS)

Cary, J.R.

1979-08-01

Nonlinear wave evolution in Vlasov plasma is analyzed using the Lie transform, a powerful mathematical tool which is applicable to Hamiltonian systems. The first part of this thesis is an exposition of the Lie transform. Dewar's general Lie transform theory is explained and is used to construct Deprit's Lie transform perturbation technique. The basic theory is illustrated by simple examples

Contributions to mathematical analysis and to numerical approximation in plasma physics

International Nuclear Information System (INIS)

Besse, N.

2009-01-01

The author's scientific works deal with numerical analysis and the simulation of the partial differential equations that intervene in the transport of charged particles and in plasma physics. In the chapters 2 and 3, a reduction of the Vlasov equation is presented, this method is based on the Liouville geometric invariants and it leads to a mathematical model named water-bag model that can be coupled with various equations of the electromagnetic field: the Poisson equation, the quasi-neutral equation or Maxwell equations. In the chapter 3 this reduction method is applied to the Vlasov gyro-kinetic equation to form the gyro-water-bag model. The mathematical analysis of this model produces interesting analytical results such as: threshold instabilities, instability growth rate, transport coefficient and non-linear turbulence mechanisms. Simulations have been performed to study turbulence in magnetized plasmas. In these plasmas occurred numerous instabilities due to the presence of high density and temperature gradients. These instabilities generate turbulence that deteriorates plasma confinement conditions required for thermonuclear fusion. The numerical calculation of turbulent thermal diffusivities is important since confinement time is determined by these transport coefficients. The chapter 4 gathers mathematical analysis issues like convergence or prior knowledge of errors concerning several high-order numerical methods used to solve Vlasov-Poisson or Vlasov-Einstein equation systems as well as the induction equation of an idealistic MHD system. The chapter 5 presents original numerical methods to solve several non-linear Vlasov equations such as Vlasov-Poisswell, Vlasov-Darwin, Vlasov-Maxwell and Vlasov-gyrokinetic that are involved either in inertial fusion or in magnetic confinement fusion

Nonlinear behavior of a monochromatic wave in a one-dimensional Vlasov plasma

International Nuclear Information System (INIS)

Shoucri, M.M.; Gagne, R.R.J.

1978-01-01

The nonlinear evolution of a monochromatic wave in a one-dimensional Vlasov plasma is studied numerically. The numerical results are carried out far enough in time for phase mixing to dominate the asymptotic state of the system. A qualitative comparison with previously reported simulations is given

Equal-Time and Equal-Space Poisson Brackets of the N -Component Coupled NLS Equation

International Nuclear Information System (INIS)

Zhou Ru-Guang; Li Pei-Yao; Gao Yuan

2017-01-01

Two Poisson brackets for the N-component coupled nonlinear Schrödinger (NLS) equation are derived by using the variantional principle. The first one is called the equal-time Poisson bracket which does not depend on time but only on the space variable. Actually it is just the usual one describing the time evolution of system in the traditional theory of integrable Hamiltonian systems. The second one is equal-space and new. It is shown that the spatial part of Lax pair with respect to the equal-time Poisson bracket and temporal part of Lax pair with respect to the equal-space Poisson bracket share the same r-matrix formulation. These properties are similar to that of the NLS equation. (paper)

Fully Electromagnetic Nonlinear Gyrokinetic Equations for Tokamak Edge Turbulence

International Nuclear Information System (INIS)

Hahm, T.S.; Wang, Lu; Madsen, J.

2008-01-01

An energy conserving set of the fully electromagnetic nonlinear gyrokinetic Vlasov equation and Maxwell's equations, which is applicable to both L-mode turbulence with large amplitude and H-mode turbulence in the presence of high E x B shear has been derived. The phase-space action variational Lie perturbation method ensures the preservation of the conservation laws of the underlying Vlasov-Maxwell system. Our generalized ordering takes ρ i θi ∼ L E ∼ L p i is the thermal ion Larmor radius and ρ θi = B/B θ ρ i ), as typically observed in the tokamak H-mode edge, with L E and L p being the radial electric field and pressure gradient lengths. We take k # perpendicular# ρ i ∼ 1 for generality, and keep the relative fluctuation amplitudes e(delta)φ/T i ∼ (delta)B/B up to the second order. Extending the electrostatic theory in the presence of high E x B shear [Hahm, Phys. Plasmas 3, 4658 (1996)], contributions of electromagnetic fluctuations to the particle charge density and current are explicitly evaluated via pull-back transformation from the gyrocenter distribution function in the gyrokinetic Maxwell's equation

Local WKB dispersion relation for the Vlasov-Maxwell equations

International Nuclear Information System (INIS)

Berk, H.L.; Dominguez, R.R.

1982-10-01

A formalism for analyzing systems of integral equations, based on the Wentzel-Kramers-Brillouin (WKB) approximation, is applied to the Vlasov-Maxwell integral equations in an arbitrary-β, spatially inhomogenous plasma model. It is shown that when treating frequencies comparable with and larger than the cyclotron frequency, relevant new terms must be accounted for to treat waves that depend upon local spatial gradients. For a specific model, the response for very short wavelength and high frequency is shown to reduce to the straight-line orbit approximation when the WKB rules are correctly followed

Scaling of the distribution function and the critical exponents near the point of a marginal stability under the Vlasov-Poisson equations

International Nuclear Information System (INIS)

Ivanov, Alexei

2000-08-01

A model system, described by the consistent Vlasov-Poisson equations under periodical boundary conditions, has been studied numerically near the point of a marginal stability. The power laws, typical for a system, undergoing a second-order phase transition, hold in a vicinity of the critical point: (i) A ∝ -θ β , β=1.907±0.006 for θ ≤ 0, where A is the saturated amplitude of the marginally-stable mode; (ii) χ ∝ θ -γ as θ → 0, γ=γ - =1.020±0.008 for θ + =0.995±0.020 for θ > 0, where χ=∂A/∂F 1 at F 1 → 0 is the susceptibility to external drive of the strain F 1 ; (iii) at θ=0 the system responds to external drive as A ∝ F 1 1/δ , and δ=1.544±0.002. θ=( 2 >- cr 2 >)/ cr 2 > is the dimensionless reduced velocity dispersion. Within the error of computation these critical exponents satisfy to equality γ=β(δ-1), known in thermodynamics as the Widom equality, which is direct consequence of scaling invariance of the Fourier components f m of the distribution function f at |θ| m (λ at t, λ av v, λ aθ θ, λ aA0 A 0 , λ aF F 1 )=λf m (t, v, θ, A 0 , F 1 ) at θ approx. = 0. On the contrary to thermodynamics these critical indices indicate to a very wide critical area. In turn, it means that critical phenomena may determine macroscopic dynamics of a large fraction of systems. (author)

Controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with delay and Poisson jumps

Directory of Open Access Journals (Sweden)

Diem Dang Huan

2015-12-01

Full Text Available The current paper is concerned with the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps in Hilbert spaces. Using the theory of a strongly continuous cosine family of bounded linear operators, stochastic analysis theory and with the help of the Banach fixed point theorem, we derive a new set of sufficient conditions for the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps. Finally, an application to the stochastic nonlinear wave equation with infinite delay and Poisson jumps is given.

Single particle dynamics of many-body systems described by Vlasov-Fokker-Planck equations

International Nuclear Information System (INIS)

Frank, T.D.

2003-01-01

Using Langevin equations we describe the random walk of single particles that belong to particle systems satisfying Vlasov-Fokker-Planck equations. In doing so, we show that Haissinski distributions of bunched particles in electron storage rings can be derived from a particle dynamics model

Fully nonlinear phenomenology of the Berk-Breizman augmentation of the Vlasov-Maxwell system

International Nuclear Information System (INIS)

Vann, R.G.L.; Dendy, R.O.; Rowlands, G.; Arber, T.D.; D'Ambrumenil, N.

2003-01-01

The Berk-Breizman augmentation of the Vlasov-Maxwell system is widely used to model self-consistent resonant excitation and damping of wave fields by evolving energetic particle populations in magnetic fusion plasmas. The key model parameters are the particle annihilation rate ν a , which drives bump-on-tail structure, and the linear wave damping rate γ d . A code, based on the piecewise parabolic method, is used to integrate the fully nonlinear Berk-Breizman system of equations across the whole (ν a ,γ d ) parameter space. The results of this code show that the system's behavior can be classified into one of four types, each of which occurs in a well-defined region of parameter space: chaotic, periodic, steady state, and damped. The corresponding evolution in (x,v) phase space is also examined

Numerical Integration of the Vlasov Equation of Two Colliding Beams

CERN Document Server

Zorzano-Mier, M P

2000-01-01

In a circular collider the motion of particles of one beam is strongly perturbed at the interaction points by the electro-magnetic field associated with the counter-rotating beam. For any two arbitrary initial particle distributions the time evolution of the two beams can be known by solving the coupled system of two Vlasov equations. This collective description is mandatory when the two beams have similar strengths, as in the case of LEP or LHC. The coherent modes excited by this beam-beam interaction can be a strong limitation for the operation of LHC. In this work, the coupled Vlasov equations of two colliding flat beams are solved numerically using a finite difference scheme. The results suggest that, for the collision of beams with equal tunes, the tune shift between the $\\sigma$- and $\\pi$- coherent dipole mode depends on the unperturbed tune $q$ because of the deformation that the so-called dynamic beta effect induces on the beam distribution. Only when the unperturbed tune $q\\rightarrow 0.25$ this tun...

Numerical study of a Vlasov equation for systems with interacting particles

Energy Technology Data Exchange (ETDEWEB)

Herrera, Dianela; Curilef, Sergio [Departamento de Física, Universidad Católica del Norte, Avenida Angamos 0610, Antofagasta (Chile)

2015-03-10

We solve numerically the Vlasov equation for the self-gravitating sheet model. We used the method introduced by Cheng and Knorr [Comput Phys 22, 330-351 (1976)]. We discuss the quasi-stationary state for some thermodynamical observables, specifically the kinetic energy, whose trend is depicted for early evolution.

From nonlinear Schroedinger hierarchy to some (2+1)-dimensional nonlinear pseudodifferential equations

International Nuclear Information System (INIS)

Yang Xiao; Du Dianlou

2010-01-01

The Poisson structure on C N xR N is introduced to give the Hamiltonian system associated with a spectral problem which yields the nonlinear Schroedinger (NLS) hierarchy. The Hamiltonian system is proven to be Liouville integrable. Some (2+1)-dimensional equations including NLS equation, Kadomtesev-Petviashvili I (KPI) equation, coupled KPI equation, and modified Kadomtesev-Petviashvili (mKP) equation, are decomposed into Hamilton flows via the NLS hierarchy. The algebraic curve, Abel-Jacobi coordinates, and Riemann-Jacobi inversion are used to obtain the algebrogeometric solutions of these equations.

Scaling of the distribution function and the critical exponents near the point of a marginal stability under the Vlasov-Poisson equations

Energy Technology Data Exchange (ETDEWEB)

Lvanov, Alexei [Theory and Computer Simulation Center, National Inst. for Fusion Science, Toki, Gifu (Japan)

2000-08-01

A model system, described by the consistent Vlasov-Poisson equations under periodical boundary conditions, has been studied numerically near the point of a marginal stability. The power laws, typical for a system, undergoing a second-order phase transition, hold in a vicinity of the critical point: (i) A {proportional_to} -{theta}{sup {beta}}, {beta}=1.907{+-}0.006 for {theta} {<=} 0, where A is the saturated amplitude of the marginally-stable mode; (ii) {chi} {proportional_to} {theta}{sup -{gamma}} as {theta} {yields} 0, {gamma}={gamma}{sub -}=1.020{+-}0.008 for {theta} < 0, and {gamma}={gamma}{sub +}=0.995{+-}0.020 for {theta} > 0, where {chi}={partial_derivative}A/{partial_derivative}F{sub 1} at F{sub 1} {yields} 0 is the susceptibility to external drive of the strain F{sub 1}; (iii) at {theta}=0 the system responds to external drive as A {proportional_to} F{sub 1}{sup 1/{delta}}, and {delta}=1.544{+-}0.002. {theta}=(-)/ is the dimensionless reduced velocity dispersion. Within the error of computation these critical exponents satisfy to equality {gamma}={beta}({delta}-1), known in thermodynamics as the Widom equality, which is direct consequence of scaling invariance of the Fourier components f{sub m} of the distribution function f at |{theta}| << 1, i.e. f{sub m}({lambda}{sup at}t, {lambda}{sup av}v, {lambda}{sup a{theta}}{theta}, {lambda}{sup aA0}A{sub 0}, {lambda}{sup aF}F{sub 1})={lambda}f{sub m}(t, v, {theta}, A{sub 0}, F{sub 1}) at {theta} approx. = 0. On the contrary to thermodynamics these critical indices indicate to a very wide critical area. In turn, it means that critical phenomena may determine macroscopic dynamics of a large fraction of systems. (author)

Non-Linear Excitation of Ion Acoustic Waves

DEFF Research Database (Denmark)

Michelsen, Poul; Hirsfield, J. L.

1974-01-01

The excitation of ion acoustic waves by nonlinear coupling of two transverse magnetic waves generated in a microwave cavity was investigated. Measurements of the wave amplitude showed good agreement with calculations based on the Vlasov equation.......The excitation of ion acoustic waves by nonlinear coupling of two transverse magnetic waves generated in a microwave cavity was investigated. Measurements of the wave amplitude showed good agreement with calculations based on the Vlasov equation....

AQUASOL: An efficient solver for the dipolar Poisson-Boltzmann-Langevin equation.

Science.gov (United States)

Koehl, Patrice; Delarue, Marc

2010-02-14

The Poisson-Boltzmann (PB) formalism is among the most popular approaches to modeling the solvation of molecules. It assumes a continuum model for water, leading to a dielectric permittivity that only depends on position in space. In contrast, the dipolar Poisson-Boltzmann-Langevin (DPBL) formalism represents the solvent as a collection of orientable dipoles with nonuniform concentration; this leads to a nonlinear permittivity function that depends both on the position and on the local electric field at that position. The differences in the assumptions underlying these two models lead to significant differences in the equations they generate. The PB equation is a second order, elliptic, nonlinear partial differential equation (PDE). Its response coefficients correspond to the dielectric permittivity and are therefore constant within each subdomain of the system considered (i.e., inside and outside of the molecules considered). While the DPBL equation is also a second order, elliptic, nonlinear PDE, its response coefficients are nonlinear functions of the electrostatic potential. Many solvers have been developed for the PB equation; to our knowledge, none of these can be directly applied to the DPBL equation. The methods they use may adapt to the difference; their implementations however are PBE specific. We adapted the PBE solver originally developed by Holst and Saied [J. Comput. Chem. 16, 337 (1995)] to the problem of solving the DPBL equation. This solver uses a truncated Newton method with a multigrid preconditioner. Numerical evidences suggest that it converges for the DPBL equation and that the convergence is superlinear. It is found however to be slow and greedy in memory requirement for problems commonly encountered in computational biology and computational chemistry. To circumvent these problems, we propose two variants, a quasi-Newton solver based on a simplified, inexact Jacobian and an iterative self-consistent solver that is based directly on the PBE

VLASOVIA 2013 - International workshop on the theory and applications of Vlasov equation - Slides of the presentations

International Nuclear Information System (INIS)

Aunai, N.; Belmont, G.; Smets, R.; Chandre, C.; Tassi, E.; Morrison, P.J.; Back, A.; Guillebon, L. de; Qin, H.; Squire, J.; Tang, W.M.; Garbet, X.; Esteve, D.; Sarazin, Y.; Abiteboul, J.; Bourdelle, C.; Dif-Pradalier, G.; Ghendrih, P.; Grandgirard, V.; Latu, G.; Smolyakov, A.; Hervieux, P.A.; Manfredi, G.; Jasiak, R.; Kraus, M.; Mora, P.; Morel, P.; Dreydemy Ghiro, F.; Berionni, V.; Gurcan, O.D.; Morrison, P.J.; Negulescu, C.; Pegoraro, F.; Bulanov, S.V.; Califano, F.; Fedeli, L.; Grassi, A.; Macchi, A.; Petri, J.; Pezzi, O.; Valentini, F.; Perrone, D.; Veltri, P.; Taccogna, F.; Minelli, P.; Thide, B.; Tamburini, F.; Throumoulopoulos, G.; Tasso, H.

2014-01-01

The Vlasov equation is used for the modelling of a wide range of phenomena occurring in natural and man-made plasmas, as well as in other many-particle systems displaying a collective behaviour. The purpose of this workshop is to bring together scientists to discuss the latest results on Vlasov theory and related applications. The topics discussed include: space plasmas, inertial confinement plasmas, magnetic confinement plasmas, quantum effects in collisionless plasmas, gravitational systems, Hamiltonian Vlasov dynamics, and computational and numerical approaches. This document gathers the slides of the presentations.

Action principles for the Vlasov equation: Four old, one new

International Nuclear Information System (INIS)

Ye, Huanchun; Morrison, P.J.

1991-01-01

Action principles for the Vlasov equation are presented. Four previously known action principles, which differ by the choice of dynamical variables, are described and the interrelationship between them discussed. A new action principle called the leaf action, which manifestly preserves the Casimir invariants and possess a single function as the dynamical variable, is presented. The relationship to the noncanonical Hamiltonian formalism is also explored. 21 refs

Quantum kinetic field theory in curved spacetime: Covariant Wigner function and Liouville-Vlasov equations

International Nuclear Information System (INIS)

Calzetta, E.; Habib, S.; Hu, B.L.

1988-01-01

We consider quantum fields in an external potential and show how, by using the Fourier transform on propagators, one can obtain the mass-shell constraint conditions and the Liouville-Vlasov equation for the Wigner distribution function. We then consider the Hadamard function G 1 (x 1 ,x 2 ) of a real, free, scalar field in curved space. We postulate a form for the Fourier transform F/sup (//sup Q//sup )/(X,k) of the propagator with respect to the difference variable x = x 1 -x 2 on a Riemann normal coordinate centered at Q. We show that F/sup (//sup Q//sup )/ is the result of applying a certain Q-dependent operator on a covariant Wigner function F. We derive from the wave equations for G 1 a covariant equation for the distribution function and show its consistency. We seek solutions to the set of Liouville-Vlasov equations for the vacuum and nonvacuum cases up to the third adiabatic order. Finally we apply this method to calculate the Hadamard function in the Einstein universe. We show that the covariant Wigner function can incorporate certain relevant global properties of the background spacetime. Covariant Wigner functions and Liouville-Vlasov equations are also derived for free fermions in curved spacetime. The method presented here can serve as a basis for constructing quantum kinetic theories in curved spacetime or for near-uniform systems under quasiequilibrium conditions. It can also be useful to the development of a transport theory of quantum fields for the investigation of grand unification and post-Planckian quantum processes in the early Universe

Evolution of the phase-space density and the Jeans scale for dark matter derived from the Vlasov-Einstein equation

International Nuclear Information System (INIS)

Piattella, O.F.; Rodrigues, D.C.; Fabris, J.C.; Pacheco, J.A. de Freitas

2013-01-01

We discuss solutions of Vlasov-Einstein equation for collisionless dark matter particles in the context of a flat Friedmann universe. We show that, after decoupling from the primordial plasma, the dark matter phase-space density indicator Q = ρ/(σ 1D 2 ) 3/2 remains constant during the expansion of the universe, prior to structure formation. This well known result is valid for non-relativistic particles and is not ''observer dependent'' as in solutions derived from the Vlasov-Poisson system. In the linear regime, the inclusion of velocity dispersion effects permits to define a physical Jeans length for collisionless matter as function of the primordial phase-space density indicator: λ J = (5π/G) 1/2 Q −1/3 ρ dm −1/6 . The comoving Jeans wavenumber at matter-radiation equality is smaller by a factor of 2-3 than the comoving wavenumber due to free-streaming, contributing to the cut-off of the density fluctuation power spectrum at the lowest scales. We discuss the physical differences between these two scales. For dark matter particles of mass equal to 200 GeV, the derived Jeans mass is 4.3 × 10 −6 M ⊙

Exact solutions for a system of nonlinear plasma fluid equations

International Nuclear Information System (INIS)

Prahovic, M.G.; Hazeltine, R.D.; Morrison, P.J.

1991-04-01

A method is presented for constructing exact solutions to a system of nonlinear plasma fluid equations that combines the physics of reduced magnetohydrodynamics and the electrostatic drift-wave description of the Charney-Hasegawa-Mima equation. The system has nonlinearities that take the form of Poisson brackets involving the fluid field variables. The method relies on modifying a class of simple equilibrium solutions, but no approximations are made. A distinguishing feature is that the original nonlinear problem is reduced to the solution of two linear partial differential equations, one fourth-order and the other first-order. The first-order equation has Hamiltonian characteristics and is easily integrated, supplying information about the general structure of solutions. 6 refs

Global well posedness of the relativistic Vlasov-Yukawa system with small data

International Nuclear Information System (INIS)

Ha, Seung-Yeal; Lee, Ho

2007-01-01

In this paper, we present an existence theory and uniform L 1 -stability estimate for classical solutions with small data to the Vlasov-Yukawa system. The Vlasov-Yukawa system corresponds to a short-range correction of the Vlasov-Poisson system appearing in plasma physics and astrophysics. For the existence and stability of classical solutions, we crucially use dispersion estimates due to the smallness of data

Vlasov treatment of coherent synchrotron radiation from arbitrary planar orbits

International Nuclear Information System (INIS)

Warnock, R.; Bassi, G.; Ellison, J.A.

2006-01-01

We study the influence of coherent synchrotron radiation (CSR) on particle bunches traveling on arbitrary planar orbits between parallel conducting plates which represent the vacuum chamber. Our goal is to follow the time evolution of the phase space distribution by solving the Vlasov-Maxwell equations in the time domain. This should provide simulations with lower numerical noise than the macro-particle method, and allow one to study such issues as emittance degradation and microbunching due to CSR in bunch compressors. The fields excited by the bunch are computed in the laboratory frame from a new formula that leads to much simpler computations than usual methods. The nonlinear Vlasov equation, formulated in the interaction picture, is integrated in the beam frame by approximating the Perron-Frobenius operator. For application to a chicane bunch compressor we take steps to deal with energy chirp

Integral propagator solvers for Vlasov-Fokker-Planck equations

International Nuclear Information System (INIS)

Donoso, J M; Rio, E del

2007-01-01

We briefly discuss the use of short-time integral propagators on solving the so-called Vlasov-Fokker-Planck equation for the dynamics of a distribution function. For this equation, the diffusion tensor is singular and the usual Gaussian representation of the short-time propagator is no longer valid. However, we prove that the path-integral approach on solving the equation is, in fact, reliable by means of our generalized propagator, which is obtained through the construction of an auxiliary solvable Fokker-Planck equation. The new representation of the grid-free advancing scheme describes the inherent cross- and self-diffusion processes, in both velocity and configuration spaces, in a natural manner, although these processes are not explicitly depicted in the differential equation. We also show that some splitting methods, as well as some finite-difference schemes, could fail in describing the aforementioned diffusion processes, governed in the whole phase space only by the velocity diffusion tensor. The short-time transition probability offers a stable and robust numerical algorithm that preserves the distribution positiveness and its norm, ensuring the smoothness of the evolving solution at any time step. (fast track communication)

On kinetic Boltzmann equations and related hydrodynamic flows with dry viscosity

Directory of Open Access Journals (Sweden)

Nikolai N. Bogoliubov (Jr.

2007-01-01

Full Text Available A two-component particle model of Boltzmann-Vlasov type kinetic equations in the form of special nonlinear integro-differential hydrodynamic systems on an infinite-dimensional functional manifold is discussed. We show that such systems are naturally connected with the nonlinear kinetic Boltzmann-Vlasov equations for some one-dimensional particle flows with pointwise interaction potential between particles. A new type of hydrodynamic two-component Benney equations is constructed and their Hamiltonian structure is analyzed.

From nonlinear Schrödinger hierarchy to some (2+1)-dimensional nonlinear pseudodifferential equations

Science.gov (United States)

Yang, Xiao; Du, Dianlou

2010-08-01

The Poisson structure on CN×RN is introduced to give the Hamiltonian system associated with a spectral problem which yields the nonlinear Schrödinger (NLS) hierarchy. The Hamiltonian system is proven to be Liouville integrable. Some (2+1)-dimensional equations including NLS equation, Kadomtesev-Petviashvili I (KPI) equation, coupled KPI equation, and modified Kadomtesev-Petviashvili (mKP) equation, are decomposed into Hamilton flows via the NLS hierarchy. The algebraic curve, Abel-Jacobi coordinates, and Riemann-Jacobi inversion are used to obtain the algebrogeometric solutions of these equations.

Vlasov simulation of the relativistic effect on the breaking of large amplitude plasma waves

International Nuclear Information System (INIS)

Xu Hui; Sheng Zhengming; Zhang Jie

2007-01-01

The influence of relativistic and thermal effects on plasma wave breaking has been studied by solving the coupled Vlasov-Poisson equations. When the relativistic effect is not considered, the wave breaking will not occur, provided the initial perturbation is less than certain value as predicted previously, and the largest amplitude of the plasma wave will decrease with the increase of the initial temperature. When the relativistic effect is considered, wave breaking always occurs during the time evolution, irrespective of the initial perturbation amplitude. Yet the smaller the initial perturbation amplitude is, the longer is the time for wave breaking to occur. With large initial perturbations, wave breaking can always occur with the without the relativistic effect. However, the results are significantly different in the two cases. The thermal effects of electrons decrease the threshold value to initial amplitude for wave breaking and large phase velocity makes the nonlinear phenomenon occur more easily. (authors)

Vlasov equation for photons and quasi-particles in a plasma

International Nuclear Information System (INIS)

Mendonca, J.T.

2014-01-01

We show that, in quite general conditions, a Vlasov equation can be derived for photons in a medium. The same is true for other quasi-particles, such as plasmons, phonons or driftons, associated with other wave modes in a plasma. The range of validity of this equation is discussed. We also discuss the Landau resonance, and its relation with photon acceleration. Exact and approximate expressions for photon and quasi-particle Landau damping are stated. Photon and quasi-particle acceleration and trapping is also discussed. Specific applications to laser-plasma interaction, and to magnetic fusion turbulence, are considered as illustrations of the general approach. (author)

Nonlinear poisson brackets geometry and quantization

CERN Document Server

Karasev, M V

2012-01-01

This book deals with two old mathematical problems. The first is the problem of constructing an analog of a Lie group for general nonlinear Poisson brackets. The second is the quantization problem for such brackets in the semiclassical approximation (which is the problem of exact quantization for the simplest classes of brackets). These problems are progressively coming to the fore in the modern theory of differential equations and quantum theory, since the approach based on constructions of algebras and Lie groups seems, in a certain sense, to be exhausted. The authors' main goal is to describe in detail the new objects that appear in the solution of these problems. Many ideas of algebra, modern differential geometry, algebraic topology, and operator theory are synthesized here. The authors prove all statements in detail, thus making the book accessible to graduate students.

A high order solver for the unbounded Poisson equation

DEFF Research Database (Denmark)

Hejlesen, Mads Mølholm; Rasmussen, Johannes Tophøj; Chatelain, Philippe

2013-01-01

. The method is extended to directly solve the derivatives of the solution to Poissonʼs equation. In this way differential operators such as the divergence or curl of the solution field can be solved to the same high order convergence without additional computational effort. The method, is applied......A high order converging Poisson solver is presented, based on the Greenʼs function solution to Poissonʼs equation subject to free-space boundary conditions. The high order convergence is achieved by formulating regularised integration kernels, analogous to a smoothing of the solution field...... and validated, however not restricted, to the equations of fluid mechanics, and can be used in many applications to solve Poissonʼs equation on a rectangular unbounded domain....

On Landau Scenario of Chaotization for Beam Distribution

International Nuclear Information System (INIS)

Parsa, Z.; Zadorozhny, V.

1999-01-01

We examine a problem in nonlinear dynamics in which both regular and chaotic motions are possible. Thus we deal with some of the fundamental theoretical problem of accelerator physics, mathematics theory of dynamical systems, and other fields of physics. The focus is on the appearance of chaos in a beam distribution. A study of the problem is based on two observations. The First observation is that using Lyapunov method and its extension we obtain solutions of partial differential equations. Using this approach we discuss the problem of finding a solution of Vlasov-Poisson equation, i.e., some stationary solution where we consider magnetic field as some disturbance with a small parameter. Thus the solution of Vlasov equation yields an asymptotic series such that the solution of Vlasov-Poisson equation is the basis solution for one. The second observation is that physical chaos is weakly limit of, well known, the Landau bifurcation's. This fact we have proved using ideas on the Nature of Turbulence

On classical solutions of the relativistic Vlasov-Klein-Gordon system

Directory of Open Access Journals (Sweden)

Michael Kunzinger

2005-01-01

Full Text Available We consider a collisionless ensemble of classical particles coupled with a Klein-Gordon field. For the resulting nonlinear system of partial differential equations, the relativistic Vlasov-Klein-Gordon system, we prove local-in-time existence of classical solutions and a continuation criterion which says that a solution can blow up only if the particle momenta become large. We also show that classical solutions are global in time in the one-dimensional case.

Continuum Vlasov Simulation in Four Phase-space Dimensions

Science.gov (United States)

Cohen, B. I.; Banks, J. W.; Berger, R. L.; Hittinger, J. A.; Brunner, S.

2010-11-01

In the VALHALLA project, we are developing scalable algorithms for the continuum solution of the Vlasov-Maxwell equations in two spatial and two velocity dimensions. We use fourth-order temporal and spatial discretizations of the conservative form of the equations and a finite-volume representation to enable adaptive mesh refinement and nonlinear oscillation control [1]. The code has been implemented with and without adaptive mesh refinement, and with electromagnetic and electrostatic field solvers. A goal is to study the efficacy of continuum Vlasov simulations in four phase-space dimensions for laser-plasma interactions. We have verified the code in examples such as the two-stream instability, the weak beam-plasma instability, Landau damping, electron plasma waves with electron trapping and nonlinear frequency shifts [2]^ extended from 1D to 2D propagation, and light wave propagation.^ We will report progress on code development, computational methods, and physics applications. This work was performed under the auspices of the U.S. DOE by LLNL under contract no. DE-AC52-07NA27344. This work was funded by the Lab. Dir. Res. and Dev. Prog. at LLNL under project tracking code 08-ERD-031. [1] J.W. Banks and J.A.F. Hittinger, to appear in IEEE Trans. Plas. Sci. (Sept., 2010). [2] G.J. Morales and T.M. O'Neil, Phys. Rev. Lett. 28,417 (1972); R. L. Dewar, Phys. Fluids 15,712 (1972).

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

International Nuclear Information System (INIS)

Fisicaro, G.; Goedecker, S.; Genovese, L.; Andreussi, O.; Marzari, N.

2016-01-01

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments.

Science.gov (United States)

Fisicaro, G; Genovese, L; Andreussi, O; Marzari, N; Goedecker, S

2016-01-07

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

Energy Technology Data Exchange (ETDEWEB)

Fisicaro, G., E-mail: giuseppe.fisicaro@unibas.ch; Goedecker, S. [Department of Physics, University of Basel, Klingelbergstrasse 82, 4056 Basel (Switzerland); Genovese, L. [University of Grenoble Alpes, CEA, INAC-SP2M, L-Sim, F-38000 Grenoble (France); Andreussi, O. [Institute of Computational Science, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6904 Lugano (Switzerland); Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, Station 12, CH-1015 Lausanne (Switzerland); Marzari, N. [Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, Station 12, CH-1015 Lausanne (Switzerland)

2016-01-07

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.

Equations of motion of test particles for solving the spin-dependent Boltzmann–Vlasov equation

Energy Technology Data Exchange (ETDEWEB)

Xia, Yin [Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800 (China); University of Chinese Academy of Science, Beijing 100049 (China); Xu, Jun, E-mail: xujun@sinap.ac.cn [Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800 (China); Li, Bao-An [Department of Physics and Astronomy, Texas A& M University-Commerce, Commerce, TX 75429-3011 (United States); Department of Applied Physics, Xi' an Jiao Tong University, Xi' an 710049 (China); Shen, Wen-Qing [Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800 (China)

2016-08-10

A consistent derivation of the equations of motion (EOMs) of test particles for solving the spin-dependent Boltzmann–Vlasov equation is presented. The resulting EOMs in phase space are similar to the canonical equations in Hamiltonian dynamics, and the EOM of spin is the same as that in the Heisenburg picture of quantum mechanics. Considering further the quantum nature of spin and choosing the direction of total angular momentum in heavy-ion reactions as a reference of measuring nucleon spin, the EOMs of spin-up and spin-down nucleons are given separately. The key elements affecting the spin dynamics in heavy-ion collisions are identified. The resulting EOMs provide a solid foundation for using the test-particle approach in studying spin dynamics in heavy-ion collisions at intermediate energies. Future comparisons of model simulations with experimental data will help to constrain the poorly known in-medium nucleon spin–orbit coupling relevant for understanding properties of rare isotopes and their astrophysical impacts.

Development of a Generalized Version of the Poisson-Nernst-Planck Equations Using the Hybrid Mixture Theory: Presentation of 2D Numerical Examples

DEFF Research Database (Denmark)

Johannesson, Björn

2010-01-01

A numerical scheme for the transient solution of generalized version of the Poisson-Nernst-Planck equations is presented. The finite element method is used to establish the coupled non-linear matrix system of equations capable of solving the present problem iteratively. The Poisson......-scale and that it includes the volume fractions of phases in its structure. The background to the Poisson-Nernst-Planck equations can by the HMT approach be described by using the postulates of mass conservation of constituents together with the Gauss’ law used together with consistent constitutive laws. The HMT theory......-Nernst-Planck equations represent a set of diffusion equations for charged species, i.e. dissolved ions. These equations are coupled to the ‘internally’ induced electrical field and to the velocity field of the fluid. The Nernst-Planck equations describing the diffusion of the ionic species and the Gauss’ law in used are...

Modified Poisson eigenfunctions for electrostatic Bernstein--Greene--Kruskal equilibria

International Nuclear Information System (INIS)

Ling, K.; Abraham-Shrauner, B.

1981-01-01

The stability of an electrostatic Bernstein--Greene--Kruskal equilibrium by Lewis and Symon's general linear stability analysis for spatially inhomogeneous Vlasov equilibria, which employs eigenfunctions and eigenvalues of the equilibrium Liouville operator and the modified Poisson operator, is considered. Analytic expressions for the Liouville eigenfuctions and eigenvalues have already been given; approximate analytic expressions for the dominant eigenfunction and eigenvalue of the modified Poisson operator are given. In the kinetic limit three methods are given: (i) the perturbation method, (ii) the Rayleigh--Ritz method, and (iii) a method based on a Hill's equation. In the fluid limit the Rayleigh--Ritz method is used. The dominant eigenfunction and eigenvalue are then substituted in the dispersion relation and the growth rate calculated. The growth rate agrees very well with previous results found by numerical simulation and by modified Poisson eigenfunctions calculated numerically

Truncated Painleve expansion: Tanh-traveling wave solutions and reduction of sine-Poisson equation to a quadrature for stationary and nonstationary three-dimensional collisionless cold plasma

International Nuclear Information System (INIS)

Ibrahim, R. S.; El-Kalaawy, O. H.

2006-01-01

The relativistic nonlinear self-consistent equations for a collisionless cold plasma with stationary ions [R. S. Ibrahim, IMA J. Appl. Math. 68, 523 (2003)] are extended to 3 and 3+1 dimensions. The resulting system of equations is reduced to the sine-Poisson equation. The truncated Painleve expansion and reduction of the partial differential equation to a quadrature problem (RQ method) are described and applied to obtain the traveling wave solutions of the sine-Poisson equation for stationary and nonstationary equations in 3 and 3+1 dimensions describing the charge-density equilibrium configuration model

The Einstein-Vlasov System/Kinetic Theory

Directory of Open Access Journals (Sweden)

Håkan Andréasson

2002-12-01

Full Text Available The main purpose of this article is to provide a guide to theorems on global properties of solutions to the Einstein-Vlasov system. This system couples Einstein's equations to a kinetic matter model. Kinetic theory has been an important field of research during several decades in which the main focus has been on nonrelativistic and special relativistic physics, i.e., to model the dynamics of neutral gases, plasmas, and Newtonian self-gravitating systems. In 1990, Rendall and Rein initiated a mathematical study of the Einstein-Vlasov system. Since then many theorems on global properties of solutions to this system have been established. The Vlasov equation describes matter phenomenologically, and it should be stressed that most of the theorems presented in this article are not presently known for other such matter models (i.e., fluid models. This paper gives introductions to kinetic theory in non-curved spacetimes and then the Einstein-Vlasov system is introduced. We believe that a good understanding of kinetic theory in non-curved spacetimes is fundamental to good comprehension of kinetic theory in general relativity.

An adaptive, implicit, conservative, 1D-2V multi-species Vlasov-Fokker-Planck multi-scale solver in planar geometry

Science.gov (United States)

Taitano, W. T.; Chacón, L.; Simakov, A. N.

2018-07-01

We consider a 1D-2V Vlasov-Fokker-Planck multi-species ionic description coupled to fluid electrons. We address temporal stiffness with implicit time stepping, suitably preconditioned. To address temperature disparity in time and space, we extend the conservative adaptive velocity-space discretization scheme proposed in [Taitano et al., J. Comput. Phys., 318, 391-420, (2016)] to a spatially inhomogeneous system. In this approach, we normalize the velocity-space coordinate to a temporally and spatially varying local characteristic speed per species. We explicitly consider the resulting inertial terms in the Vlasov equation, and derive a discrete formulation that conserves mass, momentum, and energy up to a prescribed nonlinear tolerance upon convergence. Our conservation strategy employs nonlinear constraints to enforce these properties discretely for both the Vlasov operator and the Fokker-Planck collision operator. Numerical examples of varying degrees of complexity, including shock-wave propagation, demonstrate the favorable efficiency and accuracy properties of the scheme.

Numerical solution of the Maxwell-Vlasov equations in the periodic regime. Application to the study of isotope separation by ion cyclotron resonance; Resolution numerique des equations de Maxwell-Vlasov en regime periodique. Application a l'etude de la separation isotopique par resonance cyclotron ionique

Energy Technology Data Exchange (ETDEWEB)

Omnes, P

1999-01-25

This work is dedicated to the study of the behaviour of a magnetic confined plasma that is excited by a purely sinusoidal electric current delivered by an antenna. The response of the electrons to the electromagnetic field is considered as linear,whereas the ions of the plasma are represented by a non-relativistic Vlasov equation. In order to avoid transients, the coupled Maxwell-Vlasov equations are solved in a periodic mode and in a bounded domain. An equivalent electric conductivity tensor has been defined, this tensor is a linear operator that links the electric current generated by the movement of the particles to the electromagnetic field. Theoretical considerations can assure the existence and uniqueness of a periodical solution to Vlasov equations and of a solution to Maxwell equations in harmonic mode. The system of equations is periodical and has been solved by using an iterative method. The application of this method to the simulation of a isotopic separation device based on ionic cyclotron resonance has shown that the convergence is reached in a few iterations and that the solution is valid. Furthermore a method based on a finite-volume formulation of Maxwell equations in the time domain is presented. 2 new variables are defined in order to better take into account the Gauss' law and the conservation of the magnetic flux, the new system is still hyperbolic. The parallelization of the process has been successfully realized. (A.C.)

Solving (2 + 1)-dimensional sine-Poisson equation by a modified variable separated ordinary differential equation method

International Nuclear Information System (INIS)

Ka-Lin, Su; Yuan-Xi, Xie

2010-01-01

By introducing a more general auxiliary ordinary differential equation (ODE), a modified variable separated ordinary differential equation method is presented for solving the (2 + 1)-dimensional sine-Poisson equation. As a result, many explicit and exact solutions of the (2 + 1)-dimensional sine-Poisson equation are derived in a simple manner by this technique. (general)

The Poisson equation on Klein surfaces

Directory of Open Access Journals (Sweden)

Monica Rosiu

2016-04-01

Full Text Available We obtain a formula for the solution of the Poisson equation with Dirichlet boundary condition on a region of a Klein surface. This formula reveals the symmetric character of the solution.

Multi-dimensional, fully-implicit, spectral method for the Vlasov-Maxwell equations with exact conservation laws in discrete form

Science.gov (United States)

Delzanno, G. L.

2015-11-01

A spectral method for the numerical solution of the multi-dimensional Vlasov-Maxwell equations is presented. The plasma distribution function is expanded in Fourier (for the spatial part) and Hermite (for the velocity part) basis functions, leading to a truncated system of ordinary differential equations for the expansion coefficients (moments) that is discretized with an implicit, second order accurate Crank-Nicolson time discretization. The discrete non-linear system is solved with a preconditioned Jacobian-Free Newton-Krylov method. It is shown analytically that the Fourier-Hermite method features exact conservation laws for total mass, momentum and energy in discrete form. Standard tests involving plasma waves and the whistler instability confirm the validity of the conservation laws numerically. The whistler instability test also shows that we can step over the fastest time scale in the system without incurring in numerical instabilities. Some preconditioning strategies are presented, showing that the number of linear iterations of the Krylov solver can be drastically reduced and a significant gain in performance can be obtained.

Comparison of two forms of Vlasov-type relativistic kinetic equations in hadrodynamics

International Nuclear Information System (INIS)

Mashnik, S.G.; Maino, G.

1996-01-01

A comparison of two methods in the relativistic kinetic theory of the Fermi systems is carried out assuming, as an example, the simplest σω-version of quantum hadrodynamics with allowance for strong mean meson fields. It is shown that the Vlasov-type relativistic kinetic equation (VRKE) obtained by means of the procedure of squaring at an intermediate step is responsible for unphysical features. A direct method of derivation of kinetic equations is proposed. This method does not contain such drawback and gives rise to VRKE in hydrodynamics of a non-contradictory form in which both spin degrees of freedom and states with positive and negative energies are taken into account. 17 refs

Selective Contrast Adjustment by Poisson Equation

Directory of Open Access Journals (Sweden)

Ana-Belen Petro

2013-09-01

Full Text Available Poisson Image Editing is a new technique permitting to modify the gradient vector field of an image, and then to recover an image with a gradient approaching this modified gradient field. This amounts to solve a Poisson equation, an operation which can be efficiently performed by Fast Fourier Transform (FFT. This paper describes an algorithm applying this technique, with two different variants. The first variant enhances the contrast by increasing the gradient in the dark regions of the image. This method is well adapted to images with back light or strong shadows, and reveals details in the shadows. The second variant of the same Poisson technique enhances all small gradients in the image, thus also sometimes revealing details and texture.

A discontinuous Poisson-Boltzmann equation with interfacial jump: homogenisation and residual error estimate.

Science.gov (United States)

Fellner, Klemens; Kovtunenko, Victor A

2016-01-01

A nonlinear Poisson-Boltzmann equation with inhomogeneous Robin type boundary conditions at the interface between two materials is investigated. The model describes the electrostatic potential generated by a vector of ion concentrations in a periodic multiphase medium with dilute solid particles. The key issue stems from interfacial jumps, which necessitate discontinuous solutions to the problem. Based on variational techniques, we derive the homogenisation of the discontinuous problem and establish a rigorous residual error estimate up to the first-order correction.

Nonlinear evolution equations

CERN Document Server

Uraltseva, N N

1995-01-01

This collection focuses on nonlinear problems in partial differential equations. Most of the papers are based on lectures presented at the seminar on partial differential equations and mathematical physics at St. Petersburg University. Among the topics explored are the existence and properties of solutions of various classes of nonlinear evolution equations, nonlinear imbedding theorems, bifurcations of solutions, and equations of mathematical physics (Navier-Stokes type equations and the nonlinear Schrödinger equation). The book will be useful to researchers and graduate students working in p

Fast and Accurate Poisson Denoising With Trainable Nonlinear Diffusion.

Science.gov (United States)

Feng, Wensen; Qiao, Peng; Chen, Yunjin; Wensen Feng; Peng Qiao; Yunjin Chen; Feng, Wensen; Chen, Yunjin; Qiao, Peng

2018-06-01

The degradation of the acquired signal by Poisson noise is a common problem for various imaging applications, such as medical imaging, night vision, and microscopy. Up to now, many state-of-the-art Poisson denoising techniques mainly concentrate on achieving utmost performance, with little consideration for the computation efficiency. Therefore, in this paper we aim to propose an efficient Poisson denoising model with both high computational efficiency and recovery quality. To this end, we exploit the newly developed trainable nonlinear reaction diffusion (TNRD) model which has proven an extremely fast image restoration approach with performance surpassing recent state-of-the-arts. However, the straightforward direct gradient descent employed in the original TNRD-based denoising task is not applicable in this paper. To solve this problem, we resort to the proximal gradient descent method. We retrain the model parameters, including the linear filters and influence functions by taking into account the Poisson noise statistics, and end up with a well-trained nonlinear diffusion model specialized for Poisson denoising. The trained model provides strongly competitive results against state-of-the-art approaches, meanwhile bearing the properties of simple structure and high efficiency. Furthermore, our proposed model comes along with an additional advantage, that the diffusion process is well-suited for parallel computation on graphics processing units (GPUs). For images of size , our GPU implementation takes less than 0.1 s to produce state-of-the-art Poisson denoising performance.

Solving the Fluid Pressure Poisson Equation Using Multigrid-Evaluation and Improvements.

Science.gov (United States)

Dick, Christian; Rogowsky, Marcus; Westermann, Rudiger

2016-11-01

In many numerical simulations of fluids governed by the incompressible Navier-Stokes equations, the pressure Poisson equation needs to be solved to enforce mass conservation. Multigrid solvers show excellent convergence in simple scenarios, yet they can converge slowly in domains where physically separated regions are combined at coarser scales. Moreover, existing multigrid solvers are tailored to specific discretizations of the pressure Poisson equation, and they cannot easily be adapted to other discretizations. In this paper we analyze the convergence properties of existing multigrid solvers for the pressure Poisson equation in different simulation domains, and we show how to further improve the multigrid convergence rate by using a graph-based extension to determine the coarse grid hierarchy. The proposed multigrid solver is generic in that it can be applied to different kinds of discretizations of the pressure Poisson equation, by using solely the specification of the simulation domain and pre-assembled computational stencils. We analyze the proposed solver in combination with finite difference and finite volume discretizations of the pressure Poisson equation. Our evaluations show that, despite the common assumption, multigrid schemes can exploit their potential even in the most complicated simulation scenarios, yet this behavior is obtained at the price of higher memory consumption.

The Lie-Poisson structure of integrable classical non-linear sigma models

International Nuclear Information System (INIS)

Bordemann, M.; Forger, M.; Schaeper, U.; Laartz, J.

1993-01-01

The canonical structure of classical non-linear sigma models on Riemannian symmetric spaces, which constitute the most general class of classical non-linear sigma models known to be integrable, is shown to be governed by a fundamental Poisson bracket relation that fits into the r-s-matrix formalism for non-ultralocal integrable models first discussed by Maillet. The matrices r and s are computed explicitly and, being field dependent, satisfy fundamental Poisson bracket relations of their own, which can be expressed in terms of a new numerical matrix c. It is proposed that all these Poisson brackets taken together are representation conditions for a new kind of algebra which, for this class of models, replaces the classical Yang-Baxter algebra governing the canonical structure of ultralocal models. The Poisson brackets for the transition matrices are also computed, and the notorious regularization problem associated with the definition of the Poisson brackets for the monodromy matrices is discussed. (orig.)

Nonlinear Dirac Equations

Directory of Open Access Journals (Sweden)

Wei Khim Ng

2009-02-01

Full Text Available We construct nonlinear extensions of Dirac's relativistic electron equation that preserve its other desirable properties such as locality, separability, conservation of probability and Poincaré invariance. We determine the constraints that the nonlinear term must obey and classify the resultant non-polynomial nonlinearities in a double expansion in the degree of nonlinearity and number of derivatives. We give explicit examples of such nonlinear equations, studying their discrete symmetries and other properties. Motivated by some previously suggested applications we then consider nonlinear terms that simultaneously violate Lorentz covariance and again study various explicit examples. We contrast our equations and construction procedure with others in the literature and also show that our equations are not gauge equivalent to the linear Dirac equation. Finally we outline various physical applications for these equations.

Poisson-Nernst-Planck Equations for Simulating Biomolecular Diffusion-Reaction Processes I: Finite Element Solutions.

Science.gov (United States)

Lu, Benzhuo; Holst, Michael J; McCammon, J Andrew; Zhou, Y C

2010-09-20

In this paper we developed accurate finite element methods for solving 3-D Poisson-Nernst-Planck (PNP) equations with singular permanent charges for electrodiffusion in solvated biomolecular systems. The electrostatic Poisson equation was defined in the biomolecules and in the solvent, while the Nernst-Planck equation was defined only in the solvent. We applied a stable regularization scheme to remove the singular component of the electrostatic potential induced by the permanent charges inside biomolecules, and formulated regular, well-posed PNP equations. An inexact-Newton method was used to solve the coupled nonlinear elliptic equations for the steady problems; while an Adams-Bashforth-Crank-Nicolson method was devised for time integration for the unsteady electrodiffusion. We numerically investigated the conditioning of the stiffness matrices for the finite element approximations of the two formulations of the Nernst-Planck equation, and theoretically proved that the transformed formulation is always associated with an ill-conditioned stiffness matrix. We also studied the electroneutrality of the solution and its relation with the boundary conditions on the molecular surface, and concluded that a large net charge concentration is always present near the molecular surface due to the presence of multiple species of charged particles in the solution. The numerical methods are shown to be accurate and stable by various test problems, and are applicable to real large-scale biophysical electrodiffusion problems.

A high order solver for the unbounded Poisson equation

DEFF Research Database (Denmark)

Hejlesen, Mads Mølholm; Rasmussen, Johannes Tophøj; Chatelain, Philippe

In mesh-free particle methods a high order solution to the unbounded Poisson equation is usually achieved by constructing regularised integration kernels for the Biot-Savart law. Here the singular, point particles are regularised using smoothed particles to obtain an accurate solution with an order...... of convergence consistent with the moments conserved by the applied smoothing function. In the hybrid particle-mesh method of Hockney and Eastwood (HE) the particles are interpolated onto a regular mesh where the unbounded Poisson equation is solved by a discrete non-cyclic convolution of the mesh values...... and the integration kernel. In this work we show an implementation of high order regularised integration kernels in the HE algorithm for the unbounded Poisson equation to formally achieve an arbitrary high order convergence. We further present a quantitative study of the convergence rate to give further insight...

Vlasov Treatment of Coherent Synchrotron Radiation from Arbitrary Planar Orbits

International Nuclear Information System (INIS)

Warnock, R

2004-01-01

We study the influence of coherent synchrotron radiation (CSR) on particle bunches traveling on arbitrary planar orbits between parallel conducting plates. The plates represent shielding due to the vacuum chamber. The vertical distribution of charge is an arbitrary fixed function. Our goal is to follow the time evolution of the phase space distribution by solving the Vlasov-Maxwell equations in the time domain. This provides simulations with lower numerical noise than the macroparticle method, and allows one to study such issues as emittance degradation and microbunching due to CSR in bunch compressors. The fields excited by the bunch are computed in the laboratory frame from a new formula that leads to much simpler computations than the usual retarded potentials or Lienard-Wiechert potentials. The nonlinear Vlasov equation, formulated in the interaction picture, is integrated in the beam frame by approximating the Perron-Frobenius operator. The distribution function is represented by B-splines, in a scheme preserving positivity and normalization of the distribution. For application to a chicane bunch compressor we take steps to deal with energy chirp, an initial near-perfect correlation of energy with position in the bunch

The Vlasov equation with strong magnetic field and oscillating electric field as a model for isotop resonant separation

Directory of Open Access Journals (Sweden)

Emmanuel Frenod

2002-01-01

Full Text Available We study the qualitative behavior of solutions to the Vlasov equation with strong external magnetic field and oscillating electric field. This model is relevant to the understanding of isotop resonant separation. We show that the effective equation is a kinetic equation with a memory term. This memory term involves a pseudo-differential operator whose kernel is characterized by an integral equation involving Bessel functions. The kernel is explicitly given in some particular cases.

Poisson equation for weak gravitational lensing

International Nuclear Information System (INIS)

Kling, Thomas P.; Campbell, Bryan

2008-01-01

Using the Newman and Penrose [E. T. Newman and R. Penrose, J. Math. Phys. (N.Y.) 3, 566 (1962).] spin-coefficient formalism, we examine the full Bianchi identities of general relativity in the context of gravitational lensing, where the matter and space-time curvature are projected into a lens plane perpendicular to the line of sight. From one component of the Bianchi identity, we provide a rigorous, new derivation of a Poisson equation for the projected matter density where the source term involves second derivatives of the observed weak gravitational lensing shear. We also show that the other components of the Bianchi identity reveal no new results. Numerical integration of the Poisson equation in test cases shows an accurate mass map can be constructed from the combination of a ground-based, wide-field image and a Hubble Space Telescope image of the same system

Transient finite element analysis of electric double layer using Nernst-Planck-Poisson equations with a modified Stern layer.

Science.gov (United States)

Lim, Jongil; Whitcomb, John; Boyd, James; Varghese, Julian

2007-01-01

A finite element implementation of the transient nonlinear Nernst-Planck-Poisson (NPP) and Nernst-Planck-Poisson-modified Stern (NPPMS) models is presented. The NPPMS model uses multipoint constraints to account for finite ion size, resulting in realistic ion concentrations even at high surface potential. The Poisson-Boltzmann equation is used to provide a limited check of the transient models for low surface potential and dilute bulk solutions. The effects of the surface potential and bulk molarity on the electric potential and ion concentrations as functions of space and time are studied. The ability of the models to predict realistic energy storage capacity is investigated. The predicted energy is much more sensitive to surface potential than to bulk solution molarity.

Nonlinear Gyrokinetic Theory With Polarization Drift

International Nuclear Information System (INIS)

Wang, L.; Hahm, T.S.

2010-01-01

A set of the electrostatic toroidal gyrokinetic Vlasov equation and the Poisson equation, which explicitly includes the polarization drift, is derived systematically by using Lie-transform method. The polarization drift is introduced in the gyrocenter equations of motion, and the corresponding polarization density is derived. Contrary to the wide-spread expectation, the inclusion of the polarization drift in the gyrocenter equations of motion does not affect the expression for the polarization density significantly. This is due to modification of the gyrocenter phase-space volume caused by the electrostatic potential [T. S. Hahm, Phys. Plasmas 3, 4658 (1996)].

Hamiltonian dynamics of spatially-homogeneous Vlasov-Einstein systems

International Nuclear Information System (INIS)

Okabe, Takahide; Morrison, P. J.; Friedrichsen, J. E. III; Shepley, L. C.

2011-01-01

We introduce a new matter action principle, with a wide range of applicability, for the Vlasov equation in terms of a conjugate pair of functions. Here we apply this action principle to the study of matter in Bianchi cosmological models in general relativity. The Bianchi models are spatially-homogeneous solutions to the Einstein field equations, classified by the three-dimensional Lie algebra that describes the symmetry group of the model. The Einstein equations for these models reduce to a set of coupled ordinary differential equations. The class A Bianchi models admit a Hamiltonian formulation in which the components of the metric tensor and their time derivatives yield the canonical coordinates. The evolution of anisotropy in the vacuum Bianchi models is determined by a potential due to the curvature of the model, according to its symmetry. For illustrative purposes, we examine the evolution of anisotropy in models with Vlasov matter. The Vlasov content is further simplified by the assumption of cold, counter-streaming matter, a kind of matter that is far from thermal equilibrium and is not describable by an ordinary fluid model nor other more simplistic matter models. Qualitative differences and similarities are found in the dynamics of certain vacuum class A Bianchi models and Bianchi type I models with cold, counter-streaming Vlasov-matter potentials analogous to the curvature potentials of corresponding vacuum models.

Nonlinear Delta-f Particle Simulations of Collective Effects in High-Intensity Bunched Beams

CERN Document Server

Qin, Hong; Hudson, Stuart R; Startsev, Edward

2005-01-01

The collective effects in high-intensity 3D bunched beams are described self-consistently by the nonlinear Vlasov-Maxwell equations.* The nonlinear delta-f method,** a particle simulation method for solving the nonlinear Vlasov-Maxwell equations, is being used to study the collective effects in high-intensity 3D bunched beams. The delta-f method, as a nonlinear perturbative scheme, splits the distribution function into equilibrium and perturbed parts. The perturbed distribution function is represented as a weighted summation over discrete particles, where the particle orbits are advanced by equations of motion in the focusing field and self-consistent fields, and the particle weights are advanced by the coupling between the perturbed fields and the zero-order distribution function. The nonlinear delta-f method exhibits minimal noise and accuracy problems in comparison with standard particle-in-cell simulations. A self-consistent 3D kinetic equilibrium is first established for high intensity bunched beams. The...

Diagnosing collisionless energy transfer using field-particle correlations: Vlasov-Poisson plasmas

Science.gov (United States)

Howes, Gregory G.; Klein, Kristopher G.; Li, Tak Chu

2017-02-01

Turbulence plays a key role in the conversion of the energy of large-scale fields and flows to plasma heat, impacting the macroscopic evolution of the heliosphere and other astrophysical plasma systems. Although we have long been able to make direct spacecraft measurements of all aspects of the electromagnetic field and plasma fluctuations in near-Earth space, our understanding of the physical mechanisms responsible for the damping of the turbulent fluctuations in heliospheric plasmas remains incomplete. Here we propose an innovative field-particle correlation technique that can be used to measure directly the secular energy transfer from fields to particles associated with collisionless damping of the turbulent fluctuations. Furthermore, this novel procedure yields information about the collisionless energy transfer as a function of particle velocity, providing vital new information that can help to identify the dominant collisionless mechanism governing the damping of the turbulent fluctuations. Kinetic plasma theory is used to devise the appropriate correlation to diagnose Landau damping, and the field-particle correlation technique is thoroughly illustrated using the simplified case of the Landau damping of Langmuir waves in a 1D-1V (one dimension in physical space and one dimension in velocity space) Vlasov-Poisson plasma. Generalizations necessary to apply the field-particle correlation technique to diagnose the collisionless damping of turbulent fluctuations in the solar wind are discussed, highlighting several caveats. This novel field-particle correlation technique is intended to be used as a primary analysis tool for measurements from current, upcoming and proposed spacecraft missions that are focused on the kinetic microphysics of weakly collisional heliospheric plasmas, including the Magnetospheric Multiscale (MMS), Solar Probe Plus, Solar Orbiter and Turbulence Heating ObserveR (THOR) missions.

Wigner transformation in curved space-time and the curvature correction of the Vlasov equation for semiclassical gravitating systems

International Nuclear Information System (INIS)

Winter, J.

1985-01-01

A covariant generalization of the Wigner transformation of quantum equations is proposed for gravitating many-particle systems, which modifies the Einstein-Liouville equations for the coupled gravity-matter problem by inclusion of quantum effects of the matter moving in its self-consistent classical gravitational field, in order to extend their realm of validity to higher particle densities. The corrections of the Vlasov equation (Liouville equation in one-particle phase space) are exhibited as combined effects of quantum mechanics and the curvature of space-time arranged in a semiclassical expansion in powers of h 2 , the first-order term of which is explicitly calculated. It is linear in the Riemann tensor and in its gradient; the Riemann tensor occurs in a similar position as the tensor of the Yang-Mills field strength in a corresponding Vlasov equation for systems with local gauge invariance in the purely classical limit. The performance of the Wigner transformation is based on expressing the equation of motion for the two-point function of the Klein-Gordon field, in particular the Beltrami operator, in terms of a midpoint and a distance vector covariantly defined for the two points. This implies the calculation of deviations of the geodesic between these points, the standard concept of which has to be refined to include infinitesimal variations of the second order. A differential equation for the second-order deviation is established

Numerical solution of the Maxwell-Vlasov equations in the periodic regime. Application to the study of isotope separation by ion cyclotron resonance

International Nuclear Information System (INIS)

Omnes, P.

1999-01-01

This work is dedicated to the study of the behaviour of a magnetic confined plasma that is excited by a purely sinusoidal electric current delivered by an antenna. The response of the electrons to the electromagnetic field is considered as linear, whereas the ions of the plasma are represented by a non-relativistic Vlasov equation. In order to avoid transients, the coupled Maxwell-Vlasov equations are solved in a periodic mode and in a bounded domain. An equivalent electric conductivity tensor has been defined, this tensor is a linear operator that links the electric current generated by the movement of the particles to the electromagnetic field. Theoretical considerations can assure the existence and uniqueness of a periodical solution to Vlasov equations and of a solution to Maxwell equations in harmonic mode. The system of equations is periodical and has been solved by using an iterative method. The application of this method to the simulation of a isotopic separation device based on ionic cyclotron resonance has shown that the convergence is reached in a few iterations and that the solution is valid. Furthermore a method based on a finite-volume formulation of Maxwell equations in the time domain is presented. 2 new variables are defined in order to better take into account the Gauss' law and the conservation of the magnetic flux, the new system is still hyperbolic. The parallelization of the process has been successfully realized. (A.C.)

A Stabilized Finite Element Method for Modified Poisson-Nernst-Planck Equations to Determine Ion Flow Through a Nanopore

Science.gov (United States)

Chaudhry, Jehanzeb Hameed; Comer, Jeffrey; Aksimentiev, Aleksei; Olson, Luke N.

2013-01-01

The conventional Poisson-Nernst-Planck equations do not account for the finite size of ions explicitly. This leads to solutions featuring unrealistically high ionic concentrations in the regions subject to external potentials, in particular, near highly charged surfaces. A modified form of the Poisson-Nernst-Planck equations accounts for steric effects and results in solutions with finite ion concentrations. Here, we evaluate numerical methods for solving the modified Poisson-Nernst-Planck equations by modeling electric field-driven transport of ions through a nanopore. We describe a novel, robust finite element solver that combines the applications of the Newton's method to the nonlinear Galerkin form of the equations, augmented with stabilization terms to appropriately handle the drift-diffusion processes. To make direct comparison with particle-based simulations possible, our method is specifically designed to produce solutions under periodic boundary conditions and to conserve the number of ions in the solution domain. We test our finite element solver on a set of challenging numerical experiments that include calculations of the ion distribution in a volume confined between two charged plates, calculations of the ionic current though a nanopore subject to an external electric field, and modeling the effect of a DNA molecule on the ion concentration and nanopore current. PMID:24363784

A Combined MPI-CUDA Parallel Solution of Linear and Nonlinear Poisson-Boltzmann Equation

Directory of Open Access Journals (Sweden)

José Colmenares

2014-01-01

Full Text Available The Poisson-Boltzmann equation models the electrostatic potential generated by fixed charges on a polarizable solute immersed in an ionic solution. This approach is often used in computational structural biology to estimate the electrostatic energetic component of the assembly of molecular biological systems. In the last decades, the amount of data concerning proteins and other biological macromolecules has remarkably increased. To fruitfully exploit these data, a huge computational power is needed as well as software tools capable of exploiting it. It is therefore necessary to move towards high performance computing and to develop proper parallel implementations of already existing and of novel algorithms. Nowadays, workstations can provide an amazing computational power: up to 10 TFLOPS on a single machine equipped with multiple CPUs and accelerators such as Intel Xeon Phi or GPU devices. The actual obstacle to the full exploitation of modern heterogeneous resources is efficient parallel coding and porting of software on such architectures. In this paper, we propose the implementation of a full Poisson-Boltzmann solver based on a finite-difference scheme using different and combined parallel schemes and in particular a mixed MPI-CUDA implementation. Results show great speedups when using the two schemes, achieving an 18.9x speedup using three GPUs.

An exterior Poisson solver using fast direct methods and boundary integral equations with applications to nonlinear potential flow

Science.gov (United States)

Young, D. P.; Woo, A. C.; Bussoletti, J. E.; Johnson, F. T.

1986-01-01

A general method is developed combining fast direct methods and boundary integral equation methods to solve Poisson's equation on irregular exterior regions. The method requires O(N log N) operations where N is the number of grid points. Error estimates are given that hold for regions with corners and other boundary irregularities. Computational results are given in the context of computational aerodynamics for a two-dimensional lifting airfoil. Solutions of boundary integral equations for lifting and nonlifting aerodynamic configurations using preconditioned conjugate gradient are examined for varying degrees of thinness.

Natural Poisson structures of nonlinear plasma dynamics

International Nuclear Information System (INIS)

Kaufman, A.N.

1982-01-01

Hamiltonian field theories, for models of nonlinear plasma dynamics, require a Poisson bracket structure for functionals of the field variables. These are presented, applied, and derived for several sets of field variables: coherent waves, incoherent waves, particle distributions, and multifluid electrodynamics. Parametric coupling of waves and plasma yields concise expressions for ponderomotive effects (in kinetic and fluid models) and for induced scattering. (Auth.)

Natural Poisson structures of nonlinear plasma dynamics

International Nuclear Information System (INIS)

Kaufman, A.N.

1982-06-01

Hamiltonian field theories, for models of nonlinear plasma dynamics, require a Poisson bracket structure for functionals of the field variables. These are presented, applied, and derived for several sets of field variables: coherent waves, incoherent waves, particle distributions, and multifluid electrodynamics. Parametric coupling of waves and plasma yields concise expressions for ponderomotive effects (in kinetic and fluid models) and for induced scattering

Nonlinear turbulence theory and simulation of Buneman instability

International Nuclear Information System (INIS)

Yoon, P. H.; Umeda, T.

2010-01-01

In the present paper, the weak turbulence theory for reactive instabilities, formulated in a companion paper [P. H. Yoon, Phys. Plasmas 17, 112316 (2010)], is applied to the strong electron-ion two-stream (or Buneman) instability. The self-consistent theory involves quasilinear velocity space diffusion equation for the particles and nonlinear wave kinetic equation that includes quasilinear (or induced emission) term as well as nonlinear wave-particle interaction term (or a term that represents an induced scattering off ions). We have also performed one-dimensional electrostatic Vlasov simulation in order to benchmark the theoretical analysis. Under the assumption of self-similar drifting Gaussian distribution function for the electrons it is shown that the current reduction and the accompanying electron heating as well as electric field turbulence generation can be discussed in a self-consistent manner. Upon comparison with the Vlasov simulation result it is found that quasilinear wave kinetic equation alone is insufficient to account for the final saturation amplitude. Upon including the nonlinear scattering term in the wave kinetic equation, however, we find that a qualitative agreement with the simulation is recovered. From this, we conclude that the combined quasilinear particle diffusion plus induced emission and scattering (off ions) processes adequately account for the nonlinear development of the Buneman instability.

Poisson Autoregression

DEFF Research Database (Denmark)

Fokianos, Konstantinos; Rahbek, Anders Christian; Tjøstheim, Dag

This paper considers geometric ergodicity and likelihood based inference for linear and nonlinear Poisson autoregressions. In the linear case the conditional mean is linked linearly to its past values as well as the observed values of the Poisson process. This also applies to the conditional...... variance, implying an interpretation as an integer valued GARCH process. In a nonlinear conditional Poisson model, the conditional mean is a nonlinear function of its past values and a nonlinear function of past observations. As a particular example an exponential autoregressive Poisson model for time...

Poisson Autoregression

DEFF Research Database (Denmark)

Fokianos, Konstantinos; Rahbæk, Anders; Tjøstheim, Dag

This paper considers geometric ergodicity and likelihood based inference for linear and nonlinear Poisson autoregressions. In the linear case the conditional mean is linked linearly to its past values as well as the observed values of the Poisson process. This also applies to the conditional...... variance, making an interpretation as an integer valued GARCH process possible. In a nonlinear conditional Poisson model, the conditional mean is a nonlinear function of its past values and a nonlinear function of past observations. As a particular example an exponential autoregressive Poisson model...

Antisymmetrized molecular dynamics of wave packets with stochastic incorporation of Vlasov equation

International Nuclear Information System (INIS)

Ono, Akira; Horiuchi, Hisashi.

1996-01-01

The first purpose of this report is to present an extended AMD model which can generally describe such minor branching processes by removing the restriction on the one-body distribution function. This is done not by generalizing the wave packets to arbitrary single-particle wave functions but by representing the diffused and/or deformed wave packet as an ensemble of Gaussian wave packets. In other words, stochastic displacements are given to the wave packets in phase space so that the ensemble-average of the time evolution of the one-body distribution function is essentially equivalent to the solution of Vlasov equation which does not have any restriction on the shape of wave packets. This new model is called AMD-V. Although AMD-V is equivalent to Vlasov equation in the instantaneous time evolution of the one-body distribution function for an AMD wave function, AMD-V describes the branching into channels and the fluctuation of the mean field which are caused by the spreading or the splitting of the single-particle wave function. The second purpose of this report is to show the drastic effect of this new stochastic process of wave packet splitting on the dynamics of heavy ion collisions, especially in the fragmentation mechanism. We take the 40 Ca + 40 Ca system at the incident energy 35 MeV/nucleon. It will be shown that the reproduction of data by the AMD-V calculation is surprisingly good. We will see that the effect of the wave packet diffusion is crucially important to remove the spurious binary feature of the AMD calculation and to enable the multi-fragment final state. (J.P.N.)

Electromagnetic nonlinear gyrokinetics with polarization drift

International Nuclear Information System (INIS)

Duthoit, F.-X.; Hahm, T. S.; Wang, Lu

2014-01-01

A set of new nonlinear electromagnetic gyrokinetic Vlasov equation with polarization drift and gyrokinetic Maxwell equations is systematically derived by using the Lie-transform perturbation method in toroidal geometry. For the first time, we recover the drift-kinetic expression for parallel acceleration [R. M. Kulsrud, in Basic Plasma Physics, edited by A. A. Galeev and R. N. Sudan (North-Holland, Amsterdam, 1983)] from the nonlinear gyrokinetic equations, thereby bridging a gap between the two formulations. This formalism should be useful in addressing nonlinear ion Compton scattering of intermediate-mode-number toroidal Alfvén eigenmodes for which the polarization current nonlinearity [T. S. Hahm and L. Chen, Phys. Rev. Lett. 74, 266 (1995)] and the usual finite Larmor radius effects should compete

Semiclassical limit and well-posedness of nonlinear Schrodinger-Poisson systems

Directory of Open Access Journals (Sweden)

Hailiang Li

2003-09-01

Full Text Available This paper concerns the well-posedness and semiclassical limit of nonlinear Schrodinger-Poisson systems. We show the local well-posedness and the existence of semiclassical limit of the two models for initial data with Sobolev regularity, before shocks appear in the limit system. We establish the existence of a global solution and show the time-asymptotic behavior of a classical solutions of Schrodinger-Poisson system for a fixed re-scaled Planck constant.

Robust iterative observer for source localization for Poisson equation

KAUST Repository

Majeed, Muhammad Usman

2017-01-05

Source localization problem for Poisson equation with available noisy boundary data is well known to be highly sensitive to noise. The problem is ill posed and lacks to fulfill Hadamards stability criteria for well posedness. In this work, first a robust iterative observer is presented for boundary estimation problem for Laplace equation, and then this algorithm along with the available noisy boundary data from the Poisson problem is used to localize point sources inside a rectangular domain. The algorithm is inspired from Kalman filter design, however one of the space variables is used as time-like. Numerical implementation along with simulation results is detailed towards the end.

Robust iterative observer for source localization for Poisson equation

KAUST Repository

Majeed, Muhammad Usman; Laleg-Kirati, Taous-Meriem

2017-01-01

Source localization problem for Poisson equation with available noisy boundary data is well known to be highly sensitive to noise. The problem is ill posed and lacks to fulfill Hadamards stability criteria for well posedness. In this work, first a robust iterative observer is presented for boundary estimation problem for Laplace equation, and then this algorithm along with the available noisy boundary data from the Poisson problem is used to localize point sources inside a rectangular domain. The algorithm is inspired from Kalman filter design, however one of the space variables is used as time-like. Numerical implementation along with simulation results is detailed towards the end.

Resolution of the Vlasov-Maxwell system by PIC discontinuous Galerkin method on GPU with OpenCL

Directory of Open Access Journals (Sweden)

Crestetto Anaïs

2013-01-01

Full Text Available We present an implementation of a Vlasov-Maxwell solver for multicore processors. The Vlasov equation describes the evolution of charged particles in an electromagnetic field, solution of the Maxwell equations. The Vlasov equation is solved by a Particle-In-Cell method (PIC, while the Maxwell system is computed by a Discontinuous Galerkin method. We use the OpenCL framework, which allows our code to run on multicore processors or recent Graphic Processing Units (GPU. We present several numerical applications to two-dimensional test cases.

Poisson's theorem and integrals of KdV equation

International Nuclear Information System (INIS)

Tasso, H.

1978-01-01

Using Poisson's theorem it is proved that if F = integral sub(-infinity)sup(+infinity) T(u,usub(x),...usub(n,t))dx is an invariant functional of KdV equation, then integral sub(-infinity)sup(+infinity) delta F/delta u dx integral sub(-infinity)sup(+infinity) delta T/delta u dx is also an invariant functional. In the case of a polynomial T, one finds in a simple way the known recursion ΔTr/Δu = Tsub(r-1). This note gives an example of the usefulness of Poisson's theorem. (author)

Poisson's equation in de Sitter space-time

Energy Technology Data Exchange (ETDEWEB)

Pessa, E [Rome Univ. (Italy). Ist. di Matematica

1980-11-01

Based on a suitable generalization of Poisson's equation for de Sitter space-time the form of gravitation's law in 'projective relativity' is examined; it is found that, in the interior case, a small difference with the customary Newtonian law arises. This difference, of a repulsive character, can be very important in cosmological problems.

The Einstein-Vlasov System/Kinetic Theory.

Science.gov (United States)

Andréasson, Håkan

2011-01-01

The main purpose of this article is to provide a guide to theorems on global properties of solutions to the Einstein-Vlasov system. This system couples Einstein's equations to a kinetic matter model. Kinetic theory has been an important field of research during several decades in which the main focus has been on non-relativistic and special relativistic physics, i.e., to model the dynamics of neutral gases, plasmas, and Newtonian self-gravitating systems. In 1990, Rendall and Rein initiated a mathematical study of the Einstein-Vlasov system. Since then many theorems on global properties of solutions to this system have been established. This paper gives introductions to kinetic theory in non-curved spacetimes and then the Einstein-Vlasov system is introduced. We believe that a good understanding of kinetic theory in non-curved spacetimes is fundamental to a good comprehension of kinetic theory in general relativity.

Stability of the trivial solution for linear stochastic differential equations with Poisson white noise

International Nuclear Information System (INIS)

Grigoriu, Mircea; Samorodnitsky, Gennady

2004-01-01

Two methods are considered for assessing the asymptotic stability of the trivial solution of linear stochastic differential equations driven by Poisson white noise, interpreted as the formal derivative of a compound Poisson process. The first method attempts to extend a result for diffusion processes satisfying linear stochastic differential equations to the case of linear equations with Poisson white noise. The developments for the method are based on Ito's formula for semimartingales and Lyapunov exponents. The second method is based on a geometric ergodic theorem for Markov chains providing a criterion for the asymptotic stability of the solution of linear stochastic differential equations with Poisson white noise. Two examples are presented to illustrate the use and evaluate the potential of the two methods. The examples demonstrate limitations of the first method and the generality of the second method

Structural interactions in ionic liquids linked to higher-order Poisson-Boltzmann equations

Science.gov (United States)

Blossey, R.; Maggs, A. C.; Podgornik, R.

2017-06-01

We present a derivation of generalized Poisson-Boltzmann equations starting from classical theories of binary fluid mixtures, employing an approach based on the Legendre transform as recently applied to the case of local descriptions of the fluid free energy. Under specific symmetry assumptions, and in the linearized regime, the Poisson-Boltzmann equation reduces to a phenomenological equation introduced by Bazant et al. [Phys. Rev. Lett. 106, 046102 (2011)], 10.1103/PhysRevLett.106.046102, whereby the structuring near the surface is determined by bulk coefficients.

The non-linear coupled spin 2-spin 3 Cotton equation in three dimensions

Energy Technology Data Exchange (ETDEWEB)

Linander, Hampus; Nilsson, Bengt E.W. [Department of Physics, Theoretical PhysicsChalmers University of Technology, S-412 96 Göteborg (Sweden)

2016-07-05

In the context of three-dimensional conformal higher spin theory we derive, in the frame field formulation, the full non-linear spin 3 Cotton equation coupled to spin 2. This is done by solving the corresponding Chern-Simons gauge theory system of equations, that is, using F=0 to eliminate all auxiliary fields and thus expressing the Cotton equation in terms of just the spin 3 frame field and spin 2 covariant derivatives and tensors (Schouten). In this derivation we neglect the spin 4 and higher spin sectors and approximate the star product commutator by a Poisson bracket. The resulting spin 3 Cotton equation is complicated but can be related to linearized versions in the metric formulation obtained previously by other authors. The expected symmetry (spin 3 “translation”, “Lorentz” and “dilatation”) properties are verified for Cotton and other relevant tensors but some perhaps unexpected features emerge in the process, in particular in relation to the non-linear equations. We discuss the structure of this non-linear spin 3 Cotton equation but its explicit form is only presented here, in an exact but not completely refined version, in appended files obtained by computer algebra methods. Both the frame field and metric formulations are provided.

The Integral Equation Method and the Neumann Problem for the Poisson Equation on NTA Domains

Czech Academy of Sciences Publication Activity Database

Medková, Dagmar

2009-01-01

Roč. 63, č. 21 (2009), s. 227-247 ISSN 0378-620X Institutional research plan: CEZ:AV0Z10190503 Keywords : Poisson equation * Neumann problem * integral equation method Subject RIV: BA - General Mathematics Impact factor: 0.477, year: 2009

Lie-Nambu and Lie-Poisson structures in linear and nonlinear quantum mechanics

International Nuclear Information System (INIS)

Czachor, M.

1996-01-01

Space of density matrices in quantum mechanics can be regarded as a Poisson manifold with the dynamics given by certain Lie-Poisson bracket corresponding to an infinite dimensional Lie algebra. The metric structure associated with this Lie algebra is given by a metric tensor which is not equivalent to the Cartan-Killing metric. The Lie-Poisson bracket can be written in a form involving a generalized (Lie-)Nambu bracket. This bracket can be used to generate a generalized, nonlinear and completely integrable dynamics of density matrices. (author)

Degenerate nonlinear diffusion equations

CERN Document Server

Favini, Angelo

2012-01-01

The aim of these notes is to include in a uniform presentation style several topics related to the theory of degenerate nonlinear diffusion equations, treated in the mathematical framework of evolution equations with multivalued m-accretive operators in Hilbert spaces. The problems concern nonlinear parabolic equations involving two cases of degeneracy. More precisely, one case is due to the vanishing of the time derivative coefficient and the other is provided by the vanishing of the diffusion coefficient on subsets of positive measure of the domain. From the mathematical point of view the results presented in these notes can be considered as general results in the theory of degenerate nonlinear diffusion equations. However, this work does not seek to present an exhaustive study of degenerate diffusion equations, but rather to emphasize some rigorous and efficient techniques for approaching various problems involving degenerate nonlinear diffusion equations, such as well-posedness, periodic solutions, asympt...

The forced nonlinear Schroedinger equation

International Nuclear Information System (INIS)

Kaup, D.J.; Hansen, P.J.

1985-01-01

The nonlinear Schroedinger equation describes the behaviour of a radio frequency wave in the ionosphere near the reflexion point where nonlinear processes are important. A simple model of this phenomenon leads to the forced nonlinear Schroedinger equation in terms of a nonlinear boundary value problem. A WKB analysis of the time evolution equations for the nonlinear Schroedinger equation in the inverse scattering transform formalism gives a crude order of magnitude estimation of the qualitative behaviour of the solutions. This estimation is compared with the numerical solutions. (D.Gy.)

The free energy of Maxwell-Vlasov equilibria

International Nuclear Information System (INIS)

Morrison, P.J.; Pfirsch, D.

1989-10-01

A previously derived expression for the energy of arbitrary perturbations about arbitrary Vlasov-Maxwell equilibria is transformed into a very compact form. The new form is also obtained by a canonical transformation method for solving Vlasov's equation, which is based on Lie group theory. This method is simpler than the one used before and provides better physical insight. Finally a procedure is presented for determining the existence of negative-energy modes. In this context the question of why there is an accessibility constraint for the particles, but not for the fields, is discussed. 16 refs

Generalized master equations for non-Poisson dynamics on networks.

Science.gov (United States)

Hoffmann, Till; Porter, Mason A; Lambiotte, Renaud

2012-10-01

The traditional way of studying temporal networks is to aggregate the dynamics of the edges to create a static weighted network. This implicitly assumes that the edges are governed by Poisson processes, which is not typically the case in empirical temporal networks. Accordingly, we examine the effects of non-Poisson inter-event statistics on the dynamics of edges, and we apply the concept of a generalized master equation to the study of continuous-time random walks on networks. We show that this equation reduces to the standard rate equations when the underlying process is Poissonian and that its stationary solution is determined by an effective transition matrix whose leading eigenvector is easy to calculate. We conduct numerical simulations and also derive analytical results for the stationary solution under the assumption that all edges have the same waiting-time distribution. We discuss the implications of our work for dynamical processes on temporal networks and for the construction of network diagnostics that take into account their nontrivial stochastic nature.

Plasma physics on the TI-85 calculator or Down with supercomputers

International Nuclear Information System (INIS)

Sedlacek, Z.

1998-10-01

In the Fourier transformed velocity space the Vlasov plasma oscillations may be interpreted as a wave propagation process corresponding to an imperfectly trapped (leaking) wave. The Landau damped solutions of the Vlasov-Poisson equation then become genuine Eigenmodes corresponding to complex eigenvalues. To illustrate this new interpretation we solve numerically the Fourier transformed Vlasov-Poisson equation, essentially a perturbed advective equation, on the TI-85 pocket graphics calculator. A program is described, based on the method of lines: A finite-difference scheme is utilized to discretize the transformed equation and the resulting set of ordinary differential equations is then solved in time. The user can choose from several possible finite-difference differentiation schemes differing in the total number of points and the number of downwind points. The resulting evolution of the electric field showing the Landau damped plasma oscillations is displayed on the screen of the calculator. In addition, calculation of the eigenvalues of the Fourier transformed Vlasov-Poisson operator is possible. The user can also experiment with the numerical solution of the advective equation which describes free streaming. (author)

A nonlinear equation for ionic diffusion in a strong binary electrolyte

Science.gov (United States)

Ghosal, Sandip; Chen, Zhen

2010-01-01

The problem of the one-dimensional electro-diffusion of ions in a strong binary electrolyte is considered. The mathematical description, known as the Poisson–Nernst–Planck (PNP) system, consists of a diffusion equation for each species augmented by transport owing to a self-consistent electrostatic field determined by the Poisson equation. This description is also relevant to other important problems in physics, such as electron and hole diffusion across semiconductor junctions and the diffusion of ions in plasmas. If concentrations do not vary appreciably over distances of the order of the Debye length, the Poisson equation can be replaced by the condition of local charge neutrality first introduced by Planck. It can then be shown that both species diffuse at the same rate with a common diffusivity that is intermediate between that of the slow and fast species (ambipolar diffusion). Here, we derive a more general theory by exploiting the ratio of the Debye length to a characteristic length scale as a small asymptotic parameter. It is shown that the concentration of either species may be described by a nonlinear partial differential equation that provides a better approximation than the classical linear equation for ambipolar diffusion, but reduces to it in the appropriate limit. PMID:21818176

Vlasov dynamics of periodically driven systems

Science.gov (United States)

Banerjee, Soumyadip; Shah, Kushal

2018-04-01

Analytical solutions of the Vlasov equation for periodically driven systems are of importance in several areas of plasma physics and dynamical systems and are usually approximated using ponderomotive theory. In this paper, we derive the plasma distribution function predicted by ponderomotive theory using Hamiltonian averaging theory and compare it with solutions obtained by the method of characteristics. Our results show that though ponderomotive theory is relatively much easier to use, its predictions are very restrictive and are likely to be very different from the actual distribution function of the system. We also analyse all possible initial conditions which lead to periodic solutions of the Vlasov equation for periodically driven systems and conjecture that the irreducible polynomial corresponding to the initial condition must only have squares of the spatial and momentum coordinate. The resulting distribution function for other initial conditions is aperiodic and can lead to complex relaxation processes within the plasma.

Introduction to nonlinear dispersive equations

CERN Document Server

Linares, Felipe

2015-01-01

This textbook introduces the well-posedness theory for initial-value problems of nonlinear, dispersive partial differential equations, with special focus on two key models, the Korteweg–de Vries equation and the nonlinear Schrödinger equation. A concise and self-contained treatment of background material (the Fourier transform, interpolation theory, Sobolev spaces, and the linear Schrödinger equation) prepares the reader to understand the main topics covered: the initial-value problem for the nonlinear Schrödinger equation and the generalized Korteweg–de Vries equation, properties of their solutions, and a survey of general classes of nonlinear dispersive equations of physical and mathematical significance. Each chapter ends with an expert account of recent developments and open problems, as well as exercises. The final chapter gives a detailed exposition of local well-posedness for the nonlinear Schrödinger equation, taking the reader to the forefront of recent research. The second edition of Introdu...

Cosmology in one dimension: Vlasov dynamics.

Science.gov (United States)

Manfredi, Giovanni; Rouet, Jean-Louis; Miller, Bruce; Shiozawa, Yui

2016-04-01

Numerical simulations of self-gravitating systems are generally based on N-body codes, which solve the equations of motion of a large number of interacting particles. This approach suffers from poor statistical sampling in regions of low density. In contrast, Vlasov codes, by meshing the entire phase space, can reach higher accuracy irrespective of the density. Here, we perform one-dimensional Vlasov simulations of a long-standing cosmological problem, namely, the fractal properties of an expanding Einstein-de Sitter universe in Newtonian gravity. The N-body results are confirmed for high-density regions and extended to regions of low matter density, where the N-body approach usually fails.

Turbulence in unmagnetized Vlasov plasmas

International Nuclear Information System (INIS)

Kuo, S.P.

1985-01-01

The classical technique of transformation and characteristics is employed to analyze the problem of strong turbulence in unmagnetized plasmas. The effect of resonance broadening and perturbation expansion are treated simultaneously, without time secularities. The renormalization procedure of Dupree and Tetreault is used in the transformed Vlasov equation to analyze the turbulence and to derive explicitly a diffusion equation. Analyses are extended to inhomogeneous plasmas and the relationship between the transformation and ponderomotive force is obtained. (author)

High-Order Finite-Difference Solution of the Poisson Equation Involving Complex Geometries in Embedded Meshes

Science.gov (United States)

Marques, Alexandre; Nave, Jean-Christophe; Rosales, Ruben

2011-11-01

The Poisson equation is of central importance in the description of fluid flows and other physical phenomena. In prior work, Marques, Nave, and Rosales introduced the Correction Function Method (CFM) to obtain fourth-order accurate solutions for the constant coefficient Poisson problem with prescribed jump conditions for the solution and its normal derivative across arbitrary interfaces. Here we combine this method with the ideas introduced by Mayo to solve other Poisson problems involving complex geometries. In summary, we are able to rewrite the problem as a boundary integral equation in terms of a potential distribution over the boundary or interface. The solution of this integral equation is discontinuous across the boundary or interface. Hence, after this integral equation is solved using standard techniques, the potential distribution can be used to determine the jump discontinuities. We are then able to use the CFM to solve the resulting Poisson equation with jump discontinuities. The outcome is a fourth-order accurate scheme to solve general Poisson problems which, over arbitrary geometries, has a cost that is approximately twice that of a fast Poisson solver using FFT on a rectangular geometry of the same size. Details of the method and applications will be presented.

Poisson equation in the Kohn-Sham Coulomb problem

OpenAIRE

Manby, F. R.; Knowles, Peter James

2001-01-01

We apply the Poisson equation to the quantum mechanical Coulomb problem for many-particle systems. By introducing a suitable basis set, the two-electron Coulomb integrals become simple overlaps. This offers the possibility of very rapid linear-scaling treatment of the Coulomb contribution to Kohn-Sham theory.

Exact Solutions to Nonlinear Schroedinger Equation and Higher-Order Nonlinear Schroedinger Equation

International Nuclear Information System (INIS)

Ren Ji; Ruan Hangyu

2008-01-01

We study solutions of the nonlinear Schroedinger equation (NLSE) and higher-order nonlinear Schroedinger equation (HONLSE) with variable coefficients. By considering all the higher-order effect of HONLSE as a new dependent variable, the NLSE and HONLSE can be changed into one equation. Using the generalized Lie group reduction method (GLGRM), the abundant solutions of NLSE and HONLSE are obtained

Parallel computing in plasma physics: Nonlinear instabilities

International Nuclear Information System (INIS)

Pohn, E.; Kamelander, G.; Shoucri, M.

2000-01-01

A Vlasov-Poisson-system is used for studying the time evolution of the charge-separation at a spatial one- as well as a two-dimensional plasma-edge. Ions are advanced in time using the Vlasov-equation. The whole three-dimensional velocity-space is considered leading to very time-consuming four-resp. five-dimensional fully kinetic simulations. In the 1D simulations electrons are assumed to behave adiabatic, i.e. they are Boltzmann-distributed, leading to a nonlinear Poisson-equation. In the 2D simulations a gyro-kinetic approximation is used for the electrons. The plasma is assumed to be initially neutral. The simulations are performed at an equidistant grid. A constant time-step is used for advancing the density-distribution function in time. The time-evolution of the distribution function is performed using a splitting scheme. Each dimension (x, y, υ x , υ y , υ z ) of the phase-space is advanced in time separately. The value of the distribution function for the next time is calculated from the value of an - in general - interstitial point at the present time (fractional shift). One-dimensional cubic-spline interpolation is used for calculating the interstitial function values. After the fractional shifts are performed for each dimension of the phase-space, a whole time-step for advancing the distribution function is finished. Afterwards the charge density is calculated, the Poisson-equation is solved and the electric field is calculated before the next time-step is performed. The fractional shift method sketched above was parallelized for p processors as follows. Considering first the shifts in y-direction, a proper parallelization strategy is to split the grid into p disjoint υ z -slices, which are sub-grids, each containing a different 1/p-th part of the υ z range but the whole range of all other dimensions. Each processor is responsible for performing the y-shifts on a different slice, which can be done in parallel without any communication between

Connection between hydrodynamic, water bag and Vlasov models

International Nuclear Information System (INIS)

Gros, M.; Bertrand, P.; Feix, M.R.

1978-01-01

The connection between hydrodynamic, water bag and Vlasov models is still under consideration with numerical experiments. For long wavelength, slightly non linear excitations and initial preparations such as the usual adiabatic invariant Pn -3 is space independent, the hydrodynamic model is equivalent to the water bag, and for long wavelengths a nice agreement is found with the full numerical solution of the Vlasov equation. For other initial conditions when the water bag cannot be defined, the hydrodynamic approach does not represent the correct behaviour. (author)

A System of Poisson Equations for a Nonconstant Varadhan Functional on a Finite State Space

International Nuclear Information System (INIS)

Cavazos-Cadena, Rolando; Hernandez-Hernandez, Daniel

2006-01-01

Given a discrete-time Markov chain with finite state space and a stationary transition matrix, a system of 'local' Poisson equations characterizing the (exponential) Varadhan's functional J(.) is given. The main results, which are derived for an arbitrary transition structure so that J(.) may be nonconstant, are as follows: (i) Any solution to the local Poisson equations immediately renders Varadhan's functional, and (ii) a solution of the system always exist. The proof of this latter result is constructive and suggests a method to solve the local Poisson equations

The solution of the Poisson-Boltzmann's equation for self-consistent potential of infinite, random, nonlinear and non-uniform system

International Nuclear Information System (INIS)

Rasulova, M.Yu

1998-01-01

A study has been made of a system of charged particles and inhomogeneities randomly distributed in accordance with the same law in the neighborhoods of corresponding sites of a planar crystal lattice. The existence and uniqueness of the solution of the generalized Poisson-Boltzmann's equation for the average self-consistent potential and average density of surface charges are proved. (author)

Comparison of free-streaming ELM formulae to a Vlasov simulation

Energy Technology Data Exchange (ETDEWEB)

Moulton, D., E-mail: david.moulton@cea.fr [CEA, IRFM, F-13108 Saint-Paul Lez Durance (France); Fundamenski, W. [Imperial College of Science, Technology and Medicine, London (United Kingdom); Manfredi, G. [Institut de Physique et Chimie des Matériaux, CNRS and Université de Strasbourg, BP 43, F-67034 Strasbourg (France); Hirstoaga, S. [INRIA Nancy Grand-Est and Institut de Recherche en Mathématiques Avancées, 7 rue René Descartes, F-67084 Strasbourg (France); Tskhakaya, D. [Association EURATOM-ÖAW, University of Innsbruck, A-6020 Innsbruck (Austria)

2013-07-15

The main drawbacks of the original free-streaming equations for edge localised mode transport in the scrape-off layer [W. Fundamenski, R.A. Pitts, Plasma Phys. Control Fusion 48 (2006) 109] are that the plasma potential is not accounted for and that only solutions for ion quantities are considered. In this work, the equations are modified and augmented in order to address these two issues. The new equations are benchmarked against (and justified by) a numerical simulation which solves the Vlasov equation in 1d1v. When the source function due to an edge localised mode is instantaneous, the modified free-streaming ‘impulse response’ equations agree closely with the Vlasov simulation results. When the source has a finite duration in time, the agreement worsens. However, in all cases the match is encouragingly good, thus justifying the applicability of the free-streaming approach.

DL_MG: A Parallel Multigrid Poisson and Poisson-Boltzmann Solver for Electronic Structure Calculations in Vacuum and Solution.

Science.gov (United States)

Womack, James C; Anton, Lucian; Dziedzic, Jacek; Hasnip, Phil J; Probert, Matt I J; Skylaris, Chris-Kriton

2018-03-13

The solution of the Poisson equation is a crucial step in electronic structure calculations, yielding the electrostatic potential-a key component of the quantum mechanical Hamiltonian. In recent decades, theoretical advances and increases in computer performance have made it possible to simulate the electronic structure of extended systems in complex environments. This requires the solution of more complicated variants of the Poisson equation, featuring nonhomogeneous dielectric permittivities, ionic concentrations with nonlinear dependencies, and diverse boundary conditions. The analytic solutions generally used to solve the Poisson equation in vacuum (or with homogeneous permittivity) are not applicable in these circumstances, and numerical methods must be used. In this work, we present DL_MG, a flexible, scalable, and accurate solver library, developed specifically to tackle the challenges of solving the Poisson equation in modern large-scale electronic structure calculations on parallel computers. Our solver is based on the multigrid approach and uses an iterative high-order defect correction method to improve the accuracy of solutions. Using two chemically relevant model systems, we tested the accuracy and computational performance of DL_MG when solving the generalized Poisson and Poisson-Boltzmann equations, demonstrating excellent agreement with analytic solutions and efficient scaling to ∼10 9 unknowns and 100s of CPU cores. We also applied DL_MG in actual large-scale electronic structure calculations, using the ONETEP linear-scaling electronic structure package to study a 2615 atom protein-ligand complex with routinely available computational resources. In these calculations, the overall execution time with DL_MG was not significantly greater than the time required for calculations using a conventional FFT-based solver.

Modeling Microbunching from Shot Noise Using Vlasov Solvers

International Nuclear Information System (INIS)

Venturini, Marco; Venturini, Marco; Zholents, Alexander

2008-01-01

Unlike macroparticle simulations, which are sensitive to unphysical statistical fluctuations when the number of macroparticles is smaller than the bunch population, direct methods for solving the Vlasov equation are free from sampling noise and are ideally suited for studying microbunching instabilities evolving from shot noise. We review a 2D (longitudinal dynamics) Vlasov solver we have recently developed to study the microbunching instability in the beam delivery systems for x-ray FELs and present an application to FERMI(at)Elettra. We discuss, in particular, the impact of the spreader design on microbunching

Electron-acoustic Instability Simulated By Modified Zakharov Equations

Science.gov (United States)

Jásenský, V.; Fiala, V.; Vána, O.; Trávnícek, P.; Hellinger, P.

We present non-linear equations describing processes in plasma when electron - acoustic waves are excited. These waves are present for instance in the vicinity of Earth's bow shock and in the polar ionosphere. Frequently they are excited by an elec- tron beam in a plasma with two electron populations, a cold and hot one. We derive modified Zakharov equations from kinetic theory for such a case together with numer- ical method for solving of this type of equations. Bispectral analysis is used to show which non-linear wave processes are of importance in course of the instability. Finally, we compare these results with similar simulations using Vlasov approach.

Poisson Autoregression

DEFF Research Database (Denmark)

Fokianos, Konstantinos; Rahbek, Anders Christian; Tjøstheim, Dag

2009-01-01

In this article we consider geometric ergodicity and likelihood-based inference for linear and nonlinear Poisson autoregression. In the linear case, the conditional mean is linked linearly to its past values, as well as to the observed values of the Poisson process. This also applies...... to the conditional variance, making possible interpretation as an integer-valued generalized autoregressive conditional heteroscedasticity process. In a nonlinear conditional Poisson model, the conditional mean is a nonlinear function of its past values and past observations. As a particular example, we consider...... an exponential autoregressive Poisson model for time series. Under geometric ergodicity, the maximum likelihood estimators are shown to be asymptotically Gaussian in the linear model. In addition, we provide a consistent estimator of their asymptotic covariance matrix. Our approach to verifying geometric...

Modifying Poisson equation for near-solute dielectric polarization and solvation free energy

Energy Technology Data Exchange (ETDEWEB)

Yang, Pei-Kun, E-mail: peikun@isu.edu.tw

2016-06-15

Highlights: • We modify the Poisson equation. • The dielectric polarization was calculated from the modified Poisson equation. • The solvation free energies of the solutes were calculated from the dielectric polarization. • The calculated solvation free energies were similar to those obtained from MD simulations. - Abstract: The dielectric polarization P is important for calculating the stability of protein conformation and the binding affinity of protein–protein/ligand interactions and for exploring the nonthermal effect of an external electric field on biomolecules. P was decomposed into the product of the electric dipole moment per molecule p; bulk solvent density N{sub bulk}; and relative solvent molecular density g. For a molecular solute, 4πr{sup 2}p(r) oscillates with the distance r to the solute, and g(r) has a large peak in the near-solute region, as observed in molecular dynamics (MD) simulations. Herein, the Poisson equation was modified for computing p based on the modified Gauss’s law of Maxwell’s equations, and the potential of the mean force was used for computing g. For one or two charged atoms in a water cluster, the solvation free energies of the solutes obtained by these equations were similar to those obtained from MD simulations.

A high order solver for the unbounded Poisson equation

DEFF Research Database (Denmark)

Hejlesen, Mads Mølholm; Rasmussen, Johannes Tophøj; Chatelain, Philippe

2012-01-01

This work improves upon Hockney and Eastwood's Fourier-based algorithm for the unbounded Poisson equation to formally achieve arbitrary high order of convergence without any additional computational cost. We assess the methodology on the kinematic relations between the velocity and vorticity fields....

Nonlinear theory of surface-wave--particle interactions in a cylindrical plasma

International Nuclear Information System (INIS)

Dengra, A.; Palop, J.I.F.

1994-01-01

This work is an application of the specular reflection hypothesis to the study of the nonlinear surface-wave--particle interactions in a cylindrical plasma. The model is based on nonlinear resolution of the Vlasov equation by the method of characteristics. The expression obtained for the rate of increase of kinetic energy per electron has permitted us to investigate the temporal behavior of nonlinear collisionless damping for different situations as a function of the critical parameters

Nonlinear wave equations

CERN Document Server

Li, Tatsien

2017-01-01

This book focuses on nonlinear wave equations, which are of considerable significance from both physical and theoretical perspectives. It also presents complete results on the lower bound estimates of lifespan (including the global existence), which are established for classical solutions to the Cauchy problem of nonlinear wave equations with small initial data in all possible space dimensions and with all possible integer powers of nonlinear terms. Further, the book proposes the global iteration method, which offers a unified and straightforward approach for treating these kinds of problems. Purely based on the properties of solut ions to the corresponding linear problems, the method simply applies the contraction mapping principle.

Coefficient Inverse Problem for Poisson's Equation in a Cylinder

NARCIS (Netherlands)

Solov'ev, V. V.

2011-01-01

The inverse problem of determining the coefficient on the right-hand side of Poisson's equation in a cylindrical domain is considered. The Dirichlet boundary value problem is studied. Two types of additional information (overdetermination) can be specified: (i) the trace of the solution to the

Poisson structure of the equations of ideal multispecies fluid electrodynamics

International Nuclear Information System (INIS)

Spencer, R.G.

1984-01-01

The equations of the two- (or multi-) fluid model of plasma physics are recast in Hamiltonian form, following general methods of symplectic geometry. The dynamical variables are the fields of physical interest, but are noncanonical, so that the Poisson bracket in the theory is not the standard one. However, it is a skew-symmetric bilinear form which, from the method of derivation, automatically satisfies the Jacobi identity; therefore, this noncanonical structure has all the essential properties of a canonical Poisson bracket

Numerical solution of continuous-time DSGE models under Poisson uncertainty

DEFF Research Database (Denmark)

Posch, Olaf; Trimborn, Timo

We propose a simple and powerful method for determining the transition process in continuous-time DSGE models under Poisson uncertainty numerically. The idea is to transform the system of stochastic differential equations into a system of functional differential equations of the retarded type. We...... classes of models. We illustrate the algorithm simulating both the stochastic neoclassical growth model and the Lucas model under Poisson uncertainty which is motivated by the Barro-Rietz rare disaster hypothesis. We find that, even for non-linear policy functions, the maximum (absolute) error is very...

Nonlinear screening effect in an ultrarelativistic degenerate electron-positron gas

International Nuclear Information System (INIS)

Tsintsadze, N. L.; Rasheed, A.; Shah, H. A.; Murtaza, G.

2009-01-01

Nonlinear screening process in an ultrarelativistic degenerate electron-positron gas has been investigated by deriving a generalized nonlinear Poisson equation for the electrostatic potential. In the simple one-dimensional case, the nonlinear Poisson equation leads to Debye-like (Coulomb-like) solutions at distances larger (less) than the characteristic length. When the electrostatic energy is larger than the thermal energy, this nonlinear Poisson equation converts into the relativistic Thomas-Fermi equation whose asymptotic solution in three dimensions shows that the potential field goes to zero at infinity much more slowly than the Debye potential. The possibility of the formation of a bound state in electron-positron plasma is also indicated. Further, it is investigated that the strong spatial fluctuations of the potential field may reduce the screening length and that the root mean square of this spatial fluctuating potential goes to zero for large r rather slowly as compared to the case of the Debye potential.

Nonlinear Alfvén Waves in a Vlasov Plasma

DEFF Research Database (Denmark)

Bell, T.F.

1965-01-01

Stationary solutions to the nonlinear Vlasov—Boltzmann equations are considered which represent one-dimensional electromagnetic waves in a hot magnetoplasma. These solutions appear in arbitrary reference frames as circularly polarized, sinusoidal waves of unlimited amplitude, i.e., as nonlinear...... Alfvén waves. Solutions are found implicitly by deriving a set of integral dispersion relations which link the wave characteristics with the particle distribution functions. A physical discussion is given of the way in which the Alfvén waves can trap particles, and it is shown that the presence...

Cellular solutions for the Poisson equation in extended systems

International Nuclear Information System (INIS)

Zhang, X.; Butler, W.H.; MacLaren, J.M.; van Ek, J.

1994-01-01

The Poisson equation for the electrostatic potential in a solid is solved using three different cellular techniques. The relative merits of these different approaches are discussed for two test charge densities for which an analytic solution to the Poisson equation is known. The first approach uses full-cell multiple-scattering theory and results in the famililar structure constant and multipole moment expansion. This solution is shown to be valid everywhere inside the cell, although for points outside the muffin-tin sphere but inside the cell the sums must be performed in the correct order to yield meaningful results. A modification of the multiple-scattering-theory approach yields a second method, a Green-function cellular method, which only requires the solution of a nearest-neighbor linear system of equations. A third approach, a related variational cellular method, is also derived. The variational cellular approach is shown to be the most accurate and reliable, and to have the best convergence in angular momentum of the three methods. Coulomb energies accurate to within 10 -6 hartree are easily achieved with the variational cellular approach, demonstrating the practicality of the approach in electronic structure calculations

Instability of the time splitting scheme for the one-dimensional and relativistic Vlasov-Maxwell system

CERN Document Server

Huot, F; Bertrand, P; Sonnendrücker, E; Coulaud, O

2003-01-01

The Time Splitting Scheme (TSS) has been examined within the context of the one-dimensional (1D) relativistic Vlasov-Maxwell model. In the strongly relativistic regime of the laser-plasma interaction, the TSS cannot be applied to solve the Vlasov equation. We propose a new semi-Lagrangian scheme based on a full 2D advection and study its advantages over the classical Splitting procedure. Details of the underlying integration of the Vlasov equation appear to be important in achieving accurate plasma simulations. Examples are given which are related to the relativistic modulational instability and the self-induced transparency of an ultra-intense electromagnetic pulse in the relativistic regime.

Nonlinear differential equations

Energy Technology Data Exchange (ETDEWEB)

Dresner, L.

1988-01-01

This report is the text of a graduate course on nonlinear differential equations given by the author at the University of Wisconsin-Madison during the summer of 1987. The topics covered are: direction fields of first-order differential equations; the Lie (group) theory of ordinary differential equations; similarity solutions of second-order partial differential equations; maximum principles and differential inequalities; monotone operators and iteration; complementary variational principles; and stability of numerical methods. The report should be of interest to graduate students, faculty, and practicing scientists and engineers. No prior knowledge is required beyond a good working knowledge of the calculus. The emphasis is on practical results. Most of the illustrative examples are taken from the fields of nonlinear diffusion, heat and mass transfer, applied superconductivity, and helium cryogenics.

Nonlinear differential equations

International Nuclear Information System (INIS)

Dresner, L.

1988-01-01

This report is the text of a graduate course on nonlinear differential equations given by the author at the University of Wisconsin-Madison during the summer of 1987. The topics covered are: direction fields of first-order differential equations; the Lie (group) theory of ordinary differential equations; similarity solutions of second-order partial differential equations; maximum principles and differential inequalities; monotone operators and iteration; complementary variational principles; and stability of numerical methods. The report should be of interest to graduate students, faculty, and practicing scientists and engineers. No prior knowledge is required beyond a good working knowledge of the calculus. The emphasis is on practical results. Most of the illustrative examples are taken from the fields of nonlinear diffusion, heat and mass transfer, applied superconductivity, and helium cryogenics

A generalized auxiliary equation method and its application to nonlinear Klein-Gordon and generalized nonlinear Camassa-Holm equations

International Nuclear Information System (INIS)

Yomba, Emmanuel

2008-01-01

With the aid of symbolic computation, a generalized auxiliary equation method is proposed to construct more general exact solutions to two types of NLPDEs. First, we present new family of solutions to a nonlinear Klein-Gordon equation, by using this auxiliary equation method including a new first-order nonlinear ODE with six-degree nonlinear term proposed by Sirendaoreji. Then, we apply an indirect F-function method very close to the F-expansion method to solve the generalized Camassa-Holm equation with fully nonlinear dispersion and fully nonlinear convection C(l,n,p). Taking advantage of the new first-order nonlinear ODE with six degree nonlinear term, this indirect F-function method is used to map the solutions of C(l,n,p) equations to those of that nonlinear ODE. As a result, we can successfully obtain in a unified way, many exact solutions

L2-stability of the Vlasov-Maxwell-Boltzmann system near global Maxwellians

International Nuclear Information System (INIS)

Ha, Seung-Yeal; Xiao, Qinghua; Xiong, Linjie; Zhao, Huijiang

2013-01-01

We present a L 2 -stability theory of the Vlasov-Maxwell-Boltzmann system for the two-species collisional plasma. We show that in a perturbative regime of a global Maxwellian, the L 2 -distance between two strong solutions can be controlled by that between initial data in a Lipschitz manner. Our stability result extends earlier results [Ha, S.-Y. and Xiao, Q.-H., “A revisiting to the L 2 -stability theory of the Boltzmann equation near global Maxwellians,” (submitted) and Ha, S.-Y., Yang, X.-F., and Yun, S.-B., “L 2 stability theory of the Boltzmann equation near a global Maxwellian,” Arch. Ration. Mech. Anal. 197, 657–688 (2010)] on the L 2 -stability of the Boltzmann equation to the Boltzmann equation coupled with self-consistent external forces. As a direct application of our stability result, we show that classical solutions in Duan et al. [“Optimal large-time behavior of the Vlasov-Maxwell-Boltzmann system in the whole space,” Commun. Pure Appl. Math. 24, 1497–1546 (2011)] and Guo [“The Vlasov-Maxwell-Boltzmann system near Maxwellians,” Invent. Math. 153(3), 593–630 (2003)] satisfy a uniform L 2 -stability estimate. This is the first result on the L 2 -stability of the Boltzmann equation coupled with self-consistent field equations in three dimensions

Solving Nonlinear Coupled Differential Equations

Science.gov (United States)

Mitchell, L.; David, J.

1986-01-01

Harmonic balance method developed to obtain approximate steady-state solutions for nonlinear coupled ordinary differential equations. Method usable with transfer matrices commonly used to analyze shaft systems. Solution to nonlinear equation, with periodic forcing function represented as sum of series similar to Fourier series but with form of terms suggested by equation itself.

Analytical solution of the Poisson-Nernst-Planck equations for an electrochemical system close to electroneutrality

International Nuclear Information System (INIS)

Pabst, M.

2014-01-01

Single charge densities and the potential are used to describe models of electrochemical systems. These quantities can be calculated by solving a system of time dependent nonlinear coupled partial differential equations, the Poisson-Nernst-Planck equations. Assuming small deviations from the electroneutral equilibrium, the linearized and decoupled equations are solved for a radial symmetric geometry, which represents the interface between a cell and a sensor device. The densities and the potential are expressed by Fourier-Bessels series. The system considered has a ratio between the Debye-length and its geometric dimension on the order of 10 −4 so the Fourier-Bessel series can be approximated by elementary functions. The time development of the system is characterized by two time constants, τ c and τ g . The constant τ c describes the approach to the stationary state of the total charge and the potential. τ c is several orders of magnitude smaller than the geometry-dependent constant τ g , which is on the order of 10 ms characterizing the transition to the stationary state of the single ion densities

Controlled Nonlinear Stochastic Delay Equations: Part II: Approximations and Pipe-Flow Representations

International Nuclear Information System (INIS)

Kushner, Harold J.

2012-01-01

This is the second part of a work dealing with key issues that have not been addressed in the modeling and numerical optimization of nonlinear stochastic delay systems. We consider new classes of models, such as those with nonlinear functions of several controls (such as products), each with is own delay, controlled random Poisson measure driving terms, admissions control with delayed retrials, and others. Part I was concerned with issues concerning the class of admissible controls and their approximations, since the classical definitions are inadequate for our models. This part is concerned with transportation equation representations and their approximations. Such representations of nonlinear stochastic delay models have been crucial in the development of numerical algorithms with much reduced memory and computational requirements. The representations for the new models are not obvious and are developed. They also provide a template for the adaptation of the Markov chain approximation numerical methods.

Application of the Poisson-Nernst-Planck equations to the migration test

DEFF Research Database (Denmark)

Krabbenhøft, Kristian; Krabbenhøft, Jørgen

2008-01-01

The Poisson-Nernst-Planck (PNP) equations are applied to model the migration test. A detailed analysis of the equations is presented and the effects of a number of common, simplifying assumptions are quantified. In addition, closed-form solutions for the effective chloride diffusivity based...... on the full PNP equations are derived, a number of experiments are analyzed in detail, and a new, truly accelerated migration test is proposed. Finally, we present a finite element procedure for numerical solution of the PNP equations....

A modified Poisson-Boltzmann equation applied to protein adsorption.

Science.gov (United States)

Gama, Marlon de Souza; Santos, Mirella Simões; Lima, Eduardo Rocha de Almeida; Tavares, Frederico Wanderley; Barreto, Amaro Gomes Barreto

2018-01-05

Ion-exchange chromatography has been widely used as a standard process in purification and analysis of protein, based on the electrostatic interaction between the protein and the stationary phase. Through the years, several approaches are used to improve the thermodynamic description of colloidal particle-surface interaction systems, however there are still a lot of gaps specifically when describing the behavior of protein adsorption. Here, we present an improved methodology for predicting the adsorption equilibrium constant by solving the modified Poisson-Boltzmann (PB) equation in bispherical coordinates. By including dispersion interactions between ions and protein, and between ions and surface, the modified PB equation used can describe the Hofmeister effects. We solve the modified Poisson-Boltzmann equation to calculate the protein-surface potential of mean force, treated as spherical colloid-plate system, as a function of process variables. From the potential of mean force, the Henry constants of adsorption, for different proteins and surfaces, are calculated as a function of pH, salt concentration, salt type, and temperature. The obtained Henry constants are compared with experimental data for several isotherms showing excellent agreement. We have also performed a sensitivity analysis to verify the behavior of different kind of salts and the Hofmeister effects. Copyright © 2017 Elsevier B.V. All rights reserved.

The Poisson equation at second order in relativistic cosmology

International Nuclear Information System (INIS)

Hidalgo, J.C.; Christopherson, Adam J.; Malik, Karim A.

2013-01-01

We calculate the relativistic constraint equation which relates the curvature perturbation to the matter density contrast at second order in cosmological perturbation theory. This relativistic ''second order Poisson equation'' is presented in a gauge where the hydrodynamical inhomogeneities coincide with their Newtonian counterparts exactly for a perfect fluid with constant equation of state. We use this constraint to introduce primordial non-Gaussianity in the density contrast in the framework of General Relativity. We then derive expressions that can be used as the initial conditions of N-body codes for structure formation which probe the observable signature of primordial non-Gaussianity in the statistics of the evolved matter density field

Auxiliary equation method for solving nonlinear partial differential equations

International Nuclear Information System (INIS)

Sirendaoreji,; Jiong, Sun

2003-01-01

By using the solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct several kinds of exact travelling wave solutions for some nonlinear partial differential equations. By this method some physically important nonlinear equations are investigated and new exact travelling wave solutions are explicitly obtained with the aid of symbolic computation

Steady state solution of the Poisson-Nernst-Planck equations

International Nuclear Information System (INIS)

Golovnev, A.; Trimper, S.

2010-01-01

The exact steady state solution of the Poisson-Nernst-Planck equations (PNP) is given in terms of Jacobi elliptic functions. A more tractable approximate solution is derived which can be used to compare the results with experimental observations in binary electrolytes. The breakdown of the PNP for high concentration and high applied voltage is discussed.

The multi-order envelope periodic solutions to the nonlinear Schrodinger equation and cubic nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Xiao Yafeng; Xue Haili; Zhang Hongqing

2011-01-01

Based on Jacobi elliptic function and the Lame equation, the perturbation method is applied to get the multi-order envelope periodic solutions of the nonlinear Schrodinger equation and cubic nonlinear Schrodinger equation. These multi-order envelope periodic solutions can degenerate into the different envelope solitary solutions. (authors)

A multiresolution method for solving the Poisson equation using high order regularization

DEFF Research Database (Denmark)

Hejlesen, Mads Mølholm; Walther, Jens Honore

2016-01-01

We present a novel high order multiresolution Poisson solver based on regularized Green's function solutions to obtain exact free-space boundary conditions while using fast Fourier transforms for computational efficiency. Multiresolution is a achieved through local refinement patches and regulari......We present a novel high order multiresolution Poisson solver based on regularized Green's function solutions to obtain exact free-space boundary conditions while using fast Fourier transforms for computational efficiency. Multiresolution is a achieved through local refinement patches...... and regularized Green's functions corresponding to the difference in the spatial resolution between the patches. The full solution is obtained utilizing the linearity of the Poisson equation enabling super-position of solutions. We show that the multiresolution Poisson solver produces convergence rates...

A reliable treatment for nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Khani, F.; Hamedi-Nezhad, S.; Molabahrami, A.

2007-01-01

Exp-function method is used to find a unified solution of nonlinear wave equation. Nonlinear Schroedinger equations with cubic and power law nonlinearity are selected to illustrate the effectiveness and simplicity of the method. It is shown that the Exp-function method, with the help of symbolic computation, provides a powerful mathematical tool for solving nonlinear equation

Lifting particle coordinate changes of magnetic moment type to Vlasov-Maxwell Hamiltonian dynamics

International Nuclear Information System (INIS)

Morrison, P. J.; Vittot, M.; Guillebon, L. de

2013-01-01

Techniques for coordinate changes that depend on both dependent and independent variables are developed and applied to the Maxwell-Vlasov Hamiltonian theory. Particle coordinate changes with a new velocity variable dependent on the magnetic field, with spatial coordinates unchanged, are lifted to the field theoretic level, by transforming the noncanonical Poisson bracket and Hamiltonian structure of the Vlasov-Maxwell dynamics. Several examples are given including magnetic coordinates, where the velocity is decomposed into components parallel and perpendicular to the local magnetic field, and the case of spherical velocity coordinates. An example of the lifting procedure is performed to obtain a simplified version of gyrokinetics, where the magnetic moment is used as a coordinate and the dynamics is reduced by elimination of the electric field energy in the Hamiltonian.

Stability analysis for neutral stochastic differential equation of second order driven by Poisson jumps

Science.gov (United States)

Chadha, Alka; Bora, Swaroop Nandan

2017-11-01

This paper studies the existence, uniqueness, and exponential stability in mean square for the mild solution of neutral second order stochastic partial differential equations with infinite delay and Poisson jumps. By utilizing the Banach fixed point theorem, first the existence and uniqueness of the mild solution of neutral second order stochastic differential equations is established. Then, the mean square exponential stability for the mild solution of the stochastic system with Poisson jumps is obtained with the help of an established integral inequality.

Vlasov modelling of parallel transport in a tokamak scrape-off layer

International Nuclear Information System (INIS)

Manfredi, G; Hirstoaga, S; Devaux, S

2011-01-01

A one-dimensional Vlasov-Poisson model is used to describe the parallel transport in a tokamak scrape-off layer. Thanks to a recently developed 'asymptotic-preserving' numerical scheme, it is possible to lift numerical constraints on the time step and grid spacing, which are no longer limited by, respectively, the electron plasma period and Debye length. The Vlasov approach provides a good velocity-space resolution even in regions of low density. The model is applied to the study of parallel transport during edge-localized modes, with particular emphasis on the particles and energy fluxes on the divertor plates. The numerical results are compared with analytical estimates based on a free-streaming model, with good general agreement. An interesting feature is the observation of an early electron energy flux, due to suprathermal electrons escaping the ions' attraction. In contrast, the long-time evolution is essentially quasi-neutral and dominated by the ion dynamics.

Vlasov modelling of parallel transport in a tokamak scrape-off layer

Energy Technology Data Exchange (ETDEWEB)

Manfredi, G [Institut de Physique et Chimie des Materiaux, CNRS and Universite de Strasbourg, BP 43, F-67034 Strasbourg (France); Hirstoaga, S [INRIA Nancy Grand-Est and Institut de Recherche en Mathematiques Avancees, 7 rue Rene Descartes, F-67084 Strasbourg (France); Devaux, S, E-mail: Giovanni.Manfredi@ipcms.u-strasbg.f, E-mail: hirstoaga@math.unistra.f, E-mail: Stephane.Devaux@ccfe.ac.u [JET-EFDA, Culham Science Centre, Abingdon, OX14 3DB (United Kingdom)

2011-01-15

A one-dimensional Vlasov-Poisson model is used to describe the parallel transport in a tokamak scrape-off layer. Thanks to a recently developed 'asymptotic-preserving' numerical scheme, it is possible to lift numerical constraints on the time step and grid spacing, which are no longer limited by, respectively, the electron plasma period and Debye length. The Vlasov approach provides a good velocity-space resolution even in regions of low density. The model is applied to the study of parallel transport during edge-localized modes, with particular emphasis on the particles and energy fluxes on the divertor plates. The numerical results are compared with analytical estimates based on a free-streaming model, with good general agreement. An interesting feature is the observation of an early electron energy flux, due to suprathermal electrons escaping the ions' attraction. In contrast, the long-time evolution is essentially quasi-neutral and dominated by the ion dynamics.

A study of the disintegration of highly excited nuclei with the Vlasov-Uehling-Uhlenbeck equation

International Nuclear Information System (INIS)

Vinet, L.; Gregoire, C.; Schuck, P.; Remaud, B.; Sebille, F.

1987-01-01

The disintegration of hot and/or compressed nuclei is studied using (i) the Vlasov equation (VE) with imposed spherical symmetry, (ii) the VE in three dimensions (3D) and (iii) the VE in three dimensions supplemented by the Uehling-Uhlenbeck collision term (VUU). We find that case (ii) is slightly more unstable with respect to disintegration compared to case (i) whereas (iii) tends to make nuclei more stable. In all cases the thermal energies (15-20 MeV per nucleon) needed to totally disintegrate a nucleus seem to be higher than those found in static and hydrodynamic calculation. On the contrary, compressional energy very much helps disintegration. Some comments on the introduction of fluctuations and corresponding fragmentation are added. (orig.)

Fourier–Hermite spectral representation for the Vlasov–Poisson system in the weakly collisional limit

KAUST Repository

Parker, Joseph T.; Dellar, Paul J.

2015-01-01

Copyright © Cambridge University Press 2015. We study Landau damping in the 1+1D Vlasov-Poisson system using a Fourier-Hermite spectral representation. We describe the propagation of free energy in Fourier-Hermite phase space using forwards

Lyapunov stability and poisson structure of the thermal TDHF and RPA equations

International Nuclear Information System (INIS)

Balian, R.; Veneroni, M.

1989-01-01

The thermal TDHF equation is analyzed in the Liouville representation of quantum mechanics, where the matrix elements of the single-particle (s.p) density ρ behave as classical dynamical variables. By introducing the Lie--Poisson bracket associated with the unitary group of the s.p. Hilbert space, we show that TDHF has a Hamiltonian, but non-canonical, classical form. Within this Poisson structure, either the s.p. energy or the s.p. grand potential Ω(ρ) act as a Hamilton function. The Lyapunov stability of both the TDHF and RPA equations around a HF state then follows, since the HF approximation for thermal equilibrium is determined by minimizing Ω(ρ). The RPA matrix in the Liouville space is expressed as the product of the Poisson tensor with the HF stability matrix, interpreted as a metric tensor generated by the entropy. This factorization displays the roles of the energy and entropy terms arising from Ω(ρ) in the RPA dynamics, and it helps to construct the RPA modes. Several extensions are considered. copyright 1989 Academic Press, Inc

Lyapunov stability and Poisson structure of the thermal TDHF and RPA equations

International Nuclear Information System (INIS)

Veneroni, M.; Balian, R.

1989-01-01

The thermal TDHF equation is analyzed in the Liouville representation of quantum mechanics, where the matrix elements of the single-particle (s.p.) density ρ behave as classical dynamical variables. By introducing the Lie-Poisson bracket associated with the unitary group of the s.p. Hilbert space, we show that TDHF has a hamiltonian, but non-canonical, classical form. Within this Poisson structure, either the s.p. energy or the s.p. grand potential Ω(ρ) act as a Hamilton function. The Lyapunov stability of both the TDHF and RPA equations around a HF state then follows, since the HF approximation for thermal equilibrium is determined by minimizing Ω(ρ). The RPA matrix in the Liouville space is expressed as the product of the Poisson tensor with the HF stability matrix, interpreted as a metric tensor generated by the entropy. This factorization displays the roles of the energy and entropy terms arising from Ω(ρ) in the RPA dynamics, and it helps to construct the RPA modes. Several extensions are considered

Nambu-Poisson reformulation of the finite dimensional dynamical systems

International Nuclear Information System (INIS)

Baleanu, D.; Makhaldiani, N.

1998-01-01

A system of nonlinear ordinary differential equations which in a particular case reduces to Volterra's system is introduced. We found in two simplest cases the complete sets of the integrals of motion using Nambu-Poisson reformulation of the Hamiltonian dynamics. In these cases we have solved the systems by quadratures

Simple equation method for nonlinear partial differential equations and its applications

Directory of Open Access Journals (Sweden)

Taher A. Nofal

2016-04-01

Full Text Available In this article, we focus on the exact solution of the some nonlinear partial differential equations (NLPDEs such as, Kodomtsev–Petviashvili (KP equation, the (2 + 1-dimensional breaking soliton equation and the modified generalized Vakhnenko equation by using the simple equation method. In the simple equation method the trial condition is the Bernoulli equation or the Riccati equation. It has been shown that the method provides a powerful mathematical tool for solving nonlinear wave equations in mathematical physics and engineering problems.

Numerical simulation of Vlasov equation with parallel tools

International Nuclear Information System (INIS)

Peyroux, J.

2005-11-01

This project aims to make even more powerful the resolution of Vlasov codes through the various parallelization tools (MPI, OpenMP...). A simplified test case served as a base for constructing the parallel codes for obtaining a data-processing skeleton which, thereafter, could be re-used for increasingly complex models (more than four variables of phase space). This will thus make it possible to treat more realistic situations linked, for example, to the injection of ultra short and ultra intense impulses in inertial fusion plasmas, or the study of the instability of trapped ions now taken as being responsible for the generation of turbulence in tokamak plasmas. (author)

Direct solution of the biharmonic equation on rectangular regions and the Poisson equation on irregular regions

International Nuclear Information System (INIS)

Buzbee, B.L.; Dorr, F.W.

1974-01-01

The discrete biharmonic equation on a rectangular region and the discrete Poisson equation on an irregular region can be treated as modifications to matrix problems with very special structure. It is shown how to use the direct method of matrix decomposition to formulate an effective numerical algorithm for these problems. For typical applications the operation count is O(N 3 ) for an N x N grid. Numerical comparisons with other techniques are included. (U.S.)

Direct numerical solution of Poisson's equation in cylindrical (r, z) coordinates

International Nuclear Information System (INIS)

Chao, E.H.; Paul, S.F.; Davidson, R.C.; Fine, K.S.

1997-01-01

A direct solver method is developed for solving Poisson's equation numerically for the electrostatic potential φ(r,z) in a cylindrical region (r wall , 0 wall , z) are specified, and ∂φ/∂z = 0 at the axial boundaries (z = 0, L)

Relations between nonlinear Riccati equations and other equations in fundamental physics

International Nuclear Information System (INIS)

Schuch, Dieter

2014-01-01

Many phenomena in the observable macroscopic world obey nonlinear evolution equations while the microscopic world is governed by quantum mechanics, a fundamental theory that is supposedly linear. In order to combine these two worlds in a common formalism, at least one of them must sacrifice one of its dogmas. Linearizing nonlinear dynamics would destroy the fundamental property of this theory, however, it can be shown that quantum mechanics can be reformulated in terms of nonlinear Riccati equations. In a first step, it will be shown that the information about the dynamics of quantum systems with analytical solutions can not only be obtainable from the time-dependent Schrödinger equation but equally-well from a complex Riccati equation. Comparison with supersymmetric quantum mechanics shows that even additional information can be obtained from the nonlinear formulation. Furthermore, the time-independent Schrödinger equation can also be rewritten as a complex Riccati equation for any potential. Extension of the Riccati formulation to include irreversible dissipative effects is straightforward. Via (real and complex) Riccati equations, other fields of physics can also be treated within the same formalism, e.g., statistical thermodynamics, nonlinear dynamical systems like those obeying a logistic equation as well as wave equations in classical optics, Bose- Einstein condensates and cosmological models. Finally, the link to abstract ''quantizations'' such as the Pythagorean triples and Riccati equations connected with trigonometric and hyperbolic functions will be shown

Numerical simulation of Vlasov equation with parallel tools; Simulations numeriques de l'equation de Vlasov a l'aide d'outils paralleles

Energy Technology Data Exchange (ETDEWEB)

Peyroux, J

2005-11-15

This project aims to make even more powerful the resolution of Vlasov codes through the various parallelization tools (MPI, OpenMP...). A simplified test case served as a base for constructing the parallel codes for obtaining a data-processing skeleton which, thereafter, could be re-used for increasingly complex models (more than four variables of phase space). This will thus make it possible to treat more realistic situations linked, for example, to the injection of ultra short and ultra intense impulses in inertial fusion plasmas, or the study of the instability of trapped ions now taken as being responsible for the generation of turbulence in tokamak plasmas. (author)

Solution of the linearised Vlasov equation for collisionless plasmas evolving in external fields of arbitrary spatial and time dependence: Pt. 2

International Nuclear Information System (INIS)

Skarka, V.; Coveney, P.V.

1990-01-01

We solve perturbatively the linearised Vlasov equation describing inhomogeneous collisionless plasmas evolving in time-dependent external fields. The method employs an explicitly time-dependent formalism and is facilitated by the used of diagrammatic techniques. It leads to a straightforward algorithm for computing the contribution to the solution, order by order in the external field. In the previous paper we provided the solution to first order; higher orders are described in the present paper. (author)

Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation

Science.gov (United States)

Vitanov, Nikolay K.; Dimitrova, Zlatinka I.

2018-03-01

We consider an extension of the methodology of the modified method of simplest equation to the case of use of two simplest equations. The extended methodology is applied for obtaining exact solutions of model nonlinear partial differential equations for deep water waves: the nonlinear Schrödinger equation. It is shown that the methodology works also for other equations of the nonlinear Schrödinger kind.

A New Approach and Solution Technique to Solve Time Fractional Nonlinear Reaction-Diffusion Equations

Directory of Open Access Journals (Sweden)

Inci Cilingir Sungu

2015-01-01

Full Text Available A new application of the hybrid generalized differential transform and finite difference method is proposed by solving time fractional nonlinear reaction-diffusion equations. This method is a combination of the multi-time-stepping temporal generalized differential transform and the spatial finite difference methods. The procedure first converts the time-evolutionary equations into Poisson equations which are then solved using the central difference method. The temporal differential transform method as used in the paper takes care of stability and the finite difference method on the resulting equation results in a system of diagonally dominant linear algebraic equations. The Gauss-Seidel iterative procedure then used to solve the linear system thus has assured convergence. To have optimized convergence rate, numerical experiments were done by using a combination of factors involving multi-time-stepping, spatial step size, and degree of the polynomial fit in time. It is shown that the hybrid technique is reliable, accurate, and easy to apply.

Symmetry and exact solutions of nonlinear spinor equations

International Nuclear Information System (INIS)

Fushchich, W.I.; Zhdanov, R.Z.

1989-01-01

This review is devoted to the application of algebraic-theoretical methods to the problem of constructing exact solutions of the many-dimensional nonlinear systems of partial differential equations for spinor, vector and scalar fields widely used in quantum field theory. Large classes of nonlinear spinor equations invariant under the Poincare group P(1, 3), Weyl group (i.e. Poincare group supplemented by a group of scale transformations), and the conformal group C(1, 3) are described. Ansaetze invariant under the Poincare and the Weyl groups are constructed. Using these we reduce the Poincare-invariant nonlinear Dirac equations to systems of ordinary differential equations and construct large families of exact solutions of the nonlinear Dirac-Heisenberg equation depending on arbitrary parameters and functions. In a similar way we have obtained new families of exact solutions of the nonlinear Maxwell-Dirac and Klein-Gordon-Dirac equations. The obtained solutions can be used for quantization of nonlinear equations. (orig.)

On a new series of integrable nonlinear evolution equations

International Nuclear Information System (INIS)

Ichikawa, Y.H.; Wadati, Miki; Konno, Kimiaki; Shimizu, Tohru.

1980-10-01

Recent results of our research are surveyed in this report. The derivative nonlinear Schroedinger equation for the circular polarized Alfven wave admits the spiky soliton solutions for the plane wave boundary condition. The nonlinear equation for complex amplitude associated with the carrier wave is shown to be a generalized nonlinear Schroedinger equation, having the ordinary cubic nonlinear term and the derivative of cubic nonlinear term. A generalized scheme of the inverse scattering transformation has confirmed that superposition of the A-K-N-S scheme and the K-N scheme for the component equations valids for the generalized nonlinear Schroedinger equation. Then, two types of new integrable nonlinear evolution equation have been derived from our scheme of the inverse scattering transformation. One is the type of nonlinear Schroedinger equation, while the other is the type of Korteweg-de Vries equation. Brief discussions are presented for physical phenomena, which could be accounted by the second type of the new integrable nonlinear evolution equation. Lastly, the stationary solitary wave solutions have been constructed for the integrable nonlinear evolution equation of the second type. These solutions have peculiar structure that they are singular and discrete. It is a new challenge to construct singular potentials by the inverse scattering transformation. (author)

Mathematical study and numerical simulations of bi-kinetic plasma sheaths

International Nuclear Information System (INIS)

Badsi, Mehdi

2016-01-01

This thesis focuses on the construction and the numerical simulation theoretical models of plasmas in interaction with an absorbing wall. These models are based on two species Vlasov-Poisson or Vlasov-Ampere systems in the presence of boundary conditions. The expected stationary solutions must verify the balance of the flux of charges in the orthogonal direction to the wall. This feature is called the ambi-polarity. Through the study of a non linear Poisson equation, we prove the well-posedness of 1d-1v stationary Vlasov-Poisson system, for which we determine incoming particles distributions and a wall potential that induces the ambi-polarity as well as a non negative charge density hold. We also give a quantitative estimates of the thickness of the boundary layer that develops at the wall. These results are illustrated numerically. We prove the linear stability of the electronic stationary solution for a non-stationary Vlasov-Ampere system. Finally, we study a 1d-3v stationary Vlasov-Poisson system in the presence of a constant and parallel to the wall magnetic field. We determine incoming particles distributions and a wall potential so that the ambi-polarity holds. We study a non linear Poisson equation through a non linear functional energy that admits minimizers. We established some bounds on the numerical parameters inside which, our model is relevant and we propose an interpretation of the results. (author)

Semiconductor device simulation by a new method of solving poisson, Laplace and Schrodinger equations

International Nuclear Information System (INIS)

Sharifi, M. J.; Adibi, A.

2000-01-01

In this paper, we have extended and completed our previous work, that was introducing a new method for finite differentiation. We show the applicability of the method for solving a wide variety of equations such as poisson, Laplace and Schrodinger. These equations are fundamental to the most semiconductor device simulators. In a section, we solve the Shordinger equation by this method in several cases including the problem of finding electron concentration profile in the channel of a HEMT. In another section, we solve the Poisson equation by this method, choosing the problem of SBD as an example. Finally we solve the Laplace equation in two dimensions and as an example, we focus on the VED. In this paper, we have shown that, the method can get stable and precise results in solving all of these problems. Also the programs which have been written based on this method become considerably faster, more clear, and more abstract

Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation

Directory of Open Access Journals (Sweden)

Vitanov Nikolay K.

2018-03-01

Full Text Available We consider an extension of the methodology of the modified method of simplest equation to the case of use of two simplest equations. The extended methodology is applied for obtaining exact solutions of model nonlinear partial differential equations for deep water waves: the nonlinear Schrödinger equation. It is shown that the methodology works also for other equations of the nonlinear Schrödinger kind.

Numerical study of non-ideal Vlasov-BGK plasmas

International Nuclear Information System (INIS)

Levchenko, V.D.; Sigov, Y.S.; Premuda, F.

1995-01-01

A relatively simple quasi-classical description of quantum plasmas using as first approximation the Bhatnagar-Gross-Krook (BGK) collision integral, if combined with the modern numerical simulation methods, might be effective tool of a deep study of non-ideal plasma kinetics in a variety of urgent applications as inertial confinement and cold fusion, transport and collective properties of highly condensed plasmas in liquid metals, semi- and superconductors and others. Consider one-dimensional degenerate plasma consisting of thermal electrons and thermal bosons (deuterons) in the vicinity of the equilibrium Fermi- and Bose-type distributions respectively. In the frame of our rough mixed model we solve Vlasov-BGK-Poisson eqs using simplified version of the SUR code

Fractional analysis for nonlinear electrical transmission line and nonlinear Schroedinger equations with incomplete sub-equation

Science.gov (United States)

Fendzi-Donfack, Emmanuel; Nguenang, Jean Pierre; Nana, Laurent

2018-02-01

We use the fractional complex transform with the modified Riemann-Liouville derivative operator to establish the exact and generalized solutions of two fractional partial differential equations. We determine the solutions of fractional nonlinear electrical transmission lines (NETL) and the perturbed nonlinear Schroedinger (NLS) equation with the Kerr law nonlinearity term. The solutions are obtained for the parameters in the range (0<α≤1) of the derivative operator and we found the traditional solutions for the limiting case of α =1. We show that according to the modified Riemann-Liouville derivative, the solutions found can describe physical systems with memory effect, transient effects in electrical systems and nonlinear transmission lines, and other systems such as optical fiber.

Thermoviscous Model Equations in Nonlinear Acoustics

DEFF Research Database (Denmark)

Rasmussen, Anders Rønne

Four nonlinear acoustical wave equations that apply to both perfect gasses and arbitrary fluids with a quadratic equation of state are studied. Shock and rarefaction wave solutions to the equations are studied. In order to assess the accuracy of the wave equations, their solutions are compared...... to solutions of the basic equations from which the wave equations are derived. A straightforward weakly nonlinear equation is the most accurate for shock modeling. A higher order wave equation is the most accurate for modeling of smooth disturbances. Investigations of the linear stability properties...... of solutions to the wave equations, reveal that the solutions may become unstable. Such instabilities are not found in the basic equations. Interacting shocks and standing shocks are investigated....

Relativistic extension of a charge-conservative finite element solver for time-dependent Maxwell-Vlasov equations

Science.gov (United States)

Na, D.-Y.; Moon, H.; Omelchenko, Y. A.; Teixeira, F. L.

2018-01-01

Accurate modeling of relativistic particle motion is essential for physical predictions in many problems involving vacuum electronic devices, particle accelerators, and relativistic plasmas. A local, explicit, and charge-conserving finite-element time-domain (FETD) particle-in-cell (PIC) algorithm for time-dependent (non-relativistic) Maxwell-Vlasov equations on irregular (unstructured) meshes was recently developed by Moon et al. [Comput. Phys. Commun. 194, 43 (2015); IEEE Trans. Plasma Sci. 44, 1353 (2016)]. Here, we extend this FETD-PIC algorithm to the relativistic regime by implementing and comparing three relativistic particle-pushers: (relativistic) Boris, Vay, and Higuera-Cary. We illustrate the application of the proposed relativistic FETD-PIC algorithm for the analysis of particle cyclotron motion at relativistic speeds, harmonic particle oscillation in the Lorentz-boosted frame, and relativistic Bernstein modes in magnetized charge-neutral (pair) plasmas.

Nonlinear diffusion equations

CERN Document Server

Wu Zhuo Qun; Li Hui Lai; Zhao Jun Ning

2001-01-01

Nonlinear diffusion equations, an important class of parabolic equations, come from a variety of diffusion phenomena which appear widely in nature. They are suggested as mathematical models of physical problems in many fields, such as filtration, phase transition, biochemistry and dynamics of biological groups. In many cases, the equations possess degeneracy or singularity. The appearance of degeneracy or singularity makes the study more involved and challenging. Many new ideas and methods have been developed to overcome the special difficulties caused by the degeneracy and singularity, which

Integrable dissipative nonlinear second order differential equations via factorizations and Abel equations

Energy Technology Data Exchange (ETDEWEB)

Mancas, Stefan C. [Department of Mathematics, Embry–Riddle Aeronautical University, Daytona Beach, FL 32114-3900 (United States); Rosu, Haret C., E-mail: hcr@ipicyt.edu.mx [IPICYT, Instituto Potosino de Investigacion Cientifica y Tecnologica, Apdo Postal 3-74 Tangamanga, 78231 San Luis Potosí, SLP (Mexico)

2013-09-02

We emphasize two connections, one well known and another less known, between the dissipative nonlinear second order differential equations and the Abel equations which in their first-kind form have only cubic and quadratic terms. Then, employing an old integrability criterion due to Chiellini, we introduce the corresponding integrable dissipative equations. For illustration, we present the cases of some integrable dissipative Fisher, nonlinear pendulum, and Burgers–Huxley type equations which are obtained in this way and can be of interest in applications. We also show how to obtain Abel solutions directly from the factorization of second order nonlinear equations.

Solution of the Dirichlet Problem for the Poisson's Equation in a Multidimensional Infinite Layer

Directory of Open Access Journals (Sweden)

O. D. Algazin

2015-01-01

Full Text Available The paper considers the multidimensional Poisson equation in the domain bounded by two parallel hyperplanes (in the multidimensional infinite layer. For an n-dimensional half-space method of solving boundary value problems for linear partial differential equations with constant coefficients is a Fourier transform to the variables in the boundary hyperplane. The same method can be used for an infinite layer, as is done in this paper in the case of the Dirichlet problem for the Poisson equation. For strip and infinite layer in three-dimensional space the solutions of this problem are known. And in the three-dimensional case Green's function is written as an infinite series. In this paper, the solution is obtained in the integral form and kernels of integrals are expressed in a finite form in terms of elementary functions and Bessel functions. A recurrence relation between the kernels of integrals for n-dimensional and (n + 2 -dimensional layers was obtained. In particular, is built the Green's function of the Laplace operator for the Dirichlet problem, through which the solution of the problem is recorded. Even in three-dimensional case we obtained new formula compared to the known. It is shown that the kernel of the integral representation of the solution of the Dirichlet problem for a homogeneous Poisson equation (Laplace equation is an approximate identity (δ-shaped system of functions. Therefore, if the boundary values are generalized functions of slow growth, the solution of the Dirichlet problem for the homogeneous equation (Laplace is written as a convolution of kernels with these functions.

Which solutions of the third problem for the Poisson equation are bounded?

Czech Academy of Sciences Publication Activity Database

Medková, Dagmar

-, č. 6 (2004), s. 501-510 ISSN 1085-3375 R&D Projects: GA ČR GA201/00/1515 Institutional research plan: CEZ:AV0Z1019905 Keywords : Poisson equation * Robin problem * boundedness Subject RIV: BA - General Mathematics

POSSOL, 2-D Poisson Equation Solver for Nonuniform Grid

International Nuclear Information System (INIS)

Orvis, W.J.

1988-01-01

1 - Description of program or function: POSSOL is a two-dimensional Poisson equation solver for problems with arbitrary non-uniform gridding in Cartesian coordinates. It is an adaptation of the uniform grid PWSCRT routine developed by Schwarztrauber and Sweet at the National Center for Atmospheric Research (NCAR). 2 - Method of solution: POSSOL will solve the Helmholtz equation on an arbitrary, non-uniform grid on a rectangular domain allowing only one type of boundary condition on any one side. It can also be used to handle more than one type of boundary condition on a side by means of a capacitance matrix technique. There are three types of boundary conditions that can be applied: fixed, derivative, or periodic

Longitudinal traveling waves bifurcating from Vlasov plasma equilibria

International Nuclear Information System (INIS)

Holloway, J.P.

1989-01-01

The kinetic equations governing longitudinal motion along a straight magnetic field in a multi-species collisionless plasma are investigated. A necessary condition for the existence of small amplitude spatially periodic equilibria and traveling waves near a given spatially uniform background equilibrium is derived, and the wavelengths which such solutions must approach as their amplitude decreases to zero are discussed. A sufficient condition for the existence of these small amplitude waves is also established. This is accomplished by studying the nonlinear ODE for the potential which arises when the distribution functions are represented in a BGK form; the arbitrary functions of energy that describe the BGK representation are tested as an infinite dimensional set of parameters in a bifurcation theory for the ODE. The positivity and zero current condition in the wave frame of the BGK distribution functions are maintained. The undamped small amplitude nonlinear waves so constructed can be made to satisfy the Vlasov dispersion relation exactly, but in general they need only satisfy it approximately. Numerical calculations reveal that even a thermal equilibrium electron-proton plasma with equal ion and electron temperatures will support undamped traveling waves with phase speeds greater than 1.3 times the electron velocity; the dispersion relation for this case exhibits both Langmuir and ion-acoustic branches as long wavelength limits, and shows how these branches are in fact connected by short wavelength waves of intermediate frequency. In apparent contradiction to the linear theory of Landau, these exact solutions of the kinetic equations do not damp; this contradiction is explained by observing that the linear theory is, in general, fundamentally incapable of describing undamped traveling waves

Linear superposition solutions to nonlinear wave equations

International Nuclear Information System (INIS)

Liu Yu

2012-01-01

The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K(2,2,1), the Oliver water wave equation, and the k(n, n) equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed

Decomposition of a hierarchy of nonlinear evolution equations

International Nuclear Information System (INIS)

Geng Xianguo

2003-01-01

The generalized Hamiltonian structures for a hierarchy of nonlinear evolution equations are established with the aid of the trace identity. Using the nonlinearization approach, the hierarchy of nonlinear evolution equations is decomposed into a class of new finite-dimensional Hamiltonian systems. The generating function of integrals and their generator are presented, based on which the finite-dimensional Hamiltonian systems are proved to be completely integrable in the Liouville sense. As an application, solutions for the hierarchy of nonlinear evolution equations are reduced to solving the compatible Hamiltonian systems of ordinary differential equations

Poisson-Nernst-Planck Equations for Simulating Biomolecular Diffusion-Reaction Processes II: Size Effects on Ionic Distributions and Diffusion-Reaction Rates

Science.gov (United States)

Lu, Benzhuo; Zhou, Y.C.

2011-01-01

The effects of finite particle size on electrostatics, density profiles, and diffusion have been a long existing topic in the study of ionic solution. The previous size-modified Poisson-Boltzmann and Poisson-Nernst-Planck models are revisited in this article. In contrast to many previous works that can only treat particle species with a single uniform size or two sizes, we generalize the Borukhov model to obtain a size-modified Poisson-Nernst-Planck (SMPNP) model that is able to treat nonuniform particle sizes. The numerical tractability of the model is demonstrated as well. The main contributions of this study are as follows. 1), We show that an (arbitrarily) size-modified PB model is indeed implied by the SMPNP equations under certain boundary/interface conditions, and can be reproduced through numerical solutions of the SMPNP. 2), The size effects in the SMPNP effectively reduce the densities of highly concentrated counterions around the biomolecule. 3), The SMPNP is applied to the diffusion-reaction process for the first time, to our knowledge. In the case of low substrate density near the enzyme reactive site, it is observed that the rate coefficients predicted by SMPNP model are considerably larger than those by the PNP model, suggesting both ions and substrates are subject to finite size effects. 4), An accurate finite element method and a convergent Gummel iteration are developed for the numerical solution of the completely coupled nonlinear system of SMPNP equations. PMID:21575582

A new auxiliary equation and exact travelling wave solutions of nonlinear equations

International Nuclear Information System (INIS)

Sirendaoreji

2006-01-01

A new auxiliary ordinary differential equation and its solutions are used for constructing exact travelling wave solutions of nonlinear partial differential equations in a unified way. The main idea of this method is to take full advantage of the auxiliary equation which has more new exact solutions. More new exact travelling wave solutions are obtained for the quadratic nonlinear Klein-Gordon equation, the combined KdV and mKdV equation, the sine-Gordon equation and the Whitham-Broer-Kaup equations

Quasi-exact solutions of nonlinear differential equations

OpenAIRE

Kudryashov, Nikolay A.; Kochanov, Mark B.

2014-01-01

The concept of quasi-exact solutions of nonlinear differential equations is introduced. Quasi-exact solution expands the idea of exact solution for additional values of parameters of differential equation. These solutions are approximate solutions of nonlinear differential equations but they are close to exact solutions. Quasi-exact solutions of the the Kuramoto--Sivashinsky, the Korteweg--de Vries--Burgers and the Kawahara equations are founded.

One-Dimensional Vlasov-Maxwell Equilibrium for the Force-Free Harris Sheet

International Nuclear Information System (INIS)

Harrison, Michael G.; Neukirch, Thomas

2009-01-01

In this Letter, the first nonlinear force-free Vlasov-Maxwell equilibrium is presented. One component of the equilibrium magnetic field has the same spatial structure as the Harris sheet, but whereas the Harris sheet is kept in force balance by pressure gradients, in the force-free solution presented here force balance is maintained by magnetic shear. Magnetic pressure, plasma pressure and plasma density are constant. The method used to find the equilibrium is based on the analogy of the one-dimensional Vlasov-Maxwell equilibrium problem to the motion of a pseudoparticle in a two-dimensional conservative potential. The force-free solution can be generalized to a complete family of equilibria that describe the transition between the purely pressure-balanced Harris sheet to the force-free Harris sheet

Generalized solutions of nonlinear partial differential equations

CERN Document Server

Rosinger, EE

1987-01-01

During the last few years, several fairly systematic nonlinear theories of generalized solutions of rather arbitrary nonlinear partial differential equations have emerged. The aim of this volume is to offer the reader a sufficiently detailed introduction to two of these recent nonlinear theories which have so far contributed most to the study of generalized solutions of nonlinear partial differential equations, bringing the reader to the level of ongoing research.The essence of the two nonlinear theories presented in this volume is the observation that much of the mathematics concernin

Regionally Implicit Discontinuous Galerkin Methods for Solving the Relativistic Vlasov-Maxwell System Submitted to Iowa State University

Science.gov (United States)

Guthrey, Pierson Tyler

) argument requires. The maximum stable time-step scales inversely with the highest degree in the DG polynomial approximation space and becomes progressively smaller with each added spatial dimension. In this work, we overcome this difficulty by introducing a novel time-stepping strategy: the regionally-implicit discontinuous Galerkin (RIDG) method. The RIDG is method is based on an extension of the Lax-Wendroff DG (LxW-DG) method, which previously had been shown to be equivalent (for linear constant coefficient problems) to a predictor-corrector approach, where the prediction is computed by a space-time DG method (STDG). The corrector is an explicit method that uses the space-time reconstructed solution from the predictor step. In this work, we modify the predictor to include not just local information, but also neighboring information. With this modification, we show that the stability is greatly enhanced; we show that we can remove the polynomial degree dependence of the maximum time-step and show vastly improved time-steps in multiple spatial dimensions. Upon the development of the general RIDG method, we apply it to the non-relativistic 1D1V Vlasov-Poisson equations and the relativistic 1D2V Vlasov-Maxwell equations. For each we validate the high-order method on several test cases. In the final test case, we demonstrate the ability of the method to simulate the acceleration of electrons to relativistic speeds in a simplified test case.

Hamiltonian structures of some non-linear evolution equations

International Nuclear Information System (INIS)

Tu, G.Z.

1983-06-01

The Hamiltonian structure of the O(2,1) non-linear sigma model, generalized AKNS equations, are discussed. By reducing the O(2,1) non-linear sigma model to its Hamiltonian form some new conservation laws are derived. A new hierarchy of non-linear evolution equations is proposed and shown to be generalized Hamiltonian equations with an infinite number of conservation laws. (author)

Nonlinear von Neumann equations for quantum dissipative systems

International Nuclear Information System (INIS)

Messer, J.; Baumgartner, B.

1978-01-01

For pure states nonlinear Schroedinger equations, the so-called Schroedinger-Langevin equations are well-known to model quantum dissipative systems of the Langevin type. For mixtures it is shown that these wave equations do not extend to master equations, but to corresponding nonlinear von Neumann equations. Solutions for the damped harmonic oscillator are discussed. (Auth.)

Nonlinear von Neumann equations for quantum dissipative systems

International Nuclear Information System (INIS)

Messer, J.; Baumgartner, B.

For pure states nonlinear Schroedinger equations, the so-called Schroedinger-Langevin equations are well-known to model quantum dissipative systems of the Langevin type. For mixtures it is shown that these wave equations do not extend to master equations, but to corresponding nonlinear von Neumann equations. Solutions for the damped harmonic oscillator are discussed. (Author)

Multi-diffusive nonlinear Fokker–Planck equation

International Nuclear Information System (INIS)

Ribeiro, Mauricio S; Casas, Gabriela A; Nobre, Fernando D

2017-01-01

Nonlinear Fokker–Planck equations, characterized by more than one diffusion term, have appeared recently in literature. Here, it is shown that these equations may be derived either from approximations in a master equation, or from a Langevin-type approach. An H-theorem is proven, relating these Fokker–Planck equations to an entropy composed by a sum of contributions, each of them associated with a given diffusion term. Moreover, the stationary state of the Fokker–Planck equation is shown to coincide with the equilibrium state, obtained by extremization of the entropy, in the sense that both procedures yield precisely the same equation. Due to the nonlinear character of this equation, the equilibrium probability may be obtained, in most cases, only by means of numerical approaches. Some examples are worked out, where the equilibrium probability distribution is computed for nonlinear Fokker–Planck equations presenting two diffusion terms, corresponding to an entropy characterized by a sum of two contributions. It is shown that the resulting equilibrium distribution, in general, presents a form that differs from a sum of the equilibrium distributions that maximizes each entropic contribution separately, although in some cases one may construct such a linear combination as a good approximation for the equilibrium distribution. (paper)

Conservation of energy and momentum in nonrelativistic plasmas

International Nuclear Information System (INIS)

Sugama, H.; Watanabe, T.-H.; Nunami, M.

2013-01-01

Conservation laws of energy and momentum for nonrelativistic plasmas are derived from applying Noether's theorem to the action integral for the Vlasov-Poisson-Ampère system [Sugama, Phys. Plasmas 7, 466 (2000)]. The symmetric pressure tensor is obtained from modifying the asymmetric canonical pressure tensor with using the rotational symmetry of the action integral. Differences between the resultant conservation laws and those for the Vlasov-Maxwell system including the Maxwell displacement current are clarified. These results provide a useful basis for gyrokinetic conservation laws because gyrokinetic equations are derived as an approximation of the Vlasov-Poisson-Ampère system.

Exact Solutions for Nonlinear Differential Difference Equations in Mathematical Physics

Directory of Open Access Journals (Sweden)

Khaled A. Gepreel

2013-01-01

Full Text Available We modified the truncated expansion method to construct the exact solutions for some nonlinear differential difference equations in mathematical physics via the general lattice equation, the discrete nonlinear Schrodinger with a saturable nonlinearity, the quintic discrete nonlinear Schrodinger equation, and the relativistic Toda lattice system. Also, we put a rational solitary wave function method to find the rational solitary wave solutions for some nonlinear differential difference equations. The proposed methods are more effective and powerful to obtain the exact solutions for nonlinear difference differential equations.

Exact solutions of a nonpolynomially nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Parwani, R.; Tan, H.S.

2007-01-01

A nonlinear generalisation of Schrodinger's equation had previously been obtained using information-theoretic arguments. The nonlinearities in that equation were of a nonpolynomial form, equivalent to the occurrence of higher-derivative nonlinear terms at all orders. Here we construct some exact solutions to that equation in 1+1 dimensions. On the half-line, the solutions resemble (exponentially damped) Bloch waves even though no external periodic potential is included. The solutions are nonperturbative as they do not reduce to solutions of the linear theory in the limit that the nonlinearity parameter vanishes. An intriguing feature of the solutions is their infinite degeneracy: for a given energy, there exists a very large arbitrariness in the normalisable wavefunctions. We also consider solutions to a q-deformed version of the nonlinear equation and discuss a natural discretisation implied by the nonpolynomiality. Finally, we contrast the properties of our solutions with other solutions of nonlinear Schrodinger equations in the literature and suggest some possible applications of our results in the domains of low-energy and high-energy physics

A Poisson equation formulation for pressure calculations in penalty finite element models for viscous incompressible flows

Science.gov (United States)

Sohn, J. L.; Heinrich, J. C.

1990-01-01

The calculation of pressures when the penalty-function approximation is used in finite-element solutions of laminar incompressible flows is addressed. A Poisson equation for the pressure is formulated that involves third derivatives of the velocity field. The second derivatives appearing in the weak formulation of the Poisson equation are calculated from the C0 velocity approximation using a least-squares method. The present scheme is shown to be efficient, free of spurious oscillations, and accurate. Examples of applications are given and compared with results obtained using mixed formulations.

Solitonlike solutions of the generalized discrete nonlinear Schrödinger equation

DEFF Research Database (Denmark)

Rasmussen, Kim; Henning, D.; Gabriel, H.

1996-01-01

We investigate the solution properties oi. a generalized discrete nonlinear Schrodinger equation describing a nonlinear lattice chain. The generalized equation interpolates between the integrable discrete Ablowitz-Ladik equation and the nonintegrable discrete Schrodinger equation. Special interes...... nonlinear Schrodinger equation. In this way eve are able to construct coherent solitonlike structures of profile determined by the map parameters.......We investigate the solution properties oi. a generalized discrete nonlinear Schrodinger equation describing a nonlinear lattice chain. The generalized equation interpolates between the integrable discrete Ablowitz-Ladik equation and the nonintegrable discrete Schrodinger equation. Special interest...

Controlled Nonlinear Stochastic Delay Equations: Part I: Modeling and Approximations

International Nuclear Information System (INIS)

Kushner, Harold J.

2012-01-01

This two-part paper deals with “foundational” issues that have not been previously considered in the modeling and numerical optimization of nonlinear stochastic delay systems. There are new classes of models, such as those with nonlinear functions of several controls (such as products), each with is own delay, controlled random Poisson measure driving terms, admissions control with delayed retrials, and others. There are two basic and interconnected themes for these models. The first, dealt with in this part, concerns the definition of admissible control. The classical definition of an admissible control as a nonanticipative relaxed control is inadequate for these models and needs to be extended. This is needed for the convergence proofs of numerical approximations for optimal controls as well as to have a well-defined model. It is shown that the new classes of admissible controls do not enlarge the range of the value functions, is closed (together with the associated paths) under weak convergence, and is approximatable by ordinary controls. The second theme, dealt with in Part II, concerns transportation equation representations, and their role in the development of numerical algorithms with much reduced memory and computational requirements.

Particular solutions of generalized Euler-Poisson-Darboux equation

Directory of Open Access Journals (Sweden)

Rakhila B. Seilkhanova

2015-01-01

Full Text Available In this article we consider the generalized Euler-Poisson-Darboux equation $$ {u}_{tt}+\\frac{2\\gamma }{t}{{u}_{t}}={u}_{xx}+{u}_{yy} +\\frac{2\\alpha }{x}{{u}_{x}}+\\frac{2\\beta }{y}{{u}_y},\\quad x>0,\\;y>0,\\;t>0. $$ We construct particular solutions in an explicit form expressed by the Lauricella hypergeometric function of three variables. Properties of each constructed solutions have been investigated in sections of surfaces of the characteristic cone. Precisely, we prove that found solutions have singularity $1/r$ at $r\\to 0$, where ${{r}^2}={{( x-{{x}_0}}^2}+{{( y-{{y}_0}}^2}-{{( t-{{t}_0}}^2}$.

Nonlinear hydrodynamic equations for superfluid helium in aerogel

International Nuclear Information System (INIS)

Brusov, Peter N.; Brusov, Paul P.

2003-01-01

Aerogel in superfluids is studied very intensively during last decade. The importance of these systems is connected to the fact that this allows to investigate the influence of impurities on superfluidity. We have derived for the first time nonlinear hydrodynamic equations for superfluid helium in aerogel. These equations are generalization of McKenna et al. equations for nonlinear hydrodynamics case and could be used to study sound propagation phenomena in aerogel-superfluid system, in particular--to study sound conversion phenomena. We have obtained two alternative sets of equations, one of which is a generalization of a traditional set of nonlinear hydrodynamics equations for the case of an aerogel-superfluid system and, the other one represents a la Putterman equations (equation for v→ s is replaced by equation for A→=((ρ n )/(ρσ))w→, where w→=v→ n -v→ s )

Team behaviour analysis in sports using the poisson equation

OpenAIRE

Direkoglu, Cem; O'Connor, Noel E.

2012-01-01

We propose a novel physics-based model for analysing team play- ers’ positions and movements on a sports playing field. The goal is to detect for each frame the region with the highest population of a given team’s players and the region towards which the team is moving as they press for territorial advancement, termed the region of intent. Given the positions of team players from a plan view of the playing field at any given time, we solve a particular Poisson equation to generate a smooth di...

Nonlinear wave equation in frequency domain: accurate modeling of ultrafast interaction in anisotropic nonlinear media

DEFF Research Database (Denmark)

Guo, Hairun; Zeng, Xianglong; Zhou, Binbin

2013-01-01

We interpret the purely spectral forward Maxwell equation with up to third-order induced polarizations for pulse propagation and interactions in quadratic nonlinear crystals. The interpreted equation, also named the nonlinear wave equation in the frequency domain, includes quadratic and cubic...... nonlinearities, delayed Raman effects, and anisotropic nonlinearities. The full potential of this wave equation is demonstrated by investigating simulations of solitons generated in the process of ultrafast cascaded second-harmonic generation. We show that a balance in the soliton delay can be achieved due...

Linear differential equations to solve nonlinear mechanical problems: A novel approach

OpenAIRE

Nair, C. Radhakrishnan

2004-01-01

Often a non-linear mechanical problem is formulated as a non-linear differential equation. A new method is introduced to find out new solutions of non-linear differential equations if one of the solutions of a given non-linear differential equation is known. Using the known solution of the non-linear differential equation, linear differential equations are set up. The solutions of these linear differential equations are found using standard techniques. Then the solutions of the linear differe...

A new sine-Gordon equation expansion algorithm to investigate some special nonlinear differential equations

International Nuclear Information System (INIS)

Yan Zhenya

2005-01-01

A new transformation method is developed using the general sine-Gordon travelling wave reduction equation and a generalized transformation. With the aid of symbolic computation, this method can be used to seek more types of solutions of nonlinear differential equations, which include not only the known solutions derived by some known methods but new solutions. Here we choose the double sine-Gordon equation, the Magma equation and the generalized Pochhammer-Chree (PC) equation to illustrate the method. As a result, many types of new doubly periodic solutions are obtained. Moreover when using the method to these special nonlinear differential equations, some transformations are firstly needed. The method can be also extended to other nonlinear differential equations

Symplectic and Hamiltonian structures of nonlinear evolution equations

International Nuclear Information System (INIS)

Dorfman, I.Y.

1993-01-01

A Hamiltonian structure on a finite-dimensional manifold can be introduced either by endowing it with a (pre)symplectic structure, or by describing the Poisson bracket with the help of a tensor with two upper indices named the Poisson structure. Under the assumption of nondegeneracy, the Poisson structure is nothing else than the inverse of the symplectic structure. Also in the degenerate case the distinction between the two approaches is almost insignificant, because both presymplectic and Poisson structures split into symplectic structures on leaves of appropriately chosen foliations. Hamiltonian structures that arise in the theory of evolution equations demonstrate something new in this respect: trying to operate in local terms, one is induced to develop both approaches independently. Hamiltonian operators, being the infinite-dimensional counterparts of Poisson structures, were the first to become the subject of investigations. A considerable period of time passed before the papers initiated research in the theory of symplectic operators, being the counterparts of presymplectic structures. In what follows, we focus on the main achievements in this field

Nonlinear Fokker-Planck Equations Fundamentals and Applications

CERN Document Server

Frank, Till Daniel

2005-01-01

Providing an introduction to the theory of nonlinear Fokker-Planck equations, this book discusses fundamental properties of transient and stationary solutions, emphasizing the stability analysis of stationary solutions by means of self-consistency equations, linear stability analysis, and Lyapunov's direct method. Also treated are Langevin equations and correlation functions. Nonlinear Fokker-Planck Equations addresses various phenomena such as phase transitions, multistability of systems, synchronization, anomalous diffusion, cut-off solutions, travelling-wave solutions and the emergence of power law solutions. A nonlinear Fokker-Planck perspective to quantum statistics, generalized thermodynamics, and linear nonequilibrium thermodynamics is given. Theoretical concepts are illustrated where possible by simple examples. The book also reviews several applications in the fields of condensed matter physics, the physics of porous media and liquid crystals, accelerator physics, neurophysics, social sciences, popul...

Discrete Time McKean–Vlasov Control Problem: A Dynamic Programming Approach

Energy Technology Data Exchange (ETDEWEB)

Pham, Huyên, E-mail: pham@math.univ-paris-diderot.fr; Wei, Xiaoli, E-mail: tyswxl@gmail.com [Laboratoire de Probabilités et Modèles Aléatoires, CNRS, UMR 7599, Université Paris Diderot (France)

2016-12-15

We consider the stochastic optimal control problem of nonlinear mean-field systems in discrete time. We reformulate the problem into a deterministic control problem with marginal distribution as controlled state variable, and prove that dynamic programming principle holds in its general form. We apply our method for solving explicitly the mean-variance portfolio selection and the multivariate linear-quadratic McKean–Vlasov control problem.

Discrete Time McKean–Vlasov Control Problem: A Dynamic Programming Approach

International Nuclear Information System (INIS)

Pham, Huyên; Wei, Xiaoli

2016-01-01

We consider the stochastic optimal control problem of nonlinear mean-field systems in discrete time. We reformulate the problem into a deterministic control problem with marginal distribution as controlled state variable, and prove that dynamic programming principle holds in its general form. We apply our method for solving explicitly the mean-variance portfolio selection and the multivariate linear-quadratic McKean–Vlasov control problem.

General solution of Poisson equation in three dimensions for disk-like galaxies

International Nuclear Information System (INIS)

Tong, Y.; Zheng, X.; Peng, O.

1982-01-01

The general solution of the Poisson equation is solved by means of integral transformations for Vertical BarkVertical Barr>>1 provided that the perturbed density of disk-like galaxies distributes along the radial direction according to the Hankel function. This solution can more accurately represent the outer spiral arms of disk-like galaxies

The Cauchy problem for non-linear Klein-Gordon equations

International Nuclear Information System (INIS)

Simon, J.C.H.; Taflin, E.

1993-01-01

We consider in R n+1 , n≥2, the non-linear Klein-Gordon equation. We prove for such an equation that there is neighbourhood of zero in a Hilbert space of initial conditions for which the Cauchy problem has global solutions and on which there is asymptotic completeness. The inverse of the wave operator linearizes the non-linear equation. If, moreover, the equation is manifestly Poincare covariant then the non-linear representation of the Poincare-Lie algebra, associated with the non-linear Klein-Gordon equation is integrated to a non-linear representation of the Poincare group on an invariant neighbourhood of zero in the Hilbert space. This representation is linearized by the inverse of the wave operator. The Hilbert space is, in both cases, the closure of the space of the differentiable vectors for the linear representation of the Poincare group, associated with the Klein-Gordon equation, with respect to a norm defined by the representation of the enveloping algebra. (orig.)

International Conference on Differential Equations and Nonlinear Mechanics

CERN Document Server

2001-01-01

The International Conference on Differential Equations and Nonlinear Mechanics was hosted by the University of Central Florida in Orlando from March 17-19, 1999. One of the conference days was dedicated to Professor V. Lakshmikantham in th honor of his 75 birthday. 50 well established professionals (in differential equations, nonlinear analysis, numerical analysis, and nonlinear mechanics) attended the conference from 13 countries. Twelve of the attendees delivered hour long invited talks and remaining thirty-eight presented invited forty-five minute talks. In each of these talks, the focus was on the recent developments in differential equations and nonlinear mechanics and their applications. This book consists of 29 papers based on the invited lectures, and I believe that it provides a good selection of advanced topics of current interest in differential equations and nonlinear mechanics. I am indebted to the Department of Mathematics, College of Arts and Sciences, Department of Mechanical, Materials and Ae...

Study of the nonlinear three-dimensional Debye screening in plasmas

International Nuclear Information System (INIS)

Lin Chang; Zhao Jinbao; Zhang Xiulian

2000-01-01

The nonlinear three-dimensional Debye screening in plasmas is investigated. New analytical solutions for the three-dimensional Poisson equation have been obtained for the nonlinear Debye potential for the first time. We derive exact analytical expression for the special case of the nonlinear three-dimensional Debye screening in plasmas. (orig.)

Numerical simulation of Vlasov equation with parallel tools; Simulations numeriques de l'equation de Vlasov a l'aide d'outils paralleles

Energy Technology Data Exchange (ETDEWEB)

Peyroux, J

2005-11-15

This project aims to make even more powerful the resolution of Vlasov codes through the various parallelization tools (MPI, OpenMP...). A simplified test case served as a base for constructing the parallel codes for obtaining a data-processing skeleton which, thereafter, could be re-used for increasingly complex models (more than four variables of phase space). This will thus make it possible to treat more realistic situations linked, for example, to the injection of ultra short and ultra intense impulses in inertial fusion plasmas, or the study of the instability of trapped ions now taken as being responsible for the generation of turbulence in tokamak plasmas. (author)

Nonlinear elliptic equations of the second order

CERN Document Server

Han, Qing

2016-01-01

Nonlinear elliptic differential equations are a diverse subject with important applications to the physical and social sciences and engineering. They also arise naturally in geometry. In particular, much of the progress in the area in the twentieth century was driven by geometric applications, from the Bernstein problem to the existence of Kähler-Einstein metrics. This book, designed as a textbook, provides a detailed discussion of the Dirichlet problems for quasilinear and fully nonlinear elliptic differential equations of the second order with an emphasis on mean curvature equations and on Monge-Ampère equations. It gives a user-friendly introduction to the theory of nonlinear elliptic equations with special attention given to basic results and the most important techniques. Rather than presenting the topics in their full generality, the book aims at providing self-contained, clear, and "elementary" proofs for results in important special cases. This book will serve as a valuable resource for graduate stu...

Multilevel and Multi-index Monte Carlo methods for the McKean–Vlasov equation

KAUST Repository

Haji-Ali, Abdul-Lateef

2017-09-12

We address the approximation of functionals depending on a system of particles, described by stochastic differential equations (SDEs), in the mean-field limit when the number of particles approaches infinity. This problem is equivalent to estimating the weak solution of the limiting McKean–Vlasov SDE. To that end, our approach uses systems with finite numbers of particles and a time-stepping scheme. In this case, there are two discretization parameters: the number of time steps and the number of particles. Based on these two parameters, we consider different variants of the Monte Carlo and Multilevel Monte Carlo (MLMC) methods and show that, in the best case, the optimal work complexity of MLMC, to estimate the functional in one typical setting with an error tolerance of $$\\\\mathrm {TOL}$$TOL, is when using the partitioning estimator and the Milstein time-stepping scheme. We also consider a method that uses the recent Multi-index Monte Carlo method and show an improved work complexity in the same typical setting of . Our numerical experiments are carried out on the so-called Kuramoto model, a system of coupled oscillators.

Derivation of Poisson and Nernst-Planck equations in a bath and channel from a molecular model.

Science.gov (United States)

Schuss, Z; Nadler, B; Eisenberg, R S

2001-09-01

Permeation of ions from one electrolytic solution to another, through a protein channel, is a biological process of considerable importance. Permeation occurs on a time scale of micro- to milliseconds, far longer than the femtosecond time scales of atomic motion. Direct simulations of atomic dynamics are not yet possible for such long-time scales; thus, averaging is unavoidable. The question is what and how to average. In this paper, we average a Langevin model of ionic motion in a bulk solution and protein channel. The main result is a coupled system of averaged Poisson and Nernst-Planck equations (CPNP) involving conditional and unconditional charge densities and conditional potentials. The resulting NP equations contain the averaged force on a single ion, which is the sum of two components. The first component is the gradient of a conditional electric potential that is the solution of Poisson's equation with conditional and permanent charge densities and boundary conditions of the applied voltage. The second component is the self-induced force on an ion due to surface charges induced only by that ion at dielectric interfaces. The ion induces surface polarization charge that exerts a significant force on the ion itself, not present in earlier PNP equations. The proposed CPNP system is not complete, however, because the electric potential satisfies Poisson's equation with conditional charge densities, conditioned on the location of an ion, while the NP equations contain unconditional densities. The conditional densities are closely related to the well-studied pair-correlation functions of equilibrium statistical mechanics. We examine a specific closure relation, which on the one hand replaces the conditional charge densities by the unconditional ones in the Poisson equation, and on the other hand replaces the self-induced force in the NP equation by an effective self-induced force. This effective self-induced force is nearly zero in the baths but is approximately

Fuchs indices and the first integrals of nonlinear differential equations

International Nuclear Information System (INIS)

Kudryashov, Nikolai A.

2005-01-01

New method of finding the first integrals of nonlinear differential equations in polynomial form is presented. Basic idea of our approach is to use the scaling of solution of nonlinear differential equation and to find the dimensions of arbitrary constants in the Laurent expansion of the general solution. These dimensions allows us to obtain the scalings of members for the first integrals of nonlinear differential equations. Taking the polynomials with unknown coefficients into account we present the algorithm of finding the first integrals of nonlinear differential equations in the polynomial form. Our method is applied to look for the first integrals of eight nonlinear ordinary differential equations of the fourth order. The general solution of one of the fourth order ordinary differential equations is given

On localization in the discrete nonlinear Schrödinger equation

DEFF Research Database (Denmark)

Bang, O.; Juul Rasmussen, J.; Christiansen, P.L.

1993-01-01

For some values of the grid resolution, depending on the nonlinearity, the discrete nonlinear Schrodinger equation with arbitrary power nonlinearity can be approximated by the corresponding continuum version of the equation. When the discretization becomes too coarse, the discrete equation exhibits...

Exact solutions for the higher-order nonlinear Schoerdinger equation in nonlinear optical fibres

International Nuclear Information System (INIS)

Liu Chunping

2005-01-01

First, by using the generally projective Riccati equation method, many kinds of exact solutions for the higher-order nonlinear Schoerdinger equation in nonlinear optical fibres are obtained in a unified way. Then, some relations among these solutions are revealed

Explosions in Landau Vlasov dynamics

International Nuclear Information System (INIS)

Suraud, E.; Cussol, D.; Gregoire, C.; Boilley, D.; Pi, M.; Schuck, P.; Remaud, B.; Sebille, F.

1988-01-01

A microscopic study of the quasi-fusion/explosion transition is presented in the framework of Landau-Vlasov simulations of intermediate energy heavy-ion collisions (bombarding energies between 10 and 100 MeV/A). A detailed analysis in terms of the Equation of State of the system is performed. In agreement with schematic models we find that the composite nuclear system formed in the collision does explode when it stays long enough in the mechanically unstable region (spinodal region). Quantitative estimates of the explosion threshold are given for central symmetric reactions (Ca+Ca and Ar+Ti). The effect of the nuclear matter compressibility modulus is discussed

Simultaneous multigrid techniques for nonlinear eigenvalue problems: Solutions of the nonlinear Schrödinger-Poisson eigenvalue problem in two and three dimensions

Science.gov (United States)

Costiner, Sorin; Ta'asan, Shlomo

1995-07-01

Algorithms for nonlinear eigenvalue problems (EP's) often require solving self-consistently a large number of EP's. Convergence difficulties may occur if the solution is not sought in an appropriate region, if global constraints have to be satisfied, or if close or equal eigenvalues are present. Multigrid (MG) algorithms for nonlinear problems and for EP's obtained from discretizations of partial differential EP have often been shown to be more efficient than single level algorithms. This paper presents MG techniques and a MG algorithm for nonlinear Schrödinger Poisson EP's. The algorithm overcomes the above mentioned difficulties combining the following techniques: a MG simultaneous treatment of the eigenvectors and nonlinearity, and with the global constrains; MG stable subspace continuation techniques for the treatment of nonlinearity; and a MG projection coupled with backrotations for separation of solutions. These techniques keep the solutions in an appropriate region, where the algorithm converges fast, and reduce the large number of self-consistent iterations to only a few or one MG simultaneous iteration. The MG projection makes it possible to efficiently overcome difficulties related to clusters of close and equal eigenvalues. Computational examples for the nonlinear Schrödinger-Poisson EP in two and three dimensions, presenting special computational difficulties that are due to the nonlinearity and to the equal and closely clustered eigenvalues are demonstrated. For these cases, the algorithm requires O(qN) operations for the calculation of q eigenvectors of size N and for the corresponding eigenvalues. One MG simultaneous cycle per fine level was performed. The total computational cost is equivalent to only a few Gauss-Seidel relaxations per eigenvector. An asymptotic convergence rate of 0.15 per MG cycle is attained.

Nonlinear H-infinity control, Hamiltonian systems and Hamilton-Jacobi equations

CERN Document Server

Aliyu, MDS

2011-01-01

A comprehensive overview of nonlinear Haeu control theory for both continuous-time and discrete-time systems, Nonlinear Haeu-Control, Hamiltonian Systems and Hamilton-Jacobi Equations covers topics as diverse as singular nonlinear Haeu-control, nonlinear Haeu -filtering, mixed H2/ Haeu-nonlinear control and filtering, nonlinear Haeu-almost-disturbance-decoupling, and algorithms for solving the ubiquitous Hamilton-Jacobi-Isaacs equations. The link between the subject and analytical mechanics as well as the theory of partial differential equations is also elegantly summarized in a single chapter

OSCILLATION OF NONLINEAR DELAY DIFFERENCE EQUATIONS

Institute of Scientific and Technical Information of China (English)

无

2006-01-01

This paper deals with the oscillatory properties of a class of nonlinear difference equations with several delays. Sufficient criteria in the form of infinite sum for the equations to be oscillatory are obtained.

High-Order Finite-Difference Solution of the Poisson Equation with Interface Jump Conditions II

Science.gov (United States)

Marques, Alexandre; Nave, Jean-Christophe; Rosales, Rodolfo

2010-11-01

The Poisson equation with jump discontinuities across an interface is of central importance in Computational Fluid Dynamics. In prior work, Marques, Nave, and Rosales have introduced a method to obtain fourth-order accurate solutions for the constant coefficient Poisson problem. Here we present an extension of this method to solve the variable coefficient Poisson problem to fourth-order of accuracy. The extended method is based on local smooth extrapolations of the solution field across the interface. The extrapolation procedure uses a combination of cubic Hermite interpolants and a high-order representation of the interface using the Gradient-Augmented Level-Set technique. This procedure is compatible with the use of standard discretizations for the Laplace operator, and leads to modified linear systems which have the same sparsity pattern as the standard discretizations. As a result, standard Poisson solvers can be used with only minimal modifications. Details of the method and applications will be presented.

Algorithms For Integrating Nonlinear Differential Equations

Science.gov (United States)

Freed, A. D.; Walker, K. P.

1994-01-01

Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.

Jacobi Elliptic Solutions for Nonlinear Differential Difference Equations in Mathematical Physics

Directory of Open Access Journals (Sweden)

Khaled A. Gepreel

2012-01-01

Full Text Available We put a direct new method to construct the rational Jacobi elliptic solutions for nonlinear differential difference equations which may be called the rational Jacobi elliptic functions method. We use the rational Jacobi elliptic function method to construct many new exact solutions for some nonlinear differential difference equations in mathematical physics via the lattice equation and the discrete nonlinear Schrodinger equation with a saturable nonlinearity. The proposed method is more effective and powerful to obtain the exact solutions for nonlinear differential difference equations.

Iterative observer based method for source localization problem for Poisson equation in 3D

KAUST Repository

Majeed, Muhammad Usman; Laleg-Kirati, Taous-Meriem

2017-01-01

A state-observer based method is developed to solve point source localization problem for Poisson equation in a 3D rectangular prism with available boundary data. The technique requires a weighted sum of solutions of multiple boundary data

Exact Solution of a Generalized Nonlinear Schrodinger Equation Dimer

DEFF Research Database (Denmark)

Christiansen, Peter Leth; Maniadis, P.; Tsironis, G.P.

1998-01-01

We present exact solutions for a nonlinear dimer system defined throught a discrete nonlinear Schrodinger equation that contains also an integrable Ablowitz-Ladik term. The solutions are obtained throught a transformation that maps the dimer into a double Sine-Gordon like ordinary nonlinear...... differential equation....

Advances in iterative methods for nonlinear equations

CERN Document Server

Busquier, Sonia

2016-01-01

This book focuses on the approximation of nonlinear equations using iterative methods. Nine contributions are presented on the construction and analysis of these methods, the coverage encompassing convergence, efficiency, robustness, dynamics, and applications. Many problems are stated in the form of nonlinear equations, using mathematical modeling. In particular, a wide range of problems in Applied Mathematics and in Engineering can be solved by finding the solutions to these equations. The book reveals the importance of studying convergence aspects in iterative methods and shows that selection of the most efficient and robust iterative method for a given problem is crucial to guaranteeing a good approximation. A number of sample criteria for selecting the optimal method are presented, including those regarding the order of convergence, the computational cost, and the stability, including the dynamics. This book will appeal to researchers whose field of interest is related to nonlinear problems and equations...

Exact solutions of nonlinear generalizations of the Klein Gordon and Schrodinger equations

International Nuclear Information System (INIS)

Burt, P.B.

1978-01-01

Exact solutions of sine Gordon and multiple sine Gordon equations are constructed in terms of solutions of a linear base equation, the Klein Gordon equation and also in terms of nonlinear base equations where the nonlinearity is polynomial in the dependent variable. Further, exact solutions of nonlinear generalizations of the Schrodinger equation and of additional nonlinear generalizations of the Klein Gordon equation are constructed in terms of solutions of linear base equations. Finally, solutions with spherical symmetry, of nonlinear Klein Gordon equations are given. 14 references

Hybrid (Vlasov-Fluid) simulation of ion-acoustic soliton chain formation and validity of Korteweg de-Vries model

Energy Technology Data Exchange (ETDEWEB)

Aminmansoor, F.; Abbasi, H., E-mail: abbasi@aut.ac.ir [Faculty of Energy Engineering and Physics, Amirkabir University of Technology, P.O. Box 15875-4413, Tehran (Iran, Islamic Republic of)

2015-08-15

The present paper is devoted to simulation of nonlinear disintegration of a localized perturbation into ion-acoustic solitons train in a plasma with hot electrons and cold ions. A Gaussian initial perturbation is used to model the localized perturbation. For this purpose, first, we reduce fluid system of equations to a Korteweg de-Vries equation by the following well-known assumptions. (i) On the ion-acoustic evolution time-scale, the electron velocity distribution function (EVDF) is assumed to be stationary. (ii) The calculation is restricted to small amplitude cases. Next, in order to generalize the model to finite amplitudes cases, the evolution of EVDF is included. To this end, a hybrid code is designed to simulate the case, in which electrons dynamics is governed by Vlasov equation, while cold ions dynamics is, like before, studied by the fluid equations. A comparison between the two models shows that although the fluid model is capable of demonstrating the general features of the process, to have a better insight into the relevant physics resulting from the evolution of EVDF, the use of kinetic treatment is of great importance.

The modified simplest equation method to look for exact solutions of nonlinear partial differential equations

OpenAIRE

Efimova, Olga Yu.

2010-01-01

The modification of simplest equation method to look for exact solutions of nonlinear partial differential equations is presented. Using this method we obtain exact solutions of generalized Korteweg-de Vries equation with cubic source and exact solutions of third-order Kudryashov-Sinelshchikov equation describing nonlinear waves in liquids with gas bubbles.

On the Stochastic Wave Equation with Nonlinear Damping

International Nuclear Information System (INIS)

Kim, Jong Uhn

2008-01-01

We discuss an initial boundary value problem for the stochastic wave equation with nonlinear damping. We establish the existence and uniqueness of a solution. Our method for the existence of pathwise solutions consists of regularization of the equation and data, the Galerkin approximation and an elementary measure-theoretic argument. We also prove the existence of an invariant measure when the equation has pure nonlinear damping

Poisson-Nernst-Planck equations with steric effects - non-convexity and multiple stationary solutions

Science.gov (United States)

Gavish, Nir

2018-04-01

We study the existence and stability of stationary solutions of Poisson-Nernst-Planck equations with steric effects (PNP-steric equations) with two counter-charged species. We show that within a range of parameters, steric effects give rise to multiple solutions of the corresponding stationary equation that are smooth. The PNP-steric equation, however, is found to be ill-posed at the parameter regime where multiple solutions arise. Following these findings, we introduce a novel PNP-Cahn-Hilliard model, show that it is well-posed and that it admits multiple stationary solutions that are smooth and stable. The various branches of stationary solutions and their stability are mapped utilizing bifurcation analysis and numerical continuation methods.

Green functions and Langevin equations for nonlinear diffusion equations: A comment on ‘Markov processes, Hurst exponents, and nonlinear diffusion equations’ by Bassler et al.

Science.gov (United States)

Frank, T. D.

2008-02-01

We discuss two central claims made in the study by Bassler et al. [K.E. Bassler, G.H. Gunaratne, J.L. McCauley, Physica A 369 (2006) 343]. Bassler et al. claimed that Green functions and Langevin equations cannot be defined for nonlinear diffusion equations. In addition, they claimed that nonlinear diffusion equations are linear partial differential equations disguised as nonlinear ones. We review bottom-up and top-down approaches that have been used in the literature to derive Green functions for nonlinear diffusion equations and, in doing so, show that the first claim needs to be revised. We show that the second claim as well needs to be revised. To this end, we point out similarities and differences between non-autonomous linear Fokker-Planck equations and autonomous nonlinear Fokker-Planck equations. In this context, we raise the question whether Bassler et al.’s approach to financial markets is physically plausible because it necessitates the introduction of external traders and causes. Such external entities can easily be eliminated when taking self-organization principles and concepts of nonextensive thermostatistics into account and modeling financial processes by means of nonlinear Fokker-Planck equations.

Stability analysis of a Vlasov-Wave system describing particles interacting with their environment

Science.gov (United States)

De Bièvre, Stephan; Goudon, Thierry; Vavasseur, Arthur

2018-06-01

We study a kinetic equation of the Vlasov-Wave type, which arises in the description of the behavior of a large number of particles interacting weakly with an environment, composed of an infinite collection of local vibrational degrees of freedom, modeled by wave equations. We use variational techniques to establish the existence of large families of stationary states for this system, and analyze their stability.

Generalized nonlinear Proca equation and its free-particle solutions

Energy Technology Data Exchange (ETDEWEB)

Nobre, F.D. [Centro Brasileiro de Pesquisas Fisicas and National Institute of Science and Technology for Complex Systems, Rio de Janeiro, RJ (Brazil); Plastino, A.R. [Universidad Nacional Buenos Aires-Noreoeste, CeBio y Secretaria de Investigacion, Junin (Argentina)

2016-06-15

We introduce a nonlinear extension of Proca's field theory for massive vector (spin 1) bosons. The associated relativistic nonlinear wave equation is related to recently advanced nonlinear extensions of the Schroedinger, Dirac, and Klein-Gordon equations inspired on the non-extensive generalized thermostatistics. This is a theoretical framework that has been applied in recent years to several problems in nuclear and particle physics, gravitational physics, and quantum field theory. The nonlinear Proca equation investigated here has a power-law nonlinearity characterized by a real parameter q (formally corresponding to the Tsallis entropic parameter) in such a way that the standard linear Proca wave equation is recovered in the limit q → 1. We derive the nonlinear Proca equation from a Lagrangian, which, besides the usual vectorial field Ψ{sup μ}(vector x,t), involves an additional field Φ{sup μ}(vector x,t). We obtain exact time-dependent soliton-like solutions for these fields having the form of a q-plane wave, and we show that both field equations lead to the relativistic energy-momentum relation E{sup 2} = p{sup 2}c{sup 2} + m{sup 2}c{sup 4} for all values of q. This suggests that the present nonlinear theory constitutes a new field theoretical representation of particle dynamics. In the limit of massless particles the present q-generalized Proca theory reduces to Maxwell electromagnetism, and the q-plane waves yield localized, transverse solutions of Maxwell equations. Physical consequences and possible applications are discussed. (orig.)

Nonlinear Electrostatic Wave Equations for Magnetized Plasmas

DEFF Research Database (Denmark)

Dysthe, K.B.; Mjølhus, E.; Pécseli, Hans

1984-01-01

The lowest order kinetic effects are included in the equations for nonlinear electrostatic electron waves in a magnetized plasma. The modifications of the authors' previous analysis based on a fluid model are discussed.......The lowest order kinetic effects are included in the equations for nonlinear electrostatic electron waves in a magnetized plasma. The modifications of the authors' previous analysis based on a fluid model are discussed....

The matrix nonlinear Schrodinger equation in dimension 2

DEFF Research Database (Denmark)

Zuhan, L; Pedersen, Michael

2001-01-01

In this paper we study the existence of global solutions to the Cauchy problem for the matrix nonlinear Schrodinger equation (MNLS) in 2 space dimensions. A sharp condition for the global existence is obtained for this equation. This condition is in terms of an exact stationary solution...... of a semilinear elliptic equation. In the scalar case, the MNLS reduces to the well-known cubic nonlinear Schrodinger equation for which existence of solutions has been studied by many authors. (C) 2001 Academic Press....

Global solutions of nonlinear Schrödinger equations

CERN Document Server

Bourgain, J

1999-01-01

This volume presents recent progress in the theory of nonlinear dispersive equations, primarily the nonlinear Schrödinger (NLS) equation. The Cauchy problem for defocusing NLS with critical nonlinearity is discussed. New techniques and results are described on global existence and properties of solutions with large Cauchy data. Current research in harmonic analysis around Strichartz's inequalities and its relevance to nonlinear PDE is presented. Several topics in NLS theory on bounded domains are reviewed. Using the NLS as an example, the book offers comprehensive insight on current research r

Prolongation Structure of Semi-discrete Nonlinear Evolution Equations

International Nuclear Information System (INIS)

Bai Yongqiang; Wu Ke; Zhao Weizhong; Guo Hanying

2007-01-01

Based on noncommutative differential calculus, we present a theory of prolongation structure for semi-discrete nonlinear evolution equations. As an illustrative example, a semi-discrete model of the nonlinear Schroedinger equation is discussed in terms of this theory and the corresponding Lax pairs are also given.

Noncanonical Hamiltonian methods in plasma dynamics

International Nuclear Information System (INIS)

Kaufman, A.N.

1982-01-01

A Hamiltonian approach to plasma dynamics is described. The Poisson bracket of two observables g 1 and g 2 is given by using an antisymmetric tensor J, and must satisfy the Jacobi condition. The J can be obtained by elementary tensor analysis. The evolution in time of an observable g is given in terms of the Poisson bracket and a Hamiltonian H(Z). The guiding-center description of particle motion was presented by Littlejohn. The ponderomotive drift and force, the wave-induced oscillation-center velocity, and the gyrofrequency shift are obtained. The Lie transform yields the wave-induced increment to the gyromomentum. In the coulomb model for a Vlasov system, the dynamical variable is the Vlasov distribution f(z). The Hamiltonian functional and the Poisson bracket are obtained. The coupling of f(z) to the Maxwell field appears in the Poisson bracket. The evolution equation yields the Vlasov-Maxwell system. (Kato, T.)

Blow-Up Time for Nonlinear Heat Equations with Transcendental Nonlinearity

Directory of Open Access Journals (Sweden)

Hee Chul Pak

2012-01-01

Full Text Available A blow-up time for nonlinear heat equations with transcendental nonlinearity is investigated. An upper bound of the first blow-up time is presented. It is pointed out that the upper bound of the first blow-up time depends on the support of the initial datum.

Solitary wave solutions of two-dimensional nonlinear Kadomtsev ...

Indian Academy of Sciences (India)

Aly R Seadawy

2017-09-13

Sep 13, 2017 ... We considered the two-dimensional DASWs in colli- sionless, unmagnetized cold plasma consisting of dust fluid, ions and electrons. The dynamics of DASWs is governed by the normalized fluid equations of nonlin- ear continuity (1), nonlinear motion of system (2) and. (3) and linear Poisson equation (4) as.

Handbook of Nonlinear Partial Differential Equations

CERN Document Server

Polyanin, Andrei D

2011-01-01

New to the Second Edition More than 1,000 pages with over 1,500 new first-, second-, third-, fourth-, and higher-order nonlinear equations with solutions Parabolic, hyperbolic, elliptic, and other systems of equations with solutions Some exact methods and transformations Symbolic and numerical methods for solving nonlinear PDEs with Maple(t), Mathematica(R), and MATLAB(R) Many new illustrative examples and tables A large list of references consisting of over 1,300 sources To accommodate different mathematical backgrounds, the authors avoid wherever possible the use of special terminology. They

Exact solutions to two higher order nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Xu Liping; Zhang Jinliang

2007-01-01

Using the homogeneous balance principle and F-expansion method, the exact solutions to two higher order nonlinear Schroedinger equations which describe the propagation of femtosecond pulses in nonlinear fibres are obtained with the aid of a set of subsidiary higher order ordinary differential equations (sub-equations for short)

A nonlinear bounce kinetic equation for trapped electrons

International Nuclear Information System (INIS)

Gang, F.Y.

1990-03-01

A nonlinear bounce averaged drift kinetic equation for trapped electrons is derived. This equation enables one to compute the nonlinear response of the trapped electron distribution function in terms of the field-line projection of a potential fluctuation left-angle e -inqθ φ n right-angle b . It is useful for both analytical and computational studies of the nonlinear evolution of short wavelength (n much-gt 1) trapped electron mode-driven turbulence. 7 refs

Integrable discretization s of derivative nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Tsuchida, Takayuki

2002-01-01

We propose integrable discretizations of derivative nonlinear Schroedinger (DNLS) equations such as the Kaup-Newell equation, the Chen-Lee-Liu equation and the Gerdjikov-Ivanov equation by constructing Lax pairs. The discrete DNLS systems admit the reduction of complex conjugation between two dependent variables and possess bi-Hamiltonian structure. Through transformations of variables and reductions, we obtain novel integrable discretizations of the nonlinear Schroedinger (NLS), modified KdV (mKdV), mixed NLS, matrix NLS, matrix KdV, matrix mKdV, coupled NLS, coupled Hirota, coupled Sasa-Satsuma and Burgers equations. We also discuss integrable discretizations of the sine-Gordon equation, the massive Thirring model and their generalizations. (author)

Nonlocal and nonlinear electrostatics of a dipolar Coulomb fluid.

Science.gov (United States)

Sahin, Buyukdagli; Ralf, Blossey

2014-07-16

We study a model Coulomb fluid consisting of dipolar solvent molecules of finite extent which generalizes the point-like dipolar Poisson-Boltzmann model (DPB) previously introduced by Coalson and Duncan (1996 J. Phys. Chem. 100 2612) and Abrashkin et al (2007 Phys. Rev. Lett. 99 077801). We formulate a nonlocal Poisson-Boltzmann equation (NLPB) and study both linear and nonlinear dielectric response in this model for the case of a single plane geometry. Our results shed light on the relevance of nonlocal versus nonlinear effects in continuum models of material electrostatics.

Evans functions and bifurcations of nonlinear waves of some nonlinear reaction diffusion equations

Science.gov (United States)

Zhang, Linghai

2017-10-01

The main purposes of this paper are to accomplish the existence, stability, instability and bifurcation of the nonlinear waves of the nonlinear system of reaction diffusion equations ut =uxx + α [ βH (u - θ) - u ] - w, wt = ε (u - γw) and to establish the existence, stability, instability and bifurcation of the nonlinear waves of the nonlinear scalar reaction diffusion equation ut =uxx + α [ βH (u - θ) - u ], under different conditions on the model constants. To establish the bifurcation for the system, we will study the existence and instability of a standing pulse solution if 0 1; the existence and instability of two standing wave fronts if 2 (1 + αγ) θ = αβγ and 0 traveling wave front as well as the existence and instability of a standing pulse solution if 0 traveling wave front as well as the existence and instability of an upside down standing pulse solution if 0 traveling wave back of the nonlinear scalar reaction diffusion equation ut =uxx + α [ βH (u - θ) - u ] -w0, where w0 = α (β - 2 θ) > 0 is a positive constant, if 0 motivation to study the existence, stability, instability and bifurcations of the nonlinear waves is to study the existence and stability/instability of infinitely many fast/slow multiple traveling pulse solutions of the nonlinear system of reaction diffusion equations. The existence and stability of infinitely many fast multiple traveling pulse solutions are of great interests in mathematical neuroscience.

Localization of Point Sources for Poisson Equation using State Observers

KAUST Repository

Majeed, Muhammad Usman

2016-08-09

A method based On iterative observer design is presented to solve point source localization problem for Poisson equation with riven boundary data. The procedure involves solution of multiple boundary estimation sub problems using the available Dirichlet and Neumann data from different parts of the boundary. A weighted sum of these solution profiles of sub-problems localizes point sources inside the domain. Method to compute these weights is also provided. Numerical results are presented using finite differences in a rectangular domain. (C) 2016, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

Localization of Point Sources for Poisson Equation using State Observers

KAUST Repository

Majeed, Muhammad Usman; Laleg-Kirati, Taous-Meriem

2016-01-01

A method based On iterative observer design is presented to solve point source localization problem for Poisson equation with riven boundary data. The procedure involves solution of multiple boundary estimation sub problems using the available Dirichlet and Neumann data from different parts of the boundary. A weighted sum of these solution profiles of sub-problems localizes point sources inside the domain. Method to compute these weights is also provided. Numerical results are presented using finite differences in a rectangular domain. (C) 2016, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

Bäcklund transformation of fractional Riccati equation and its applications to nonlinear fractional partial differential equations

International Nuclear Information System (INIS)

Lu, Bin

2012-01-01

In this Letter, the fractional derivatives in the sense of modified Riemann–Liouville derivative and the Bäcklund transformation of fractional Riccati equation are employed for constructing the exact solutions of nonlinear fractional partial differential equations. The power of this manageable method is presented by applying it to several examples. This approach can also be applied to other nonlinear fractional differential equations. -- Highlights: ► Backlund transformation of fractional Riccati equation is presented. ► A new method for solving nonlinear fractional differential equations is proposed. ► Three important fractional differential equations are solved successfully. ► Some new exact solutions of the fractional differential equations are obtained.

Spurious Solutions Of Nonlinear Differential Equations

Science.gov (United States)

Yee, H. C.; Sweby, P. K.; Griffiths, D. F.

1992-01-01

Report utilizes nonlinear-dynamics approach to investigate possible sources of errors and slow convergence and non-convergence of steady-state numerical solutions when using time-dependent approach for problems containing nonlinear source terms. Emphasizes implications for development of algorithms in CFD and computational sciences in general. Main fundamental conclusion of study is that qualitative features of nonlinear differential equations cannot be adequately represented by finite-difference method and vice versa.

Stochastic foundations of undulatory transport phenomena: generalized Poisson-Kac processes—part III extensions and applications to kinetic theory and transport

Science.gov (United States)

Giona, Massimiliano; Brasiello, Antonio; Crescitelli, Silvestro

2017-08-01

This third part extends the theory of Generalized Poisson-Kac (GPK) processes to nonlinear stochastic models and to a continuum of states. Nonlinearity is treated in two ways: (i) as a dependence of the parameters (intensity of the stochastic velocity, transition rates) of the stochastic perturbation on the state variable, similarly to the case of nonlinear Langevin equations, and (ii) as the dependence of the stochastic microdynamic equations of motion on the statistical description of the process itself (nonlinear Fokker-Planck-Kac models). Several numerical and physical examples illustrate the theory. Gathering nonlinearity and a continuum of states, GPK theory provides a stochastic derivation of the nonlinear Boltzmann equation, furnishing a positive answer to the Kac’s program in kinetic theory. The transition from stochastic microdynamics to transport theory within the framework of the GPK paradigm is also addressed.

Fourth-Order Conservative Vlasov-Maxwell Solver for Cartesian and Cylindrical Phase Space Coordinates

Science.gov (United States)

Vogman, Genia

Plasmas are made up of charged particles whose short-range and long-range interactions give rise to complex behavior that can be difficult to fully characterize experimentally. One of the most complete theoretical descriptions of a plasma is that of kinetic theory, which treats each particle species as a probability distribution function in a six-dimensional position-velocity phase space. Drawing on statistical mechanics, these distribution functions mathematically represent a system of interacting particles without tracking individual ions and electrons. The evolution of the distribution function(s) is governed by the Boltzmann equation coupled to Maxwell's equations, which together describe the dynamics of the plasma and the associated electromagnetic fields. When collisions can be neglected, the Boltzmann equation is reduced to the Vlasov equation. High-fidelity simulation of the rich physics in even a subset of the full six-dimensional phase space calls for low-noise high-accuracy numerical methods. To that end, this dissertation investigates a fourth-order finite-volume discretization of the Vlasov-Maxwell equation system, and addresses some of the fundamental challenges associated with applying these types of computationally intensive enhanced-accuracy numerical methods to phase space simulations. The governing equations of kinetic theory are described in detail, and their conservation-law weak form is derived for Cartesian and cylindrical phase space coordinates. This formulation is well known when it comes to Cartesian geometries, as it is used in finite-volume and finite-element discretizations to guarantee local conservation for numerical solutions. By contrast, the conservation-law weak form of the Vlasov equation in cylindrical phase space coordinates is largely unexplored, and to the author's knowledge has never previously been solved numerically. Thereby the methods described in this dissertation for simulating plasmas in cylindrical phase space

Yang-Mills-Vlasov system in the temporal gauge

International Nuclear Information System (INIS)

Choquet-Bruhat, Y.; Noutchegueme, N.

1991-01-01

We prove a local in time existence theorem of a solution of the Cauchy problem for the Yang-Mills-Vlasov integrodifferential system. Such equations govern the evolution of plasmas, for instance of quarks and gluons (quagmas), where non abelian gauge fields and Yang-Mills charges replace the usual electromagnetic field and electric charge. We work with the temporal gauge and use functional spaces with appropriate weight on the momenta, but no fall off is required in the space direction [fr

MINPACK-1, Subroutine Library for Nonlinear Equation System

International Nuclear Information System (INIS)

Garbow, Burton S.

1984-01-01

1 - Description of problem or function: MINPACK1 is a package of FORTRAN subprograms for the numerical solution of systems of non- linear equations and nonlinear least-squares problems. The individual programs are: Identification/Description: - CHKDER: Check gradients for consistency with functions, - DOGLEG: Determine combination of Gauss-Newton and gradient directions, - DPMPAR: Provide double precision machine parameters, - ENORM: Calculate Euclidean norm of vector, - FDJAC1: Calculate difference approximation to Jacobian (nonlinear equations), - FDJAC2: Calculate difference approximation to Jacobian (least squares), - HYBRD: Solve system of nonlinear equations (approximate Jacobian), - HYBRD1: Easy-to-use driver for HYBRD, - HYBRJ: Solve system of nonlinear equations (analytic Jacobian), - HYBRJ1: Easy-to-use driver for HYBRJ, - LMDER: Solve nonlinear least squares problem (analytic Jacobian), - LMDER1: Easy-to-use driver for LMDER, - LMDIF: Solve nonlinear least squares problem (approximate Jacobian), - LMDIF1: Easy-to-use driver for LMDIF, - LMPAR: Determine Levenberg-Marquardt parameter - LMSTR: Solve nonlinear least squares problem (analytic Jacobian, storage conserving), - LMSTR1: Easy-to-use driver for LMSTR, - QFORM: Accumulate orthogonal matrix from QR factorization QRFAC Compute QR factorization of rectangular matrix, - QRSOLV: Complete solution of least squares problem, - RWUPDT: Update QR factorization after row addition, - R1MPYQ: Apply orthogonal transformations from QR factorization, - R1UPDT: Update QR factorization after rank-1 addition, - SPMPAR: Provide single precision machine parameters. 4. Method of solution - MINPACK1 uses the modified Powell hybrid method and the Levenberg-Marquardt algorithm

Analysis of the Poisson-Nernst-Planck equation in a ball for modeling the Voltage-Current relation in neurobiological microdomains

Science.gov (United States)

Cartailler, J.; Schuss, Z.; Holcman, D.

2017-01-01

The electro-diffusion of ions is often described by the Poisson-Nernst-Planck (PNP) equations, which couple nonlinearly the charge concentration and the electric potential. This model is used, among others, to describe the motion of ions in neuronal micro-compartments. It remains at this time an open question how to determine the relaxation and the steady state distribution of voltage when an initial charge of ions is injected into a domain bounded by an impermeable dielectric membrane. The purpose of this paper is to construct an asymptotic approximation to the solution of the stationary PNP equations in a d-dimensional ball (d = 1 , 2 , 3) in the limit of large total charge. In this geometry the PNP system reduces to the Liouville-Gelfand-Bratú (LGB) equation, with the difference that the boundary condition is Neumann, not Dirichlet, and there is a minus sign in the exponent of the exponential term. The entire boundary is impermeable to ions and the electric field satisfies the compatibility condition of Poisson's equation. These differences replace attraction by repulsion in the LGB equation, thus completely changing the solution. We find that the voltage is maximal in the center and decreases toward the boundary. We also find that the potential drop between the center and the surface increases logarithmically in the total number of charges and not linearly, as in classical capacitance theory. This logarithmic singularity is obtained for d = 3 from an asymptotic argument and cannot be derived from the analysis of the phase portrait. These results are used to derive the relation between the outward current and the voltage in a dendritic spine, which is idealized as a dielectric sphere connected smoothly to the nerve axon by a narrow neck. This is a fundamental microdomain involved in neuronal communication. We compute the escape rate of an ion from the steady density in a ball, which models a neuronal spine head, to a small absorbing window in the sphere. We

Nonlinear elliptic equations and nonassociative algebras

CERN Document Server

Nadirashvili, Nikolai; Vlăduţ, Serge

2014-01-01

This book presents applications of noncommutative and nonassociative algebras to constructing unusual (nonclassical and singular) solutions to fully nonlinear elliptic partial differential equations of second order. The methods described in the book are used to solve a longstanding problem of the existence of truly weak, nonsmooth viscosity solutions. Moreover, the authors provide an almost complete description of homogeneous solutions to fully nonlinear elliptic equations. It is shown that even in the very restricted setting of "Hessian equations", depending only on the eigenvalues of the Hessian, these equations admit homogeneous solutions of all orders compatible with known regularity for viscosity solutions provided the space dimension is five or larger. To the contrary, in dimension four or less the situation is completely different, and our results suggest strongly that there are no nonclassical homogeneous solutions at all in dimensions three and four. Thus this book gives a complete list of dimensions...

Discrete coupled derivative nonlinear Schroedinger equations and their quasi-periodic solutions

International Nuclear Information System (INIS)

Geng Xianguo; Su Ting

2007-01-01

A hierarchy of nonlinear differential-difference equations associated with a discrete isospectral problem is proposed, in which a typical differential-difference equation is a discrete coupled derivative nonlinear Schroedinger equation. With the help of the nonlinearization of the Lax pairs, the hierarchy of nonlinear differential-difference equations is decomposed into a new integrable symplectic map and a class of finite-dimensional integrable Hamiltonian systems. Based on the theory of algebraic curve, the Abel-Jacobi coordinates are introduced to straighten out the corresponding flows, from which quasi-periodic solutions for these differential-difference equations are obtained resorting to the Riemann-theta functions. Moreover, a (2+1)-dimensional discrete coupled derivative nonlinear Schroedinger equation is proposed and its quasi-periodic solutions are derived

Nonlinear quantum fluid equations for a finite temperature Fermi plasma

International Nuclear Information System (INIS)

Eliasson, Bengt; Shukla, Padma K

2008-01-01

Nonlinear quantum electron fluid equations are derived, taking into account the moments of the Wigner equation and by using the Fermi-Dirac equilibrium distribution for electrons with an arbitrary temperature. A simplified formalism with the assumptions of incompressibility of the distribution function is used to close the moments in velocity space. The nonlinear quantum diffraction effects into the fluid equations are incorporated. In the high-temperature limit, we retain the nonlinear fluid equations for a dense hot plasma and in the low-temperature limit, we retain the correct fluid equations for a fully degenerate plasma

Construction of Nodal Bubbling Solutions for the Weighted Sinh-Poisson Equation

Directory of Open Access Journals (Sweden)

Yibin Zhang

2013-01-01

Full Text Available We consider the weighted sinh-Poisson equation in , on , where is a small parameter, , and is a unit ball in . By a constructive way, we prove that for any positive integer , there exists a nodal bubbling solution which concentrates at the origin and the other -points , , such that as , , where and is an odd integer with , or is an even integer. The same techniques lead also to a more general result on general domains.

Positive Solutions for Coupled Nonlinear Fractional Differential Equations

Directory of Open Access Journals (Sweden)

Wenning Liu

2014-01-01

Full Text Available We consider the existence of positive solutions for a coupled system of nonlinear fractional differential equations with integral boundary values. Assume the nonlinear term is superlinear in one equation and sublinear in the other equation. By constructing two cones K1, K2 and computing the fixed point index in product cone K1×K2, we obtain that the system has a pair of positive solutions. It is remarkable that it is established on the Cartesian product of two cones, in which the feature of two equations can be opposite.

Superposition of elliptic functions as solutions for a large number of nonlinear equations

International Nuclear Information System (INIS)

Khare, Avinash; Saxena, Avadh

2014-01-01

For a large number of nonlinear equations, both discrete and continuum, we demonstrate a kind of linear superposition. We show that whenever a nonlinear equation admits solutions in terms of both Jacobi elliptic functions cn(x, m) and dn(x, m) with modulus m, then it also admits solutions in terms of their sum as well as difference. We have checked this in the case of several nonlinear equations such as the nonlinear Schrödinger equation, MKdV, a mixed KdV-MKdV system, a mixed quadratic-cubic nonlinear Schrödinger equation, the Ablowitz-Ladik equation, the saturable nonlinear Schrödinger equation, λϕ 4 , the discrete MKdV as well as for several coupled field equations. Further, for a large number of nonlinear equations, we show that whenever a nonlinear equation admits a periodic solution in terms of dn 2 (x, m), it also admits solutions in terms of dn 2 (x,m)±√(m) cn (x,m) dn (x,m), even though cn(x, m)dn(x, m) is not a solution of these nonlinear equations. Finally, we also obtain superposed solutions of various forms for several coupled nonlinear equations

Exact analytical solutions for nonlinear reaction-diffusion equations

International Nuclear Information System (INIS)

Liu Chunping

2003-01-01

By using a direct method via the computer algebraic system of Mathematica, some exact analytical solutions to a class of nonlinear reaction-diffusion equations are presented in closed form. Subsequently, the hyperbolic function solutions and the triangular function solutions of the coupled nonlinear reaction-diffusion equations are obtained in a unified way

Incompressible spectral-element method: Derivation of equations

Science.gov (United States)

Deanna, Russell G.

1993-01-01

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Induced scattering due to nonlinear Landau and cyclotron damping of electromagnetic and electrostatic waves in a magnetized plasma

International Nuclear Information System (INIS)

Sugaya, Reiji

1989-01-01

General expressions of the matrix elements for nonlinear wave-particle scattering (nonlinear Landau and cyclotron damping) of electromagnetic and electrostatic waves in a homogeneous magnetized plasma are derived from the Vlasov-Maxwell equations. The kinetic wave equations obtained for electromagnetic waves are expressed by four-order tensors in the rotating and cartesian coordinates. No restrictions are imposed on the propagation angle to a uniform magnetic field, the Larmor radius, the frequencies, or the wave numbers. By electrostatic approximation of the dielectric tensor and the matrix elements the kinetic wave equations can be applied to the case in which two scattering waves are electrostatic or they are partially electrostatic. Further, the matrix elements in the limit of parallel or perpendicular propagation to the magnetic field are given. (author)

Analytical exact solution of the non-linear Schroedinger equation

International Nuclear Information System (INIS)

Martins, Alisson Xavier; Rocha Filho, Tarcisio Marciano da

2011-01-01

Full text: In this work we present how to classify and obtain analytical solutions of the Schroedinger equation with a generic non-linearity in 1+1 dimensions. Our approach is based on the determination of Lie symmetry transformation mapping solutions into solutions, and non-classical symmetry transformations, mapping a given solution into itself. From these symmetries it is then possible to reduce the equation to a system of ordinary differential equations which can then be solved using standard methods. The generic non-linearity is handled by considering it as an additional unknown in the determining equations for the symmetry transformations. This results in an over-determined system of non-linear partial differential equations. Its solution can then be determined in some cases by reducing it to the so called involutive (triangular) form, and then solved. This reduction is very tedious and can only performed using a computer algebra system. Once the determining system is solved, we obtain the explicit form for the non-linearity admitting a Lie or non-classical symmetry. The analytical solutions are then derived by solving the reduced ordinary differential equations. The non-linear determining system for the non-classical symmetry transformations and Lie symmetry generators are obtaining using the computer algebra package SADE (symmetry analysis of differential equations), developed at our group. (author)

Progress on a Vlasov Treatment of Coherent Synchrotron Radiation from Arbitrary Planar Orbits

CERN Document Server

Bassi, Gabriele; Warnock, Robert L

2005-01-01

We study the influence of coherent synchrotron radiation (CSR) on particle bunches traveling on arbitrary planar orbits between parallel conducting plates (shielding). The time evolution of the phase space distribution is determined by solving the Vlasov-Maxwell equations in the time domain. This provides lower numerical noise than the macroparticle method, and allows the study of emittance degradation and microbunching in bunch compressors. We calculate the fields excited by the bunch in the lab frame using a formula simpler than that based on retarded potentials.* We have developed an algorithm for solving the Vlasov equation in the beam frame using arc length as the independent variable and our method of local characteristics (discretized Perron-Frobenius operator).We integrate in the interaction picture in the hope that we can adopt a fixed grid. The distribution function will be represented by B-splines, in a scheme preserving positivity and normalization of the distribution. The transformation between l...

Renormalization theory of stationary homogeneous strong turbulence in a collisionless plasma

International Nuclear Information System (INIS)

Zhang, Y.Z.

1984-01-01

A renormalization procedure for the perturbation expansion of the Vlasov-Poisson equation is presented to describe stationary homogeneous turbulence. By using the diagramatic scheme the theory is shown to be renormalizable to any order. The expressions for the renormalized propagator, the renormalized dielectric function, and the intrinsically incoherent source are given. The renormalization leads to a complete separation of the fluctuating distribution function f/sub k/ into two parts, the coherent part, which is proved to represent the dielectric effect of the medium, and the intrinsically incoherent part, which represents the effect of nonlinear source. The turbulent collisional operator in the transport equation is proved equal to GAMMA 0 , the frequency broadening when k = 0

A modified SOR method for the Poisson equation in unsteady free-surface flow calculations.

NARCIS (Netherlands)

Botta, E.F.F.; Ellenbroek, Marcellinus Hermannus Maria

1985-01-01

Convergence difficulties that sometimes occur if the successive overrelaxation (SOR) method is applied to the Poisson equation on a region with irregular free boundaries are analyzed. It is shown that these difficulties are related to the treatment of the free boundaries and caused by the appearance

Nontrivial Periodic Solutions for Nonlinear Second-Order Difference Equations

Directory of Open Access Journals (Sweden)

Tieshan He

2011-01-01

Full Text Available This paper is concerned with the existence of nontrivial periodic solutions and positive periodic solutions to a nonlinear second-order difference equation. Under some conditions concerning the first positive eigenvalue of the linear equation corresponding to the nonlinear second-order equation, we establish the existence results by using the topological degree and fixed point index theories.

Nonlinear streak computation using boundary region equations

Energy Technology Data Exchange (ETDEWEB)

Martin, J A; Martel, C, E-mail: juanangel.martin@upm.es, E-mail: carlos.martel@upm.es [Depto. de Fundamentos Matematicos, E.T.S.I Aeronauticos, Universidad Politecnica de Madrid, Plaza Cardenal Cisneros 3, 28040 Madrid (Spain)

2012-08-01

The boundary region equations (BREs) are applied for the simulation of the nonlinear evolution of a spanwise periodic array of streaks in a flat plate boundary layer. The well-known BRE formulation is obtained from the complete Navier-Stokes equations in the high Reynolds number limit, and provides the correct asymptotic description of three-dimensional boundary layer streaks. In this paper, a fast and robust streamwise marching scheme is introduced to perform their numerical integration. Typical streak computations present in the literature correspond to linear streaks or to small-amplitude nonlinear streaks computed using direct numerical simulation (DNS) or the nonlinear parabolized stability equations (PSEs). We use the BREs to numerically compute high-amplitude streaks, a method which requires much lower computational effort than DNS and does not have the consistency and convergence problems of the PSE. It is found that the flow configuration changes substantially as the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, which end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results. (paper)

Spline Collocation Method for Nonlinear Multi-Term Fractional Differential Equation

OpenAIRE

Choe, Hui-Chol; Kang, Yong-Suk

2013-01-01

We study an approximation method to solve nonlinear multi-term fractional differential equations with initial conditions or boundary conditions. First, we transform the nonlinear multi-term fractional differential equations with initial conditions and boundary conditions to nonlinear fractional integral equations and consider the relations between them. We present a Spline Collocation Method and prove the existence, uniqueness and convergence of approximate solution as well as error estimatio...

Solutions to nonlinear Schrodinger equations for special initial data

Directory of Open Access Journals (Sweden)

Takeshi Wada

2015-11-01

Full Text Available This article concerns the solvability of the nonlinear Schrodinger equation with gauge invariant power nonlinear term in one space dimension. The well-posedness of this equation is known only for $H^s$ with $s\\ge 0$. Under some assumptions on the nonlinearity, this paper shows that this equation is uniquely solvable for special but typical initial data, namely the linear combinations of $\\delta(x$ and p.v. (1/x, which belong to $H^{-1/2-0}$. The proof in this article allows $L^2$-perturbations on the initial data.

Finite difference method and algebraic polynomial interpolation for numerically solving Poisson's equation over arbitrary domains

Directory of Open Access Journals (Sweden)

Tsugio Fukuchi

2014-06-01

Full Text Available The finite difference method (FDM based on Cartesian coordinate systems can be applied to numerical analyses over any complex domain. A complex domain is usually taken to mean that the geometry of an immersed body in a fluid is complex; here, it means simply an analytical domain of arbitrary configuration. In such an approach, we do not need to treat the outer and inner boundaries differently in numerical calculations; both are treated in the same way. Using a method that adopts algebraic polynomial interpolations in the calculation around near-wall elements, all the calculations over irregular domains reduce to those over regular domains. Discretization of the space differential in the FDM is usually derived using the Taylor series expansion; however, if we use the polynomial interpolation systematically, exceptional advantages are gained in deriving high-order differences. In using the polynomial interpolations, we can numerically solve the Poisson equation freely over any complex domain. Only a particular type of partial differential equation, Poisson's equations, is treated; however, the arguments put forward have wider generality in numerical calculations using the FDM.

Relativistic electron beam acceleration by cascading nonlinear Landau damping of electromagnetic waves in a plasma

International Nuclear Information System (INIS)

Sugaya, R.; Ue, A.; Maehara, T.; Sugawa, M.

1996-01-01

Acceleration and heating of a relativistic electron beam by cascading nonlinear Landau damping involving three or four intense electromagnetic waves in a plasma are studied theoretically based on kinetic wave equations and transport equations derived from relativistic Vlasov endash Maxwell equations. Three or four electromagnetic waves excite successively two or three nonresonant beat-wave-driven relativistic electron plasma waves with a phase velocity near the speed of light [v p =c(1-γ -2 p ) 1/2 , γ p =ω/ω pe ]. Three beat waves interact nonlinearly with the electron beam and accelerate it to a highly relativistic energy γ p m e c 2 more effectively than by the usual nonlinear Landau damping of two electromagnetic waves. It is proved that the electron beam can be accelerated to more highly relativistic energy in the plasma whose electron density decreases temporally with an appropriate rate because of the temporal increase of γ p . copyright 1996 American Institute of Physics

An inverse source problem of the Poisson equation with Cauchy data

Directory of Open Access Journals (Sweden)

Ji-Chuan Liu

2017-05-01

Full Text Available In this article, we study an inverse source problem of the Poisson equation with Cauchy data. We want to find iterative algorithms to detect the hidden source within a body from measurements on the boundary. Our goal is to reconstruct the location, the size and the shape of the hidden source. This problem is ill-posed, regularization techniques should be employed to obtain the regularized solution. Numerical examples show that our proposed algorithms are valid and effective.

Numerical Simulations of Self-Focused Pulses Using the Nonlinear Maxwell Equations

Science.gov (United States)

Goorjian, Peter M.; Silberberg, Yaron; Kwak, Dochan (Technical Monitor)

1994-01-01

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations. Abstract of a proposed paper for presentation at the meeting NONLINEAR OPTICS: Materials, Fundamentals, and Applications, Hyatt Regency Waikaloa, Waikaloa, Hawaii, July 24-29, 1994, Cosponsored by IEEE/Lasers and Electro-Optics Society and Optical Society of America

Poisson Stochastic Process and Basic Schauder and Sobolev Estimates in the Theory of Parabolic Equations

Science.gov (United States)

Krylov, N. V.; Priola, E.

2017-09-01

We show, among other things, how knowing Schauder or Sobolev-space estimates for the one-dimensional heat equation allows one to derive their multidimensional analogs for equations with coefficients depending only on the time variable with the same constants as in the case of the one-dimensional heat equation. The method is quite general and is based on using the Poisson stochastic process. It also applies to equations involving non-local operators. It looks like no other methods are available at this time and it is a very challenging problem to find a purely analytical approach to proving such results.

Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations

Energy Technology Data Exchange (ETDEWEB)

Qin, Hong; Liu, Jian; Xiao, Jianyuan; Zhang, Ruili; He, Yang; Wang, Yulei; Sun, Yajuan; Burby, Joshua W.; Ellison, Leland; Zhou, Yao

2015-12-14

Particle-in-cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretising its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root searching or global matrix inversion, enabling applications of the proposed method to very large-scale plasma simulations with many, e.g. 10(9), degrees of freedom. The long-term accuracy and fidelity of the algorithm enables us to numerically confirm Mouhot and Villani's theory and conjecture on nonlinear Landau damping over several orders of magnitude using the PIC method, and to calculate the nonlinear evolution of the reflectivity during the mode conversion process from extraordinary waves to Bernstein waves.

Generalized Kapchinskij-Vladimirskij Distribution and Envelope Equation for High-intensity Beams in a Coupled Transverse Focusing Lattice

International Nuclear Information System (INIS)

Qin, Hong; Chung, Moses; Davidson, Ronald C.

2009-01-01

In an uncoupled lattice, the Kapchinskij-Vladimirskij (KV) distribution function first analyzed in 1959 is the only known exact solution of the nonlinear Vlasov-Maxwell equations for high- intensity beams including self-fields in a self-consistent manner. The KV solution is generalized here to high-intensity beams in a coupled transverse lattice using the recently developed generalized Courant-Snyder invariant for coupled transverse dynamics. This solution projects to a rotating, pulsating elliptical beam in transverse configuration space, determined by the generalized matrix envelope equation.

Universality in an information-theoretic motivated nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Parwani, R; Tabia, G

2007-01-01

Using perturbative methods, we analyse a nonlinear generalization of Schrodinger's equation that had previously been obtained through information-theoretic arguments. We obtain analytical expressions for the leading correction, in terms of the nonlinearity scale, to the energy eigenvalues of the linear Schrodinger equation in the presence of an external potential and observe some generic features. In one space dimension these are (i) for nodeless ground states, the energy shifts are subleading in the nonlinearity parameter compared to the shifts for the excited states; (ii) the shifts for the excited states are due predominantly to contribution from the nodes of the unperturbed wavefunctions, and (iii) the energy shifts for excited states are positive for small values of a regulating parameter and negative at large values, vanishing at a universal critical value that is not manifest in the equation. Some of these features hold true for higher dimensional problems. We also study two exactly solved nonlinear Schrodinger equations so as to contrast our observations. Finally, we comment on the possible significance of our results if the nonlinearity is physically realized

Lectures on nonlinear evolution equations initial value problems

CERN Document Server

Racke, Reinhard

2015-01-01

This book mainly serves as an elementary, self-contained introduction to several important aspects of the theory of global solutions to initial value problems for nonlinear evolution equations. The book employs the classical method of continuation of local solutions with the help of a priori estimates obtained for small data. The existence and uniqueness of small, smooth solutions that are defined for all values of the time parameter are investigated. Moreover, the asymptotic behavior of the solutions is described as time tends to infinity. The methods for nonlinear wave equations are discussed in detail. Other examples include the equations of elasticity, heat equations, the equations of thermoelasticity, Schrödinger equations, Klein-Gordon equations, Maxwell equations and plate equations. To emphasize the importance of studying the conditions under which small data problems offer global solutions, some blow-up results are briefly described. Moreover, the prospects for corresponding initial-boundary value p...

Integrability in Dynamical Systems: Florida Workshop in Nonlinear Astronomy, 3rd, University of Florida, Gainesville, Oct. 1, 2, 1987, Proceedings

International Nuclear Information System (INIS)

Buchler, J.R.; Ipser, J.R.; Williams, C.A.

1988-01-01

Recent advances in theoretical celestial mechanics are examined in reviews and reports. Topics addressed include resonant integrable models of galaxies, new integrable systems, Painleve expansions for integrable and nonintegrable ordinary differential equations, and particle-simulation solutions of the Vlasov equation in general relativity. Consideration is given to repulsive and attractive double-bubble space-times, the integrability of magnetic-confinement systems, Hannay's angle and Berry's phase in the classical adiabatic motion of charged particles, the integrability of the nonlinear wave equations, normalization in the face of integrability, and simplifications toward the integrability of perturbed Keplerian systems

Statistical shape analysis using 3D Poisson equation--A quantitatively validated approach.

Science.gov (United States)

Gao, Yi; Bouix, Sylvain

2016-05-01

Statistical shape analysis has been an important area of research with applications in biology, anatomy, neuroscience, agriculture, paleontology, etc. Unfortunately, the proposed methods are rarely quantitatively evaluated, and as shown in recent studies, when they are evaluated, significant discrepancies exist in their outputs. In this work, we concentrate on the problem of finding the consistent location of deformation between two population of shapes. We propose a new shape analysis algorithm along with a framework to perform a quantitative evaluation of its performance. Specifically, the algorithm constructs a Signed Poisson Map (SPoM) by solving two Poisson equations on the volumetric shapes of arbitrary topology, and statistical analysis is then carried out on the SPoMs. The method is quantitatively evaluated on synthetic shapes and applied on real shape data sets in brain structures. Copyright © 2016 Elsevier B.V. All rights reserved.

Polynomial solutions of nonlinear integral equations

International Nuclear Information System (INIS)

Dominici, Diego

2009-01-01

We analyze the polynomial solutions of a nonlinear integral equation, generalizing the work of Bender and Ben-Naim (2007 J. Phys. A: Math. Theor. 40 F9, 2008 J. Nonlinear Math. Phys. 15 (Suppl. 3) 73). We show that, in some cases, an orthogonal solution exists and we give its general form in terms of kernel polynomials

Polynomial solutions of nonlinear integral equations

Energy Technology Data Exchange (ETDEWEB)

Dominici, Diego [Department of Mathematics, State University of New York at New Paltz, 1 Hawk Dr. Suite 9, New Paltz, NY 12561-2443 (United States)], E-mail: dominicd@newpaltz.edu

2009-05-22

We analyze the polynomial solutions of a nonlinear integral equation, generalizing the work of Bender and Ben-Naim (2007 J. Phys. A: Math. Theor. 40 F9, 2008 J. Nonlinear Math. Phys. 15 (Suppl. 3) 73). We show that, in some cases, an orthogonal solution exists and we give its general form in terms of kernel polynomials.

The energy-momentum tensor for the linearized Maxwell-Vlasov and kinetic guiding center theories

International Nuclear Information System (INIS)

Pfirsch, D.; Morrison, P.J.; Texas Univ., Austin

1990-02-01

A modified Hamilton-Jacobi formalism is introduced as a tool to obtain the energy-momentum and angular-momentum tensors for any kind of nonlinear or linearized Maxwell-collisionless kinetic theories. The emphasis is on linearized theories, for which these tensors are derived for the first time. The kinetic theories treated - which need not be the same for all particle species in a plasma - are the Vlasov and kinetic guiding center theories. The Hamiltonian for the guiding center motion is taken in the form resulting from Dirac's constraint theory for non-standard Lagrangian systems. As an example of the Maxwell-kinetic guiding center theory, the second-order energy for a perturbed homogeneous magnetized plasma is calculated with initially vanishing field perturbations. The expression obtained is compared with the corresponding one of Maxwell-Vlasov theory. (orig.)

The energy-momentum tensor for the linearized Maxwell-Vlasov and kinetic guiding center theories

International Nuclear Information System (INIS)

Pfirsch, D.; Morrison, P.J.

1990-02-01

A modified Hamilton-Jacobi formalism is introduced as a tool to obtain the energy-momentum and angular-momentum tensors for any king of nonlinear or linearized Maxwell-collisionless kinetic theories. The emphasis is on linearized theories, for which these tensors are derived for the first time. The kinetic theories treated --- which need not be the same for all particle species in a plasma --- are the Vlasov and kinetic guiding center theories. The Hamiltonian for the guiding center motion is taken in the form resulting from Dirac's constraint theory for non-standard Lagrangian systems. As an example of the Maxwell-kinetic guiding center theory, the second-order energy for a perturbed homogeneous magnetized plasma is calculated with initially vanishing field perturbations. The expression obtained is compared with the corresponding one of Maxwell-Vlasov theory. 11 refs

Construction of alternative Hamiltonian structures for field equations

Energy Technology Data Exchange (ETDEWEB)

Herrera, Mauricio [Departamento de Fisica, Facultad de Ciencias Fisicas y Matematicas, Universidad de Chile, Santiago (Chile); Hojman, Sergio A. [Departamento de Fisica, Facultad de Ciencias, Universidad de Chile, Santiago (Chile); Facultad de Educacion, Universidad Nacional Andres Bello, Santiago (Chile); Centro de Recursos Educativos Avanzados, CREA, Santiago (Chile)

2001-08-10

We use symmetry vectors of nonlinear field equations to build alternative Hamiltonian structures. We construct such structures even for equations which are usually believed to be non-Hamiltonian such as heat, Burger and potential Burger equations. We improve on a previous version of the approach using recursion operators to increase the rank of the Poisson bracket matrices. Cole-Hopf and Miura-type transformations allow the mapping of these structures from one equation to another. (author)

Steric effects in the dynamics of electrolytes at large applied voltages. II. Modified Poisson-Nernst-Planck equations.

Science.gov (United States)

Kilic, Mustafa Sabri; Bazant, Martin Z; Ajdari, Armand

2007-02-01

In situations involving large potentials or surface charges, the Poisson-Boltzman (PB) equation has shortcomings because it neglects ion-ion interactions and steric effects. This has been widely recognized by the electrochemistry community, leading to the development of various alternative models resulting in different sets "modified PB equations," which have had at least qualitative success in predicting equilibrium ion distributions. On the other hand, the literature is scarce in terms of descriptions of concentration dynamics in these regimes. Here, adapting strategies developed to modify the PB equation, we propose a simple modification of the widely used Poisson-Nernst-Planck (PNP) equations for ionic transport, which at least qualitatively accounts for steric effects. We analyze numerical solutions of these modified PNP equations on the model problem of the charging of a simple electrolyte cell, and compare the outcome to that of the standard PNP equations. Finally, we repeat the asymptotic analysis of Bazant, Thornton, and Ajdari [Phys. Rev. E 70, 021506 (2004)] for this new system of equations to further document the interest and limits of validity of the simpler equivalent electrical circuit models introduced in Part I [Kilic, Bazant, and Ajdari, Phys. Rev. E 75, 021502 (2007)] for such problems.

Entropy and convexity for nonlinear partial differential equations.

Science.gov (United States)

Ball, John M; Chen, Gui-Qiang G

2013-12-28

Partial differential equations are ubiquitous in almost all applications of mathematics, where they provide a natural mathematical description of many phenomena involving change in physical, chemical, biological and social processes. The concept of entropy originated in thermodynamics and statistical physics during the nineteenth century to describe the heat exchanges that occur in the thermal processes in a thermodynamic system, while the original notion of convexity is for sets and functions in mathematics. Since then, entropy and convexity have become two of the most important concepts in mathematics. In particular, nonlinear methods via entropy and convexity have been playing an increasingly important role in the analysis of nonlinear partial differential equations in recent decades. This opening article of the Theme Issue is intended to provide an introduction to entropy, convexity and related nonlinear methods for the analysis of nonlinear partial differential equations. We also provide a brief discussion about the content and contributions of the papers that make up this Theme Issue.

Poisson hierarchy of discrete strings

International Nuclear Information System (INIS)

Ioannidou, Theodora; Niemi, Antti J.

2016-01-01

The Poisson geometry of a discrete string in three dimensional Euclidean space is investigated. For this the Frenet frames are converted into a spinorial representation, the discrete spinor Frenet equation is interpreted in terms of a transfer matrix formalism, and Poisson brackets are introduced in terms of the spinor components. The construction is then generalised, in a self-similar manner, into an infinite hierarchy of Poisson algebras. As an example, the classical Virasoro (Witt) algebra that determines reparametrisation diffeomorphism along a continuous string, is identified as a particular sub-algebra, in the hierarchy of the discrete string Poisson algebra. - Highlights: • Witt (classical Virasoro) algebra is derived in the case of discrete string. • Infinite dimensional hierarchy of Poisson bracket algebras is constructed for discrete strings. • Spinor representation of discrete Frenet equations is developed.

Poisson hierarchy of discrete strings

Energy Technology Data Exchange (ETDEWEB)

Ioannidou, Theodora, E-mail: ti3@auth.gr [Faculty of Civil Engineering, School of Engineering, Aristotle University of Thessaloniki, 54249, Thessaloniki (Greece); Niemi, Antti J., E-mail: Antti.Niemi@physics.uu.se [Department of Physics and Astronomy, Uppsala University, P.O. Box 803, S-75108, Uppsala (Sweden); Laboratoire de Mathematiques et Physique Theorique CNRS UMR 6083, Fédération Denis Poisson, Université de Tours, Parc de Grandmont, F37200, Tours (France); Department of Physics, Beijing Institute of Technology, Haidian District, Beijing 100081 (China)

2016-01-28

The Poisson geometry of a discrete string in three dimensional Euclidean space is investigated. For this the Frenet frames are converted into a spinorial representation, the discrete spinor Frenet equation is interpreted in terms of a transfer matrix formalism, and Poisson brackets are introduced in terms of the spinor components. The construction is then generalised, in a self-similar manner, into an infinite hierarchy of Poisson algebras. As an example, the classical Virasoro (Witt) algebra that determines reparametrisation diffeomorphism along a continuous string, is identified as a particular sub-algebra, in the hierarchy of the discrete string Poisson algebra. - Highlights: • Witt (classical Virasoro) algebra is derived in the case of discrete string. • Infinite dimensional hierarchy of Poisson bracket algebras is constructed for discrete strings. • Spinor representation of discrete Frenet equations is developed.

Analysis of an Nth-order nonlinear differential-delay equation

Science.gov (United States)

Vallée, Réal; Marriott, Christopher

1989-01-01

The problem of a nonlinear dynamical system with delay and an overall response time which is distributed among N individual components is analyzed. Such a system can generally be modeled by an Nth-order nonlinear differential delay equation. A linear-stability analysis as well as a numerical simulation of that equation are performed and a comparison is made with the experimental results. Finally, a parallel is established between the first-order differential equation with delay and the Nth-order differential equation without delay.

Instability and dynamics of two nonlinearly coupled intense laser beams in a quantum plasma

International Nuclear Information System (INIS)

Wang Yunliang; Shukla, P. K.; Eliasson, B.

2013-01-01

We consider nonlinear interactions between two relativistically strong laser beams and a quantum plasma composed of degenerate electron fluids and immobile ions. The collective behavior of degenerate electrons is modeled by quantum hydrodynamic equations composed of the electron continuity, quantum electron momentum (QEM) equation, as well as the Poisson and Maxwell equations. The QEM equation accounts the quantum statistical electron pressure, the quantum electron recoil due to electron tunneling through the quantum Bohm potential, electron-exchange, and electron-correlation effects caused by electron spin, and relativistic ponderomotive forces (RPFs) of two circularly polarized electromagnetic (CPEM) beams. The dynamics of the latter are governed by nonlinear wave equations that include nonlinear currents arising from the relativistic electron mass increase in the CPEM wave fields, as well as from the beating of the electron quiver velocity and electron density variations reinforced by the RPFs of the two CPEM waves. Furthermore, nonlinear electron density variations associated with the driven (by the RPFs) quantum electron plasma oscillations obey a coupled nonlinear Schrödinger and Poisson equations. The nonlinearly coupled equations for our purposes are then used to obtain a general dispersion relation (GDR) for studying the parametric instabilities and the localization of CPEM wave packets in a quantum plasma. Numerical analyses of the GDR reveal that the growth rate of a fastest growing parametrically unstable mode is in agreement with the result that has been deduced from numerical simulations of the governing nonlinear equations. Explicit numerical results for two-dimensional (2D) localized CPEM wave packets at nanoscales are also presented. Possible applications of our investigation to intense laser-solid density compressed plasma experiments are highlighted.

The nonlinear Schrödinger equation singular solutions and optical collapse

CERN Document Server

Fibich, Gadi

2015-01-01

This book is an interdisciplinary introduction to optical collapse of laser beams, which is modelled by singular (blow-up) solutions of the nonlinear Schrödinger equation. With great care and detail, it develops the subject including the mathematical and physical background and the history of the subject. It combines rigorous analysis, asymptotic analysis, informal arguments, numerical simulations, physical modelling, and physical experiments. It repeatedly emphasizes the relations between these approaches, and the intuition behind the results. The Nonlinear Schrödinger Equation will be useful to graduate students and researchers in applied mathematics who are interested in singular solutions of partial differential equations, nonlinear optics and nonlinear waves, and to graduate students and researchers in physics and engineering who are interested in nonlinear optics and Bose-Einstein condensates. It can be used for courses on partial differential equations, nonlinear waves, and nonlinear optics. Gadi Fib...

Relative and Absolute Error Control in a Finite-Difference Method Solution of Poisson's Equation

Science.gov (United States)

Prentice, J. S. C.

2012-01-01

An algorithm for error control (absolute and relative) in the five-point finite-difference method applied to Poisson's equation is described. The algorithm is based on discretization of the domain of the problem by means of three rectilinear grids, each of different resolution. We discuss some hardware limitations associated with the algorithm,…

Separation of variables for the nonlinear wave equation in polar coordinates

International Nuclear Information System (INIS)

Shermenev, Alexander

2004-01-01

Some classical types of nonlinear wave motion in polar coordinates are studied within quadratic approximation. When the nonlinear quadratic terms in the wave equation are arbitrary, the usual perturbation techniques used in polar coordinates leads to overdetermined systems of linear algebraic equations for the unknown coefficients. However, we show that these overdetermined systems are compatible with the special case of the nonlinear shallow water equation and express explicitly the coefficients of the first two harmonics as polynomials of the Bessel functions of radius and of the trigonometric functions of angle. It gives a series of solutions to the nonlinear shallow water equation that are periodic in time and found with the same accuracy as the equation is derived

Nonlinear scalar field equations. Pt. 1

International Nuclear Information System (INIS)

Berestycki, H.; Lions, P.L.

1983-01-01

This paper as well as a subsequent one is concerned with the existence of nontrivial solutions for some semi-linear elliptic equations in Rsup(N). Such problems are motivated in particular by the search for certain kinds of solitary waves (stationary states) in nonlinear equations of the Klein-Gordon or Schroedinger type. (orig./HSI)

Solitary wave for a nonintegrable discrete nonlinear Schrödinger equation in nonlinear optical waveguide arrays

Science.gov (United States)

Ma, Li-Yuan; Ji, Jia-Liang; Xu, Zong-Wei; Zhu, Zuo-Nong

2018-03-01

We study a nonintegrable discrete nonlinear Schrödinger (dNLS) equation with the term of nonlinear nearest-neighbor interaction occurred in nonlinear optical waveguide arrays. By using discrete Fourier transformation, we obtain numerical approximations of stationary and travelling solitary wave solutions of the nonintegrable dNLS equation. The analysis of stability of stationary solitary waves is performed. It is shown that the nonlinear nearest-neighbor interaction term has great influence on the form of solitary wave. The shape of solitary wave is important in the electric field propagating. If we neglect the nonlinear nearest-neighbor interaction term, much important information in the electric field propagating may be missed. Our numerical simulation also demonstrates the difference of chaos phenomenon between the nonintegrable dNLS equation with nonlinear nearest-neighbor interaction and another nonintegrable dNLS equation without the term. Project supported by the National Natural Science Foundation of China (Grant Nos. 11671255 and 11701510), the Ministry of Economy and Competitiveness of Spain (Grant No. MTM2016-80276-P (AEI/FEDER, EU)), and the China Postdoctoral Science Foundation (Grant No. 2017M621964).

A high order multi-resolution solver for the Poisson equation with application to vortex methods

DEFF Research Database (Denmark)

Hejlesen, Mads Mølholm; Spietz, Henrik Juul; Walther, Jens Honore

A high order method is presented for solving the Poisson equation subject to mixed free-space and periodic boundary conditions by using fast Fourier transforms (FFT). The high order convergence is achieved by deriving mollified Green’s functions from a high order regularization function which...

On the invariant measure for the nonlinear Schroedinger equation

International Nuclear Information System (INIS)

Zhidkov, P.R.

1991-01-01

The invariant measure for the nonlinear Schroedinger equation is constructed. In fact, it is assumed that the nonlinearity in the equation is semilinear. The main aim of the paper is the explanation of the Fermi - Past - Ulam phenomenon. Poincare theorem gives the answer to this question. 17 refs

Non-linear wave equations:Mathematical techniques

International Nuclear Information System (INIS)

1978-01-01

An account of certain well-established mathematical methods, which prove useful to deal with non-linear partial differential equations is presented. Within the strict framework of Functional Analysis, it describes Semigroup Techniques in Banach Spaces as well as variational approaches towards critical points. Detailed proofs are given of the existence of local and global solutions of the Cauchy problem and of the stability of stationary solutions. The formal approach based upon invariance under Lie transformations deserves attention due to its wide range of applicability, even if the explicit solutions thus obtained do not allow for a deep analysis of the equations. A compre ensive introduction to the inverse scattering approach and to the solution concept for certain non-linear equations of physical interest are also presented. A detailed discussion is made about certain convergence and stability problems which arise in importance need not be emphasized. (author) [es

Regarding on the exact solutions for the nonlinear fractional differential equations

Directory of Open Access Journals (Sweden)

Kaplan Melike

2016-01-01

Full Text Available In this work, we have considered the modified simple equation (MSE method for obtaining exact solutions of nonlinear fractional-order differential equations. The space-time fractional equal width (EW and the modified equal width (mEW equation are considered for illustrating the effectiveness of the algorithm. It has been observed that all exact solutions obtained in this paper verify the nonlinear ordinary differential equations which was obtained from nonlinear fractional-order differential equations under the terms of wave transformation relationship. The obtained results are shown graphically.

Inverse Boundary Value Problem for Non-linear Hyperbolic Partial Differential Equations

OpenAIRE

Nakamura, Gen; Vashisth, Manmohan

2017-01-01

In this article we are concerned with an inverse boundary value problem for a non-linear wave equation of divergence form with space dimension $n\\geq 3$. This non-linear wave equation has a trivial solution, i.e. zero solution. By linearizing this equation at the trivial solution, we have the usual linear isotropic wave equation with the speed $\\sqrt{\\gamma(x)}$ at each point $x$ in a given spacial domain. For any small solution $u=u(t,x)$ of this non-linear equation, we have the linear isotr...

Quasi-periodic solutions of nonlinear beam equations with quintic quasi-periodic nonlinearities

Directory of Open Access Journals (Sweden)

Qiuju Tuo

2015-01-01

Full Text Available In this article, we consider the one-dimensional nonlinear beam equations with quasi-periodic quintic nonlinearities $$ u_{tt}+u_{xxxx}+(B+ \\varepsilon\\phi(tu^5=0 $$ under periodic boundary conditions, where B is a positive constant, $\\varepsilon$ is a small positive parameter, $\\phi(t$ is a real analytic quasi-periodic function in t with frequency vector $\\omega=(\\omega_1,\\omega_2,\\dots,\\omega_m$. It is proved that the above equation admits many quasi-periodic solutions by KAM theory and partial Birkhoff normal form.

A vectorized Poisson solver over a spherical shell and its application to the quasi-geostrophic omega-equation

Science.gov (United States)

Mullenmeister, Paul

1988-01-01

The quasi-geostrophic omega-equation in flux form is developed as an example of a Poisson problem over a spherical shell. Solutions of this equation are obtained by applying a two-parameter Chebyshev solver in vector layout for CDC 200 series computers. The performance of this vectorized algorithm greatly exceeds the performance of its scalar analog. The algorithm generates solutions of the omega-equation which are compared with the omega fields calculated with the aid of the mass continuity equation.

A Generalized FDM for solving the Poisson's Equation on 3D Irregular Domains

Directory of Open Access Journals (Sweden)

J. Izadian

2014-01-01

Full Text Available In this paper a new method for solving the Poisson's equation with Dirichlet conditions on irregular domains is presented. For this purpose a generalized finite differences method is applied for numerical differentiation on irregular meshes. Three examples on cylindrical and spherical domains are considered. The numerical results are compared with analytical solution. These results show the performance and efficiency of the proposed method.

On the solution of the nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Zayed, E.M.E.; Zedan, Hassan A.

2003-01-01

In this paper we study the nonlinear Schrodinger equation with respect to the unknown function S(x,t). New dimensional reduction and exact solution for a nonlinear Schrodinger equation are presented and a complete group classification is given with respect to the function S(x,t). Moreover, specializing the potential function S(x,t), new classes of invariant solution and group classification are obtained in the cases of physical interest

A Fortran program (RELAX3D) to solve the 3 dimensional Poisson (Laplace) equation

International Nuclear Information System (INIS)

Houtman, H.; Kost, C.J.

1983-09-01

RELAX3D is an efficient, user friendly, interactive FORTRAN program which solves the Poisson (Laplace) equation Λ 2 =p for a general 3 dimensional geometry consisting of Dirichlet and Neumann boundaries approximated to lie on a regular 3 dimensional mesh. The finite difference equations at these nodes are solved using a successive point-iterative over-relaxation method. A menu of commands, supplemented by HELP facility, controls the dynamic loading of the subroutine describing the problem case, the iterations to converge to a solution, and the contour plotting of any desired slices, etc

A new sub-equation method applied to obtain exact travelling wave solutions of some complex nonlinear equations

International Nuclear Information System (INIS)

Zhang Huiqun

2009-01-01

By using a new coupled Riccati equations, a direct algebraic method, which was applied to obtain exact travelling wave solutions of some complex nonlinear equations, is improved. And the exact travelling wave solutions of the complex KdV equation, Boussinesq equation and Klein-Gordon equation are investigated using the improved method. The method presented in this paper can also be applied to construct exact travelling wave solutions for other nonlinear complex equations.

Sensitivity theory for general non-linear algebraic equations with constraints

International Nuclear Information System (INIS)

Oblow, E.M.

1977-04-01

Sensitivity theory has been developed to a high state of sophistication for applications involving solutions of the linear Boltzmann equation or approximations to it. The success of this theory in the field of radiation transport has prompted study of possible extensions of the method to more general systems of non-linear equations. Initial work in the U.S. and in Europe on the reactor fuel cycle shows that the sensitivity methodology works equally well for those non-linear problems studied to date. The general non-linear theory for algebraic equations is summarized and applied to a class of problems whose solutions are characterized by constrained extrema. Such equations form the basis of much work on energy systems modelling and the econometrics of power production and distribution. It is valuable to have a sensitivity theory available for these problem areas since it is difficult to repeatedly solve complex non-linear equations to find out the effects of alternative input assumptions or the uncertainties associated with predictions of system behavior. The sensitivity theory for a linear system of algebraic equations with constraints which can be solved using linear programming techniques is discussed. The role of the constraints in simplifying the problem so that sensitivity methodology can be applied is highlighted. The general non-linear method is summarized and applied to a non-linear programming problem in particular. Conclusions are drawn in about the applicability of the method for practical problems

Soliton solution for nonlinear partial differential equations by cosine-function method

International Nuclear Information System (INIS)

Ali, A.H.A.; Soliman, A.A.; Raslan, K.R.

2007-01-01

In this Letter, we established a traveling wave solution by using Cosine-function algorithm for nonlinear partial differential equations. The method is used to obtain the exact solutions for five different types of nonlinear partial differential equations such as, general equal width wave equation (GEWE), general regularized long wave equation (GRLW), general Korteweg-de Vries equation (GKdV), general improved Korteweg-de Vries equation (GIKdV), and Coupled equal width wave equations (CEWE), which are the important soliton equations

Multi-symplectic Runge-Kutta methods for nonlinear Dirac equations

International Nuclear Information System (INIS)

Hong Jialin; Li Chun

2006-01-01

In this paper, we consider the multi-symplectic Runge-Kutta (MSRK) methods applied to the nonlinear Dirac equation in relativistic quantum physics, based on a discovery of the multi-symplecticity of the equation. In particular, the conservation of energy, momentum and charge under MSRK discretizations is investigated by means of numerical experiments and numerical comparisons with non-MSRK methods. Numerical experiments presented reveal that MSRK methods applied to the nonlinear Dirac equation preserve exactly conservation laws of charge and momentum, and conserve the energy conservation in the corresponding numerical accuracy to the method utilized. It is verified numerically that MSRK methods are stable and convergent with respect to the conservation laws of energy, momentum and charge, and MSRK methods preserve not only the inner geometric structure of the equation, but also some crucial conservative properties in quantum physics. A remarkable advantage of MSRK methods applied to the nonlinear Dirac equation is the precise preservation of charge conservation law

Iterative solution of the semiconductor device equations

Energy Technology Data Exchange (ETDEWEB)

Bova, S.W.; Carey, G.F. [Univ. of Texas, Austin, TX (United States)

1996-12-31

Most semiconductor device models can be described by a nonlinear Poisson equation for the electrostatic potential coupled to a system of convection-reaction-diffusion equations for the transport of charge and energy. These equations are typically solved in a decoupled fashion and e.g. Newton`s method is used to obtain the resulting sequences of linear systems. The Poisson problem leads to a symmetric, positive definite system which we solve iteratively using conjugate gradient. The transport equations lead to nonsymmetric, indefinite systems, thereby complicating the selection of an appropriate iterative method. Moreover, their solutions exhibit steep layers and are subject to numerical oscillations and instabilities if standard Galerkin-type discretization strategies are used. In the present study, we use an upwind finite element technique for the transport equations. We also evaluate the performance of different iterative methods for the transport equations and investigate various preconditioners for a few generalized gradient methods. Numerical examples are given for a representative two-dimensional depletion MOSFET.

Nonlinear differential equations with exact solutions expressed via the Weierstrass function

NARCIS (Netherlands)

Kudryashov, NA

2004-01-01

A new problem is studied, that is to find nonlinear differential equations with special solutions expressed via the Weierstrass function. A method is discussed to construct nonlinear ordinary differential equations with exact solutions. The main step of our method is the assumption that nonlinear

The numerical dynamic for highly nonlinear partial differential equations

Science.gov (United States)

Lafon, A.; Yee, H. C.

1992-01-01

Problems associated with the numerical computation of highly nonlinear equations in computational fluid dynamics are set forth and analyzed in terms of the potential ranges of spurious behaviors. A reaction-convection equation with a nonlinear source term is employed to evaluate the effects related to spatial and temporal discretizations. The discretization of the source term is described according to several methods, and the various techniques are shown to have a significant effect on the stability of the spurious solutions. Traditional linearized stability analyses cannot provide the level of confidence required for accurate fluid dynamics computations, and the incorporation of nonlinear analysis is proposed. Nonlinear analysis based on nonlinear dynamical systems complements the conventional linear approach and is valuable in the analysis of hypersonic aerodynamics and combustion phenomena.

Nonlinear electrostatic wave equations for magnetized plasmas - II

DEFF Research Database (Denmark)

Dysthe, K. B.; Mjølhus, E.; Pécseli, H. L.

1985-01-01

For pt.I see ibid., vol.26, p.443-7 (1984). The problem of extending the high frequency part of the Zakharov equations for nonlinear electrostatic waves to magnetized plasmas, is considered. Weak electromagnetic and thermal effects are retained on an equal footing. Direction dependent (electrosta......For pt.I see ibid., vol.26, p.443-7 (1984). The problem of extending the high frequency part of the Zakharov equations for nonlinear electrostatic waves to magnetized plasmas, is considered. Weak electromagnetic and thermal effects are retained on an equal footing. Direction dependent...... (electrostatic) cut-off implies that various cases must be considered separately, leading to equations with rather different properties. Various equations encountered previously in the literature are recovered as limiting cases....

Exact solutions for nonlinear evolution equations using Exp-function method

International Nuclear Information System (INIS)

Bekir, Ahmet; Boz, Ahmet

2008-01-01

In this Letter, the Exp-function method is used to construct solitary and soliton solutions of nonlinear evolution equations. The Klein-Gordon, Burger-Fisher and Sharma-Tasso-Olver equations are chosen to illustrate the effectiveness of the method. The method is straightforward and concise, and its applications are promising. The Exp-function method presents a wider applicability for handling nonlinear wave equations

GHM method for obtaining rationalsolutions of nonlinear differential equations.

Science.gov (United States)

Vazquez-Leal, Hector; Sarmiento-Reyes, Arturo

2015-01-01

In this paper, we propose the application of the general homotopy method (GHM) to obtain rational solutions of nonlinear differential equations. It delivers a high precision representation of the nonlinear differential equation using a few linear algebraic terms. In order to assess the benefits of this proposal, three nonlinear problems are solved and compared against other semi-analytic methods or numerical methods. The obtained results show that GHM is a powerful tool, capable to generate highly accurate rational solutions. AMS subject classification 34L30.

Explicit Solutions for Generalized (2+1)-Dimensional Nonlinear Zakharov-Kuznetsov Equation

International Nuclear Information System (INIS)

Sun Yuhuai; Ma Zhimin; Li Yan

2010-01-01

The exact solutions of the generalized (2+1)-dimensional nonlinear Zakharov-Kuznetsov (Z-K) equation are explored by the method of the improved generalized auxiliary differential equation. Many explicit analytic solutions of the Z-K equation are obtained. The methods used to solve the Z-K equation can be employed in further work to establish new solutions for other nonlinear partial differential equations. (general)

Geometric discretization of the multidimensional Dirac delta distribution - Application to the Poisson equation with singular source terms

Science.gov (United States)

Egan, Raphael; Gibou, Frédéric

2017-10-01

We present a discretization method for the multidimensional Dirac distribution. We show its applicability in the context of integration problems, and for discretizing Dirac-distributed source terms in Poisson equations with constant or variable diffusion coefficients. The discretization is cell-based and can thus be applied in a straightforward fashion to Quadtree/Octree grids. The method produces second-order accurate results for integration. Superlinear convergence is observed when it is used to model Dirac-distributed source terms in Poisson equations: the observed order of convergence is 2 or slightly smaller. The method is consistent with the discretization of Dirac delta distribution for codimension one surfaces presented in [1,2]. We present Quadtree/Octree construction procedures to preserve convergence and present various numerical examples, including multi-scale problems that are intractable with uniform grids.

A Study of Schrödinger–Type Equations Appearing in Bohmian Mechanics and in the Theory of Bose–Einstein Condensates

KAUST Repository

Sierra Nunez, Jesus Alfredo

2018-05-16

The Schrödinger equations have had a profound impact on a wide range of fields of modern science, including quantum mechanics, superfluidity, geometrical optics, Bose-Einstein condensates, and the analysis of dispersive phenomena in the theory of PDE. The main purpose of this thesis is to explore two Schrödinger-type equations appearing in the so-called Bohmian formulation of quantum mechanics and in the study of exciton-polariton condensates. For the first topic, the linear Schrödinger equation is the starting point in the formulation of a phase-space model proposed in [1] for the Bohmian interpretation of quantum mechanics. We analyze this model, a nonlinear Vlasov-type equation, as a Hamiltonian system defined on an appropriate Poisson manifold built on Wasserstein spaces, the aim being to establish its existence theory. For this purpose, we employ results from the theory of PDE, optimal transportation, differential geometry and algebraic topology. The second topic of the thesis is the study of a nonlinear Schrödinger equation, called the complex Gross-Pitaevskii equation, appearing in the context of Bose-Einstein condensation of exciton-polaritons. This model can be roughly described as a driven-damped Gross-Pitaevskii equation which shares some similarities with the complex Ginzburg-Landau equation. The difficulties in the analysis of this equation stem from the fact that, unlike the complex Ginzburg-Landau equation, the complex Gross-Pitaevskii equation does not include a viscous dissipation term. Our approach to this equation will be in the framework of numerical computations, using two main tools: collocation methods and numerical continuation for the stationary solutions and a time-splitting spectral method for the dynamics. After performing a linear stability analysis on the computed stationary solutions, we are led to postulate the existence of radially symmetric stationary ground state solutions only for certain values of the parameters in the

Weierstrass Elliptic Function Solutions to Nonlinear Evolution Equations

International Nuclear Information System (INIS)

Yu Jianping; Sun Yongli

2008-01-01

This paper is based on the relations between projection Riccati equations and Weierstrass elliptic equation, combined with the Groebner bases in the symbolic computation. Then the novel method for constructing the Weierstrass elliptic solutions to the nonlinear evolution equations is given by using the above relations

Multiphase Weakly Nonlinear Geometric Optics for Schrödinger Equations

KAUST Repository

Carles, Ré mi; Dumas, Eric; Sparber, Christof

2010-01-01

We describe and rigorously justify the nonlinear interaction of highly oscillatory waves in nonlinear Schrödinger equations, posed on Euclidean space or on the torus. Our scaling corresponds to a weakly nonlinear regime where the nonlinearity affects the leading order amplitude of the solution, but does not alter the rapid oscillations. We consider initial states which are superpositions of slowly modulated plane waves, and use the framework of Wiener algebras. A detailed analysis of the corresponding nonlinear wave mixing phenomena is given, including a geometric interpretation of the resonance structure for cubic nonlinearities. As an application, we recover and extend some instability results for the nonlinear Schrödinger equation on the torus in negative order Sobolev spaces. © 2010 Society for Industrial and Applied Mathematics.

Discrete maximum principle for Poisson equation with mixed boundary conditions solved by hp-FEM

Czech Academy of Sciences Publication Activity Database

Vejchodský, Tomáš; Šolín, P.

2009-01-01

Roč. 1, č. 2 (2009), s. 201-214 ISSN 2070-0733 R&D Projects: GA AV ČR IAA100760702; GA ČR(CZ) GA102/07/0496; GA ČR GA102/05/0629 Institutional research plan: CEZ:AV0Z10190503 Keywords : discrete maximum principle * hp-FEM * Poisson equation * mixed boundary conditions Subject RIV: BA - General Mathematics

Quantum hydrodynamics and nonlinear differential equations for degenerate Fermi gas

International Nuclear Information System (INIS)

Bettelheim, Eldad; Abanov, Alexander G; Wiegmann, Paul B

2008-01-01

We present new nonlinear differential equations for spacetime correlation functions of Fermi gas in one spatial dimension. The correlation functions we consider describe non-stationary processes out of equilibrium. The equations we obtain are integrable equations. They generalize known nonlinear differential equations for correlation functions at equilibrium [1-4] and provide vital tools for studying non-equilibrium dynamics of electronic systems. The method we developed is based only on Wick's theorem and the hydrodynamic description of the Fermi gas. Differential equations appear directly in bilinear form. (fast track communication)

Lie-transformed action principle for classical plasma dynamics

International Nuclear Information System (INIS)

Kaufman, A.N.

1984-06-01

The Lie transform for a single particle in a wave is embedded in a Lagrangian action principle for self-consistent plasma dynamics. Variation of the action then yields the Vlasov equation for the oscillation-center distribution, the ray equations and amplitude transport for the wave, and the Poisson equation for the quasistatic field

Collapse in a forced three-dimensional nonlinear Schrodinger equation

DEFF Research Database (Denmark)

Lushnikov, P.M.; Saffman, M.

2000-01-01

We derive sufficient conditions for the occurrence of collapse in a forced three-dimensional nonlinear Schrodinger equation without dissipation. Numerical studies continue the results to the case of finite dissipation.......We derive sufficient conditions for the occurrence of collapse in a forced three-dimensional nonlinear Schrodinger equation without dissipation. Numerical studies continue the results to the case of finite dissipation....

Nonlinear localized dust acoustic waves in a charge varying dusty plasma with nonthermal ions

International Nuclear Information System (INIS)

Tribeche, Mouloud; Amour, Rabia

2007-01-01

A numerical investigation is presented to show the existence, formation, and possible realization of large-amplitude dust acoustic (DA) solitary waves in a charge varying dusty plasma with nonthermal ions. These nonlinear localized structures are self-consistent solutions of the collisionless Vlasov equation with a population of fast particles. The spatial patterns of the variable charge DA solitary wave are significantly modified by the nonthermal effects. The results complement and provide new insights into previously published results on this problem

SMPBS: Web server for computing biomolecular electrostatics using finite element solvers of size modified Poisson-Boltzmann equation.

Science.gov (United States)

Xie, Yang; Ying, Jinyong; Xie, Dexuan

2017-03-30

SMPBS (Size Modified Poisson-Boltzmann Solvers) is a web server for computing biomolecular electrostatics using finite element solvers of the size modified Poisson-Boltzmann equation (SMPBE). SMPBE not only reflects ionic size effects but also includes the classic Poisson-Boltzmann equation (PBE) as a special case. Thus, its web server is expected to have a broader range of applications than a PBE web server. SMPBS is designed with a dynamic, mobile-friendly user interface, and features easily accessible help text, asynchronous data submission, and an interactive, hardware-accelerated molecular visualization viewer based on the 3Dmol.js library. In particular, the viewer allows computed electrostatics to be directly mapped onto an irregular triangular mesh of a molecular surface. Due to this functionality and the fast SMPBE finite element solvers, the web server is very efficient in the calculation and visualization of electrostatics. In addition, SMPBE is reconstructed using a new objective electrostatic free energy, clearly showing that the electrostatics and ionic concentrations predicted by SMPBE are optimal in the sense of minimizing the objective electrostatic free energy. SMPBS is available at the URL: smpbs.math.uwm.edu © 2017 Wiley Periodicals, Inc. © 2017 Wiley Periodicals, Inc.

Non-holonomic dynamics and Poisson geometry

International Nuclear Information System (INIS)

Borisov, A V; Mamaev, I S; Tsiganov, A V

2014-01-01

This is a survey of basic facts presently known about non-linear Poisson structures in the analysis of integrable systems in non-holonomic mechanics. It is shown that by using the theory of Poisson deformations it is possible to reduce various non-holonomic systems to dynamical systems on well-understood phase spaces equipped with linear Lie-Poisson brackets. As a result, not only can different non-holonomic systems be compared, but also fairly advanced methods of Poisson geometry and topology can be used for investigating them. Bibliography: 95 titles

A procedure to construct exact solutions of nonlinear fractional differential equations.

Science.gov (United States)

Güner, Özkan; Cevikel, Adem C

2014-01-01

We use the fractional transformation to convert the nonlinear partial fractional differential equations with the nonlinear ordinary differential equations. The Exp-function method is extended to solve fractional partial differential equations in the sense of the modified Riemann-Liouville derivative. We apply the Exp-function method to the time fractional Sharma-Tasso-Olver equation, the space fractional Burgers equation, and the time fractional fmKdV equation. As a result, we obtain some new exact solutions.

Numerical treatments for solving nonlinear mixed integral equation

Directory of Open Access Journals (Sweden)

M.A. Abdou

2016-12-01

Full Text Available We consider a mixed type of nonlinear integral equation (MNLIE of the second kind in the space C[0,T]×L2(Ω,T<1. The Volterra integral terms (VITs are considered in time with continuous kernels, while the Fredholm integral term (FIT is considered in position with singular general kernel. Using the quadratic method and separation of variables method, we obtain a nonlinear system of Fredholm integral equations (NLSFIEs with singular kernel. A Toeplitz matrix method, in each case, is then used to obtain a nonlinear algebraic system. Numerical results are calculated when the kernels take a logarithmic form or Carleman function. Moreover, the error estimates, in each case, are then computed.

Maintaining the stability of nonlinear differential equations by the enhancement of HPM

International Nuclear Information System (INIS)

Hosein Nia, S.H.; Ranjbar, A.N.; Ganji, D.D.; Soltani, H.; Ghasemi, J.

2008-01-01

Homotopy perturbation method is an effective method to find a solution of a nonlinear differential equation. In this method, a nonlinear complex differential equation is transformed to a series of linear and nonlinear parts, almost simpler differential equations. These sets of equations are then solved iteratively. Finally, a linear series of the solutions completes the answer if the convergence is maintained. In this Letter, the need for stability verification is shown through some examples. Consequently, HPM is enhanced by a preliminary assumption. The idea is to keep the inherent stability of nonlinear dynamic, even the selected linear part is not

Lattice Boltzmann model for high-order nonlinear partial differential equations.

Science.gov (United States)

Chai, Zhenhua; He, Nanzhong; Guo, Zhaoli; Shi, Baochang

2018-01-01

In this paper, a general lattice Boltzmann (LB) model is proposed for the high-order nonlinear partial differential equation with the form ∂_{t}ϕ+∑_{k=1}^{m}α_{k}∂_{x}^{k}Π_{k}(ϕ)=0 (1≤k≤m≤6), α_{k} are constant coefficients, Π_{k}(ϕ) are some known differential functions of ϕ. As some special cases of the high-order nonlinear partial differential equation, the classical (m)KdV equation, KdV-Burgers equation, K(n,n)-Burgers equation, Kuramoto-Sivashinsky equation, and Kawahara equation can be solved by the present LB model. Compared to the available LB models, the most distinct characteristic of the present model is to introduce some suitable auxiliary moments such that the correct moments of equilibrium distribution function can be achieved. In addition, we also conducted a detailed Chapman-Enskog analysis, and found that the high-order nonlinear partial differential equation can be correctly recovered from the proposed LB model. Finally, a large number of simulations are performed, and it is found that the numerical results agree with the analytical solutions, and usually the present model is also more accurate than the existing LB models [H. Lai and C. Ma, Sci. China Ser. G 52, 1053 (2009)1672-179910.1007/s11433-009-0149-3; H. Lai and C. Ma, Phys. A (Amsterdam) 388, 1405 (2009)PHYADX0378-437110.1016/j.physa.2009.01.005] for high-order nonlinear partial differential equations.

Lattice Boltzmann model for high-order nonlinear partial differential equations

Science.gov (United States)

Chai, Zhenhua; He, Nanzhong; Guo, Zhaoli; Shi, Baochang

2018-01-01

In this paper, a general lattice Boltzmann (LB) model is proposed for the high-order nonlinear partial differential equation with the form ∂tϕ +∑k=1mαk∂xkΠk(ϕ ) =0 (1 ≤k ≤m ≤6 ), αk are constant coefficients, Πk(ϕ ) are some known differential functions of ϕ . As some special cases of the high-order nonlinear partial differential equation, the classical (m)KdV equation, KdV-Burgers equation, K (n ,n ) -Burgers equation, Kuramoto-Sivashinsky equation, and Kawahara equation can be solved by the present LB model. Compared to the available LB models, the most distinct characteristic of the present model is to introduce some suitable auxiliary moments such that the correct moments of equilibrium distribution function can be achieved. In addition, we also conducted a detailed Chapman-Enskog analysis, and found that the high-order nonlinear partial differential equation can be correctly recovered from the proposed LB model. Finally, a large number of simulations are performed, and it is found that the numerical results agree with the analytical solutions, and usually the present model is also more accurate than the existing LB models [H. Lai and C. Ma, Sci. China Ser. G 52, 1053 (2009), 10.1007/s11433-009-0149-3; H. Lai and C. Ma, Phys. A (Amsterdam) 388, 1405 (2009), 10.1016/j.physa.2009.01.005] for high-order nonlinear partial differential equations.

Evolutionary algorithm based heuristic scheme for nonlinear heat transfer equations.

Science.gov (United States)

Ullah, Azmat; Malik, Suheel Abdullah; Alimgeer, Khurram Saleem

2018-01-01

In this paper, a hybrid heuristic scheme based on two different basis functions i.e. Log Sigmoid and Bernstein Polynomial with unknown parameters is used for solving the nonlinear heat transfer equations efficiently. The proposed technique transforms the given nonlinear ordinary differential equation into an equivalent global error minimization problem. Trial solution for the given nonlinear differential equation is formulated using a fitness function with unknown parameters. The proposed hybrid scheme of Genetic Algorithm (GA) with Interior Point Algorithm (IPA) is opted to solve the minimization problem and to achieve the optimal values of unknown parameters. The effectiveness of the proposed scheme is validated by solving nonlinear heat transfer equations. The results obtained by the proposed scheme are compared and found in sharp agreement with both the exact solution and solution obtained by Haar Wavelet-Quasilinearization technique which witnesses the effectiveness and viability of the suggested scheme. Moreover, the statistical analysis is also conducted for investigating the stability and reliability of the presented scheme.

Evolutionary algorithm based heuristic scheme for nonlinear heat transfer equations.

Directory of Open Access Journals (Sweden)

Azmat Ullah

Full Text Available In this paper, a hybrid heuristic scheme based on two different basis functions i.e. Log Sigmoid and Bernstein Polynomial with unknown parameters is used for solving the nonlinear heat transfer equations efficiently. The proposed technique transforms the given nonlinear ordinary differential equation into an equivalent global error minimization problem. Trial solution for the given nonlinear differential equation is formulated using a fitness function with unknown parameters. The proposed hybrid scheme of Genetic Algorithm (GA with Interior Point Algorithm (IPA is opted to solve the minimization problem and to achieve the optimal values of unknown parameters. The effectiveness of the proposed scheme is validated by solving nonlinear heat transfer equations. The results obtained by the proposed scheme are compared and found in sharp agreement with both the exact solution and solution obtained by Haar Wavelet-Quasilinearization technique which witnesses the effectiveness and viability of the suggested scheme. Moreover, the statistical analysis is also conducted for investigating the stability and reliability of the presented scheme.

Modified harmonic balance method for the solution of nonlinear jerk equations

Science.gov (United States)

Rahman, M. Saifur; Hasan, A. S. M. Z.

2018-03-01

In this paper, a second approximate solution of nonlinear jerk equations (third order differential equation) can be obtained by using modified harmonic balance method. The method is simpler and easier to carry out the solution of nonlinear differential equations due to less number of nonlinear equations are required to solve than the classical harmonic balance method. The results obtained from this method are compared with those obtained from the other existing analytical methods that are available in the literature and the numerical method. The solution shows a good agreement with the numerical solution as well as the analytical methods of the available literature.

Simplified Model of Nonlinear Landau Damping

International Nuclear Information System (INIS)

Yampolsky, N.A.; Fisch, N.J.

2009-01-01

The nonlinear interaction of a plasma wave with resonant electrons results in a plateau in the electron distribution function close to the phase velocity of the plasma wave. As a result, Landau damping of the plasma wave vanishes and the resonant frequency of the plasma wave downshifts. However, this simple picture is invalid when the external driving force changes the plasma wave fast enough so that the plateau cannot be fully developed. A new model to describe amplification of the plasma wave including the saturation of Landau damping and the nonlinear frequency shift is proposed. The proposed model takes into account the change of the plasma wave amplitude and describes saturation of the Landau damping rate in terms of a single fluid equation, which simplifies the description of the inherently kinetic nature of Landau damping. A proposed fluid model, incorporating these simplifications, is verified numerically using a kinetic Vlasov code.

Damped nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Nicholson, D.R.; Goldman, M.V.

1976-01-01

High frequency electrostatic plasma oscillations described by the nonlinear Schrodinger equation in the presence of damping, collisional or Landau, are considered. At early times, Landau damping of an initial soliton profile results in a broader, but smaller amplitude soliton, while collisional damping reduces the soliton size everywhere; soliton speeds at early times are unchanged by either kind of damping. For collisional damping, soliton speeds are unchanged for all time

Newtonian hydrodynamic equations with relativistic pressure and velocity

Energy Technology Data Exchange (ETDEWEB)

Hwang, Jai-chan [Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Daegu 702-701 (Korea, Republic of); Noh, Hyerim [Korea Astronomy and Space Science Institute, Daejeon 305-348 (Korea, Republic of); Fabris, Júlio; Piattella, Oliver F.; Zimdahl, Winfried, E-mail: jchan@knu.ac.kr, E-mail: hr@kasi.re.kr, E-mail: fabris@pq.cnpq.br, E-mail: oliver.piattella@pq.cnpq.br, E-mail: winfried.zimdahl@pq.cnpq.br [Departamento de Fisica, Universidade Federal do Espirito Santo, Vitória (Brazil)

2016-07-01

We present a new approximation to include fully general relativistic pressure and velocity in Newtonian hydrodynamics. The energy conservation, momentum conservation and two Poisson's equations are consistently derived from Einstein's gravity in the zero-shear gauge assuming weak gravity and action-at-a-distance limit. The equations show proper special relativity limit in the absence of gravity. Our approximation is complementary to the post-Newtonian approximation and the equations are valid in fully nonlinear situations.

On implicit abstract neutral nonlinear differential equations

Energy Technology Data Exchange (ETDEWEB)

Hernández, Eduardo, E-mail: lalohm@ffclrp.usp.br [Universidade de São Paulo, Departamento de Computação e Matemática, Faculdade de Filosofia Ciências e Letras de Ribeirão Preto (Brazil); O’Regan, Donal, E-mail: donal.oregan@nuigalway.ie [National University of Ireland, School of Mathematics, Statistics and Applied Mathematics (Ireland)

2016-04-15

In this paper we continue our developments in Hernández and O’Regan (J Funct Anal 261:3457–3481, 2011) on the existence of solutions for abstract neutral differential equations. In particular we extend the results in Hernández and O’Regan (J Funct Anal 261:3457–3481, 2011) for the case of implicit nonlinear neutral equations and we focus on applications to partial “nonlinear” neutral differential equations. Some applications involving partial neutral differential equations are presented.

Numerical study of fractional nonlinear Schrodinger equations

KAUST Repository

Klein, Christian

2014-10-08

Using a Fourier spectral method, we provide a detailed numerical investigation of dispersive Schrödinger-type equations involving a fractional Laplacian in an one-dimensional case. By an appropriate choice of the dispersive exponent, both mass and energy sub- and supercritical regimes can be identified. This allows us to study the possibility of finite time blow-up versus global existence, the nature of the blow-up, the stability and instability of nonlinear ground states and the long-time dynamics of solutions. The latter is also studied in a semiclassical setting. Moreover, we numerically construct ground state solutions of the fractional nonlinear Schrödinger equation.

Nonlocal and nonlinear dispersion in a nonlinear Schrodinger-type equation: exotic solitons and short-wavelength instabilities

DEFF Research Database (Denmark)

Oster, Michael; Gaididei, Yuri B.; Johansson, Magnus

2004-01-01

We study the continuum limit of a nonlinear Schrodinger lattice model with both on-site and inter-site nonlinearities, describing weakly coupled optical waveguides or Bose-Einstein condensates. The resulting continuum nonlinear Schrodinger-type equation includes both nonlocal and nonlinear...

New classical r-matrices from integrable non-linear sigma-models

International Nuclear Information System (INIS)

Laartz, J.; Bordemann, M.; Forger, M.; Schaper, U.

1993-01-01

Non-linear sigma models on Riemannian symmetric spaces constitute the most general class of classical non-linear sigma models which are known to be integrable. Using the current algebra structure of these models their canonical structure is analyzed and it is shown that their non-ultralocal fundamental Poisson bracket relation is governed by a field dependent non antisymmetric r-matrix obeying a dynamical Yang Baxter equation. The fundamental Poisson bracket relations and the r-matrix are derived explicitly and a new kind of algebra is found that is supposed to replace the classical Yang Baxter algebra governing the canonical structure of ultralocal models. (Author) 9 refs

General form of the Euler-Poisson-Darboux equation and application of the transmutation method

Directory of Open Access Journals (Sweden)

Elina L. Shishkina

2017-07-01

Full Text Available In this article, we find solution representations in the compact integral form to the Cauchy problem for a general form of the Euler-Poisson-Darboux equation with Bessel operators via generalized translation and spherical mean operators for all values of the parameter k, including also not studying before exceptional odd negative values. We use a Hankel transform method to prove results in a unified way. Under additional conditions we prove that a distributional solution is a classical one too. A transmutation property for connected generalized spherical mean is proved and importance of applying transmutation methods for differential equations with Bessel operators is emphasized. The paper also contains a short historical introduction on differential equations with Bessel operators and a rather detailed reference list of monographs and papers on mathematical theory and applications of this class of differential equations.

Spectral transform and solvability of nonlinear evolution equations

International Nuclear Information System (INIS)

Degasperis, A.

1979-01-01

These lectures deal with an exciting development of the last decade, namely the resolving method based on the spectral transform which can be considered as an extension of the Fourier analysis to nonlinear evolution equations. Since many important physical phenomena are modeled by nonlinear partial wave equations this method is certainly a major breakthrough in mathematical physics. We follow the approach, introduced by Calogero, which generalizes the usual Wronskian relations for solutions of a Sturm-Liouville problem. Its application to the multichannel Schroedinger problem will be the subject of these lectures. We will focus upon dynamical systems described at time t by a multicomponent field depending on one space coordinate only. After recalling the Fourier technique for linear evolution equations we introduce the spectral transform method taking the integral equations of potential scattering as an example. The second part contains all the basic functional relationships between the fields and their spectral transforms as derived from the Wronskian approach. In the third part we discuss a particular class of solutions of nonlinear evolution equations, solitons, which are considered by many physicists as a first step towards an elementary particle theory, because of their particle-like behaviour. The effect of the polarization time-dependence on the motion of the soliton is studied by means of the corresponding spectral transform, leading to new concepts such as the 'boomeron' and the 'trappon'. The rich dynamic structure is illustrated by a brief report on the main results of boomeron-boomeron and boomeron-trappon collisions. In the final section we discuss further results concerning important properties of the solutions of basic nonlinear equations. We introduce the Baecklund transform for the special case of scalar fields and demonstrate how it can be used to generate multisoliton solutions and how the conservation laws are obtained. (HJ)

Non-linear partial differential equations an algebraic view of generalized solutions

CERN Document Server

Rosinger, Elemer E

1990-01-01

A massive transition of interest from solving linear partial differential equations to solving nonlinear ones has taken place during the last two or three decades. The availability of better computers has often made numerical experimentations progress faster than the theoretical understanding of nonlinear partial differential equations. The three most important nonlinear phenomena observed so far both experimentally and numerically, and studied theoretically in connection with such equations have been the solitons, shock waves and turbulence or chaotical processes. In many ways, these phenomen

Reduction of the equation for lower hybrid waves in a plasma to a nonlinear Schroedinger equation

Science.gov (United States)

Karney, C. F. F.

1977-01-01

Equations describing the nonlinear propagation of waves in an anisotropic plasma are rarely exactly soluble. However it is often possible to make approximations that reduce the exact equations into a simpler equation. The use of MACSYMA to make such approximations, and so reduce the equation describing lower hybrid waves into the nonlinear Schrodinger equation which is soluble by the inverse scattering method is demonstrated. MACSYMA is used at several stages in the calculation only because there is a natural division between calculations that are easiest done by hand, and those that are easiest done by machine.

KINETIC-J: A computational kernel for solving the linearized Vlasov equation applied to calculations of the kinetic, configuration space plasma current for time harmonic wave electric fields

Science.gov (United States)

Green, David L.; Berry, Lee A.; Simpson, Adam B.; Younkin, Timothy R.

2018-04-01

We present the KINETIC-J code, a computational kernel for evaluating the linearized Vlasov equation with application to calculating the kinetic plasma response (current) to an applied time harmonic wave electric field. This code addresses the need for a configuration space evaluation of the plasma current to enable kinetic full-wave solvers for waves in hot plasmas to move beyond the limitations of the traditional Fourier spectral methods. We benchmark the kernel via comparison with the standard k →-space forms of the hot plasma conductivity tensor.

The application of He's exp-function method to a nonlinear differential-difference equation

International Nuclear Information System (INIS)

Dai Chaoqing; Cen Xu; Wu Shengsheng

2009-01-01

This paper applies He's exp-function method, which was originally proposed to find new exact travelling wave solutions of nonlinear partial differential equations (NPDEs) or coupled nonlinear partial differential equations (CNPDEs), to a nonlinear differential-difference equation, and some new travelling wave solutions are obtained.

Nonlinear Fredholm Integral Equation of the Second Kind with Quadrature Methods

Directory of Open Access Journals (Sweden)

M. Jafari Emamzadeh

2010-06-01

Full Text Available In this paper, a numerical method for solving the nonlinear Fredholm integral equation is presented. We intend to approximate the solution of this equation by quadrature methods and by doing so, we solve the nonlinear Fredholm integral equation more accurately. Several examples are given at the end of this paper

On the integrability of the generalized Fisher-type nonlinear diffusion equations

International Nuclear Information System (INIS)

Wang Dengshan; Zhang Zhifei

2009-01-01

In this paper, the geometric integrability and Lax integrability of the generalized Fisher-type nonlinear diffusion equations with modified diffusion in (1+1) and (2+1) dimensions are studied by the pseudo-spherical surface geometry method and prolongation technique. It is shown that the (1+1)-dimensional Fisher-type nonlinear diffusion equation is geometrically integrable in the sense of describing a pseudo-spherical surface of constant curvature -1 only for m = 2, and the generalized Fisher-type nonlinear diffusion equations in (1+1) and (2+1) dimensions are Lax integrable only for m = 2. This paper extends the results in Bindu et al 2001 (J. Phys. A: Math. Gen. 34 L689) and further provides the integrability information of (1+1)- and (2+1)-dimensional Fisher-type nonlinear diffusion equations for m = 2

The Volterra's integral equation theory for accelerator single-freedom nonlinear components

International Nuclear Information System (INIS)

Wang Sheng; Xie Xi

1996-01-01

The Volterra's integral equation equivalent to the dynamic equation of accelerator single-freedom nonlinear components is given, starting from which the transport operator of accelerator single-freedom nonlinear components and its inverse transport operator are obtained. Therefore, another algorithm for the expert system of the beam transport operator of accelerator single-freedom nonlinear components is developed

Analytic Approximation of the Solutions of Stochastic Differential Delay Equations with Poisson Jump and Markovian Switching

Directory of Open Access Journals (Sweden)

Hua Yang

2012-01-01

Full Text Available We are concerned with the stochastic differential delay equations with Poisson jump and Markovian switching (SDDEsPJMSs. Most SDDEsPJMSs cannot be solved explicitly as stochastic differential equations. Therefore, numerical solutions have become an important issue in the study of SDDEsPJMSs. The key contribution of this paper is to investigate the strong convergence between the true solutions and the numerical solutions to SDDEsPJMSs when the drift and diffusion coefficients are Taylor approximations.

Study of the heavy ions (Au+Au at 150 AMeV) collisions with the FOPI detector. Comparison with the Landau-Vlasov model; Etude des collisions d`ions lourds AU+AU a 150 A.MeV avec le detecteur FOPI. Comparaison avec le modele de Landau-Vlasov

Energy Technology Data Exchange (ETDEWEB)

Boussange, S

1995-09-15

In this thesis, heavy ions (Au+Au) collisions experiments are made at 150 AMeV.In the first part, a general study of the nuclear matter equation is presented. Then the used Landau-Vlasov theoretical model is describe. The third part presents the FOPI experience and the details of how to obtain this theoretical predictions (filter, cuts, corrections, possible centrality selections).At the end, experimental results and comparisons with the Landau-Vlasov model are presented. (TEC). 105 refs., 96 figs., 14 tabs.

Current interactions from the one-form sector of nonlinear higher-spin equations

Science.gov (United States)

Gelfond, O. A.; Vasiliev, M. A.

2018-06-01

The form of higher-spin current interactions in the sector of one-forms is derived from the nonlinear higher-spin equations in AdS4. Quadratic corrections to higher-spin equations are shown to be independent of the phase of the parameter η = exp ⁡ iφ in the full nonlinear higher-spin equations. The current deformation resulting from the nonlinear higher-spin equations is represented in the canonical form with the minimal number of space-time derivatives. The non-zero spin-dependent coupling constants of the resulting currents are determined in terms of the higher-spin coupling constant η η bar . Our results confirm the conjecture that (anti-)self-dual nonlinear higher-spin equations result from the full system at (η = 0) η bar = 0.

Existence and asymptotic behavior of solutions for nonlinear Schrödinger-Poisson systems with steep potential well.

Science.gov (United States)

Du, Miao; Tian, Lixin; Wang, Jun; Zhang, Fubao

2016-03-01

In this paper, we are concerned with a class of Schrödinger-Poisson systems with the asymptotically linear or asymptotically 3-linear nonlinearity. Under some suitable assumptions on V , K , a , and f , we prove the existence, nonexistence, and asymptotic behavior of solutions via variational methods. In particular, the potential V is allowed to be sign-changing for the asymptotically linear case.

Equating TIMSS Mathematics Subtests with Nonlinear Equating Methods Using NEAT Design: Circle-Arc Equating Approaches

Science.gov (United States)

Ozdemir, Burhanettin

2017-01-01

The purpose of this study is to equate Trends in International Mathematics and Science Study (TIMSS) mathematics subtest scores obtained from TIMSS 2011 to scores obtained from TIMSS 2007 form with different nonlinear observed score equating methods under Non-Equivalent Anchor Test (NEAT) design where common items are used to link two or more test…

Finding all solutions of nonlinear equations using the dual simplex method

Science.gov (United States)

Yamamura, Kiyotaka; Fujioka, Tsuyoshi

2003-03-01

Recently, an efficient algorithm has been proposed for finding all solutions of systems of nonlinear equations using linear programming. This algorithm is based on a simple test (termed the LP test) for nonexistence of a solution to a system of nonlinear equations using the dual simplex method. In this letter, an improved version of the LP test algorithm is proposed. By numerical examples, it is shown that the proposed algorithm could find all solutions of a system of 300 nonlinear equations in practical computation time.

Nonlinear evolution of f(R) cosmologies. I. Methodology

International Nuclear Information System (INIS)

Oyaizu, Hiroaki

2008-01-01

We introduce the method and the implementation of a cosmological simulation of a class of metric-variation f(R) models that accelerate the cosmological expansion without a cosmological constant and evade solar-system bounds of small-field deviations to general relativity. Such simulations are shown to reduce to solving a nonlinear Poisson equation for the scalar degree of freedom introduced by the f(R) modifications. We detail the method to efficiently solve the nonlinear Poisson equation by using a Newton-Gauss-Seidel relaxation scheme coupled with the multigrid method to accelerate the convergence. The simulations are shown to satisfy tests comparing the simulated outcome to analytical solutions for simple situations, and the dynamics of the simulations are tested with orbital and Zeldovich collapse tests. Finally, we present several static and dynamical simulations using realistic cosmological parameters to highlight the differences between standard physics and f(R) physics. In general, we find that the f(R) modifications result in stronger gravitational attraction that enhances the dark matter power spectrum by ∼20% for large but observationally allowed f(R) modifications. A more detailed study of the nonlinear f(R) effects on the power spectrum are presented in a companion paper.

Differential constraints and exact solutions of nonlinear diffusion equations

International Nuclear Information System (INIS)

Kaptsov, Oleg V; Verevkin, Igor V

2003-01-01

The differential constraints are applied to obtain explicit solutions of nonlinear diffusion equations. Certain linear determining equations with parameters are used to find such differential constraints. They generalize the determining equations used in the search for classical Lie symmetries

Allowable graphs of the nonlinear Schrödinger equation and their ...

Indian Academy of Sciences (India)

Bich Nguyen

2017-11-20

Nov 20, 2017 ... Non-linear Schrödinger equation; graphs; characteristic polynomial; .... Allowable graphs of the NLS and their applications. 795 ...... nonlinear Schroödinger equation, J. Algebra Appl. 16 (2017) 37 pp., https://doi.org/10.1142/.

Partially integrable nonlinear equations with one higher symmetry

International Nuclear Information System (INIS)

Mikhailov, A V; Novikov, V S; Wang, J P

2005-01-01

In this letter, we present a family of second order in time nonlinear partial differential equations, which have only one higher symmetry. These equations are not integrable, but have a solution depending on one arbitrary function. (letter to the editor)