Near integrability of kink lattice with higher order interactions

Science.gov (United States)

Jiang, Yun-Guo; Liu, Jia-Zhen; He, Song

2017-11-01

We make use of Manton’s analytical method to investigate the force between kinks and anti-kinks at large distances in 1+1 dimensional field theory. The related potential has infinite order corrections of exponential pattern, and the coefficients for each order are determined. These coefficients can also be obtained by solving the equation of the fluctuations around the vacuum. At the lowest order, the kink lattice represents the Toda lattice. With higher order correction terms, the kink lattice can represent one kind of generic Toda lattice. With only two sites, the kink lattice is classically integrable. If the number of sites of the lattice is larger than two, the kink lattice is not integrable but is a near integrable system. We make use of Flaschka’s variables to study the Lax pair of the kink lattice. These Flaschka’s variables have interesting algebraic relations and non-integrability can be manifested. We also discuss the higher Hamiltonians for the deformed open Toda lattice, which has a similar result to the ordinary deformed Toda. Supported by Shandong Provincial Natural Science Foundation (ZR2014AQ007), National Natural Science Foundation of China (11403015, U1531105), S. He is supported by Max-Planck fellowship in Germany and National Natural Science Foundation of China (11305235)

A lattice based solution of the collisional Boltzmann equation with applications to microchannel flows

International Nuclear Information System (INIS)

Green, B I; Vedula, Prakash

2013-01-01

An alternative approach for solution of the collisional Boltzmann equation for a lattice architecture is presented. In the proposed method, termed the collisional lattice Boltzmann method (cLBM), the effects of spatial transport are accounted for via a streaming operator, using a lattice framework, and the effects of detailed collisional interactions are accounted for using the full collision operator of the Boltzmann equation. The latter feature is in contrast to the conventional lattice Boltzmann methods (LBMs) where collisional interactions are modeled via simple equilibrium based relaxation models (e.g. BGK). The underlying distribution function is represented using weights and fixed velocity abscissas according to the lattice structure. These weights are evolved based on constraints on the evolution of generalized moments of velocity according to the collisional Boltzmann equation. It can be shown that the collision integral can be reduced to a summation of elementary integrals, which can be analytically evaluated. The proposed method is validated using studies of canonical microchannel Couette and Poiseuille flows (both body force and pressure driven) and the results are found to be in good agreement with those obtained from conventional LBMs and experiments where available. Unlike conventional LBMs, the proposed method does not involve any equilibrium based approximations and hence can be useful for simulation of highly nonequilibrium flows (for a range of Knudsen numbers) using a lattice framework. (paper)

An integrable coupling system of lattice hierarchy and its continuous limits

International Nuclear Information System (INIS)

Yu Fajun; Li Li

2009-01-01

In [E.G. Fan, Phys. Lett. A 372 (2008) 6368], Fan present a lattice hierarchy and its continuous limits. In this Letter, we extend this method, by introducing a complex discrete spectral problem, a coupling lattice hierarchy is derived. It is shown that a new sequence of combinations of complex lattice spectral problem converges to the integrable coupling couplings of soliton equation hierarchy, which has the integrable coupling system of AKNS hierarchy as a continuous limit.

Integrable light-cone lattice discretizations from the universal R-matrix

International Nuclear Information System (INIS)

Meneghelli, C.

2015-04-01

Our goal is to develop a more general scheme for constructing integrable lattice regularisations of integrable quantum field theories. Considering the affine Toda theories as examples, we show how to construct such lattice regularisations using the representation theory of quantum affine algebras. This requires us to clarify in particular the relations between the light-cone approach to integrable lattice models and the representation theory of quantum affine algebras. Both are found to be related in a very natural way, suggesting a general scheme for the construction of generalised Baxter Q-operators. One of the main difficulties we need to deal with is coming from the infinite-dimensionality of the relevant families of representations. It is handled by means of suitable renormalisation prescriptions defining what may be called the modular double of quantum affine algebras. This framework allows us to give a representation-theoretic proof of finite-difference equations generalising the Baxter equation.

Continuous limits for an integrable coupling system of Toda equation hierarchy

International Nuclear Information System (INIS)

Li Li; Yu Fajun

2009-01-01

In this Letter, we present an integrable coupling system of lattice hierarchy and its continuous limits by using of Lie algebra sl(4). By introducing a complex discrete spectral problem, the integrable coupling system of Toda lattice hierarchy is derived. It is shown that a new complex lattice spectral problem converges to the integrable couplings of discrete soliton equation hierarchy, which has the integrable coupling system of C-KdV hierarchy as a new kind of continuous limit.

Continuous limits for an integrable coupling system of Toda equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

Li Li [College of Maths and Systematic Science, Shenyang Normal University, Shenyang 110034 (China); Yu Fajun, E-mail: yfajun@163.co [College of Maths and Systematic Science, Shenyang Normal University, Shenyang 110034 (China)

2009-09-21

In this Letter, we present an integrable coupling system of lattice hierarchy and its continuous limits by using of Lie algebra sl(4). By introducing a complex discrete spectral problem, the integrable coupling system of Toda lattice hierarchy is derived. It is shown that a new complex lattice spectral problem converges to the integrable couplings of discrete soliton equation hierarchy, which has the integrable coupling system of C-KdV hierarchy as a new kind of continuous limit.

BCS @ 50: derivation of gap equations in different lattice geometries

International Nuclear Information System (INIS)

Saurabh Basu

2007-07-01

We rigorously derive BCS gap equations for a square, triangular and a honeycomb lattice using a two-dimensional t-J model. The gap equations in all the three lattice geometries look usual, with band indices appearing and a minor modification in the separable pair potential for the (two band) honeycomb lattice. In each case, the gap equation is solved (self consistently with the number equation) at low densities assuming singlet pairing. (author)

PREFACE: Symmetries and Integrability of Difference Equations

Science.gov (United States)

Doliwa, Adam; Korhonen, Risto; Lafortune, Stéphane

2007-10-01

to integrability. The first section contains a paper by T Hamamoto and K Kajiwara on hypergeometric solutions to the q-Painlevé equation of type A4(1). Discrete geometry. In this category there are three papers. J Cielinski offers a geometric definition and a spectral approach on pseudospherical surfaces on time scales, while A Doliwa considers generalized isothermic lattices. The paper by U Pinkall, B Springborn and S Weiss mann is concerned with a new doubly discrete analogue of smoke ring flow and the real time simulation of fluid flow. Integrable systems in statistical physics. Under this heading there is a paper by R J Baxter on corner transfer matrices in statistical mechanics, and a paper by S Boukraa, S Hassani, J-M Maillard, B M McCoy, J-A Weil and N Zenine where the authors consider Fuchs-Painlevé elliptic representation of the Painlevé VI equation. KP lattices and differential-difference hierarchies. In this section we have seven articles. C R Gilson, J J C Nimmo and Y Ohta consider quasideterminant solutions of a non-Abelian Hirota-Miwa equation, while B Grammaticos, A Ramani, V Papageorgiou, J Satsuma and R Willox discuss the construction of lump-like solutions of the Hirota-Miwa equation. J Hietarinta and C Viallet analyze the factorization process for lattice maps searching for integrable cases, the paper by X-B Hu and G-F Yu is concerned with integrable discretizations of the (2+1)-dimensional sinh-Gordon equation, and K Kajiwara, M Mazzocco and Y Ohta consider the Hankel determinant formula of the tau-functions of the Toda equation. Finally, V G Papageorgiou and A G Tongas study Yang-Baxter maps and multi-field integrable lattice equations, and H-Y Wang, X-B Hu and H-W Tam consider the two-dimensional Leznov lattice equation with self-consistent sources. Quantum integrable systems. This category contains a paper on q-extended eigenvectors of the integral and finite Fourier transforms by N M Atakishiyev, J P Rueda and K B Wolf, and an article by S

On discrete 2D integrable equations of higher order

International Nuclear Information System (INIS)

Adler, V E; Postnikov, V V

2014-01-01

We study two-dimensional discrete integrable equations of order 1 with respect to one independent variable and m with respect to another one. A generalization of the multidimensional consistency property is proposed for this type of equations. The examples are related to the Bäcklund–Darboux transformations for the lattice equations of Bogoyavlensky type. (paper)

New series of 3 D lattice integrable models

International Nuclear Information System (INIS)

Mangazeev, V.V.; Sergeev, S.M.; Stroganov, Yu.G.

1993-01-01

In this paper we present a new series of 3-dimensional integrable lattice models with N colors. The weight functions of the models satisfy modified tetrahedron equations with N states and give a commuting family of two-layer transfer-matrices. The dependence on the spectral parameters corresponds to the static limit of the modified tetrahedron equations and weights are parameterized in terms of elliptic functions. The models contain two free parameters: elliptic modulus and additional parameter η. 12 refs

Quiver gauge theories and integrable lattice models

International Nuclear Information System (INIS)

Yagi, Junya

2015-01-01

We discuss connections between certain classes of supersymmetric quiver gauge theories and integrable lattice models from the point of view of topological quantum field theories (TQFTs). The relevant classes include 4d N=1 theories known as brane box and brane tilling models, 3d N=2 and 2d N=(2,2) theories obtained from them by compactification, and 2d N=(0,2) theories closely related to these theories. We argue that their supersymmetric indices carry structures of TQFTs equipped with line operators, and as a consequence, are equal to the partition functions of lattice models. The integrability of these models follows from the existence of extra dimension in the TQFTs, which emerges after the theories are embedded in M-theory. The Yang-Baxter equation expresses the invariance of supersymmetric indices under Seiberg duality and its lower-dimensional analogs.

A Class of Two-Component Adler—Bobenko—Suris Lattice Equations

International Nuclear Information System (INIS)

Fu Wei; Zhang Da-Jun; Zhou Ru-Guang

2014-01-01

We study a class of two-component forms of the famous list of the Adler—Bobenko—Suris lattice equations. The obtained two-component lattice equations are still consistent around the cube and they admit solutions with ‘jumping properties’ between two levels. (general)

The multiscale expansions of difference equations in the small lattice spacing regime, and a vicinity and integrability test: I

International Nuclear Information System (INIS)

Santini, Paolo Maria

2010-01-01

We propose an algorithmic procedure (i) to study the 'distance' between an integrable PDE and any discretization of it, in the small lattice spacing ε regime, and, at the same time, (ii) to test the (asymptotic) integrability properties of such discretization. This method should provide, in particular, useful and concrete information on how good is any numerical scheme used to integrate a given integrable PDE. The procedure, illustrated on a fairly general ten-parameter family of discretizations of the nonlinear Schroedinger equation, consists of the following three steps: (i) the construction of the continuous multiscale expansion of a generic solution of the discrete system at all orders in ε, following Degasperis et al (1997 Physica D 100 187-211); (ii) the application, to such an expansion, of the Degasperis-Procesi (DP) integrability test (Degasperis A and Procesi M 1999 Asymptotic integrability Symmetry and Perturbation Theory, SPT98, ed A Degasperis and G Gaeta (Singapore: World Scientific) pp 23-37; Degasperis A 2001 Multiscale expansion and integrability of dispersive wave equations Lectures given at the Euro Summer School: 'What is integrability?' (Isaac Newton Institute, Cambridge, UK, 13-24 August); Integrability (Lecture Notes in Physics vol 767) ed A Mikhailov (Berlin: Springer)), to test the asymptotic integrability properties of the discrete system and its 'distance' from its continuous limit; (iii) the use of the main output of the DP test to construct infinitely many approximate symmetries and constants of motion of the discrete system, through novel and simple formulas.

New Positive and Negative Hierarchies of Integrable Differential-Difference Equations and Conservation Laws

International Nuclear Information System (INIS)

Li Xinyue; Zhao Qiulan

2009-01-01

Two hierarchies of nonlinear integrable positive and negative lattice equations are derived from a discrete spectral problem. The two lattice hierarchies are proved to have discrete zero curvature representations associated with a discrete spectral problem, which also shows that the positive and negative hierarchies correspond to positive and negative power expansions of Lax operators with respect to the spectral parameter, respectively. Moreover, the integrable lattice models in the positive hierarchy are of polynomial type, and the integrable lattice models in the negative hierarchy are of rational type. Further, we construct infinite conservation laws about the positive hierarchy.

Commutativity of the source generation procedure and integrable semi-discretizations: the two-dimensional Leznov lattice

International Nuclear Information System (INIS)

Hu Juan; Yu Guofu; Tam, Hon-Wah

2012-01-01

The source generation procedure (SGP) is applied to a y-directional discrete version and an x-directional discrete version of the Leznov lattice. Consequently, a y-discrete Leznov lattice equation with self-consistent sources (y-discrete Leznov ESCS) and an x-discrete Leznov ESCS are presented. Also utilizing the SGP, a new type of Leznov lattice equation with self-consistent sources (new Leznov ESCS) is derived. It is interesting that the two semi-discrete Leznov ESCS produced constitute a y-discretization for the Leznov ESCS given by Wang et al (2007 J. Phys. A: Math. Theor. 40 12691) and an x-discretization for the new Leznov ESCS, respectively. This means that the commutativity of SGP and integrable semi-discretizations is valid for the two-dimensional Leznov lattice equation. (paper)

The multiscale expansions of difference equations in the small lattice spacing regime, and a vicinity and integrability test: I

Science.gov (United States)

Santini, Paolo Maria

2010-01-01

We propose an algorithmic procedure (i) to study the 'distance' between an integrable PDE and any discretization of it, in the small lattice spacing epsilon regime, and, at the same time, (ii) to test the (asymptotic) integrability properties of such discretization. This method should provide, in particular, useful and concrete information on how good is any numerical scheme used to integrate a given integrable PDE. The procedure, illustrated on a fairly general ten-parameter family of discretizations of the nonlinear Schrödinger equation, consists of the following three steps: (i) the construction of the continuous multiscale expansion of a generic solution of the discrete system at all orders in epsilon, following Degasperis et al (1997 Physica D 100 187-211) (ii) the application, to such an expansion, of the Degasperis-Procesi (DP) integrability test (Degasperis A and Procesi M 1999 Asymptotic integrability Symmetry and Perturbation Theory, SPT98, ed A Degasperis and G Gaeta (Singapore: World Scientific) pp 23-37 Degasperis A 2001 Multiscale expansion and integrability of dispersive wave equations Lectures given at the Euro Summer School: 'What is integrability?' (Isaac Newton Institute, Cambridge, UK, 13-24 August); Integrability (Lecture Notes in Physics vol 767) ed A Mikhailov (Berlin: Springer)), to test the asymptotic integrability properties of the discrete system and its 'distance' from its continuous limit; (iii) the use of the main output of the DP test to construct infinitely many approximate symmetries and constants of motion of the discrete system, through novel and simple formulas.

An integrable coupling family of Merola-Ragnisco-Tu lattice systems, its Hamiltonian structure and related nonisospectral integrable lattice family

Energy Technology Data Exchange (ETDEWEB)

Xu Xixiang, E-mail: xu_xixiang@hotmail.co [College of Science, Shandong University of Science and Technology, Qingdao, 266510 (China)

2010-01-04

An integrable coupling family of Merola-Ragnisco-Tu lattice systems is derived from a four-by-four matrix spectral problem. The Hamiltonian structure of the resulting integrable coupling family is established by the discrete variational identity. Each lattice system in the resulting integrable coupling family is proved to be integrable discrete Hamiltonian system in Liouville sense. Ultimately, a nonisospectral integrable lattice family associated with the resulting integrable lattice family is constructed through discrete zero curvature representation.

An integrable coupling family of Merola-Ragnisco-Tu lattice systems, its Hamiltonian structure and related nonisospectral integrable lattice family

International Nuclear Information System (INIS)

Xu Xixiang

2010-01-01

An integrable coupling family of Merola-Ragnisco-Tu lattice systems is derived from a four-by-four matrix spectral problem. The Hamiltonian structure of the resulting integrable coupling family is established by the discrete variational identity. Each lattice system in the resulting integrable coupling family is proved to be integrable discrete Hamiltonian system in Liouville sense. Ultimately, a nonisospectral integrable lattice family associated with the resulting integrable lattice family is constructed through discrete zero curvature representation.

Infinitely many conservation laws for two integrable lattice hierarchies associated with a new discrete Schroedinger spectral problem

International Nuclear Information System (INIS)

Zhu, Zuo-nong; Tam, Hon-Wah; Ding, Qing

2003-01-01

In this Letter, by means of considering matrix form of a new Schroedinger discrete spectral operator equation, and constructing opportune time evolution equations, and using discrete zero curvature representation, two discrete integrable lattice hierarchies proposed by Boiti et al. [J. Phys. A: Math. Gen. 36 (2003) 139] are re-derived. From the matrix Lax representations, we demonstrate the existence of infinitely many conservation laws for the two lattice hierarchies and give the corresponding conserved densities and the associated fluxes by means of formulae. Thus their integrability is further confirmed. Specially we obtain the infinitely many conservation laws for a new discrete version of the KdV equation. A connection between the conservation laws of the discrete KdV equation and the ones of the KdV equation is discussed by two examples

Thermal equation of state for lattice Boltzmann gases

International Nuclear Information System (INIS)

Zheng, Ran

2009-01-01

The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar–Gross–Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar–Gross–Krook equation (LBGKE) induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics

Constructing New Discrete Integrable Coupling System for Soliton Equation by Kronecker Product

International Nuclear Information System (INIS)

Yu Fajun; Zhang Hongqing

2008-01-01

It is shown that the Kronecker product can be applied to constructing new discrete integrable coupling system of soliton equation hierarchy in this paper. A direct application to the fractional cubic Volterra lattice spectral problem leads to a novel integrable coupling system of soliton equation hierarchy. It is also indicated that the study of discrete integrable couplings by using the Kronecker product is an efficient and straightforward method. This method can be used generally

Symmetries and conserved quantities of discrete wave equation associated with the Ablowitz—Ladik—Lattice system

International Nuclear Information System (INIS)

Fu Jing-Li; He Yu-Fang; Hong Fang-Yu; Song Duan; Fu Hao

2013-01-01

In this paper, we present a new method to obtain the Lie symmetries and conserved quantities of the discrete wave equation with the Ablowitz—Ladik—Lattice equations. Firstly, the wave equation is transformed into a simple difference equation with the Ablowitz—Ladik—Lattice method. Secondly, according to the invariance of the discrete wave equation and the Ablowitz—Ladik—Lattice equations under infinitesimal transformation of dependent and independent variables, we derive the discrete determining equation and the discrete restricted equations. Thirdly, a series of the discrete analogs of conserved quantities, the discrete analogs of Lie groups, and the characteristic equations are obtained for the wave equation. Finally, we study a model of a biological macromolecule chain of mechanical behaviors, the Lie symmetry theory of discrete wave equation with the Ablowitz—Ladik—Lattice method is verified. (general)

New continual analogs of two-dimensional Toda lattices related with nonlinear integro-differential equations

International Nuclear Information System (INIS)

Savel'ev, M.V.

1988-01-01

Continual ''extensions'' of two-dimensional Toda lattices are proposed. They are described by integro-differential equations, generally speaking, with singular kernels, depending on new (third) variable. The problem of their integrability on the corresponding class of the initial discrete system solutions is discussed. The latter takes place, in particular, for the kernel coinciding with the causal function

The lattice Boltzmann model for the second-order Benjamin–Ono equations

International Nuclear Information System (INIS)

Lai, Huilin; Ma, Changfeng

2010-01-01

In this paper, in order to extend the lattice Boltzmann method to deal with more complicated nonlinear equations, we propose a 1D lattice Boltzmann scheme with an amending function for the second-order (1 + 1)-dimensional Benjamin–Ono equation. With the Taylor expansion and the Chapman–Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The equilibrium distribution function and the amending function are obtained. Numerical simulations are carried out for the 'good' Boussinesq equation and the 'bad' one to validate the proposed model. It is found that the numerical results agree well with the analytical solutions. The present model can be used to solve more kinds of nonlinear partial differential equations

Iterative resonance self-shielding methods using resonance integral table in heterogeneous transport lattice calculations

International Nuclear Information System (INIS)

Hong, Ser Gi; Kim, Kang-Seog

2011-01-01

This paper describes the iteration methods using resonance integral tables to estimate the effective resonance cross sections in heterogeneous transport lattice calculations. Basically, these methods have been devised to reduce an effort to convert resonance integral table into subgroup data to be used in the physical subgroup method. Since these methods do not use subgroup data but only use resonance integral tables directly, these methods do not include an error in converting resonance integral into subgroup data. The effective resonance cross sections are estimated iteratively for each resonance nuclide through the heterogeneous fixed source calculations for the whole problem domain to obtain the background cross sections. These methods have been implemented in the transport lattice code KARMA which uses the method of characteristics (MOC) to solve the transport equation. The computational results show that these iteration methods are quite promising in the practical transport lattice calculations.

Axisymmetric multiphase lattice Boltzmann method for generic equations of state

NARCIS (Netherlands)

Reijers, S.A.; Gelderblom, H.; Toschi, F.

2016-01-01

We present an axisymmetric lattice Boltzmann model based on the Kupershtokh et al. multiphase model that is capable of solving liquid–gas density ratios up to 103. Appropriate source terms are added to the lattice Boltzmann evolution equation to fully recover the axisymmetric multiphase conservation

Coupled double-distribution-function lattice Boltzmann method for the compressible Navier-Stokes equations.

Science.gov (United States)

Li, Q; He, Y L; Wang, Y; Tao, W Q

2007-11-01

A coupled double-distribution-function lattice Boltzmann method is developed for the compressible Navier-Stokes equations. Different from existing thermal lattice Boltzmann methods, this method can recover the compressible Navier-Stokes equations with a flexible specific-heat ratio and Prandtl number. In the method, a density distribution function based on a multispeed lattice is used to recover the compressible continuity and momentum equations, while the compressible energy equation is recovered by an energy distribution function. The energy distribution function is then coupled to the density distribution function via the thermal equation of state. In order to obtain an adjustable specific-heat ratio, a constant related to the specific-heat ratio is introduced into the equilibrium energy distribution function. Two different coupled double-distribution-function lattice Boltzmann models are also proposed in the paper. Numerical simulations are performed for the Riemann problem, the double-Mach-reflection problem, and the Couette flow with a range of specific-heat ratios and Prandtl numbers. The numerical results are found to be in excellent agreement with analytical and/or other solutions.

CLUB - a multigroup integral transport theory code for lattice calculations of PHWR cells

International Nuclear Information System (INIS)

Krishnani, P.D.

1992-01-01

The computer code CLUB has been developed to calculate lattice parameters as a function of burnup for a pressurised heavy water reactor (PHWR) lattice cell containing fuel in the form of cluster. It solves the multigroup integral transport equation by the method based on combination of small scale collision probability (CP) method and large scale interface current technique. The calculations are performed by using WIMS 69 group cross section library or its condensed versions of 27 or 28 group libraries. It can also compute Keff from the given geometrical buckling in the input using multigroup diffusion theory in fundamental mode. The first order differential burnup equations can be solved by either Trapezoidal rule or Runge-Kutta method. (author). 17 refs., 2 figs

Integral equations

CERN Document Server

Moiseiwitsch, B L

2005-01-01

Two distinct but related approaches hold the solutions to many mathematical problems--the forms of expression known as differential and integral equations. The method employed by the integral equation approach specifically includes the boundary conditions, which confers a valuable advantage. In addition, the integral equation approach leads naturally to the solution of the problem--under suitable conditions--in the form of an infinite series.Geared toward upper-level undergraduate students, this text focuses chiefly upon linear integral equations. It begins with a straightforward account, acco

Integrable discretizations for the short-wave model of the Camassa-Holm equation

International Nuclear Information System (INIS)

Feng Baofeng; Maruno, Ken-ichi; Ohta, Yasuhiro

2010-01-01

The link between the short-wave model of the Camassa-Holm equation (SCHE) and bilinear equations of the two-dimensional Toda lattice equation is clarified. The parametric form of the N-cuspon solution of the SCHE in Casorati determinant is then given. Based on the above finding, integrable semi-discrete and full-discrete analogues of the SCHE are constructed. The determinant solutions of both semi-discrete and fully discrete analogues of the SCHE are also presented.

Entire solutions for bistable lattice differential equations with obstacles

CERN Document Server

Hoffman, Aaron; Vleck, E S Van

2018-01-01

The authors consider scalar lattice differential equations posed on square lattices in two space dimensions. Under certain natural conditions they show that wave-like solutions exist when obstacles (characterized by "holes") are present in the lattice. Their work generalizes to the discrete spatial setting the results obtained in Berestycki, Hamel, and Matuno (2009) for the propagation of waves around obstacles in continuous spatial domains. The analysis hinges upon the development of sub and super-solutions for a class of discrete bistable reaction-diffusion problems and on a generalization of a classical result due to Aronson and Weinberger that concerns the spreading of localized disturbances.

Group integration for lattice gauge theory at large and at small coupling

International Nuclear Information System (INIS)

Brower, R.C.; Nauenberg, M.

1981-01-01

We consider the fundamental SU(N) invariant integrals encountered in Wilson's lattice QCD with an eye to analytical results for N → infinite and approximations for small g 2 at fixed N. We develop a new semiclassical technique starting from the Schwinger-Dyson equations cast in differential form to give an exact solution to the single-link integral for N → infinite. The third-order phase transition discovered by Gross and Witten for two-dimensional QCD occurs here for any dimension. Alternatively we parametrize directly the integral over the Haar measure and obtain approximate results for SU(N) using stationary phase at small g 2 . Remarkably the single-loop correction gives the exact answer at N = infinite. We show that the naive lattice string of Weingarten is obtained from N → infinite QCD in the limit of dimensions d → infinite. We discuss applications of our techniques to the 1/N expansion. (orig.)

Multidimensional singular integrals and integral equations

CERN Document Server

Mikhlin, Solomon Grigorievich; Stark, M; Ulam, S

1965-01-01

Multidimensional Singular Integrals and Integral Equations presents the results of the theory of multidimensional singular integrals and of equations containing such integrals. Emphasis is on singular integrals taken over Euclidean space or in the closed manifold of Liapounov and equations containing such integrals. This volume is comprised of eight chapters and begins with an overview of some theorems on linear equations in Banach spaces, followed by a discussion on the simplest properties of multidimensional singular integrals. Subsequent chapters deal with compounding of singular integrals

Manipulation and quantification of microtubule lattice integrity

Directory of Open Access Journals (Sweden)

Taylor A. Reid

2017-08-01

Full Text Available Microtubules are structural polymers that participate in a wide range of cellular functions. The addition and loss of tubulin subunits allows the microtubule to grow and shorten, as well as to develop and repair defects and gaps in its cylindrical lattice. These lattice defects act to modulate the interactions of microtubules with molecular motors and other microtubule-associated proteins. Therefore, tools to control and measure microtubule lattice structure will be invaluable for developing a quantitative understanding of how the structural state of the microtubule lattice may regulate its interactions with other proteins. In this work, we manipulated the lattice integrity of in vitro microtubules to create pools of microtubules with common nucleotide states, but with variations in structural states. We then developed a series of novel semi-automated analysis tools for both fluorescence and electron microscopy experiments to quantify the type and severity of alterations in microtubule lattice integrity. These techniques will enable new investigations that explore the role of microtubule lattice structure in interactions with microtubule-associated proteins.

Wave transmission in nonlinear lattices

International Nuclear Information System (INIS)

Hennig, D.; Tsironis, G.P.

1999-01-01

The interplay of nonlinearity with lattice discreteness leads to phenomena and propagation properties quite distinct from those appearing in continuous nonlinear systems. For a large variety of condensed matter and optics applications the continuous wave approximation is not appropriate. In the present review we discuss wave transmission properties in one dimensional nonlinear lattices. Our paradigmatic equations are discrete nonlinear Schroedinger equations and their study is done through a dynamical systems approach. We focus on stationary wave properties and utilize well known results from the theory of dynamical systems to investigate various aspects of wave transmission and wave localization. We analyze in detail the more general dynamical system corresponding to the equation that interpolates between the non-integrable discrete nonlinear Schroedinger equation and the integrable Albowitz-Ladik equation. We utilize this analysis in a nonlinear Kronig-Penney model and investigate transmission and band modification properties. We discuss the modifications that are effected through an electric field and the nonlinear Wannier-Stark localization effects that are induced. Several applications are described, such as polarons in one dimensional lattices, semiconductor superlattices and one dimensional nonlinear photonic band gap systems. (Copyright (c) 1999 Elsevier Science B.V., Amsterdam. All rights reserved.)

On reductions of the discrete Kadomtsev-Petviashvili-type equations

Science.gov (United States)

Fu, Wei; Nijhoff, Frank W.

2017-12-01

The reduction by restricting the spectral parameters k and k\\prime on a generic algebraic curve of degree N is performed for the discrete AKP, BKP and CKP equations, respectively. A variety of two-dimensional discrete integrable systems possessing a more general solution structure arise from the reduction, and in each case a unified formula for the generic positive integer N≥slant 2 is given to express the corresponding reduced integrable lattice equations. The obtained extended two-dimensional lattice models give rise to many important integrable partial difference equations as special degenerations. Some new integrable lattice models such as the discrete Sawada-Kotera, Kaup-Kupershmidt and Hirota-Satsuma equations in extended form are given as examples within the framework.

On some classes of breather lattice solutions to the sinh-Gordon equation

International Nuclear Information System (INIS)

Fu Zuntao; Liu Shikuo

2007-01-01

In this paper, dependent and independent variable transformations are introduced to solve the sinh-Gordon equation by using the knowledge of the elliptic equation and Jacobian elliptic functions. It is shown that different kinds of solutions can be obtained to the sinh-Gordon equation, including breather lattice solutions and periodic wave solutions. (orig.)

Lie Symmetry Analysis of the Inhomogeneous Toda Lattice Equation via Semi-Discrete Exterior Calculus

International Nuclear Information System (INIS)

Liu Jiang; Wang Deng-Shan; Yin Yan-Bin

2017-01-01

In this work, the Lie point symmetries of the inhomogeneous Toda lattice equation are obtained by semi-discrete exterior calculus, which is a semi-discrete version of Harrison and Estabrook’s geometric approach. A four-dimensional Lie algebra and its one-, two- and three-dimensional subalgebras are given. Two similarity reductions of the inhomogeneous Toda lattice equation are obtained by using the symmetry vectors. (paper)

A Greenian approach to the solution of the Schroedinger equation for periodic lattice potentials

International Nuclear Information System (INIS)

Minelli, T.A.

1976-01-01

A modified structural Green's function (MSGF), exploiting all the information contained in the previously solved Schroedinger equation for the electron interacting with a single lattice site, has been introduced and used in order to obtain, from a Dyson-type equation, a kernel whose poles and residues give the E-vs.-k relation and, respectively, the Bloch functions. Such a formulation suggests an alternative technique for the approximate solution of the KKR equations. The MSGF formalism has been also used in order to determine the structure constants of a one-dimensional lattice in a general representation

Lattice quantum phase space and Yang-Baxter equation

International Nuclear Information System (INIS)

Djemai, A.E.F.

1995-04-01

In this work, we show that it is possible to construct the quantum group which preserves the quantum symplectic structure introduced in the context of the matrix Hamiltonian formalism. We also study the braiding existing behind the lattice quantum phase space, and present another type of non-trivial solution to the resulting Yang-Baxter equation. (author). 20 refs, 1 fig

On Traveling Waves in Lattices: The Case of Riccati Lattices

Science.gov (United States)

Dimitrova, Zlatinka

2012-09-01

The method of simplest equation is applied for analysis of a class of lattices described by differential-difference equations that admit traveling-wave solutions constructed on the basis of the solution of the Riccati equation. We denote such lattices as Riccati lattices. We search for Riccati lattices within two classes of lattices: generalized Lotka-Volterra lattices and generalized Holling lattices. We show that from the class of generalized Lotka-Volterra lattices only the Wadati lattice belongs to the class of Riccati lattices. Opposite to this many lattices from the Holling class are Riccati lattices. We construct exact traveling wave solutions on the basis of the solution of Riccati equation for three members of the class of generalized Holling lattices.

Critique of the Brownian approximation to the generalized Langevin equation in lattice dynamics

International Nuclear Information System (INIS)

Diestler, D.J.; Riley, M.E.

1985-01-01

We consider the classical motion of a harmonic lattice in which only those atoms in a certain subset of the lattice (primary zone) may interact with an external force. The formally exact generalized Langevin equation (GLE) for the primary zone is an appropriate description of the dynamics. We examine a previously proposed Brownian, or frictional damping, approximation that reduces the GLE to a set of coupled ordinary Langevin equations for the primary atoms. It is shown that the solution of these equations can contain undamped motion if there is more than one atom in the primary zone. Such motion is explicitly demonstrated for a model that has been used to describe energy transfer in atom--surface collisions. The inability of the standard Brownian approximation to yield an acceptable, physically meaningful result for primary zones comprising more than one atom suggests that the Brownian approximation may introduce other spurious dynamical effects. Further work on damping of correlated motion in lattices is needed

Handbook of integral equations

CERN Document Server

Polyanin, Andrei D

2008-01-01

This handbook contains over 2,500 integral equations with solutions as well as analytical and numerical methods for solving linear and nonlinear equations. It explores Volterra, Fredholm, WienerHopf, Hammerstein, Uryson, and other equations that arise in mathematics, physics, engineering, the sciences, and economics. This second edition includes new chapters on mixed multidimensional equations and methods of integral equations for ODEs and PDEs, along with over 400 new equations with exact solutions. With many examples added for illustrative purposes, it presents new material on Volterra, Fredholm, singular, hypersingular, dual, and nonlinear integral equations, integral transforms, and special functions.

The Asymptotic Expansion of Lattice Loop Integrals Around the Continuum Limit

International Nuclear Information System (INIS)

Becher, Thomas G

2002-01-01

We present a method of computing any one-loop integral in lattice perturbation theory by systematically expanding around its continuum limit. At any order in the expansion in the lattice spacing, the result can be written as a sum of continuum loop integrals in analytic regularization and a few genuine lattice integrals (''master integrals''). These lattice master integrals are independent of external momenta and masses and can be computed numerically. At the one-loop level, there are four master integrals in a theory with only bosonic fields, seven in HQET and sixteen in QED or QCD with Wilson fermions

A two-dimensional lattice equation as an extension of the Heideman-Hogan recurrence

Science.gov (United States)

Kamiya, Ryo; Kanki, Masataka; Mase, Takafumi; Tokihiro, Tetsuji

2018-03-01

We consider a two dimensional extension of the so-called linearizable mappings. In particular, we start from the Heideman-Hogan recurrence, which is known as one of the linearizable Somos-like recurrences, and introduce one of its two dimensional extensions. The two dimensional lattice equation we present is linearizable in both directions, and has the Laurent and the coprimeness properties. Moreover, its reduction produces a generalized family of the Heideman-Hogan recurrence. Higher order examples of two dimensional linearizable lattice equations related to the Dana Scott recurrence are also discussed.

An Integrable Approximation for the Fermi Pasta Ulam Lattice

Science.gov (United States)

Rink, Bob

This contribution presents a review of results obtained from computations of approximate equations of motion for the Fermi-Pasta-Ulam lattice. These approximate equations are obtained as a finite-dimensional Birkhoff normal form. It turns out that in many cases, the Birkhoff normal form is suitable for application of the KAM theorem. In particular, this proves Nishida's 1971 conjecture stating that almost all low-energetic motions of the anharmonic Fermi-Pasta-Ulam lattice with fixed endpoints are quasi-periodic. The proof is based on the formal Birkhoff normal form computations of Nishida, the KAM theorem and discrete symmetry considerations.

Factorization method for difference equations of hypergeometric type on nonuniform lattices

Energy Technology Data Exchange (ETDEWEB)

Alvarez-Nodarse, R. [Departamento de Analisis Matematico, Universidad de Sevilla, Sevilla (Spain); Instituto Carlos I de Fisica Teorica y Computacional, Universidad de Granada, Granada (Spain); Costas-Santos, R.S. [Departamento de Analisis Matematico, Universidad de Sevilla, Sevilla (Spain)

2001-07-13

We study the factorization of the hypergeometric-type difference equation of Nikiforov and Uvarov on nonuniform lattices. An explicit form of the raising and lowering operators is derived and some relevant examples are given. (author)

Determination of space-energy distribution of resonance neutrons in reactor lattice cell and calculation of resonance integrals

International Nuclear Information System (INIS)

Zmijarevic, I.

1980-01-01

Space-energy distribution of resonance neutrons in reactor lattice cell was determined by solving the Boltzmann equation by spherical harmonics method applying P-3 approximation. Computer code SPLET used for these calculations is described. Resonance absorption and calculation of resonance integrals are described as well. Effective resonance integral values for U-238 resonance at 6.7 Ev are calculated for heavy water reactor cell with metal, oxide and carbide fuel elements

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...

A lattice Boltzmann model for the Burgers-Fisher equation.

Science.gov (United States)

Zhang, Jianying; Yan, Guangwu

2010-06-01

A lattice Boltzmann model is developed for the one- and two-dimensional Burgers-Fisher equation based on the method of the higher-order moment of equilibrium distribution functions and a series of partial differential equations in different time scales. In order to obtain the two-dimensional Burgers-Fisher equation, vector sigma(j) has been used. And in order to overcome the drawbacks of "error rebound," a new assumption of additional distribution is presented, where two additional terms, in first order and second order separately, are used. Comparisons with the results obtained by other methods reveal that the numerical solutions obtained by the proposed method converge to exact solutions. The model under new assumption gives better results than that with second order assumption. (c) 2010 American Institute of Physics.

Integration of Chandrasekhar's integral equation

International Nuclear Information System (INIS)

Tanaka, Tasuku

2003-01-01

We solve Chandrasekhar's integration equation for radiative transfer in the plane-parallel atmosphere by iterative integration. The primary thrust in radiative transfer has been to solve the forward problem, i.e., to evaluate the radiance, given the optical thickness and the scattering phase function. In the area of satellite remote sensing, our problem is the inverse problem: to retrieve the surface reflectance and the optical thickness of the atmosphere from the radiance measured by satellites. In order to retrieve the optical thickness and the surface reflectance from the radiance at the top-of-the atmosphere (TOA), we should express the radiance at TOA 'explicitly' in the optical thickness and the surface reflectance. Chandrasekhar formalized radiative transfer in the plane-parallel atmosphere in a simultaneous integral equation, and he obtained the second approximation. Since then no higher approximation has been reported. In this paper, we obtain the third approximation of the scattering function. We integrate functions derived from the second approximation in the integral interval from 1 to ∞ of the inverse of the cos of zenith angles. We can obtain the indefinite integral rather easily in the form of a series expansion. However, the integrals at the upper limit, ∞, are not yet known to us. We can assess the converged values of those series expansions at ∞ through calculus. For integration, we choose coupling pairs to avoid unnecessary terms in the outcome of integral and discover that the simultaneous integral equation can be deduced to the mere integral equation. Through algebraic calculation, we obtain the third approximation as a polynomial of the third degree in the atmospheric optical thickness

Near-Integrability of Low-Dimensional Periodic Klein-Gordon Lattices

Directory of Open Access Journals (Sweden)

Ognyan Christov

2018-01-01

Full Text Available The low-dimensional periodic Klein-Gordon lattices are studied for integrability. We prove that the periodic lattice with two particles and certain nonlinear potential is nonintegrable. However, in the cases of up to six particles, we prove that their Birkhoff-Gustavson normal forms are integrable, which allows us to apply KAM theory in most cases.

Z flux-line lattices and self-dual equations in the standard model

International Nuclear Information System (INIS)

Bimonte, G.; Lozano, G.

1994-04-01

We derive gauge covariant self-dual equations for the SU(2) x U(1) y theory of electroweak interactions and show that they admit solutions describing a periodic lattice of Z-strings. (author). 14 refs

Lattice quantum chromodynamics equation of state: A better ...

Indian Academy of Sciences (India)

Lattice gauge theory; quantum chromodynamics; finite temperature field theory. ... to a previously underappreciated feature of the plasma phase – that it is far from being a ... setting P = 0 just below Tc and the numerical integration errors. ...... for different temperatures, both above and below Tc. We draw attention to the.

Rational solutions of the discrete time Toda lattice and the alternate discrete Painleve II equation

International Nuclear Information System (INIS)

Common, Alan K; Hone, Andrew N W

2008-01-01

The Yablonskii-Vorob'ev polynomials y n (t), which are defined by a second-order bilinear differential-difference equation, provide rational solutions of the Toda lattice. They are also polynomial tau-functions for the rational solutions of the second Painleve equation (P II ). Here we define two-variable polynomials Y n (t, h) on a lattice with spacing h, by considering rational solutions of the discrete time Toda lattice as introduced by Suris. These polynomials are shown to have many properties that are analogous to those of the Yablonskii-Vorob'ev polynomials, to which they reduce when h = 0. They also provide rational solutions for a particular discretization of P II , namely the so-called alternate discrete P II , and this connection leads to an expression in terms of the Umemura polynomials for the third Painleve equation (P III ). It is shown that the Baecklund transformation for the alternate discrete Painleve equation is a symplectic map, and the shift in time is also symplectic. Finally we present a Lax pair for the alternate discrete P II , which recovers Jimbo and Miwa's Lax pair for P II in the continuum limit h → 0

Mobile localization in nonlinear Schroedinger lattices

International Nuclear Information System (INIS)

Gomez-Gardenes, J.; Falo, F.; Floria, L.M.

2004-01-01

Using continuation methods from the integrable Ablowitz-Ladik lattice, we have studied the structure of numerically exact mobile discrete breathers in the standard discrete nonlinear Schroedinger equation. We show that, away from that integrable limit, the mobile pulse is dressed by a background of resonant plane waves with wavevectors given by a certain selection rule. This background is seen to be essential for supporting mobile localization in the absence of integrability. We show how the variations of the localized pulse energy during its motion are balanced by the interaction with this background, allowing the localization mobility along the lattice

Lie Symmetry Analysis of the Inhomogeneous Toda Lattice Equation via Semi-Discrete Exterior Calculus

Science.gov (United States)

Liu, Jiang; Wang, Deng-Shan; Yin, Yan-Bin

2017-06-01

In this work, the Lie point symmetries of the inhomogeneous Toda lattice equation are obtained by semi-discrete exterior calculus, which is a semi-discrete version of Harrison and Estabrook’s geometric approach. A four-dimensional Lie algebra and its one-, two- and three-dimensional subalgebras are given. Two similarity reductions of the inhomogeneous Toda lattice equation are obtained by using the symmetry vectors. Supported by National Natural Science Foundation of China under Grant Nos. 11375030, 11472315, and Department of Science and Technology of Henan Province under Grant No. 162300410223 and Beijing Finance Funds of Natural Science Program for Excellent Talents under Grant No. 2014000026833ZK19

Neutron streaming in fast reactors: Supplementary considerations concerning the perpendicular leakage in slab lattices and detailed equations for the cylindrical channel leakage

International Nuclear Information System (INIS)

Fischer, E.A.

1981-01-01

Basic problems of the definition and calculation of the perpendicular cell diffusion coefficient in slab lattices of fast critical assemblies are discussed in terms of an integral transport formalism. Numerical estimates of correction terms are given which are usually neglected. In addition, a detailed derivation of an equation for the diffusion-coefficient in a cylindrical low-density channel is given

Generalized modification in the lattice Bhatnagar-Gross-Krook model for incompressible Navier-Stokes equations and convection-diffusion equations.

Science.gov (United States)

Yang, Xuguang; Shi, Baochang; Chai, Zhenhua

2014-07-01

In this paper, two modified lattice Boltzmann Bhatnagar-Gross-Krook (LBGK) models for incompressible Navier-Stokes equations and convection-diffusion equations are proposed via the addition of correction terms in the evolution equations. Utilizing this modification, the value of the dimensionless relaxation time in the LBGK model can be kept in a proper range, and thus the stability of the LBGK model can be improved. Although some gradient operators are included in the correction terms, they can be computed efficiently using local computational schemes such that the present LBGK models still retain the intrinsic parallelism characteristic of the lattice Boltzmann method. Numerical studies of the steady Poiseuille flow and unsteady Womersley flow show that the modified LBGK model has a second-order convergence rate in space, and the compressibility effect in the common LBGK model can be eliminated. In addition, to test the stability of the present models, we also performed some simulations of the natural convection in a square cavity, and we found that the results agree well with those reported in the previous work, even at a very high Rayleigh number (Ra = 10(12)).

Integral equations and their applications

CERN Document Server

Rahman, M

2007-01-01

For many years, the subject of functional equations has held a prominent place in the attention of mathematicians. In more recent years this attention has been directed to a particular kind of functional equation, an integral equation, wherein the unknown function occurs under the integral sign. The study of this kind of equation is sometimes referred to as the inversion of a definite integral. While scientists and engineers can already choose from a number of books on integral equations, this new book encompasses recent developments including some preliminary backgrounds of formulations of integral equations governing the physical situation of the problems. It also contains elegant analytical and numerical methods, and an important topic of the variational principles. Primarily intended for senior undergraduate students and first year postgraduate students of engineering and science courses, students of mathematical and physical sciences will also find many sections of direct relevance. The book contains eig...

Integrability and Poisson Structures of Three Dimensional Dynamical Systems and Equations of Hydrodynamic Type

Science.gov (United States)

Gumral, Hasan

Poisson structure of completely integrable 3 dimensional dynamical systems can be defined in terms of an integrable 1-form. We take advantage of this fact and use the theory of foliations in discussing the geometrical structure underlying complete and partial integrability. We show that the Halphen system can be formulated in terms of a flat SL(2,R)-valued connection and belongs to a non-trivial Godbillon-Vey class. On the other hand, for the Euler top and a special case of 3-species Lotka-Volterra equations which are contained in the Halphen system as limiting cases, this structure degenerates into the form of globally integrable bi-Hamiltonian structures. The globally integrable bi-Hamiltonian case is a linear and the sl_2 structure is a quadratic unfolding of an integrable 1-form in 3 + 1 dimensions. We complete the discussion of the Hamiltonian structure of 2-component equations of hydrodynamic type by presenting the Hamiltonian operators for Euler's equation and a continuum limit of Toda lattice. We present further infinite sequences of conserved quantities for shallow water equations and show that their generalizations by Kodama admit bi-Hamiltonian structure. We present a simple way of constructing the second Hamiltonian operators for N-component equations admitting some scaling properties. The Kodama reduction of the dispersionless-Boussinesq equations and the Lax reduction of the Benney moment equations are shown to be equivalent by a symmetry transformation. They can be cast into the form of a triplet of conservation laws which enable us to recognize a non-trivial scaling symmetry. The resulting bi-Hamiltonian structure generates three infinite sequences of conserved densities.

Dynamics and elastic interactions of the discrete multi-dark soliton solutions for the Kaup-Newell lattice equation

Science.gov (United States)

Liu, Nan; Wen, Xiao-Yong

2018-03-01

Under consideration in this paper is the Kaup-Newell (KN) lattice equation which is an integrable discretization of the KN equation. Infinitely, many conservation laws and discrete N-fold Darboux transformation (DT) for this system are constructed and established based on its Lax representation. Via the resulting N-fold DT, the discrete multi-dark soliton solutions in terms of determinants are derived from non-vanishing background. Propagation and elastic interaction structures of such solitons are shown graphically. Overtaking interaction phenomena between/among the two, three and four solitons are discussed. Numerical simulations are used to explore their dynamical behaviors of such multi-dark solitons. Numerical results show that their evolutions are stable against a small noise. Results in this paper might be helpful for understanding the propagation of nonlinear Alfvén waves in plasmas.

Remarks on an equation common to Weyl's gauge field, Yang-Mills field and Toda lattice

International Nuclear Information System (INIS)

Nishioka, M.

1984-01-01

In this letter a remark is presented on an equation of a gauge-invariant Weyl's gauge field and it is shown that the equation is common to Yang's approach to the self-duality condition for SU 2 gauge field and the simplest Toda lattice

Elliptic Painlevé equations from next-nearest-neighbor translations on the E_8^{(1)} lattice

Science.gov (United States)

Joshi, Nalini; Nakazono, Nobutaka

2017-07-01

The well known elliptic discrete Painlevé equation of Sakai is constructed by a standard translation on the E_8(1) lattice, given by nearest neighbor vectors. In this paper, we give a new elliptic discrete Painlevé equation obtained by translations along next-nearest-neighbor vectors. This equation is a generic (8-parameter) version of a 2-parameter elliptic difference equation found by reduction from Adler’s partial difference equation, the so-called Q4 equation. We also provide a projective reduction of the well known equation of Sakai.

The anisotropic Ising correlations as elliptic integrals: duality and differential equations

International Nuclear Information System (INIS)

McCoy, B M; Maillard, J-M

2016-01-01

We present the reduction of the correlation functions of the Ising model on the anisotropic square lattice to complete elliptic integrals of the first, second and third kind, the extension of Kramers–Wannier duality to anisotropic correlation functions, and the linear differential equations for these anisotropic correlations. More precisely, we show that the anisotropic correlation functions are homogeneous polynomials of the complete elliptic integrals of the first, second and third kind. We give the exact dual transformation matching the correlation functions and the dual correlation functions. We show that the linear differential operators annihilating the general two-point correlation functions are factorized in a very simple way, in operators of decreasing orders. (paper)

A new lattice Boltzmann equation to simulate density-driven convection of carbon dioxide

KAUST Repository

Allen, Rebecca

2013-01-01

The storage of CO2 in fluid-filled geological formations has been carried out for more than a decade in locations around the world. After CO2 has been injected into the aquifer and has moved laterally under the aquifer\\'s cap-rock, density-driven convection becomes an important transport process to model. However, the challenge lies in simulating this transport process accurately with high spatial resolution and low CPU cost. This issue can be addressed by using the lattice Boltzmann equation (LBE) to formulate a model for a similar scenario when a solute diffuses into a fluid and density differences lead to convective mixing. The LBE is a promising alternative to the traditional methods of computational fluid dynamics. Rather than discretizing the system of partial differential equations of classical continuum mechanics directly, the LBE is derived from a velocity-space truncation of the Boltzmann equation of classical kinetic theory. We propose an extension to the LBE, which can accurately predict the transport of dissolved CO2 in water, as a step towards fluid-filled porous media simulations. This is achieved by coupling two LBEs, one for the fluid flow and one for the convection and diffusion of CO2. Unlike existing lattice Boltzmann equations for porous media flow, our model is derived from a system of moment equations and a Crank-Nicolson discretization of the velocity-truncated Boltzmann equation. The forcing terms are updated locally without the need for additional central difference approximation. Therefore our model preserves all the computational advantages of the single-phase lattice Boltzmann equation and is formally second-order accurate in both space and time. Our new model also features a novel implementation of boundary conditions, which is simple to implement and does not suffer from the grid-dependent error that is present in the standard "bounce-back" condition. The significance of using the LBE in this work lies in the ability to efficiently

A lattice hierarchy and its continuous limits

International Nuclear Information System (INIS)

Fan Engui

2008-01-01

By introducing a discrete spectral problem, we derive a lattice hierarchy which is integrable in Liouville's sense and possesses a multi-Hamiltonian structure. It is show that the discrete spectral problem converges to the well-known AKNS spectral problem under a certain continuous limit. In particular, we construct a sequence of equations in the lattice hierarchy which approximates the AKNS hierarchy as a continuous limit

A new lattice Boltzmann equation to simulate density-driven convection of carbon dioxide

KAUST Repository

Allen, Rebecca; Reis, Tim; Sun, Shuyu

2013-01-01

-driven convection becomes an important transport process to model. However, the challenge lies in simulating this transport process accurately with high spatial resolution and low CPU cost. This issue can be addressed by using the lattice Boltzmann equation (LBE

Integral equation for Coulomb problem

International Nuclear Information System (INIS)

Sasakawa, T.

1986-01-01

For short range potentials an inhomogeneous (homogeneous) Lippmann-Schwinger integral equation of the Fredholm type yields the wave function of scattering (bound) state. For the Coulomb potential, this statement is no more valid. It has been felt difficult to express the Coulomb wave function in a form of an integral equation with the Coulomb potential as the perturbation. In the present paper, the author shows that an inhomogeneous integral equation of a Volterra type with the Coulomb potential as the perturbation can be constructed both for the scattering and the bound states. The equation yielding the binding energy is given in an integral form. The present treatment is easily extended to the coupled Coulomb problems

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.

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.

Analytical solution to the hybrid diffusion-transport equation

International Nuclear Information System (INIS)

Nanneh, M.M.; Williams, M.M.R.

1986-01-01

A special integral equation was derived in previous work using a hybrid diffusion-transport theory method for calculating the flux distribution in slab lattices. In this paper an analytical solution of this equation has been carried out on a finite reactor lattice. The analytical results of disadvantage factors are shown to be accurate in comparison with the numerical results and accurate transport theory calculations. (author)

On integrability of the Killing equation

Science.gov (United States)

Houri, Tsuyoshi; Tomoda, Kentaro; Yasui, Yukinori

2018-04-01

Killing tensor fields have been thought of as describing the hidden symmetry of space(-time) since they are in one-to-one correspondence with polynomial first integrals of geodesic equations. Since many problems in classical mechanics can be formulated as geodesic problems in curved space and spacetime, solving the defining equation for Killing tensor fields (the Killing equation) is a powerful way to integrate equations of motion. Thus it has been desirable to formulate the integrability conditions of the Killing equation, which serve to determine the number of linearly independent solutions and also to restrict the possible forms of solutions tightly. In this paper, we show the prolongation for the Killing equation in a manner that uses Young symmetrizers. Using the prolonged equations, we provide the integrability conditions explicitly.

Loop equations and bootstrap methods in the lattice

Directory of Open Access Journals (Sweden)

Peter D. Anderson

2017-08-01

Full Text Available Pure gauge theories can be formulated in terms of Wilson Loops by means of the loop equation. In the large-N limit this equation closes in the expectation value of single loops. In particular, using the lattice as a regulator, it becomes a well defined equation for a discrete set of loops. In this paper we study different numerical approaches to solving this equation. Previous ideas gave good results in the strong coupling region. Here we propose an alternative method based on the observation that certain matrices ρˆ of Wilson loop expectation values are positive definite. They also have unit trace (ρˆ⪰0,Trρˆ=1, in fact they can be defined as reduced density matrices in the space of open loops after tracing over color indices and can be used to define an entropy associated with the loss of information due to such trace SWL=−Tr[ρˆln⁡ρˆ]. The condition that such matrices are positive definite allows us to study the weak coupling region which is relevant for the continuum limit. In the exactly solvable case of two dimensions this approach gives very good results by considering just a few loops. In four dimensions it gives good results in the weak coupling region and therefore is complementary to the strong coupling expansion. We compare the results with standard Monte Carlo simulations.

An Implicit Scheme of Lattice Boltzmann Method for Sine-Gordon Equation

International Nuclear Information System (INIS)

Hui-Lin, Lai; Chang-Feng, Ma

2008-01-01

We establish an implicit scheme of lattice Boltzmann method for simulating the sine-Gordon equation, which can be transformed into the explicit one, so the computation of the scheme is simple. Moreover, the parameter θ of the implicit scheme is independent of the relaxation time, which makes the model more flexible. The numerical results show that this method is very effective. (fundamental areas of phenomenology (including applications))

Linear integral equations and soliton systems

International Nuclear Information System (INIS)

Quispel, G.R.W.

1983-01-01

A study is presented of classical integrable dynamical systems in one temporal and one spatial dimension. The direct linearizations are given of several nonlinear partial differential equations, for example the Korteweg-de Vries equation, the modified Korteweg-de Vries equation, the sine-Gordon equation, the nonlinear Schroedinger equation, and the equation of motion for the isotropic Heisenberg spin chain; the author also discusses several relations between these equations. The Baecklund transformations of these partial differential equations are treated on the basis of a singular transformation of the measure (or equivalently of the plane-wave factor) occurring in the corresponding linear integral equations, and the Baecklund transformations are used to derive the direct linearization of a chain of so-called modified partial differential equations. Finally it is shown that the singular linear integral equations lead in a natural way to the direct linearizations of various nonlinear difference-difference equations. (Auth.)

Geophysical interpretation using integral equations

CERN Document Server

Eskola, L

1992-01-01

Along with the general development of numerical methods in pure and applied to apply integral equations to geophysical modelling has sciences, the ability improved considerably within the last thirty years or so. This is due to the successful derivation of integral equations that are applicable to the modelling of complex structures, and efficient numerical algorithms for their solution. A significant stimulus for this development has been the advent of fast digital computers. The purpose of this book is to give an idea of the principles by which boundary-value problems describing geophysical models can be converted into integral equations. The end results are the integral formulas and integral equations that form the theoretical framework for practical applications. The details of mathematical analysis have been kept to a minimum. Numerical algorithms are discussed only in connection with some illustrative examples involving well-documented numerical modelling results. The reader is assu­ med to have a back...

Integral equation methods for electromagnetics

CERN Document Server

Volakis, John

2012-01-01

This text/reference is a detailed look at the development and use of integral equation methods for electromagnetic analysis, specifically for antennas and radar scattering. Developers and practitioners will appreciate the broad-based approach to understanding and utilizing integral equation methods and the unique coverage of historical developments that led to the current state-of-the-art. In contrast to existing books, Integral Equation Methods for Electromagnetics lays the groundwork in the initial chapters so students and basic users can solve simple problems and work their way up to the mo

Isotropic quantum walks on lattices and the Weyl equation

Science.gov (United States)

D'Ariano, Giacomo Mauro; Erba, Marco; Perinotti, Paolo

2017-12-01

We present a thorough classification of the isotropic quantum walks on lattices of dimension d =1 ,2 ,3 with a coin system of dimension s =2 . For d =3 there exist two isotropic walks, namely, the Weyl quantum walks presented in the work of D'Ariano and Perinotti [G. M. D'Ariano and P. Perinotti, Phys. Rev. A 90, 062106 (2014), 10.1103/PhysRevA.90.062106], resulting in the derivation of the Weyl equation from informational principles. The present analysis, via a crucial use of isotropy, is significantly shorter and avoids a superfluous technical assumption, making the result completely general.

Gauge theories and integrable lattice models

International Nuclear Information System (INIS)

Witten, E.

1989-01-01

Investigations of new knot polynomials discovered in the last few years have shown them to be intimately connected with soluble models of two dimensional lattice statistical mechanics. In this paper, these results, which in time may illuminate the whole question of why integrable lattice models exist, are reconsidered from the point of view of three dimensional gauge theory. Expectation values of Wilson lines in three dimensional Chern-Simons gauge theories can be computed by evaluating the partition functions of certain lattice models on finite graphs obtained by projecting the Wilson lines to the plane. The models in question - previously considered in both the knot theory and statistical mechanics literature - are IRF models in which the local Boltzmann weights are the matrix elements of braiding matrices in rational conformal field theories. These matrix elements, in turn, can be represented in three dimensional gauge theory in terms of the expectation value of a certain tetrahedral configuration of Wilson lines. This representation makes manifest a surprising symmetry of the braiding matrix elements in conformal field theory. (orig.)

Nonlinear integrodifferential equations as discrete systems

Science.gov (United States)

Tamizhmani, K. M.; Satsuma, J.; Grammaticos, B.; Ramani, A.

1999-06-01

We analyse a class of integrodifferential equations of the `intermediate long wave' (ILW) type. We show that these equations can be formally interpreted as discrete, differential-difference systems. This allows us to link equations of this type with previous results of ours involving differential-delay equations and, on the basis of this, propose new integrable equations of ILW type. Finally, we extend this approach to pure difference equations and propose ILW forms for the discrete lattice KdV equation.

Differential-difference equations associated with the fractional Lax operators

Energy Technology Data Exchange (ETDEWEB)

Adler, V E [LD Landau Institute for Theoretical Physics, 1A Ak. Semenov, Chernogolovka 142432 (Russian Federation); Postnikov, V V, E-mail: adler@itp.ac.ru, E-mail: postnikofvv@mail.ru [Sochi Branch of Peoples' Friendship University of Russia, 32 Kuibyshev str., 354000 Sochi (Russian Federation)

2011-10-14

We study integrable hierarchies associated with spectral problems of the form P{psi} = {lambda}Q{psi}, where P and Q are difference operators. The corresponding nonlinear differential-difference equations can be viewed as inhomogeneous generalizations of the Bogoyavlensky-type lattices. While the latter turn into the Korteweg-de Vries equation under the continuous limit, the lattices under consideration provide discrete analogs of the Sawada-Kotera and Kaup-Kupershmidt equations. The r-matrix formulation and several of the simplest explicit solutions are presented. (paper)

PREFACE: Symmetries and integrability of difference equations Symmetries and integrability of difference equations

Science.gov (United States)

Levi, Decio; Olver, Peter; Thomova, Zora; Winternitz, Pavel

2009-11-01

The concept of integrability was introduced in classical mechanics in the 19th century for finite dimensional continuous Hamiltonian systems. It was extended to certain classes of nonlinear differential equations in the second half of the 20th century with the discovery of the inverse scattering transform and the birth of soliton theory. Also at the end of the 19th century Lie group theory was invented as a powerful tool for obtaining exact analytical solutions of large classes of differential equations. Together, Lie group theory and integrability theory in its most general sense provide the main tools for solving nonlinear differential equations. Like differential equations, difference equations play an important role in physics and other sciences. They occur very naturally in the description of phenomena that are genuinely discrete. Indeed, they may actually be more fundamental than differential equations if space-time is actually discrete at very short distances. On the other hand, even when treating continuous phenomena described by differential equations it is very often necessary to resort to numerical methods. This involves a discretization of the differential equation, i.e. a replacement of the differential equation by a difference one. Given the well developed and understood techniques of symmetry and integrability for differential equations a natural question to ask is whether it is possible to develop similar techniques for difference equations. The aim is, on one hand, to obtain powerful methods for solving `integrable' difference equations and to establish practical integrability criteria, telling us when the methods are applicable. On the other hand, Lie group methods can be adapted to solve difference equations analytically. Finally, integrability and symmetry methods can be combined with numerical methods to obtain improved numerical solutions of differential equations. The origin of the SIDE meetings goes back to the early 1990s and the first

On vector analogs of the modified Volterra lattice

Energy Technology Data Exchange (ETDEWEB)

Adler, V E; Postnikov, V V [L D Landau Institute for Theoretical Physics, 1a Semenov pr, 142432 Chernogolovka (Russian Federation); Sochi Branch of Peoples' Friendship University of Russia, 32 Kuibyshev str, 354000 Sochi (Russian Federation)], E-mail: adler@itp.ac.ru, E-mail: postnikovvv@rambler.ru

2008-11-14

The zero curvature representations, Baecklund transformations, nonlinear superposition principle and the simplest explicit solutions of soliton and breather type are presented for two vector generalizations of modified Volterra lattice. The relations with some other integrable equations are established.

Completely integrable operator evolution equations. II

International Nuclear Information System (INIS)

Chudnovsky, D.V.

1979-01-01

The author continues the investigation of operator classical completely integrable systems. The main attention is devoted to the stationary operator non-linear Schroedinger equation. It is shown that this equation can be used for separation of variables for a large class of completely integrable equations. (Auth.)

Mode instability in one-dimensional anharmonic lattices: Variational equation approach

Science.gov (United States)

Yoshimura, K.

1999-03-01

The stability of normal mode oscillations has been studied in detail under the single-mode excitation condition for the Fermi-Pasta-Ulam-β lattice. Numerical experiments indicate that the mode stability depends strongly on k/N, where k is the wave number of the initially excited mode and N is the number of degrees of freedom in the system. It has been found that this feature does not change when N increases. We propose an average variational equation - approximate version of the variational equation - as a theoretical tool to facilitate a linear stability analysis. It is shown that this strong k/N dependence of the mode stability can be explained from the view point of the linear stability of the relevant orbits. We introduce a low-dimensional approximation of the average variational equation, which approximately describes the time evolution of variations in four normal mode amplitudes. The linear stability analysis based on this four-mode approximation demonstrates that the parametric instability mechanism plays a crucial role in the strong k/N dependence of the mode stability.

Resonance shielding in thermal reactor lattices

International Nuclear Information System (INIS)

Rothenstein, W.; Taviv, E.; Aminpour, M.

1982-01-01

The theoretical foundations of a new methodology for the accurate treatment of resonance absorption in thermal reactor lattice analysis are presented. This methodology is based on the solution of the point-energy transport equation in its integral or integro-differential form for a heterogeneous lattice using detailed resonance cross-section profiles. The methodology is applied to LWR benchmark analysis, with emphasis on temperature dependence of resonance absorption during fuel depletion, spatial and mutual self-shielding, integral parameter analysis and treatment of cluster geometry. The capabilities of the OZMA code, which implements the new methodology are discussed. These capabilities provide a means against which simpler and more rapid resonance absorption algorithms can be checked. (author)

A lattice Boltzmann model with an amending function for simulating nonlinear partial differential equations

International Nuclear Information System (INIS)

Lin-Jie, Chen; Chang-Feng, Ma

2010-01-01

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form u t + αuu x + βu n u x + γu xx + δu xxx + ζu xxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman–Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions. (general)

Calculation of the critical buckling of a lattice based on the integral form of the transport equation

International Nuclear Information System (INIS)

Benoist, P.

1990-06-01

The migration area, which relates the buckling to the multiplication factor, can be calculated by means of the Deniz formula. This formula involves the direct and adjoint angular fluxes. It is shown in this note that it is possible, using the integral form of the transport equation, to establish an equivalent formula in which only angle-integrated quantities appear. This formulation is more suitable for the calculation by the collision probably method [fr

Integrable Seven-Point Discrete Equations and Second-Order Evolution Chains

Science.gov (United States)

Adler, V. E.

2018-04-01

We consider differential-difference equations defining continuous symmetries for discrete equations on a triangular lattice. We show that a certain combination of continuous flows can be represented as a secondorder scalar evolution chain. We illustrate the general construction with a set of examples including an analogue of the elliptic Yamilov chain.

3D Relativistic Hydrodynamic Computations Using Lattice-QCD-Inspired Equations of State

International Nuclear Information System (INIS)

Hama, Yogiro; Andrade, Rone P.G.; Grassi, Frederique; Socolowski, Otavio; Kodama, Takeshi; Tavares, Bernardo; Padula, Sandra S.

2006-01-01

In this communication, we report results of three-dimensional hydrodynamic computations, by using equations of state with a critical end point as suggested by the lattice QCD. Some of the results are an increase of the multiplicity in the mid-rapidity region and a larger elliptic-flow parameter v 2 . We discuss also the effcts of the initial-condition fluctuations and the continuous emission

3D Relativistic Hydrodynamic Computations Using Lattice-QCD-Inspired Equations of State

Energy Technology Data Exchange (ETDEWEB)

Hama, Yogiro [Instituto de Fisica, Universidade de Sao Paulo (Brazil); Andrade, Rone P.G. [Instituto de Fisica, Universidade de Sao Paulo (Brazil); Grassi, Frederique [Instituto de Fisica, Universidade de Sao Paulo (Brazil); Socolowski, Otavio [Instituto Tecnologico da Aeronautica (Brazil); Kodama, Takeshi [Instituto de Fisica, Universidade Federal do Rio de Janeiro (Brazil); Tavares, Bernardo [Instituto de Fisica, Universidade Federal do Rio de Janeiro (Brazil); Padula, Sandra S. [Instituto de Fisica Teorica, Universidade Estadual Paulista (Brazil)

2006-08-07

In this communication, we report results of three-dimensional hydrodynamic computations, by using equations of state with a critical end point as suggested by the lattice QCD. Some of the results are an increase of the multiplicity in the mid-rapidity region and a larger elliptic-flow parameter v{sub 2}. We discuss also the effcts of the initial-condition fluctuations and the continuous emission.

Laguerre-Freud equations for the recurrence coefficients of the Laguerre-Hahn orthogonal polynomials on special nonuniform lattices

International Nuclear Information System (INIS)

Bangerezako, Gaspard; Foupouagnigni, Mama

2003-10-01

We give an algorithmic derivation of the Laguerre-Freud equations for the recurrence coefficients β n and γ n of the Laguerre-Hahn orthogonal polynomials on special nonuniform lattices. This algorithm is the most general one since it is valid for the Laguerre-Hahn orthogonal polynomials of any class k, on the special nonuniform lattices including the continuous (limiting cases), linear, q-linear and the q-nonlinear ones. Moreover, the algorithm allows to deduce an upper bound for the order of the equations in β n and γ n , which is respectively 2 k + 2 and 2 k + 3 when k is even, or 2 k + 3 and 2 k + 2 when k is odd. Finally, as applications, we discuss explicitly these equations for k = 1 in the continuous and linear cases, and k = 2 in the continuous symmetric one. (author)

An Eulerian description of the streaming process in the lattice Boltzmann equation

CERN Document Server

Lee Tae Hun

2003-01-01

This paper presents a novel strategy for solving discrete Boltzmann equation (DBE) for simulation of fluid flows. This strategy splits the solution procedure into streaming and collision steps as in the lattice Boltzmann equation (LBE) method. The streaming step can then be carried out by solving pure linear advection equations in an Eulerian framework. This offers two significant advantages over previous methods. First, the relationship between the relaxation parameter and the discretization of the collision term developed from the LBE method is directly applicable to the DBE method. The resulting DBE collision step remains local and poses no constraint on time step. Second, decoupling of the advection step from the collision step facilitates implicit discretization of the advection equation on arbitrary meshes. An implicit unstructured DBE method is constructed based on this strategy and is evaluated using several test cases of flow over a backward-facing step, lid-driven cavity flow, and flow past a circul...

Integral Equation Methods for Electromagnetic and Elastic Waves

CERN Document Server

Chew, Weng; Hu, Bin

2008-01-01

Integral Equation Methods for Electromagnetic and Elastic Waves is an outgrowth of several years of work. There have been no recent books on integral equation methods. There are books written on integral equations, but either they have been around for a while, or they were written by mathematicians. Much of the knowledge in integral equation methods still resides in journal papers. With this book, important relevant knowledge for integral equations are consolidated in one place and researchers need only read the pertinent chapters in this book to gain important knowledge needed for integral eq

Integrable relativistic Toda type lattice hierarchies, associated coupling systems and the Darboux transformation

International Nuclear Information System (INIS)

Yang Hongxiang; Xu Xixiang; Sun Yepeng; Ding Haiyong

2006-01-01

Starting from a discrete isospectral problem, integrable positive and negative relativistic Toda type lattice hierarchies are derived. The two lattice hierarchies are proven to have discrete zero-curvature representations associated with a discrete spectral problem, and the positive and negative lattice hierarchies correspond to positive and negative power expansions of Lax operators with respect to the spectral parameter, respectively. The integrable positive and negative coupling systems of the resulting hierarchies are constructed through enlarging Lax pairs. In addition, with the help of gauge transformations of spectral problems, a Darboux transformation is established for the relativistic Toda type lattice. As an application, an exact solution is explicitly presented

Numerical simulations of a family of the coupled viscous Burgers, equation using the lattice Boltzmann method

International Nuclear Information System (INIS)

He, Y B; Tang, X H

2016-01-01

In this paper, in order to extend the lattice Boltzmann method (LBM) to deal with more nonlinear systems, a one-dimensional and five-velocity lattice Boltzmann scheme with an amending function for a family of the coupled viscous Burgers’ equation (CVBE) is proposed. With the Taylor and Chapman–Enskog expansion, a family of the CVBE is recovered correctly from the lattice Boltzmann equation through selecting the equilibrium distribution functions and amending functions properly. The method is applied to some test examples with an analytical solution. The results are compared with those obtained by the finite difference method (FDM); it is shown that the numerical solutions agree well with the analytical solutions and the errors obtained by the present method are smaller than the FDM. Furthermore, some problems without analytical solutions are numerically studied by the present method and the FDM. The results show that the numerical solutions of the LBM are in good agreement with those obtained by the FDM, which can validate the effectiveness and stability of the LBM. (paper: classical statistical mechanics, equilibrium and non-equilibrium)

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.

An integral transform of the Salpeter equation

International Nuclear Information System (INIS)

Krolikowski, W.

1980-03-01

We find a new form of relativistic wave equation for two spin-1/2 particles, which arises by an integral transformation (in the position space) of the wave function in the Salpeter equation. The non-locality involved in this transformation is extended practically over the Compton wavelength of the lighter of two particles. In the case of equal masses the new equation assumes the form of the Breit equation with an effective integral interaction. In the one-body limit it reduces to the Dirac equation also with an effective integral interaction. (author)

A 2+1 non-isospectral discrete integrable system and its discrete integrable coupling system

International Nuclear Information System (INIS)

Yu Fajun; Zhang Hongqing

2006-01-01

In this Letter by considering a (2+1)-dimensional discrete non-isospectral linear problem, a new (2+1)-dimensional integrable lattice hierarchy is constructed. It shows that generalization of the Blaszak-Marciniak lattice hierarchy can be obtained as a reduction. Then an extended algebraic system X-bar of X is presented, from which the integrable coupling system of the (2+1)-dimensional discrete non-isospectral Blaszak-Marciniak lattice equations are obtained

Integration rules for scattering equations

International Nuclear Information System (INIS)

Baadsgaard, Christian; Bjerrum-Bohr, N.E.J.; Bourjaily, Jacob L.; Damgaard, Poul H.

2015-01-01

As described by Cachazo, He and Yuan, scattering amplitudes in many quantum field theories can be represented as integrals that are fully localized on solutions to the so-called scattering equations. Because the number of solutions to the scattering equations grows quite rapidly, the contour of integration involves contributions from many isolated components. In this paper, we provide a simple, combinatorial rule that immediately provides the result of integration against the scattering equation constraints for any Möbius-invariant integrand involving only simple poles. These rules have a simple diagrammatic interpretation that makes the evaluation of any such integrand immediate. Finally, we explain how these rules are related to the computation of amplitudes in the field theory limit of string theory.

The IR sector of QCD: lattice versus Schwinger-Dyson equations

International Nuclear Information System (INIS)

Binosi, Daniele

2010-01-01

Important information about the infrared dynamics of QCD is encoded in the behavior of its (of-shell) Green's functions, most notably the gluon and the ghost propagators. Due to recent improvements in the quality of lattice data and the truncation schemes employed for the Schwinger-Dyson equations we have now reached a point where the interplay between these two non-perturbative tools can be most fruitful. In this talk several of the above points will be reviewed, with particular emphasis on the implications for the ghost sector, the non-perturbative effective charge of QCD, and the Kugo-Ojima function.

Completely integrable operator evolutionary equations

International Nuclear Information System (INIS)

Chudnovsky, D.V.

1979-01-01

The authors present natural generalizations of classical completely integrable equations where the functions are replaced by arbitrary operators. Among these equations are the non-linear Schroedinger, the Korteweg-de Vries, and the modified KdV equations. The Lax representation and the Baecklund transformations are presented. (Auth.)

Scattering integral equations and four nucleon problem

International Nuclear Information System (INIS)

Narodetskii, I.M.

1980-01-01

Existing results from the application of integral equation technique to the four-nucleon bound states and scattering are reviewed. The first numerical calculations of the four-body integral equations have been done ten years ago. Yet, it is still widely believed that these equations are too complicated to solve numerically. The purpose of this review is to provide a clear and elementary introduction in the integral equation method and to demonstrate its usefulness in physical applications. The presentation is based on the quasiparticle approach. This permits a simple interpretation of the equations in terms of quasiparticle scattering. The mathematical basis for the quasiparticle approach is the Hilbert-Schmidt method of the Fredholm integral equation theory. The first part of this review contains a detailed discussion of the Hilbert-Schmidt expansion as applied to the 2-particle amplitudes and to the kernel of the four-body equations. The second part contains the discussion of the four-body quasiparticle equations and of the resed forullts obtain bound states and scattering

Hysteresis and Phase Transitions in a Lattice Regularization of an Ill-Posed Forward-Backward Diffusion Equation

Science.gov (United States)

Helmers, Michael; Herrmann, Michael

2018-03-01

We consider a lattice regularization for an ill-posed diffusion equation with a trilinear constitutive law and study the dynamics of phase interfaces in the parabolic scaling limit. Our main result guarantees for a certain class of single-interface initial data that the lattice solutions satisfy asymptotically a free boundary problem with a hysteretic Stefan condition. The key challenge in the proof is to control the microscopic fluctuations that are inevitably produced by the backward diffusion when a particle passes the spinodal region.

Comment on the consistency of truncated nonlinear integral equation based theories of freezing

International Nuclear Information System (INIS)

Cerjan, C.; Bagchi, B.; Rice, S.A.

1985-01-01

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim--Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions

Asymptotic integration of differential and difference equations

CERN Document Server

Bodine, Sigrun

2015-01-01

This book presents the theory of asymptotic integration for both linear differential and difference equations. This type of asymptotic analysis is based on some fundamental principles by Norman Levinson. While he applied them to a special class of differential equations, subsequent work has shown that the same principles lead to asymptotic results for much wider classes of differential and also difference equations. After discussing asymptotic integration in a unified approach, this book studies how the application of these methods provides several new insights and frequent improvements to results found in earlier literature. It then continues with a brief introduction to the relatively new field of asymptotic integration for dynamic equations on time scales. Asymptotic Integration of Differential and Difference Equations is a self-contained and clearly structured presentation of some of the most important results in asymptotic integration and the techniques used in this field. It will appeal to researchers i...

A space-time lattice version of scalar electrodynamics

International Nuclear Information System (INIS)

Kijowski, J.; Thielmann, A.

1993-10-01

A Minkowski-lattice version of quantum scalar electrodynamics is constructed. Quantum field is consequently described in a gauge-independent way, i.e. the algebra of quantum observables of the theory is generated by gauge-invariant operators assigned to zero-, one-, and two-dimensional elements of the lattice. The operators satisfy canonical commutation relations. Field dynamics is formulated in terms of difference equations imposed on the field operators. The dynamics is obtained from a discrete version of the path-integral. (author). 19 refs

Transformation properties of the integrable evolution equations

International Nuclear Information System (INIS)

Konopelchenko, B.G.

1981-01-01

Group-theoretical properties of partial differential equations integrable by the inverse scattering transform method are discussed. It is shown that nonlinear transformations typical to integrable equations (symmetry groups, Baecklund-transformations) and these equations themselves are contained in a certain universal nonlinear transformation group. (orig.)

Application of wavelets to singular integral scattering equations

International Nuclear Information System (INIS)

Kessler, B.M.; Payne, G.L.; Polyzou, W.N.

2004-01-01

The use of orthonormal wavelet basis functions for solving singular integral scattering equations is investigated. It is shown that these basis functions lead to sparse matrix equations which can be solved by iterative techniques. The scaling properties of wavelets are used to derive an efficient method for evaluating the singular integrals. The accuracy and efficiency of the wavelet transforms are demonstrated by solving the two-body T-matrix equation without partial wave projection. The resulting matrix equation which is characteristic of multiparticle integral scattering equations is found to provide an efficient method for obtaining accurate approximate solutions to the integral equation. These results indicate that wavelet transforms may provide a useful tool for studying few-body systems

RICM, Resonance Absorption in Multi-Region Slab or Square or Hexagonal Lattice

International Nuclear Information System (INIS)

Mizuta, H.; Aoyama, K.; Fukai, Y.

1968-01-01

1 - Nature of physical problem solved: Calculates the resonance absorption integral of resonant isotope in a multi-region lattice using the first flight collision probability. The lattice configurations considered are a slab lattice, a square or hexagonal lattice and a cylindricalized lattice with isotropic or perfect reflecting boundary condition. Cases for an isolated rod or plate and homogeneous system can also be treated. 2 - Method of solution: Slowing down of neutrons by each isotope in each region is solved by either exact numerical integration of the slowing down equation or narrow - or wide-resonance approximation. Breit-Wigner's single level formula is used for the resonance cross section and Porter-Thomas distribution of neutron width is taken into account in the unresolved region. 3 - Restrictions on the complexity of the problem: Maximum number of regions: 5; Maximum Number of groups: 100

Elimination of spurious lattice fermion solutions and noncompact lattice QCD

Energy Technology Data Exchange (ETDEWEB)

Lee, T.D.

1997-09-22

It is well known that the Dirac equation on a discrete hyper-cubic lattice in D dimension has 2{sup D} degenerate solutions. The usual method of removing these spurious solutions encounters difficulties with chiral symmetry when the lattice spacing l {ne} 0, as exemplified by the persistent problem of the pion mass. On the other hand, we recall that in any crystal in nature, all the electrons do move in a lattice and satisfy the Dirac equation; yet there is not a single physical result that has ever been entangled with a spurious fermion solution. Therefore it should not be difficult to eliminate these unphysical elements. On a discrete lattice, particle hop from point to point, whereas in a real crystal the lattice structure in embedded in a continuum and electrons move continuously from lattice cell to lattice cell. In a discrete system, the lattice functions are defined only on individual points (or links as in the case of gauge fields). However, in a crystal the electron state vector is represented by the Bloch wave functions which are continuous functions in {rvec {gamma}}, and herein lies one of the essential differences.

Involving the Navier-Stokes equations in the derivation of boundary conditions for the lattice Boltzmann method.

Science.gov (United States)

Verschaeve, Joris C G

2011-06-13

By means of the continuity equation of the incompressible Navier-Stokes equations, additional physical arguments for the derivation of a formulation of the no-slip boundary condition for the lattice Boltzmann method for straight walls at rest are obtained. This leads to a boundary condition that is second-order accurate with respect to the grid spacing and conserves mass. In addition, the boundary condition is stable for relaxation frequencies close to two.

Darboux invariants of integrable equations with variable spectral parameters

International Nuclear Information System (INIS)

Shin, H J

2008-01-01

The Darboux transformation for integrable equations with variable spectral parameters is introduced. Darboux invariant quantities are calculated, which are used in constructing the Lax pair of integrable equations. This approach serves as a systematic method for constructing inhomogeneous integrable equations and their soliton solutions. The structure functions of variable spectral parameters determine the integrability and nonlinear coupling terms. Three cases of integrable equations are treated as examples of this approach

Discrete systems and integrability

CERN Document Server

Hietarinta, J; Nijhoff, F W

2016-01-01

This first introductory text to discrete integrable systems introduces key notions of integrability from the vantage point of discrete systems, also making connections with the continuous theory where relevant. While treating the material at an elementary level, the book also highlights many recent developments. Topics include: Darboux and Bäcklund transformations; difference equations and special functions; multidimensional consistency of integrable lattice equations; associated linear problems (Lax pairs); connections with Padé approximants and convergence algorithms; singularities and geometry; Hirota's bilinear formalism for lattices; intriguing properties of discrete Painlevé equations; and the novel theory of Lagrangian multiforms. The book builds the material in an organic way, emphasizing interconnections between the various approaches, while the exposition is mostly done through explicit computations on key examples. Written by respected experts in the field, the numerous exercises and the thoroug...

Microscopic theory for coupled atomistic magnetization and lattice dynamics

Science.gov (United States)

Fransson, J.; Thonig, D.; Bessarab, P. F.; Bhattacharjee, S.; Hellsvik, J.; Nordström, L.

2017-12-01

A coupled atomistic spin and lattice dynamics approach is developed which merges the dynamics of these two degrees of freedom into a single set of coupled equations of motion. The underlying microscopic model comprises local exchange interactions between the electron spin and magnetic moment and the local couplings between the electronic charge and lattice displacements. An effective action for the spin and lattice variables is constructed in which the interactions among the spin and lattice components are determined by the underlying electronic structure. In this way, expressions are obtained for the electronically mediated couplings between the spin and lattice degrees of freedom, besides the well known interatomic force constants and spin-spin interactions. These former susceptibilities provide an atomistic ab initio description for the coupled spin and lattice dynamics. It is important to notice that this theory is strictly bilinear in the spin and lattice variables and provides a minimal model for the coupled dynamics of these subsystems and that the two subsystems are treated on the same footing. Questions concerning time-reversal and inversion symmetry are rigorously addressed and it is shown how these aspects are absorbed in the tensor structure of the interaction fields. By means of these results regarding the spin-lattice coupling, simple explanations of ionic dimerization in double-antiferromagnetic materials, as well as charge density waves induced by a nonuniform spin structure, are given. In the final parts, coupled equations of motion for the combined spin and lattice dynamics are constructed, which subsequently can be reduced to a form which is analogous to the Landau-Lifshitz-Gilbert equations for spin dynamics and a damped driven mechanical oscillator for the ionic motion. It is important to notice, however, that these equations comprise contributions that couple these descriptions into one unified formulation. Finally, Kubo-like expressions for

New multidimensional partially integrable generalization of S-integrable N-wave equation

International Nuclear Information System (INIS)

Zenchuk, A. I.

2007-01-01

This paper develops a modification of the dressing method based on the inhomogeneous linear integral equation with integral operator having nonempty kernel. The method allows one to construct the systems of multidimensional partial differential equations having differential polynomial structure in any dimension n. The associated solution space is not full, although it is parametrized by certain number of arbitrary functions of (n-1) variables. We consider four-dimensional generalization of the classical (2+1)-dimensional S-integrable N-wave equation as an example

Integrable discretizations of the short pulse equation

International Nuclear Information System (INIS)

Feng Baofeng; Maruno, Ken-ichi; Ohta, Yasuhiro

2010-01-01

In this paper, we propose integrable semi-discrete and full-discrete analogues of the short pulse (SP) equation. The key construction is the bilinear form and determinant structure of solutions of the SP equation. We also give the determinant formulas of N-soliton solutions of the semi-discrete and full-discrete analogues of the SP equations, from which the multi-loop and multi-breather solutions can be generated. In the continuous limit, the full-discrete SP equation converges to the semi-discrete SP equation, and then to the continuous SP equation. Based on the semi-discrete SP equation, an integrable numerical scheme, i.e. a self-adaptive moving mesh scheme, is proposed and used for the numerical computation of the short pulse equation.

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)

Topological horseshoes in travelling waves of discretized nonlinear wave equations

International Nuclear Information System (INIS)

Chen, Yi-Chiuan; Chen, Shyan-Shiou; Yuan, Juan-Ming

2014-01-01

Applying the concept of anti-integrable limit to coupled map lattices originated from space-time discretized nonlinear wave equations, we show that there exist topological horseshoes in the phase space formed by the initial states of travelling wave solutions. In particular, the coupled map lattices display spatio-temporal chaos on the horseshoes

Topological horseshoes in travelling waves of discretized nonlinear wave equations

Energy Technology Data Exchange (ETDEWEB)

Chen, Yi-Chiuan, E-mail: YCChen@math.sinica.edu.tw [Institute of Mathematics, Academia Sinica, Taipei 10617, Taiwan (China); Chen, Shyan-Shiou, E-mail: sschen@ntnu.edu.tw [Department of Mathematics, National Taiwan Normal University, Taipei 11677, Taiwan (China); Yuan, Juan-Ming, E-mail: jmyuan@pu.edu.tw [Department of Financial and Computational Mathematics, Providence University, Shalu, Taichung 43301, Taiwan (China)

2014-04-15

Applying the concept of anti-integrable limit to coupled map lattices originated from space-time discretized nonlinear wave equations, we show that there exist topological horseshoes in the phase space formed by the initial states of travelling wave solutions. In particular, the coupled map lattices display spatio-temporal chaos on the horseshoes.

Integrable peakon equations with cubic nonlinearity

International Nuclear Information System (INIS)

Hone, Andrew N W; Wang, J P

2008-01-01

We present a new integrable partial differential equation found by Vladimir Novikov. Like the Camassa-Holm and Degasperis-Procesi equations, this new equation admits peaked soliton (peakon) solutions, but it has nonlinear terms that are cubic, rather than quadratic. We give a matrix Lax pair for V Novikov's equation, and show how it is related by a reciprocal transformation to a negative flow in the Sawada-Kotera hierarchy. Infinitely many conserved quantities are found, as well as a bi-Hamiltonian structure. The latter is used to obtain the Hamiltonian form of the finite-dimensional system for the interaction of N peakons, and the two-body dynamics (N = 2) is explicitly integrated. Finally, all of this is compared with some analogous results for another cubic peakon equation derived by Zhijun Qiao. (fast track communication)

Dynamical barrier for the formation of solitary waves in discrete lattices

International Nuclear Information System (INIS)

Kevrekidis, P.G.; Espinola-Rocha, J.A.; Drossinos, Y.; Stefanov, A.

2008-01-01

We consider the problem of the existence of a dynamical barrier of 'mass' that needs to be excited on a lattice site to lead to the formation and subsequent persistence of localized modes for a nonlinear Schroedinger lattice. We contrast the existence of a dynamical barrier with its absence in the static theory of localized modes in one spatial dimension. We suggest an energetic criterion that provides a sufficient, but not necessary, condition on the amplitude of a single-site initial condition required to form a solitary wave. We show that this effect is not one-dimensional by considering its two-dimensional analog. The existence of a sufficient condition for the excitation of localized modes in the non-integrable, discrete, nonlinear Schroedinger equation is compared to the dynamics of excitations in the integrable, both discrete and continuum, version of the nonlinear Schroedinger equation

Picard-Fuchs equations of dimensionally regulated Feynman integrals

Energy Technology Data Exchange (ETDEWEB)

Zayadeh, Raphael

2013-12-15

This thesis is devoted to studying differential equations of Feynman integrals. A Feynman integral depends on a dimension D. For integer values of D it can be written as a projective integral, which is called the Feynman parameter prescription. A major complication arises from the fact that for some values of D the integral can diverge. This problem is solved within dimensional regularization by continuing the integral as a meromorphic function on the complex plane and replacing the ill-defined quantity by a Laurent series in a dimensional regularization parameter. All terms in such a Laurent expansion are periods in the sense of Kontsevich and Zagier. We describe a new method to compute differential equations of Feynman integrals. So far, the standard has been to use integration-by-parts (IBP) identities to obtain coupled systems of linear differential equations for the master integrals. Our method is based on the theory of Picard-Fuchs equations. In the case we are interested in, that of projective and quasiprojective families, a Picard-Fuchs equation can be computed by means of the Griffiths-Dwork reduction. We describe a method that is designed for fixed integer dimension. After a suitable integer shift of dimension we obtain a period of a family of hypersurfaces, hence a Picard-Fuchs equation. This equation is inhomogeneous because the domain of integration has a boundary and we only obtain a relative cycle. As a second step we shift back the dimension using Tarasov's generalized dimensional recurrence relations. Furthermore, we describe a method to directly compute the differential equation for general D without shifting the dimension. This is based on the Griffiths-Dwork reduction. The success of this method depends on the ability to solve large systems of linear equations. We give examples of two and three-loop graphs. Tarasov classifies two-loop two-point functions and we give differential equations for these. For us the most interesting example is

Picard-Fuchs equations of dimensionally regulated Feynman integrals

International Nuclear Information System (INIS)

Zayadeh, Raphael

2013-12-01

This thesis is devoted to studying differential equations of Feynman integrals. A Feynman integral depends on a dimension D. For integer values of D it can be written as a projective integral, which is called the Feynman parameter prescription. A major complication arises from the fact that for some values of D the integral can diverge. This problem is solved within dimensional regularization by continuing the integral as a meromorphic function on the complex plane and replacing the ill-defined quantity by a Laurent series in a dimensional regularization parameter. All terms in such a Laurent expansion are periods in the sense of Kontsevich and Zagier. We describe a new method to compute differential equations of Feynman integrals. So far, the standard has been to use integration-by-parts (IBP) identities to obtain coupled systems of linear differential equations for the master integrals. Our method is based on the theory of Picard-Fuchs equations. In the case we are interested in, that of projective and quasiprojective families, a Picard-Fuchs equation can be computed by means of the Griffiths-Dwork reduction. We describe a method that is designed for fixed integer dimension. After a suitable integer shift of dimension we obtain a period of a family of hypersurfaces, hence a Picard-Fuchs equation. This equation is inhomogeneous because the domain of integration has a boundary and we only obtain a relative cycle. As a second step we shift back the dimension using Tarasov's generalized dimensional recurrence relations. Furthermore, we describe a method to directly compute the differential equation for general D without shifting the dimension. This is based on the Griffiths-Dwork reduction. The success of this method depends on the ability to solve large systems of linear equations. We give examples of two and three-loop graphs. Tarasov classifies two-loop two-point functions and we give differential equations for these. For us the most interesting example is the two

Measurement of the effect of the lattice pitch on the effective resonance integral of natural uranium

International Nuclear Information System (INIS)

Krcevinac, S.; Takac, S.

1966-04-01

The present work investigates the dependence of the effective resonance integral on the lattice pitch. Theoretically, the dependence is determined starting from Wigner's rational approximation in which the lattice is characterized by the effective ratio (s/M). Later this allows correlation between the lattice resonance integral and the resonance integral of the single rod (the rod in infinite medium). Using two approximations for Dancoff's factor we give the measured functional dependence of the effective resonance integral on the effective (s/M) ratio. The activation method and the differential technique of measuring absorption distribution in U-238 are used to determine the resonance integral experimentally. Since the effective lattice pitch cannot be defined with certain reliability, due to the use of cadmium in determining the cadmium ratio in the fuel rod, besides other perturbing effects, the method of comparing thermal activations of U-238 and a suitable thermal detector are used to determine the cadmium ratio

Measurement of the effect of the lattice pitch on the effective resonance integral of natural uranium

Energy Technology Data Exchange (ETDEWEB)

Krcevinac, S; Takac, S [Institute of Nuclear Sciences Boris Kidric, Vinca, Beograd (Yugoslavia)

1966-04-15

The present work investigates the dependence of the effective resonance integral on the lattice pitch. Theoretically, the dependence is determined starting from Wigner's rational approximation in which the lattice is characterized by the effective ratio (s/M). Later this allows correlation between the lattice resonance integral and the resonance integral of the single rod (the rod in infinite medium). Using two approximations for Dancoff's factor we give the measured functional dependence of the effective resonance integral on the effective (s/M) ratio. The activation method and the differential technique of measuring absorption distribution in U-238 are used to determine the resonance integral experimentally. Since the effective lattice pitch cannot be defined with certain reliability, due to the use of cadmium in determining the cadmium ratio in the fuel rod, besides other perturbing effects, the method of comparing thermal activations of U-238 and a suitable thermal detector are used to determine the cadmium ratio.

High order Fuchsian equations for the square lattice Ising model: χ-tilde(5)

International Nuclear Information System (INIS)

Bostan, A; Boukraa, S; Guttmann, A J; Jensen, I; Hassani, S; Zenine, N; Maillard, J-M

2009-01-01

We consider the Fuchsian linear differential equation obtained (modulo a prime) for χ-tilde (5) , the five-particle contribution to the susceptibility of the square lattice Ising model. We show that one can understand the factorization of the corresponding linear differential operator from calculations using just a single prime. A particular linear combination of χ-tilde (1) and χ-tilde (3) can be removed from χ-tilde (5) and the resulting series is annihilated by a high order globally nilpotent linear ODE. The corresponding (minimal order) linear differential operator, of order 29, splits into factors of small orders. A fifth-order linear differential operator occurs as the left-most factor of the 'depleted' differential operator and it is shown to be equivalent to the symmetric fourth power of L E , the linear differential operator corresponding to the elliptic integral E. This result generalizes what we have found for the lower order terms χ-tilde (3) and χ-tilde (4) . We conjecture that a linear differential operator equivalent to a symmetric (n - 1) th power of L E occurs as a left-most factor in the minimal order linear differential operators for all χ-tilde (n) 's

The two-wave X-ray field calculated by means of integral-equation methods

International Nuclear Information System (INIS)

Bremer, J.

1984-01-01

The problem of calculating the two-wave X-ray field on the basis of the Takagi-Taupin equations is discussed for the general case of curved lattice planes. A two-dimensional integral equation which incorporates the nature of the incoming radiation, the form of the crystal/vacuum boundary, and the curvature of the structure, is deduced. Analytical solutions for the symmetrical Laue case with incoming plane waves are obtained directly for perfect crystals by means of iteration. The same method permits a simple derivation of the narrow-wave Laue and Bragg cases. Modulated wave fronts are discussed, and it is shown that a cut-off in the width of an incoming plane wave leads to lateral oscillations which are superimposed on the Pendelloesung fringes. Bragg and Laue shadow fields are obtained. The influence of a non-zero kernel is discussed and a numerical procedure for calculating wave amplitudes in curved crystals is presented. (Auth.)

Exactly and completely integrable nonlinear dynamical systems

International Nuclear Information System (INIS)

Leznov, A.N.; Savel'ev, M.V.

1987-01-01

The survey is devoted to a consitent exposition of the group-algebraic methods for the integration of systems of nonlinear partial differential equations possessing a nontrivial internal symmetry algebra. Samples of exactly and completely integrable wave and evolution equations are considered in detail, including generalized (periodic and finite nonperiodic Toda lattice, nonlinear Schroedinger, Korteweg-de Vries, Lotka-Volterra equations, etc.) For exactly integrable systems the general solutions of the Cauchy and Goursat problems are given in an explicit form, while for completely integrable systems an effective method for the construction of their soliton solutions is developed. Application of the developed methods to a differential geometry problem of classification of the integrable manifolds embeddings is discussed. For exactly integrable systems the supersymmetric extensions are constructed. By the example of the generalized Toda lattice a quantization scheme is developed. It includes an explicit derivation of the corresponding Heisenberg operators and their desription in terms of the quantum algebras of the Hopf type. Among multidimensional systems the four-dimensional self-dual Yang-Mills equations are investigated most attentively with a goal of constructing their general solutions

Counting master integrals. Integration by parts vs. functional equations

International Nuclear Information System (INIS)

Kniehl, Bernd A.; Tarasov, Oleg V.

2016-01-01

We illustrate the usefulness of functional equations in establishing relationships between master integrals under the integration-by-parts reduction procedure by considering a certain two-loop propagator-type diagram as an example.

Symplectic maps for accelerator lattices

International Nuclear Information System (INIS)

Warnock, R.L.; Ruth, R.; Gabella, W.

1988-05-01

We describe a method for numerical construction of a symplectic map for particle propagation in a general accelerator lattice. The generating function of the map is obtained by integrating the Hamilton-Jacobi equation as an initial-value problem on a finite time interval. Given the generating function, the map is put in explicit form by means of a Fourier inversion technique. We give an example which suggests that the method has promise. 9 refs., 9 figs

Dynamical barrier for the formation of solitary waves in discrete lattices

Energy Technology Data Exchange (ETDEWEB)

Kevrekidis, P.G. [Department of Mathematics and Statistics, University of Massachusetts, Amherst, MA 01003 (United States)], E-mail: kevrekid@math.umass.edu; Espinola-Rocha, J.A. [Department of Mathematics and Statistics, University of Massachusetts, Amherst, MA 01003 (United States); Drossinos, Y. [European Commission, Joint Research Centre, I-21020 Ispra (Vatican City State, Holy See,) (Italy); School of Mechanical and Systems Engineering, University of Newcastle upon Tyne, Newcastle upon Tyne NE1 7RU (United Kingdom); Stefanov, A. [Department of Mathematics, University of Kansas, 1460 Jayhawk Blvd., Lawrence, KS 66045-7523 (United States)

2008-03-24

We consider the problem of the existence of a dynamical barrier of 'mass' that needs to be excited on a lattice site to lead to the formation and subsequent persistence of localized modes for a nonlinear Schroedinger lattice. We contrast the existence of a dynamical barrier with its absence in the static theory of localized modes in one spatial dimension. We suggest an energetic criterion that provides a sufficient, but not necessary, condition on the amplitude of a single-site initial condition required to form a solitary wave. We show that this effect is not one-dimensional by considering its two-dimensional analog. The existence of a sufficient condition for the excitation of localized modes in the non-integrable, discrete, nonlinear Schroedinger equation is compared to the dynamics of excitations in the integrable, both discrete and continuum, version of the nonlinear Schroedinger equation.

Integrable coupling system of fractional soliton equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

Yu Fajun, E-mail: yfajun@163.co [College of Maths and Systematic Science, Shenyang Normal University, Shenyang 110034 (China)

2009-10-05

In this Letter, we consider the derivatives and integrals of fractional order and present a class of the integrable coupling system of the fractional order soliton equations. The fractional order coupled Boussinesq and KdV equations are the special cases of this class. Furthermore, the fractional AKNS soliton equation hierarchy is obtained.

Invariant relations in Boussinesq-type equations

International Nuclear Information System (INIS)

Meletlidou, Efi; Pouget, Joeel; Maugin, Gerard; Aifantis, Elias

2004-01-01

A wide class of partial differential equations have at least three conservation laws that remain invariant for certain solutions of them and especially for solitary wave solutions. These conservation laws can be considered as the energy, pseudomomentum and mass integrals of these solutions. We investigate the invariant relation between the energy and the pseudomomentum for solitary waves in two Boussinesq-type equations that come from the theory of elasticity and lattice models

Lattice Integrable Systems of the Haldane-Shastry Type

CERN Document Server

Polychronakos, A P

1993-01-01

We present a new lattice integrable system in one dimension of the Haldane-Shastry type. It consists of spins positioned at the static equilibrium positions of particles in a corresponding classical Calogero system and interacting through an exchange term with strength inversely proportional to the square of their distance. We achieve this by viewing the Haldane-Shastry system as a high-interaction limit of the Sutherland system of particles with internal degrees of freedom and identifying the same limit in a corresponding Calogero system. The commuting integrals of motion of this system are found using the exchange operator formalism.

Inequalities for differential and integral equations

CERN Document Server

Ames, William F

1997-01-01

Inequalities for Differential and Integral Equations has long been needed; it contains material which is hard to find in other books. Written by a major contributor to the field, this comprehensive resource contains many inequalities which have only recently appeared in the literature and which can be used as powerful tools in the development of applications in the theory of new classes of differential and integral equations. For researchers working in this area, it will be a valuable source of reference and inspiration. It could also be used as the text for an advanced graduate course.Key Features* Covers a variety of linear and nonlinear inequalities which find widespread applications in the theory of various classes of differential and integral equations* Contains many inequalities which have only recently appeared in literature and cannot yet be found in other books* Provides a valuable reference to engineers and graduate students

Decoupled scheme based on the Hermite expansion to construct lattice Boltzmann models for the compressible Navier-Stokes equations with arbitrary specific heat ratio.

Science.gov (United States)

Hu, Kainan; Zhang, Hongwu; Geng, Shaojuan

2016-10-01

A decoupled scheme based on the Hermite expansion to construct lattice Boltzmann models for the compressible Navier-Stokes equations with arbitrary specific heat ratio is proposed. The local equilibrium distribution function including the rotational velocity of particle is decoupled into two parts, i.e., the local equilibrium distribution function of the translational velocity of particle and that of the rotational velocity of particle. From these two local equilibrium functions, two lattice Boltzmann models are derived via the Hermite expansion, namely one is in relation to the translational velocity and the other is connected with the rotational velocity. Accordingly, the distribution function is also decoupled. After this, the evolution equation is decoupled into the evolution equation of the translational velocity and that of the rotational velocity. The two evolution equations evolve separately. The lattice Boltzmann models used in the scheme proposed by this work are constructed via the Hermite expansion, so it is easy to construct new schemes of higher-order accuracy. To validate the proposed scheme, a one-dimensional shock tube simulation is performed. The numerical results agree with the analytical solutions very well.

A transverse lattice QCD model for mesons

Energy Technology Data Exchange (ETDEWEB)

Patel, Apoorva D.; Ratabole, Raghunath

2004-03-01

QCD is analysed with two light-front continuum dimensions and two transverse lattice dimensions. In the limit of large number of colours and strong transverse gauge coupling, the contributions of light-front and transverse directions factorise in the dynamics, and the theory can be analytically solved in a closed form. An integral equation is obtained, describing the properties of mesons, which generalises the 't Hooft equation by including spin degrees of freedom. The meson spectrum, light-front wavefunctions and form factors can be obtained by solving this equation numerically. These results would be a good starting point to model QCD observables which only weakly depend on transverse directions, e.g. deep inelastic scattering structure functions.

Continuous surface force based lattice Boltzmann equation method for simulating thermocapillary flow

International Nuclear Information System (INIS)

Zheng, Lin; Zheng, Song; Zhai, Qinglan

2016-01-01

In this paper, we extend a lattice Boltzmann equation (LBE) with continuous surface force (CSF) to simulate thermocapillary flows. The model is designed on our previous CSF LBE for athermal two phase flow, in which the interfacial tension forces and the Marangoni stresses as the results of the interface interactions between different phases are described by a conception of CSF. In this model, the sharp interfaces between different phases are separated by a narrow transition layers, and the kinetics and morphology evolution of phase separation would be characterized by an order parameter via Cahn–Hilliard equation which is solved in the frame work of LBE. The scalar convection–diffusion equation for temperature field is resolved by thermal LBE. The models are validated by thermal two layered Poiseuille flow, and two superimposed planar fluids at negligibly small Reynolds and Marangoni numbers for the thermocapillary driven convection, which have analytical solutions for the velocity and temperature. Then thermocapillary migration of two/three dimensional deformable droplet are simulated. Numerical results show that the predictions of present LBE agreed with the analytical solution/other numerical results. - Highlights: • A CSF LBE to thermocapillary flows. • Thermal layered Poiseuille flows. • Thermocapillary migration.

Continuous surface force based lattice Boltzmann equation method for simulating thermocapillary flow

Energy Technology Data Exchange (ETDEWEB)

Zheng, Lin, E-mail: lz@njust.edu.cn [School of Energy and Power Engineering, Nanjing University of Science and Technology, Nanjing 210094 (China); Zheng, Song [School of Mathematics and Statistics, Zhejiang University of Finance and Economics, Hangzhou 310018 (China); Zhai, Qinglan [School of Economics Management and Law, Chaohu University, Chaohu 238000 (China)

2016-02-05

In this paper, we extend a lattice Boltzmann equation (LBE) with continuous surface force (CSF) to simulate thermocapillary flows. The model is designed on our previous CSF LBE for athermal two phase flow, in which the interfacial tension forces and the Marangoni stresses as the results of the interface interactions between different phases are described by a conception of CSF. In this model, the sharp interfaces between different phases are separated by a narrow transition layers, and the kinetics and morphology evolution of phase separation would be characterized by an order parameter via Cahn–Hilliard equation which is solved in the frame work of LBE. The scalar convection–diffusion equation for temperature field is resolved by thermal LBE. The models are validated by thermal two layered Poiseuille flow, and two superimposed planar fluids at negligibly small Reynolds and Marangoni numbers for the thermocapillary driven convection, which have analytical solutions for the velocity and temperature. Then thermocapillary migration of two/three dimensional deformable droplet are simulated. Numerical results show that the predictions of present LBE agreed with the analytical solution/other numerical results. - Highlights: • A CSF LBE to thermocapillary flows. • Thermal layered Poiseuille flows. • Thermocapillary migration.

Bounded solutions for fuzzy differential and integral equations

Energy Technology Data Exchange (ETDEWEB)

Nieto, Juan J. [Departamento de Analisis Matematico Facultad de Matematicas Universidad de Santiago de Compostela, 15782 (Spain)] e-mail: amnieto@usc.es; Rodriguez-Lopez, Rosana [Departamento de Analisis Matematico Facultad de Matematicas Universidad de Santiago de Compostela, 15782 (Spain)] e-mail: amrosana@usc.es

2006-03-01

We find sufficient conditions for the boundness of every solution of first-order fuzzy differential equations as well as certain fuzzy integral equations. Our results are based on several theorems concerning crisp differential and integral inequalities.

Aliasing modes in the lattice Schwinger model

International Nuclear Information System (INIS)

Campos, Rafael G.; Tututi, Eduardo S.

2007-01-01

We study the Schwinger model on a lattice consisting of zeros of the Hermite polynomials that incorporates a lattice derivative and a discrete Fourier transform with many properties. Such a lattice produces a Klein-Gordon equation for the boson field and the exact value of the mass in the asymptotic limit if the boundaries are not taken into account. On the contrary, if the lattice is considered with boundaries new modes appear due to aliasing effects. In the continuum limit, however, this lattice yields also a Klein-Gordon equation with a reduced mass

The quantum group, Harper equation and structure of Bloch eigenstates on a honeycomb lattice

International Nuclear Information System (INIS)

Eliashvili, M; Tsitsishvili, G; Japaridze, G I

2012-01-01

The tight-binding model of quantum particles on a honeycomb lattice is investigated in the presence of a homogeneous magnetic field. Provided the magnetic flux per unit hexagon is a rational of the elementary flux, the one-particle Hamiltonian is expressed in terms of the generators of the quantum group U q (sl 2 ). Employing the functional representation of the quantum group U q (sl 2 ), the Harper equation is rewritten as a system of two coupled functional equations in the complex plane. For the special values of quasi-momentum, the entangled system admits solutions in terms of polynomials. The system is shown to exhibit a certain symmetry allowing us to resolve the entanglement, and a basic single equation determining the eigenvalues and eigenstates (polynomials) is obtained. Equations specifying the locations of the roots of polynomials in the complex plane are found. Employing numerical analysis, the roots of polynomials corresponding to different eigenstates are solved and diagrams exhibiting the ordered structure of one-particle eigenstates are depicted. (paper)

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.

Integral equations for four identical particles in angular momentum representation

International Nuclear Information System (INIS)

Kharchenko, V.F.; Shadchin, S.A.

1975-01-01

In integral equations of motion for a system of four identical spinless particles with central pair interactions, transition is realized from the representation of relative Jacobi momenta to the representation of their moduli and relative angular moments. As a result, the variables associated with the rotation of the system as a whole are separated in the equations. The integral equations of motion for four particles are reduced to the form of an infinite system of three-demensional integral equations. The four-particle kinematic factors contained in integral kernels are expressed in terms of three-particle type kinematic factors. In the case of separable two-particle interaction, the equations of motion for four particles have the form of an infinite system of two-dimensional integral equations

Resonance and web structure in discrete soliton systems: the two-dimensional Toda lattice and its fully discrete and ultra-discrete analogues

International Nuclear Information System (INIS)

Maruno, Ken-ichi; Biondini, Gino

2004-01-01

We present a class of solutions of the two-dimensional Toda lattice equation, its fully discrete analogue and its ultra-discrete limit. These solutions demonstrate the existence of soliton resonance and web-like structure in discrete integrable systems such as differential-difference equations, difference equations and cellular automata (ultra-discrete equations)

Microcanonical ensemble formulation of lattice gauge theory

International Nuclear Information System (INIS)

Callaway, D.J.E.; Rahman, A.

1982-01-01

A new formulation of lattice gauge theory without explicit path integrals or sums is obtained by using the microcanonical ensemble of statistical mechanics. Expectation values in the new formalism are calculated by solving a large set of coupled, nonlinear, ordinary differential equations. The average plaquette for compact electrodynamics calculated in this fashion agrees with standard Monte Carlo results. Possible advantages of the microcanonical method in applications to fermionic systems are discussed

Feynman integrals and difference equations

International Nuclear Information System (INIS)

Moch, S.; Schneider, C.

2007-09-01

We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically over difference fields, the so-called ΠΣ * -fields. We test the implementation of the Mathematica package Sigma on examples from recent higher order perturbative calculations in Quantum Chromodynamics. (orig.)

Factorization of the hypergeometric-type difference equation on non-uniform lattices: dynamical algebra

Energy Technology Data Exchange (ETDEWEB)

Alvarez-Nodarse, R [Departamento de Analisis Matematico, Universidad de Sevilla, Apdo. 1160, E-41080 Sevilla (Spain); Atakishiyev, N M [Instituto de Matematicas, UNAM, Apartado Postal 273-3, CP 62210 Cuernavaca, Morelos, Mexico (Germany); Costas-Santos, R S [Departamento de Matematicas, EPS, Universidad Carlos III de Madrid, Ave. Universidad 30, E-28911, Leganes, Madrid (Spain)

2005-01-07

We argue that one can factorize the difference equation of hypergeometric type on non-uniform lattices in the general case. It is shown that in the most cases of q-linear spectrum of the eigenvalues, this directly leads to the dynamical symmetry algebra su{sub q}(1, 1), whose generators are explicitly constructed in terms of the difference operators, obtained in the process of factorization. Thus all models with the q-linear spectrum (some of them, but not all, previously considered in a number of publications) can be treated in a unified form.

Variational Integrals of a Class of Nonhomogeneous -Harmonic Equations

Directory of Open Access Journals (Sweden)

Guanfeng Li

2014-01-01

Full Text Available We introduce a class of variational integrals whose Euler equations are nonhomogeneous -harmonic equations. We investigate the relationship between the minimization problem and the Euler equation and give a simple proof of the existence of some nonhomogeneous -harmonic equations by applying direct methods of the calculus of variations. Besides, we establish some interesting results on variational integrals.

Two-dimensional nonlinear equations of supersymmetric gauge theories

International Nuclear Information System (INIS)

Savel'ev, M.V.

1985-01-01

Supersymmetric generalization of two-dimensional nonlinear dynamical equations of gauge theories is presented. The nontrivial dynamics of a physical system in the supersymmetry and supergravity theories for (2+2)-dimensions is described by the integrable embeddings of Vsub(2/2) superspace into the flat enveloping superspace Rsub(N/M), supplied with the structure of a Lie superalgebra. An equation is derived which describes a supersymmetric generalization of the two-dimensional Toda lattice. It contains both super-Liouville and Sinh-Gordon equations

Direct numerical solution of the Ornstein-Zernike integral equation and spatial distribution of water around hydrophobic molecules

Science.gov (United States)

Ikeguchi, Mitsunori; Doi, Junta

1995-09-01

The Ornstein-Zernike integral equation (OZ equation) has been used to evaluate the distribution function of solvents around solutes, but its numerical solution is difficult for molecules with a complicated shape. This paper proposes a numerical method to directly solve the OZ equation by introducing the 3D lattice. The method employs no approximation the reference interaction site model (RISM) equation employed. The method enables one to obtain the spatial distribution of spherical solvents around solutes with an arbitrary shape. Numerical accuracy is sufficient when the grid-spacing is less than 0.5 Å for solvent water. The spatial water distribution around a propane molecule is demonstrated as an example of a nonspherical hydrophobic molecule using iso-value surfaces. The water model proposed by Pratt and Chandler is used. The distribution agrees with the molecular dynamics simulation. The distribution increases offshore molecular concavities. The spatial distribution of water around 5α-cholest-2-ene (C27H46) is visualized using computer graphics techniques and a similar trend is observed.

Method of mechanical quadratures for solving singular integral equations of various types

Science.gov (United States)

Sahakyan, A. V.; Amirjanyan, H. A.

2018-04-01

The method of mechanical quadratures is proposed as a common approach intended for solving the integral equations defined on finite intervals and containing Cauchy-type singular integrals. This method can be used to solve singular integral equations of the first and second kind, equations with generalized kernel, weakly singular equations, and integro-differential equations. The quadrature rules for several different integrals represented through the same coefficients are presented. This allows one to reduce the integral equations containing integrals of different types to a system of linear algebraic equations.

Feynman integrals and difference equations

Energy Technology Data Exchange (ETDEWEB)

Moch, S. [Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany); Schneider, C. [Johannes Kepler Univ., Linz (Austria). Research Inst. for Symbolic Computation

2007-09-15

We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically over difference fields, the so-called {pi}{sigma}{sup *}-fields. We test the implementation of the Mathematica package Sigma on examples from recent higher order perturbative calculations in Quantum Chromodynamics. (orig.)

Integral equation hierarchy for continuum percolation

International Nuclear Information System (INIS)

Given, J.A.

1988-01-01

In this thesis a projection operator technique is presented that yields hierarchies of integral equations satisfied exactly by the n-point connectedness functions in a continuum version of the site-bond percolation problem. The n-point connectedness functions carry the same structural information for a percolation problem as then-point correlation functions do for a thermal problem. This method extends the Potts model mapping of Fortuin and Kastelyn to the continuum by exploiting an s-state generalization of the Widom-Rowlinson model, a continuum model for phase separation. The projection operator technique is used to produce an integral equation hierarchy for percolation similar to the Born-Green heirarchy. The Kirkwood superposition approximation (SA) is extended to percolation in order to close this hierarchy and yield a nonlinear integral equation for the two-point connectedness function. The fact that this function, in the SA, is the analytic continuation to negative density of the two-point correlation function in a corresponding thermal problem is discussed. The BGY equation for percolation is solved numerically, both by an expansion in powers of the density, and by an iterative technique due to Kirkwood. It is argued both analytically and numerically, that the BYG equation for percolation, unlike its thermal counterpart, shows non-classical critical behavior, with η = 1 and γ = 0.05 ± .1. Finally a sequence of refinements to the superposition approximations based in the theory of fluids by Rice and Lekner is discussed

Generalized isothermic lattices

International Nuclear Information System (INIS)

Doliwa, Adam

2007-01-01

We study multi-dimensional quadrilateral lattices satisfying simultaneously two integrable constraints: a quadratic constraint and the projective Moutard constraint. When the lattice is two dimensional and the quadric under consideration is the Moebius sphere one obtains, after the stereographic projection, the discrete isothermic surfaces defined by Bobenko and Pinkall by an algebraic constraint imposed on the (complex) cross-ratio of the circular lattice. We derive the analogous condition for our generalized isothermic lattices using Steiner's projective structure of conics, and we present basic geometric constructions which encode integrability of the lattice. In particular, we introduce the Darboux transformation of the generalized isothermic lattice and we derive the corresponding Bianchi permutability principle. Finally, we study two-dimensional generalized isothermic lattices, in particular geometry of their initial boundary value problem

Master equations and the theory of stochastic path integrals

Science.gov (United States)

Weber, Markus F.; Frey, Erwin

2017-04-01

This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. Since the 1930s, master equations have served as a fundamental tool to understand the role of fluctuations in complex biological, chemical, and physical systems. Despite their simple appearance, analyses of master equations most often rely on low-noise approximations such as the Kramers-Moyal or the system size expansion, or require ad-hoc closure schemes for the derivation of low-order moment equations. We focus on numerical and analytical methods going beyond the low-noise limit and provide a unified framework for the study of master equations. After deriving the forward and backward master equations from the Chapman-Kolmogorov equation, we show how the two master equations can be cast into either of four linear partial differential equations (PDEs). Three of these PDEs are discussed in detail. The first PDE governs the time evolution of a generalized probability generating function whose basis depends on the stochastic process under consideration. Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE. The second PDE is novel and is obeyed by a distribution that is marginalized over an initial state. It proves useful for the computation of mean extinction times. The third PDE describes the time evolution of a ‘generating functional’, which generalizes the so-called Poisson representation. Subsequently, the solutions of the PDEs are expressed in terms of two path integrals: a ‘forward’ and a ‘backward’ path integral. Combined with inverse transformations, one obtains two distinct path integral representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Moreover, we show how a well-known path integral representation of averaged observables can be recovered from

Master equations and the theory of stochastic path integrals.

Science.gov (United States)

Weber, Markus F; Frey, Erwin

2017-04-01

This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. Since the 1930s, master equations have served as a fundamental tool to understand the role of fluctuations in complex biological, chemical, and physical systems. Despite their simple appearance, analyses of master equations most often rely on low-noise approximations such as the Kramers-Moyal or the system size expansion, or require ad-hoc closure schemes for the derivation of low-order moment equations. We focus on numerical and analytical methods going beyond the low-noise limit and provide a unified framework for the study of master equations. After deriving the forward and backward master equations from the Chapman-Kolmogorov equation, we show how the two master equations can be cast into either of four linear partial differential equations (PDEs). Three of these PDEs are discussed in detail. The first PDE governs the time evolution of a generalized probability generating function whose basis depends on the stochastic process under consideration. Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE. The second PDE is novel and is obeyed by a distribution that is marginalized over an initial state. It proves useful for the computation of mean extinction times. The third PDE describes the time evolution of a 'generating functional', which generalizes the so-called Poisson representation. Subsequently, the solutions of the PDEs are expressed in terms of two path integrals: a 'forward' and a 'backward' path integral. Combined with inverse transformations, one obtains two distinct path integral representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Moreover, we show how a well-known path integral representation of averaged observables can be recovered from them. Upon

Properties of one-dimensional anharmonic lattice solitons

Science.gov (United States)

Szeftel, Jacob; Laurent-Gengoux, Pascal; Ilisca, Ernest; Hebbache, Mohamed

2000-12-01

The existence of bell- and kink-shaped solitons moving at constant velocity while keeping a permanent profile is studied in infinite periodic monoatomic chains of arbitrary anharmonicity by taking advantage of the equation of motion being integrable with respect to solitons. A second-order, non-linear differential equation involving advanced and retarded terms must be solved, which is done by implementing a scheme based on the finite element and Newton's methods. If the potential has a harmonic limit, the asymptotic time-decay behaves exponentially and there is a dispersion relation between propagation velocity and decay time. Inversely if the potential has no harmonic limit, the asymptotic regime shows up either as a power-law or faster than exponential. Excellent agreement is achieved with Toda's model. Illustrative examples are also given for the Fermi-Pasta-Ulam and sine-Gordon potentials. Owing to integrability an effective one-body potential is worked out in each case. Lattice and continuum solitons differ markedly from one another as regards the amplitude versus propagation velocity relationship and the asymptotic time behavior. The relevance of the linear stability analysis when applied to solitons propagating in an infinite crystal is questioned. The reasons preventing solitons from arising in a diatomic lattice are discussed.

Abel integral equations analysis and applications

CERN Document Server

Gorenflo, Rudolf

1991-01-01

In many fields of application of mathematics, progress is crucially dependent on the good flow of information between (i) theoretical mathematicians looking for applications, (ii) mathematicians working in applications in need of theory, and (iii) scientists and engineers applying mathematical models and methods. The intention of this book is to stimulate this flow of information. In the first three chapters (accessible to third year students of mathematics and physics and to mathematically interested engineers) applications of Abel integral equations are surveyed broadly including determination of potentials, stereology, seismic travel times, spectroscopy, optical fibres. In subsequent chapters (requiring some background in functional analysis) mapping properties of Abel integral operators and their relation to other integral transforms in various function spaces are investi- gated, questions of existence and uniqueness of solutions of linear and nonlinear Abel integral equations are treated, and for equatio...

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 a Volterra Stieltjes integral equation

Directory of Open Access Journals (Sweden)

P. T. Vaz

1990-01-01

Full Text Available The paper deals with a study of linear Volterra integral equations involving Lebesgue-Stieltjes integrals in two independent variables. The authors prove an existence theorem using the Banach fixed-point principle. An explicit example is also considered.

Transmission problem for the Laplace equation and the integral equation method

Czech Academy of Sciences Publication Activity Database

Medková, Dagmar

2012-01-01

Roč. 387, č. 2 (2012), s. 837-843 ISSN 0022-247X Institutional research plan: CEZ:AV0Z10190503 Keywords : transmission problem * Laplace equation * boundary integral equation Subject RIV: BA - General Mathematics Impact factor: 1.050, year: 2012 http://www.sciencedirect.com/science/article/pii/S0022247X11008985

Numerov iteration method for second order integral-differential equation

International Nuclear Information System (INIS)

Zeng Fanan; Zhang Jiaju; Zhao Xuan

1987-01-01

In this paper, Numerov iterative method for second order integral-differential equation and system of equations are constructed. Numerical examples show that this method is better than direct method (Gauss elimination method) in CPU time and memoy requireing. Therefore, this method is an efficient method for solving integral-differential equation in nuclear physics

An integrable semi-discretization of the Boussinesq equation

International Nuclear Information System (INIS)

Zhang, Yingnan; Tian, Lixin

2016-01-01

Highlights: • A new integrable semi-discretization of the Boussinesq equation is present. • A Bäcklund transformation and a Lax pair for the differential-difference system is derived by using Hirota's bilinear method. • The soliton solutions of 'good' Boussinesq equation and numerical algorithms are investigated. - Abstract: In this paper, we present an integrable semi-discretization of the Boussinesq equation. Different from other discrete analogues, we discretize the ‘time’ variable and get an integrable differential-difference system. Under a standard limitation, the differential-difference system converges to the continuous Boussinesq equation such that the discrete system can be used to design numerical algorithms. Using Hirota's bilinear method, we find a Bäcklund transformation and a Lax pair of the differential-difference system. For the case of ‘good’ Boussinesq equation, we investigate the soliton solutions of its discrete analogue and design numerical algorithms. We find an effective way to reduce the phase shift caused by the discretization. The numerical results coincide with our analysis.

ON ASYMTOTIC APPROXIMATIONS OF FIRST INTEGRALS FOR DIFFERENTIAL AND DIFFERENCE EQUATIONS

Directory of Open Access Journals (Sweden)

W.T. van Horssen

2007-04-01

Full Text Available In this paper the concept of integrating factors for differential equations and the concept of invariance factors for difference equations to obtain first integrals or invariants will be presented. It will be shown that all integrating factors have to satisfya system of partial differential equations, and that all invariance factors have to satisfy a functional equation. In the period 1997-2001 a perturbation method based on integrating vectors was developed to approximate first integrals for systems of ordinary differential equations. This perturbation method will be reviewed shortly. Also in the paper the first results in the development of a perturbation method for difference equations based on invariance factors will be presented.

Construction and analysis of lattice Boltzmann methods applied to a 1D convection-diffusion equation

International Nuclear Information System (INIS)

Dellacherie, Stephane

2014-01-01

To solve the 1D (linear) convection-diffusion equation, we construct and we analyze two LBM schemes built on the D1Q2 lattice. We obtain these LBM schemes by showing that the 1D convection-diffusion equation is the fluid limit of a discrete velocity kinetic system. Then, we show in the periodic case that these LBM schemes are equivalent to a finite difference type scheme named LFCCDF scheme. This allows us, firstly, to prove the convergence in L∞ of these schemes, and to obtain discrete maximum principles for any time step in the case of the 1D diffusion equation with different boundary conditions. Secondly, this allows us to obtain most of these results for the Du Fort-Frankel scheme for a particular choice of the first iterate. We also underline that these LBM schemes can be applied to the (linear) advection equation and we obtain a stability result in L∞ under a classical CFL condition. Moreover, by proposing a probabilistic interpretation of these LBM schemes, we also obtain Monte-Carlo algorithms which approach the 1D (linear) diffusion equation. At last, we present numerical applications justifying these results. (authors)

Minimally coupled N-particle scattering integral equations

International Nuclear Information System (INIS)

Kowalski, K.L.

1977-01-01

A concise formalism is developed which permits the efficient representation and generalization of several known techniques for deriving connected-kernel N-particle scattering integral equations. The methods of Kouri, Levin, and Tobocman and Bencze and Redish which lead to minimally coupled integral equations are of special interest. The introduction of channel coupling arrays is characterized in a general manner and the common base of this technique and that of the so-called channel coupling scheme is clarified. It is found that in the Bencze-Redish formalism a particular coupling array has a crucial function but one different from that of the arrays employed by Kouri, Levin, and Tobocman. The apparent dependence of the proof of the minimality of the Bencze-Redish integral equations upon the form of the inhomogeneous term in these equations is eliminated. This is achieved by an investigation of the full (nonminimal) Bencze-Redish kernel. It is shown that the second power of this operator is connected, a result which is needed for the full applicability of the Bencze-Redish formalism. This is used to establish the relationship between the existence of solutions to the homogeneous form of the minimal equations and eigenvalues of the full Bencze-Redish kernel

An integral equation arising in two group neutron transport theory

International Nuclear Information System (INIS)

Cassell, J S; Williams, M M R

2003-01-01

An integral equation describing the fuel distribution necessary to maintain a flat flux in a nuclear reactor in two group transport theory is reduced to the solution of a singular integral equation. The formalism developed enables the physical aspects of the problem to be better understood and its relationship with the corresponding diffusion theory model is highlighted. The integral equation is solved by reducing it to a non-singular Fredholm equation which is then evaluated numerically

A power counting theorem for Feynman integrals on the lattice

International Nuclear Information System (INIS)

Reisz, T.

1988-01-01

A convergence theorem is proved, which states sufficient conditions for the existence of the continuum limit for a wide class of Feynman integrals on a space-time lattice. A new kind of a UV-divergence degree is introduced, which allows the formulation of the theorem in terms of power counting conditions. (orig.)

ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY

OpenAIRE

Enrique Gonzalo Reyes Garcia

2004-01-01

ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY Equations in partial derivatives appeared in the 18th century as essential tools for the analytic study of physical models and, later, they proved to be fundamental for the progress of mathematics. For example, fundamental results of modern differential geometry are based on deep theorems on differential equations. Reciprocally, it is possible to study differential equations through geometrical means just like it was done by o...

Conservative, unconditionally stable discretization methods for Hamiltonian equations, applied to wave motion in lattice equations modeling protein molecules

Science.gov (United States)

LeMesurier, Brenton

2012-01-01

A new approach is described for generating exactly energy-momentum conserving time discretizations for a wide class of Hamiltonian systems of DEs with quadratic momenta, including mechanical systems with central forces; it is well-suited in particular to the large systems that arise in both spatial discretizations of nonlinear wave equations and lattice equations such as the Davydov System modeling energetic pulse propagation in protein molecules. The method is unconditionally stable, making it well-suited to equations of broadly “Discrete NLS form”, including many arising in nonlinear optics. Key features of the resulting discretizations are exact conservation of both the Hamiltonian and quadratic conserved quantities related to continuous linear symmetries, preservation of time reversal symmetry, unconditional stability, and respecting the linearity of certain terms. The last feature allows a simple, efficient iterative solution of the resulting nonlinear algebraic systems that retain unconditional stability, avoiding the need for full Newton-type solvers. One distinction from earlier work on conservative discretizations is a new and more straightforward nearly canonical procedure for constructing the discretizations, based on a “discrete gradient calculus with product rule” that mimics the essential properties of partial derivatives. This numerical method is then used to study the Davydov system, revealing that previously conjectured continuum limit approximations by NLS do not hold, but that sech-like pulses related to NLS solitons can nevertheless sometimes arise.

The Dirac-Kaehler equation and fermions on the lattice

International Nuclear Information System (INIS)

Becher, P.

1982-05-01

The geometrical description of spinor fields by E. Kaehler is used to formulate a consistent lattice approximation of fermions. The relation to free simple Dirac fields as well as to Susskind's description of lattice fermions is clarified. The first steps towards a quantized interacting theory are given. The correspondence between the calculus of differential forms and concepts of algebraic topology is shown to be a useful method for a completely analogous treatment of the problems in the continuum and on the lattice. (orig.)

Fractional Quantum Field Theory: From Lattice to Continuum

Directory of Open Access Journals (Sweden)

Vasily E. Tarasov

2014-01-01

Full Text Available An approach to formulate fractional field theories on unbounded lattice space-time is suggested. A fractional-order analog of the lattice quantum field theories is considered. Lattice analogs of the fractional-order 4-dimensional differential operators are proposed. We prove that continuum limit of the suggested lattice field theory gives a fractional field theory for the continuum 4-dimensional space-time. The fractional field equations, which are derived from equations for lattice space-time with long-range properties of power-law type, contain the Riesz type derivatives on noninteger orders with respect to space-time coordinates.

A Generalized Analytic Operator-Valued Function Space Integral and a Related Integral Equation

International Nuclear Information System (INIS)

Chang, K.S.; Kim, B.S.; Park, C.H.; Ryu, K.S.

2003-01-01

We introduce a generalized Wiener measure associated with a Gaussian Markov process and define a generalized analytic operator-valued function space integral as a bounded linear operator from L p into L p-ci r cumflexprime (1< p ≤ 2) by the analytic continuation of the generalized Wiener integral. We prove the existence of the integral for certain functionals which involve some Borel measures. Also we show that the generalized analytic operator-valued function space integral satisfies an integral equation related to the generalized Schroedinger equation. The resulting theorems extend the theory of operator-valued function space integrals substantially and previous theorems about these integrals are generalized by our results

An integrable semi-discretization of the Boussinesq equation

Energy Technology Data Exchange (ETDEWEB)

Zhang, Yingnan, E-mail: ynzhang@njnu.edu.cn [Jiangsu Key Laboratory for NSLSCS, School of Mathematical Sciences, Nanjing Normal University, Nanjing, Jiangsu (China); Tian, Lixin, E-mail: tianlixin@njnu.edu.cn [Jiangsu Key Laboratory for NSLSCS, School of Mathematical Sciences, Nanjing Normal University, Nanjing, Jiangsu (China); Nonlinear Scientific Research Center, Jiangsu University, Zhenjiang, Jiangsu (China)

2016-10-23

Highlights: • A new integrable semi-discretization of the Boussinesq equation is present. • A Bäcklund transformation and a Lax pair for the differential-difference system is derived by using Hirota's bilinear method. • The soliton solutions of 'good' Boussinesq equation and numerical algorithms are investigated. - Abstract: In this paper, we present an integrable semi-discretization of the Boussinesq equation. Different from other discrete analogues, we discretize the ‘time’ variable and get an integrable differential-difference system. Under a standard limitation, the differential-difference system converges to the continuous Boussinesq equation such that the discrete system can be used to design numerical algorithms. Using Hirota's bilinear method, we find a Bäcklund transformation and a Lax pair of the differential-difference system. For the case of ‘good’ Boussinesq equation, we investigate the soliton solutions of its discrete analogue and design numerical algorithms. We find an effective way to reduce the phase shift caused by the discretization. The numerical results coincide with our analysis.

Normal and adjoint integral and integrodifferential neutron transport equations. Pt. 2

International Nuclear Information System (INIS)

Velarde, G.

1976-01-01

Using the simplifying hypotheses of the integrodifferential Boltzmann equations of neutron transport, given in JEN 334 report, several integral equations, and theirs adjoint ones, are obtained. Relations between the different normal and adjoint eigenfunctions are established and, in particular, proceeding from the integrodifferential Boltzmann equation it's found out the relation between the solutions of the adjoint equation of its integral one, and the solutions of the integral equation of its adjoint one (author)

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)

Relationships between lattice energies of inorganic ionic solids

Science.gov (United States)

Kaya, Savaş

2018-06-01

Lattice energy, which is a measure of the stabilities of inorganic ionic solids, is the energy required to decompose a solid into its constituent independent gaseous ions. In the present work, the relationships between lattice energies of many diatomic and triatomic inorganic ionic solids are revealed and a simple rule that can be used for the prediction of the lattice energies of inorganic ionic solids is introduced. According to this rule, the lattice energy of an AB molecule can be predicted with the help of the lattice energies of AX, BY and XY molecules in agreement with the experimental data. This rule is valid for not only diatomic molecules but also triatomic molecules. The lattice energy equations proposed in this rule provides compatible results with previously published lattice energy equations by Jenkins, Kaya, Born-Lande, Born-Mayer, Kapustinskii and Reddy. For a large set of tested molecules, calculated percent standard deviation values considering experimental data and the results of the equations proposed in this work are in general between %1-2%.

Solutions of the lattice sine–Gordon equation and the solitons of its cellular automaton

International Nuclear Information System (INIS)

Willox, R; Ramani, A; Grammaticos, B

2014-01-01

We analyse the solutions of the cellular automaton sine–Gordon equation and link them to solutions of the discrete, lattice, sine–Gordon. We show that while the ultradiscretizable, positive definite, solutions of the latter behave dispersively, certain parts of these dispersive waves nonetheless survive in the ultradiscrete limit, giving rise to the solutions of the cellular automaton. We examine the ultradiscrete solutions in the case of a generalized cellular automaton in which the dependent variable can assume non-integer values and we show that the collision of two solitary waves is inelastic, leading to the creation of a ‘bridge’ of constant height that links two outgoing structures. Based on the ultradiscrete form of the sine–Gordon equation we explain the appearance of this bridging region and we describe its interaction with a solitary wave. (paper)

On solvability of some quadratic functional-integral equation in Banach algebra

International Nuclear Information System (INIS)

Darwish, M.A.

2007-08-01

Using the technique of a suitable measure of non-compactness in Banach algebra, we prove an existence theorem for some functional-integral equations which contain, as particular cases, a lot of integral and functional-integral equations that arise in many branches of nonlinear analysis and its applications. Also, the famous Chandrasekhar's integral equation is considered as a special case. (author)

Integrable boundary conditions and modified Lax equations

International Nuclear Information System (INIS)

Avan, Jean; Doikou, Anastasia

2008-01-01

We consider integrable boundary conditions for both discrete and continuum classical integrable models. Local integrals of motion generated by the corresponding 'transfer' matrices give rise to time evolution equations for the initial Lax operator. We systematically identify the modified Lax pairs for both discrete and continuum boundary integrable models, depending on the classical r-matrix and the boundary matrix

Integral equation for inhomogeneous condensed bosons generalizing the Gross-Pitaevskii differential equation

International Nuclear Information System (INIS)

Angilella, G.G.N.; Pucci, R.; March, N.H.

2004-01-01

We give here the derivation of a Gross-Pitaevskii-type equation for inhomogeneous condensed bosons. Instead of the original Gross-Pitaevskii differential equation, we obtain an integral equation that implies less restrictive assumptions than are made in the very recent study of Pieri and Strinati [Phys. Rev. Lett. 91, 030401 (2003)]. In particular, the Thomas-Fermi approximation and the restriction to small spatial variations of the order parameter invoked in their study are avoided

Numerical solution of boundary-integral equations for molecular electrostatics.

Science.gov (United States)

Bardhan, Jaydeep P

2009-03-07

Numerous molecular processes, such as ion permeation through channel proteins, are governed by relatively small changes in energetics. As a result, theoretical investigations of these processes require accurate numerical methods. In the present paper, we evaluate the accuracy of two approaches to simulating boundary-integral equations for continuum models of the electrostatics of solvation. The analysis emphasizes boundary-element method simulations of the integral-equation formulation known as the apparent-surface-charge (ASC) method or polarizable-continuum model (PCM). In many numerical implementations of the ASC/PCM model, one forces the integral equation to be satisfied exactly at a set of discrete points on the boundary. We demonstrate in this paper that this approach to discretization, known as point collocation, is significantly less accurate than an alternative approach known as qualocation. Furthermore, the qualocation method offers this improvement in accuracy without increasing simulation time. Numerical examples demonstrate that electrostatic part of the solvation free energy, when calculated using the collocation and qualocation methods, can differ significantly; for a polypeptide, the answers can differ by as much as 10 kcal/mol (approximately 4% of the total electrostatic contribution to solvation). The applicability of the qualocation discretization to other integral-equation formulations is also discussed, and two equivalences between integral-equation methods are derived.

Irreducibility and co-primeness as an integrability criterion for discrete equations

International Nuclear Information System (INIS)

Kanki, Masataka; Mada, Jun; Mase, Takafumi; Tokihiro, Tetsuji

2014-01-01

We study the Laurent property, the irreducibility and co-primeness of discrete integrable and non-integrable equations. First we study a discrete integrable equation related to the Somos-4 sequence, and also a non-integrable equation as a comparison. We prove that the conditions of irreducibility and co-primeness hold only in the integrable case. Next, we generalize our previous results on the singularities of the discrete Korteweg–de Vries (dKdV) equation. In our previous paper (Kanki et al 2014 J. Phys. A: Math. Theor. 47 065201) we described the singularity confinement test (one of the integrability criteria) using the Laurent property, and the irreducibility, and co-primeness of the terms in the bilinear dKdV equation, in which we only considered simplified boundary conditions. This restriction was needed to obtain simple (monomial) relations between the bilinear form and the nonlinear form of the dKdV equation. In this paper, we prove the co-primeness of the terms in the nonlinear dKdV equation for general initial conditions and boundary conditions, by using the localization of Laurent rings and the interchange of the axes. We assert that co-primeness of the terms can be used as a new integrability criterion, which is a mathematical re-interpretation of the confinement of singularities in the case of discrete equations. (paper)

Set-Valued Stochastic Equation with Set-Valued Square Integrable Martingale

Directory of Open Access Journals (Sweden)

Li Jun-Gang

2017-01-01

Full Text Available In this paper, we shall introduce the stochastic integral of a stochastic process with respect to set-valued square integrable martingale. Then we shall give the Aumann integral measurable theorem, and give the set-valued stochastic Lebesgue integral and set-valued square integrable martingale integral equation. The existence and uniqueness of solution to set-valued stochastic integral equation are proved. The discussion will be useful in optimal control and mathematical finance in psychological factors.

Lattice Boltzmann model for numerical relativity.

Science.gov (United States)

Ilseven, E; Mendoza, M

2016-02-01

In the Z4 formulation, Einstein equations are written as a set of flux conservative first-order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for numerical relativity and validate it with well-established tests, also known as "apples with apples." Furthermore, we find that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improve. Finally, in order to show the potential of our approach, a linear scaling law for parallelization with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.

Nonlinear Stability of ρ-Functional Equations in Latticetic Random Banach Lattice Spaces

Directory of Open Access Journals (Sweden)

Mohammad Maleki V.

2018-02-01

Full Text Available In this paper, we prove the generalized nonlinear stability of the first and second of the following ρ -functional equations, G ( | a | Δ A * | b | Δ B * G ( | a | Δ A * * | b | − G ( | a | Δ B * * G ( | b | = ρ ( 2 G | a | Δ A * | b | 2 Δ B * G | a | Δ A * * | b | 2 − G ( | a | Δ B * * G ( | b | , and 2 G | a | Δ A * | b | 2 Δ B * G | a | Δ A * * | b | 2 − G ( | a | Δ B * * G ( | b | = ρ G ( | a | Δ A * | b | Δ B * G ( | a | Δ A * * | b | − G ( | a | Δ B * * G ( | b | in latticetic random Banach lattice spaces, where ρ is a fixed real or complex number with ρ ≠ 1 .

Babenko’s Approach to Abel’s Integral Equations

Directory of Open Access Journals (Sweden)

Chenkuan Li

2018-03-01

Full Text Available The goal of this paper is to investigate the following Abel’s integral equation of the second kind: y ( t + λ Γ ( α ∫ 0 t ( t − τ α − 1 y ( τ d τ = f ( t , ( t > 0 and its variants by fractional calculus. Applying Babenko’s approach and fractional integrals, we provide a general method for solving Abel’s integral equation and others with a demonstration of different types of examples by showing convergence of series. In particular, we extend this equation to a distributional space for any arbitrary α ∈ R by fractional operations of generalized functions for the first time and obtain several new and interesting results that cannot be realized in the classical sense or by the Laplace transform.

Symbolic-Numeric Integration of the Dynamical Cosserat Equations

KAUST Repository

Lyakhov, Dmitry A.

2017-08-29

We devise a symbolic-numeric approach to the integration of the dynamical part of the Cosserat equations, a system of nonlinear partial differential equations describing the mechanical behavior of slender structures, like fibers and rods. This is based on our previous results on the construction of a closed form general solution to the kinematic part of the Cosserat system. Our approach combines methods of numerical exponential integration and symbolic integration of the intermediate system of nonlinear ordinary differential equations describing the dynamics of one of the arbitrary vector-functions in the general solution of the kinematic part in terms of the module of the twist vector-function. We present an experimental comparison with the well-established generalized \\\\alpha -method illustrating the computational efficiency of our approach for problems in structural mechanics.

Symbolic-Numeric Integration of the Dynamical Cosserat Equations

KAUST Repository

Lyakhov, Dmitry A.; Gerdt, Vladimir P.; Weber, Andreas G.; Michels, Dominik L.

2017-01-01

We devise a symbolic-numeric approach to the integration of the dynamical part of the Cosserat equations, a system of nonlinear partial differential equations describing the mechanical behavior of slender structures, like fibers and rods. This is based on our previous results on the construction of a closed form general solution to the kinematic part of the Cosserat system. Our approach combines methods of numerical exponential integration and symbolic integration of the intermediate system of nonlinear ordinary differential equations describing the dynamics of one of the arbitrary vector-functions in the general solution of the kinematic part in terms of the module of the twist vector-function. We present an experimental comparison with the well-established generalized \\alpha -method illustrating the computational efficiency of our approach for problems in structural mechanics.

Measurement of the effect of the lattice pitch on the effective resonance integral of natural uranium

Energy Technology Data Exchange (ETDEWEB)

Krcevinac, S; Takac, S [Institut za nuklearne nauke ' Boris Kidric' , Vinca, Belgrade (Yugoslavia)

1966-07-01

The analytical theory of resonance absorption, as well as the numerical Monte Carlo method, allows calculation of the resonance integral. However, it is based on specific approximations so it may be used accurately enough in a limited number of cases. Likewise, insufficiently accurate know ledge of the basic nuclear constants (e.g. resonance parameters, etc.) used as input data in analytical calculation, leads to inaccurate determination of the resonance integral.. Therefore, experimental determination of the effective resonance integral is still indispensable. In some cases the experimental results are used as the exclusive source of information, or as the basis for the semiempirical technique of calculation, and in others as a check of the new theoretical procedures. There are several experimental methods of direct determination of the resonance integral: the activation method (1,2), the reactor oscillator and the danger coefficient method. Indirectly, using the results of critical experiments, it is possible to determine correlated values of the effective resonance integral. The present work investigates the dependence of the effective resonance integral on the lattice pitch. Theoretically, the dependence is determined starting from Wigner's rational approximation in which the lattice is characterized by the effective ratio (S/M) {sup X}. Later this allows correlation between the lattice resonance integral and the resonance integral of the single rod (the rod in infinite medium). Using two approximations for Dancoff's factor we give the measured functional dependence of the effective resonance integral on the effective (S/M){sup X} ratio. To determine the resonance integral experimentally we used the activation method and the differential technique of measuring absorption distribution in U{sup 2}38. Since, because of the use of cadmium in determining the cadmium ratio in the fuel rod the effective lattice pitch cannot be defined with certain reliability, besides

Integrability of the Gross-Pitaevskii equation with Feshbach resonance management

International Nuclear Information System (INIS)

Zhao Dun; Luo Honggang; Chai Huayue

2008-01-01

In this Letter we study the integrability of a class of Gross-Pitaevskii equations managed by Feshbach resonance in an expulsive parabolic external potential. By using WTC test, we find a condition under which the Gross-Pitaevskii equation is completely integrable. Under the present model, this integrability condition is completely consistent with that proposed by Serkin, Hasegawa, and Belyaeva [V.N. Serkin, A. Hasegawa, T.L. Belyaeva, Phys. Rev. Lett. 98 (2007) 074102]. Furthermore, this integrability can also be explicitly shown by a transformation, which can convert the Gross-Pitaevskii equation into the well-known standard nonlinear Schroedinger equation. By this transformation, each exact solution of the standard nonlinear Schroedinger equation can be converted into that of the Gross-Pitaevskii equation, which builds a systematical connection between the canonical solitons and the so-called nonautonomous ones. The finding of this transformation has a significant contribution to understanding the essential properties of the nonautonomous solitons and the dynamics of the Bose-Einstein condensates by using the Feshbach resonance technique

Comments on the integrability of the loop-space chiral equations

International Nuclear Information System (INIS)

Gu, C.; Wang, L.L.C.

1980-01-01

A demonstration is given how the ordinary space chiral equations provide the existence conditions for the infinite number of conserved currents and how these currents are related to the so-called inverse-scattering equations, whose integrability is provided by the original chiral equations. Loop-space chiral equations are introduced. The integrability conditions of the non-local currents in two possible different situations are discussed. In the first case, the generating functions are functionals of the loop alone. The integrability conditions are not satisfied and higher order conserved non-local currents do not exist. In the second case, the generating functions are functionals of the loop as well as a parameter the integrability conditions at a restricted point of the parameter are satisfied, however there is an infinite fold of arbitrariness. It indicates that additional guiding principles are needed in addition to the original loop-space chiral equation in order to uniquely determine the infinite conserved non-local currents as functionals of the loop and the parameter

First integrals of the axisymmetric shape equation of lipid membranes

Science.gov (United States)

Zhang, Yi-Heng; McDargh, Zachary; Tu, Zhan-Chun

2018-03-01

The shape equation of lipid membranes is a fourth-order partial differential equation. Under the axisymmetric condition, this equation was transformed into a second-order ordinary differential equation (ODE) by Zheng and Liu (Phys. Rev. E 48 2856 (1993)). Here we try to further reduce this second-order ODE to a first-order ODE. First, we invert the usual process of variational calculus, that is, we construct a Lagrangian for which the ODE is the corresponding Euler–Lagrange equation. Then, we seek symmetries of this Lagrangian according to the Noether theorem. Under a certain restriction on Lie groups of the shape equation, we find that the first integral only exists when the shape equation is identical to the Willmore equation, in which case the symmetry leading to the first integral is scale invariance. We also obtain the mechanical interpretation of the first integral by using the membrane stress tensor. Project supported by the National Natural Science Foundation of China (Grant No. 11274046) and the National Science Foundation of the United States (Grant No. 1515007).

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.

Towards the simplest hydrodynamic lattice-gas model.

Science.gov (United States)

Boghosian, Bruce M; Love, Peter J; Meyer, David A

2002-03-15

It has been known since 1986 that it is possible to construct simple lattice-gas cellular automata whose hydrodynamics are governed by the Navier-Stokes equations in two dimensions. The simplest such model heretofore known has six bits of state per site on a triangular lattice. In this work, we demonstrate that it is possible to construct a model with only five bits of state per site on a Kagome lattice. Moreover, the model has a simple, deterministic set of collision rules and is easily implemented on a computer. In this work, we derive the equilibrium distribution function for this lattice-gas automaton and carry out the Chapman-Enskog analysis to determine the form of the Navier-Stokes equations.

The finite temperature QCD phase transition and the thermodynamic equation of state. An investigation employing lattice QCD with Nf=2 twisted mass quarks

International Nuclear Information System (INIS)

Burger, Florian

2012-01-01

In this thesis we report about an investigation of the finite temperature crossover/phase transition of quantum chromodynamics and the evaluation of the thermodynamic equation of state. To this end the lattice method and the Wilson twisted mass discretisation of the quark action are used. This formulation is known to have an automatic improvement of lattice artifacts and thus an improved continuum limit behaviour. This work presents first robust results using this action for the non-vanishing temperature case. We investigate the chiral limit of the two flavour phase transition with several small values of the pion mass in order to address the open question of the order of the transition in the limit of vanishing quark mass. For the currently simulated pion masses in the range of 300 to 700 MeV we present evidence that the finite temperature transition is a crossover transition rather than a genuine phase transition. The chiral limit is investigated by comparing the scaling of the observed crossover temperature with the mass including several possible scenarios. Complementary to this approach the chiral condensate as the order parameter for the spontaneous breaking of chiral symmetry is analysed in comparison with the O(4) universal scaling function which characterises a second order transition. With respect to thermodynamics the equation of state is obtained from the trace anomaly employing the temperature integral method which provides the pressure and energy density in the crossover region. The continuum limit of the trace anomaly is studied by considering several values of N τ and the tree-level correction technique.

Kwong-Wong-type integral equation on time scales

Directory of Open Access Journals (Sweden)

Baoguo Jia

2011-09-01

Full Text Available Consider the second-order nonlinear dynamic equation $$ [r(tx^Delta(ho(t]^Delta+p(tf(x(t=0, $$ where $p(t$ is the backward jump operator. We obtain a Kwong-Wong-type integral equation, that is: If $x(t$ is a nonoscillatory solution of the above equation on $[T_0,infty$, then the integral equation $$ frac{r^sigma(tx^Delta(t}{f(x^sigma(t} =P^sigma(t+int^infty_{sigma(t}frac{r^sigma(s [int^1_0f'(x_h(sdh][x^Delta(s]^2}{f(x(s f(x^sigma(s}Delta s $$ is satisfied for $tgeq T_0$, where $P^sigma(t=int^infty_{sigma(t}p(sDelta s$, and $x_h(s=x(s+hmu(sx^Delta(s$. As an application, we show that the superlinear dynamic equation $$ [r(tx^{Delta}(ho(t]^Delta+p(tf(x(t=0, $$ is oscillatory, under certain conditions.

Complete integrability of the difference evolution equations

International Nuclear Information System (INIS)

Gerdjikov, V.S.; Ivanov, M.I.; Kulish, P.P.

1980-01-01

The class of exactly solvable nonlinear difference evolution equations (DEE) related to the discrete analog of the one-dimensional Dirac problem L is studied. For this starting from L we construct a special linear non-local operator Λ and obtain the expansions of w and σ 3 deltaw over its eigenfunctions, w being the potential in L. This allows us to obtain compact expressions for the integrals of motion and to prove that these DEE are completely integrable Hamiltonian systems. Moreover, it is shown that there exists a hierarchy of Hamiltonian structures, generated by Λ, and the action-angle variables are explicity calculated. As particular cases the difference analog of the non-linear Schroedinger equation and the modified Korteweg-de-Vries equation are considered. The quantization of these Hamiltonian system through the use of the quantum inverse scattering method is briefly discussed [ru

Gauge theories on a small lattice

International Nuclear Information System (INIS)

Robson, D.; Webber, D.M.

1980-01-01

We present exact solutions to U(1), SU(2), and SU(3) lattice gauge theories on a Kogut-Susskind lattice consisting of a single plaquette. We demonstrate precise equivalence between the U(1) theory and the harmonic oscillator on an infinite one-dimensional lattice, and between the SU(N) theory and an N-fermion Schroedinger equation. (orig.)

SLAROM, Neutron Flux Distribution and Spectra in Lattice Cell

International Nuclear Information System (INIS)

Nakagawa, M.; Tsuchihashi, K.

2002-01-01

1 - Description of program or function: SLAROM solves the neutron integral transport equations to determine flux distribution and spectra in a lattice and calculates cell averaged effective cross sections. 2 - Method of solution: Collision probability method for cell calculation and 1D diffusion for core calculation. 3 - Restrictions on the complexity of the problem: Variable dimensions are used throughout the program so that computer core requirements depend on a variety of program parameters

Low-energy scattering on the lattice

International Nuclear Information System (INIS)

Bour Bour, Shahin

2014-01-01

In this thesis we present precision benchmark calculations for two-component fermions in the unitarity limit using an ab initio method, namely Hamiltonian lattice formalism. We calculate the ground state energy for unpolarized four particles (Fermi gas) in a periodic cube as a fraction of the ground state energy of the non-interacting system for two independent representations of the lattice Hamiltonians. We obtain the values 0.211(2) and 0.210(2). These results are in full agreement with the Euclidean lattice and fixed-node diffusion Monte Carlo calculations. We also give an expression for the energy corrections to the binding energy of a bound state in a moving frame. These corrections contain information about the mass and number of the constituents and are topological in origin and will have a broad applications to the lattice calculations of nucleons, nuclei, hadronic molecules and cold atoms. As one of its applications we use this expression and determine the low-energy parameters for the fermion dimer elastic scattering in shallow binding limit. For our lattice calculations we use Luescher's finite volume method. From the lattice calculations we find κa fd =1.174(9) and κr fd =-0.029(13), where κ represents the binding momentum of dimer and a fd (r fd ) denotes the scattering length (effective-range). These results are confirmed by the continuum calculations using the Skorniakov-Ter-Martirosian integral equation which gives 1.17907(1) and -0.0383(3) for the scattering length and effective range, respectively.

Lattice gas cellular automata and lattice Boltzmann models an introduction

CERN Document Server

Wolf-Gladrow, Dieter A

2000-01-01

Lattice-gas cellular automata (LGCA) and lattice Boltzmann models (LBM) are relatively new and promising methods for the numerical solution of nonlinear partial differential equations. The book provides an introduction for graduate students and researchers. Working knowledge of calculus is required and experience in PDEs and fluid dynamics is recommended. Some peculiarities of cellular automata are outlined in Chapter 2. The properties of various LGCA and special coding techniques are discussed in Chapter 3. Concepts from statistical mechanics (Chapter 4) provide the necessary theoretical background for LGCA and LBM. The properties of lattice Boltzmann models and a method for their construction are presented in Chapter 5.

Initial states in integrable quantum field theory quenches from an integral equation hierarchy

Directory of Open Access Journals (Sweden)

D.X. Horváth

2016-01-01

Full Text Available We consider the problem of determining the initial state of integrable quantum field theory quenches in terms of the post-quench eigenstates. The corresponding overlaps are a fundamental input to most exact methods to treat integrable quantum quenches. We construct and examine an infinite integral equation hierarchy based on the form factor bootstrap, proposed earlier as a set of conditions determining the overlaps. Using quenches of the mass and interaction in Sinh-Gordon theory as a concrete example, we present theoretical arguments that the state has the squeezed coherent form expected for integrable quenches, and supporting an Ansatz for the solution of the hierarchy. Moreover we also develop an iterative method to solve numerically the lowest equation of the hierarchy. The iterative solution along with extensive numerical checks performed using the next equation of the hierarchy provides a strong numerical evidence that the proposed Ansatz gives a very good approximation for the solution.

Initial states in integrable quantum field theory quenches from an integral equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

Horváth, D.X., E-mail: esoxluciuslinne@gmail.com [MTA-BME “Momentum” Statistical Field Theory Research Group, Budafoki út 8, 1111 Budapest (Hungary); Department of Theoretical Physics, Budapest University of Technology and Economics, Budafoki út 8, 1111 Budapest (Hungary); Sotiriadis, S., E-mail: sotiriad@sissa.it [SISSA and INFN, Via Bonomea 265, 34136 Trieste (Italy); Takács, G., E-mail: takacsg@eik.bme.hu [MTA-BME “Momentum” Statistical Field Theory Research Group, Budafoki út 8, 1111 Budapest (Hungary); Department of Theoretical Physics, Budapest University of Technology and Economics, Budafoki út 8, 1111 Budapest (Hungary)

2016-01-15

We consider the problem of determining the initial state of integrable quantum field theory quenches in terms of the post-quench eigenstates. The corresponding overlaps are a fundamental input to most exact methods to treat integrable quantum quenches. We construct and examine an infinite integral equation hierarchy based on the form factor bootstrap, proposed earlier as a set of conditions determining the overlaps. Using quenches of the mass and interaction in Sinh-Gordon theory as a concrete example, we present theoretical arguments that the state has the squeezed coherent form expected for integrable quenches, and supporting an Ansatz for the solution of the hierarchy. Moreover we also develop an iterative method to solve numerically the lowest equation of the hierarchy. The iterative solution along with extensive numerical checks performed using the next equation of the hierarchy provides a strong numerical evidence that the proposed Ansatz gives a very good approximation for the solution.

Coupling Integrable Couplings of an Equation Hierarchy

International Nuclear Information System (INIS)

Wang Hui; Xia Tie-Cheng

2013-01-01

Based on a kind of Lie algebra G proposed by Zhang, one isospectral problem is designed. Under the framework of zero curvature equation, a new kind of integrable coupling of an equation hierarchy is generated using the methods proposed by Ma and Gao. With the help of variational identity, we get the Hamiltonian structure of the hierarchy. (general)

Multi-component bi-Hamiltonian Dirac integrable equations

Energy Technology Data Exchange (ETDEWEB)

Ma Wenxiu [Department of Mathematics and Statistics, University of South Florida, Tampa, FL 33620-5700 (United States)], E-mail: mawx@math.usf.edu

2009-01-15

A specific matrix iso-spectral problem of arbitrary order is introduced and an associated hierarchy of multi-component Dirac integrable equations is constructed within the framework of zero curvature equations. The bi-Hamiltonian structure of the obtained Dirac hierarchy is presented be means of the variational trace identity. Two examples in the cases of lower order are computed.

Quadratic algebras in the noncommutative integration method of wave equation

International Nuclear Information System (INIS)

Varaksin, O.L.

1995-01-01

The paper deals with the investigation of applications of the method of noncommutative integration of linear differential equations by partial derivatives. Nontrivial example was taken for integration of three-dimensions wave equation with the use of non-Abelian quadratic algebras

Lattice gauge theory using parallel processors

International Nuclear Information System (INIS)

Lee, T.D.; Chou, K.C.; Zichichi, A.

1987-01-01

The book's contents include: Lattice Gauge Theory Lectures: Introduction and Current Fermion Simulations; Monte Carlo Algorithms for Lattice Gauge Theory; Specialized Computers for Lattice Gauge Theory; Lattice Gauge Theory at Finite Temperature: A Monte Carlo Study; Computational Method - An Elementary Introduction to the Langevin Equation, Present Status of Numerical Quantum Chromodynamics; Random Lattice Field Theory; The GF11 Processor and Compiler; and The APE Computer and First Physics Results; Columbia Supercomputer Project: Parallel Supercomputer for Lattice QCD; Statistical and Systematic Errors in Numerical Simulations; Monte Carlo Simulation for LGT and Programming Techniques on the Columbia Supercomputer; Food for Thought: Five Lectures on Lattice Gauge Theory

Analytic solution of integral equations for molecular fluids

International Nuclear Information System (INIS)

Cummings, P.T.

1984-01-01

We review some recent progress in the analytic solution of integral equations for molecular fluids. The site-site Ornstein-Zernike (SSOZ) equation with approximate closures appropriate to homonuclear diatomic fluids both with and without attractive dispersion-like interactions has recently been solved in closed form analytically. In this paper, the close relationship between the SSOZ equation for homonuclear dumbells and the usual Ornstein-Zernike (OZ) equation for atomic fluids is carefully elucidated. This relationship is a key motivation for the analytic solutions of the SSOZ equation that have been obtained to date. (author)

Differential equations and integrable models: the SU(3) case

International Nuclear Information System (INIS)

Dorey, Patrick; Tateo, Roberto

2000-01-01

We exhibit a relationship between the massless a 2 (2) integrable quantum field theory and a certain third-order ordinary differential equation, thereby extending a recent result connecting the massless sine-Gordon model to the Schroedinger equation. This forms part of a more general correspondence involving A 2 -related Bethe ansatz systems and third-order differential equations. A non-linear integral equation for the generalised spectral problem is derived, and some numerical checks are performed. Duality properties are discussed, and a simple variant of the non-linear equation is suggested as a candidate to describe the finite volume ground state energies of minimal conformal field theories perturbed by the operators phi 12 , phi 21 and phi 15 . This is checked against previous results obtained using the thermodynamic Bethe ansatz

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)

Manifestly gauge invariant discretizations of the Schrödinger equation

International Nuclear Information System (INIS)

Halvorsen, Tore Gunnar; Kvaal, Simen

2012-01-01

Grid-based discretizations of the time dependent Schrödinger equation coupled to an external magnetic field are converted to manifest gauge invariant discretizations. This is done using generalizations of ideas used in classical lattice gauge theory, and the process defined is applicable to a large class of discretized differential operators. In particular, popular discretizations such as pseudospectral discretizations using the fast Fourier transform can be transformed to gauge invariant schemes. Also generic gauge invariant versions of generic time integration methods are considered, enabling completely gauge invariant calculations of the time dependent Schrödinger equation. Numerical examples illuminating the differences between a gauge invariant discretization and conventional discretization procedures are also presented. -- Highlights: ► We investigate the Schrödinger equation coupled to an external magnetic field. ► Any grid-based discretization is made trivially gauge invariant. ► An extension of classical lattice gauge theory.

High Efficiency Optical MEMS by the Integration of Photonic Lattices with Surface MEMS

Energy Technology Data Exchange (ETDEWEB)

FLEMING, JAMES G.; LIN, SHAWN-YU; MANI, SEETHAMBAL S.; RODGERS, M. STEVEN; DAGEL, DARYL J.

2002-11-01

This report outlines our work on the integration of high efficiency photonic lattice structures with MEMS (MicroElectroMechanical Systems). The simplest of these structures were based on 1-D mirror structures. These were integrated into a variety of devices, movable mirrors, switchable cavities and finally into Bragg fiber structures which enable the control of light in at least 2 dimensions. Of these devices, the most complex were the Bragg fibers. Bragg fibers consist of hollow tubes in which light is guided in a low index media (air) and confined by surrounding Bragg mirror stacks. In this work, structures with internal diameters from 5 to 30 microns have been fabricated and much larger structures should also be possible. We have demonstrated the fabrication of these structures with short wavelength band edges ranging from 400 to 1600nm. There may be potential applications for such structures in the fields of integrated optics and BioMEMS. We have also looked at the possibility of waveguiding in 3 dimensions by integrating defects into 3-dimensional photonic lattice structures. Eventually it may be possible to tune such structures by mechanically modulating the defects.

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

Numerical integration of asymptotic solutions of ordinary differential equations

Science.gov (United States)

Thurston, Gaylen A.

1989-01-01

Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.

Analysis of Leaky Modes in Photonic Crystal Fibers Using the Surface Integral Equation Method

Directory of Open Access Journals (Sweden)

Jung-Sheng Chiang

2018-04-01

Full Text Available A fully vectorial algorithm based on the surface integral equation method for the modelling of leaky modes in photonic crystal fibers (PCFs by solely solving the complex propagation constants of characteristic equations is presented. It can be used for calculations of the complex effective index and confinement losses of photonic crystal fibers. As complex root examination is the key technique in the solution, the new algorithm which possesses this technique can be used to solve the leaky modes of photonic crystal fibers. The leaky modes of solid-core PCFs with a hexagonal lattice of circular air-holes are reported and discussed. The simulation results indicate how the confinement loss by the imaginary part of the effective index changes with air-hole size, the number of rings of air-holes, and wavelength. Confinement loss reductions can be realized by increasing the air-hole size and the number of air-holes. The results show that the confinement loss rises with wavelength, implying that the light leaks more easily for longer wavelengths; meanwhile, the losses are decreased significantly as the air-hole size d/Λ is increased.

YM2: Continuum expectations, lattice convergence, and lassos

International Nuclear Information System (INIS)

Driver, B.K.

1989-01-01

The two dimensional Yang-Mills theory (YM 2 ) is analyzed in both the continuum and the lattice. In the complete axial gauge the continuum theory may be defined in terms of a Lie algebra valued white noise, and parallel translation may be defined by stochastic differential equations. This machinery is used to compute the expectations of gauge invariant functions of the parallel translation operators along a collection of curves C. The expectation values are expressed as finite dimensional integrals with densities that are products of the heat kernel on the structure group. The time parameters of the heat kernels are determined by the areas enclosed by the collection C, and the arguments are determined by the crossing topologies of the curves in C. The expectations for the Wilson lattice models have a similar structure, and from this it follows that in the limit of small lattice spacing the lattice expectations converge to the continuum expectations. It is also shown that the lasso variables advocated by L. Gross exist and are sufficient to generate all the measurable functions on the YM 2 -measure space. (orig.)

Energy preserving integration of bi-Hamiltonian partial differential equations

NARCIS (Netherlands)

Karasozen, B.; Simsek, G.

2013-01-01

The energy preserving average vector field (AVF) integrator is applied to evolutionary partial differential equations (PDEs) in bi-Hamiltonian form with nonconstant Poisson structures. Numerical results for the Korteweg de Vries (KdV) equation and for the Ito type coupled KdV equation confirm the

Simplifying Differential Equations for Multiscale Feynman Integrals beyond Multiple Polylogarithms.

Science.gov (United States)

Adams, Luise; Chaubey, Ekta; Weinzierl, Stefan

2017-04-07

In this Letter we exploit factorization properties of Picard-Fuchs operators to decouple differential equations for multiscale Feynman integrals. The algorithm reduces the differential equations to blocks of the size of the order of the irreducible factors of the Picard-Fuchs operator. As a side product, our method can be used to easily convert the differential equations for Feynman integrals which evaluate to multiple polylogarithms to an ϵ form.

Discrete Painlevé equations: an integrability paradigm

International Nuclear Information System (INIS)

Grammaticos, B; Ramani, A

2014-01-01

In this paper we present a review of results on discrete Painlevé equations. We begin with an introduction which serves as a refresher on the continuous Painlevé equations. Next, in the first, main part of the paper, we introduce the discrete Painlevé equations, the various methods for their derivation, and their properties as well as their classification scheme. Along the way we present a brief summary of the two major discrete integrability detectors and of Quispel–Roberts–Thompson mapping, which plays a primordial role in the derivation of discrete Painlevé equations. The second part of the paper is more technical and focuses on the presentation of new results on what are called asymmetric discrete Painlevé equations. (comment)

Lectures on differential equations for Feynman integrals

International Nuclear Information System (INIS)

Henn, Johannes M

2015-01-01

Over the last year significant progress was made in the understanding of the computation of Feynman integrals using differential equations (DE). These lectures give a review of these developments, while not assuming any prior knowledge of the subject. After an introduction to DE for Feynman integrals, we point out how they can be simplified using algorithms available in the mathematical literature. We discuss how this is related to a recent conjecture for a canonical form of the equations. We also discuss a complementary approach that is based on properties of the space–time loop integrands, and explain how the ideas of leading singularities and d-log representations can be used to find an optimal basis for the DE. Finally, as an application of these ideas we show how single-scale integrals can be bootstrapped using the Drinfeld associator of a DE. (topical review)

Bargmann Symmetry Constraint for a Family of Liouville Integrable Differential-Difference Equations

International Nuclear Information System (INIS)

Xu Xixiang

2012-01-01

A family of integrable differential-difference equations is derived from a new matrix spectral problem. The Hamiltonian forms of obtained differential-difference equations are constructed. The Liouville integrability for the obtained integrable family is proved. Then, Bargmann symmetry constraint of the obtained integrable family is presented by binary nonliearization method of Lax pairs and adjoint Lax pairs. Under this Bargmann symmetry constraints, an integrable symplectic map and a sequences of completely integrable finite-dimensional Hamiltonian systems in Liouville sense are worked out, and every integrable differential-difference equations in the obtained family is factored by the integrable symplectic map and a completely integrable finite-dimensional Hamiltonian system. (general)

Integrable semi-discretizations of the reduced Ostrovsky equation

International Nuclear Information System (INIS)

Feng, Bao-Feng; Maruno, Ken-ichi; Ohta, Yasuhiro

2015-01-01

Based on our previous work on the reduced Ostrovsky equation (J. Phys. A: Math. Theor. 45 355203), we construct its integrable semi-discretizations. Since the reduced Ostrovsky equation admits two alternative representations, one being its original form, the other the differentiated form (the short wave limit of the Degasperis–Procesi equation) two semi-discrete analogues of the reduced Ostrovsky equation are constructed possessing the same N-loop soliton solution. The relationship between these two versions of semi-discretizations is also clarified. (paper)

Block-pulse functions approach to numerical solution of Abel’s integral equation

Directory of Open Access Journals (Sweden)

Monireh Nosrati Sahlan

2015-12-01

Full Text Available This study aims to present a computational method for solving Abel’s integral equation of the second kind. The introduced method is based on the use of Block-pulse functions (BPFs via collocation method. Abel’s integral equations as singular Volterra integral equations are hard and heavy in computation, but because of the properties of BPFs, as is reported in examples, this method is more efficient and more accurate than some other methods for solving this class of integral equations. On the other hand, the benefit of this method is low cost of computing operations. The applied method transforms the singular integral equation into triangular linear algebraic system that can be solved easily. An error analysis is worked out and applications are demonstrated through illustrative examples.

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

Soliton solutions in a diatomic lattice system

International Nuclear Information System (INIS)

Yajima, Nobuo; Satsuma, Junkichi.

1979-04-01

A continuum limit is considered for a diatomic lattice system with a cubic nonlinearity. A long wave equation describing the interaction of acoustic and optical modes is obtained. It reduces, in certain approximations, to equations having coupled wave solutions. The solutions exhibit trapping of an optical mode by an acoustic soliton. The form of the trapped optical wave depends on the mass ratio of adjacent particles in the diatomic lattice. (author)

Integration of differential equations by the pseudo-linear (PL) approximation

International Nuclear Information System (INIS)

Bonalumi, Riccardo A.

1998-01-01

A new method of integrating differential equations was originated with the technique of approximately calculating the integrals called the pseudo-linear (PL) procedure: this method is A-stable. This article contains the following examples: 1st order ordinary differential equations (ODEs), 2nd order linear ODEs, stiff system of ODEs (neutron kinetics), one-dimensional parabolic (diffusion) partial differential equations. In this latter case, this PL method coincides with the Crank-Nicholson method

Matter-wave two-dimensional solitons in crossed linear and nonlinear optical lattices

International Nuclear Information System (INIS)

Luz, H. L. F. da; Gammal, A.; Abdullaev, F. Kh.; Salerno, M.; Tomio, Lauro

2010-01-01

The existence of multidimensional matter-wave solitons in a crossed optical lattice (OL) with a linear optical lattice (LOL) in the x direction and a nonlinear optical lattice (NOL) in the y direction, where the NOL can be generated by a periodic spatial modulation of the scattering length using an optically induced Feshbach resonance is demonstrated. In particular, we show that such crossed LOLs and NOLs allow for stabilizing two-dimensional solitons against decay or collapse for both attractive and repulsive interactions. The solutions for the soliton stability are investigated analytically, by using a multi-Gaussian variational approach, with the Vakhitov-Kolokolov necessary criterion for stability; and numerically, by using the relaxation method and direct numerical time integrations of the Gross-Pitaevskii equation. Very good agreement of the results corresponding to both treatments is observed.

Matter-wave two-dimensional solitons in crossed linear and nonlinear optical lattices

Science.gov (United States)

da Luz, H. L. F.; Abdullaev, F. Kh.; Gammal, A.; Salerno, M.; Tomio, Lauro

2010-10-01

The existence of multidimensional matter-wave solitons in a crossed optical lattice (OL) with a linear optical lattice (LOL) in the x direction and a nonlinear optical lattice (NOL) in the y direction, where the NOL can be generated by a periodic spatial modulation of the scattering length using an optically induced Feshbach resonance is demonstrated. In particular, we show that such crossed LOLs and NOLs allow for stabilizing two-dimensional solitons against decay or collapse for both attractive and repulsive interactions. The solutions for the soliton stability are investigated analytically, by using a multi-Gaussian variational approach, with the Vakhitov-Kolokolov necessary criterion for stability; and numerically, by using the relaxation method and direct numerical time integrations of the Gross-Pitaevskii equation. Very good agreement of the results corresponding to both treatments is observed.

On the complete integrability of the discrete Nahm equations

International Nuclear Information System (INIS)

Murray, M.K.

2000-01-01

The discrete Nahm equations, a system of matrix valued difference equations, arose in the work of Braam and Austin on half-integral mass hyperbolic monopoles. We show that the discrete Nahm equations are completely integrable in a natural sense: to any solution we can associate a spectral curve and a holomorphic line-bundle over the spectral curve, such that the discrete-time DN evolution corresponds to walking in the Jacobian of the spectral curve in a straight line through the line-bundle with steps of a fixed size. Some of the implications for hyperbolic monopoles are also discussed. (orig.)

Splines and their reciprocal-bases in volume-integral equations

International Nuclear Information System (INIS)

Sabbagh, H.A.

1993-01-01

The authors briefly outline the use of higher-order splines and their reciprocal-bases in discretizing the volume-integral equations of electromagnetics. The discretization is carried out by means of the method of moments, in which the expansion functions are the higher-order splines, and the testing functions are the corresponding reciprocal-basis functions. These functions satisfy an orthogonality condition with respect to the spline expansion functions. Thus, the method is not Galerkin, but the structure of the resulting equations is quite regular, nevertheless. The theory is applied to the volume-integral equations for the unknown current density, or unknown electric field, within a scattering body, and to the equations for eddy-current nondestructive evaluation. Numerical techniques for computing the matrix elements are also given

Evaluating four-loop conformal Feynman integrals by D-dimensional differential equations

Science.gov (United States)

Eden, Burkhard; Smirnov, Vladimir A.

2016-10-01

We evaluate a four-loop conformal integral, i.e. an integral over four four-dimensional coordinates, by turning to its dimensionally regularized version and applying differential equations for the set of the corresponding 213 master integrals. To solve these linear differential equations we follow the strategy suggested by Henn and switch to a uniformly transcendental basis of master integrals. We find a solution to these equations up to weight eight in terms of multiple polylogarithms. Further, we present an analytical result for the given four-loop conformal integral considered in four-dimensional space-time in terms of single-valued harmonic polylogarithms. As a by-product, we obtain analytical results for all the other 212 master integrals within dimensional regularization, i.e. considered in D dimensions.

Evaluating four-loop conformal Feynman integrals by D-dimensional differential equations

Energy Technology Data Exchange (ETDEWEB)

Eden, Burkhard [Institut für Mathematik und Physik, Humboldt-Universität zu Berlin,Zum großen Windkanal 6, 12489 Berlin (Germany); Smirnov, Vladimir A. [Skobeltsyn Institute of Nuclear Physics, Moscow State University,119992 Moscow (Russian Federation)

2016-10-21

We evaluate a four-loop conformal integral, i.e. an integral over four four-dimensional coordinates, by turning to its dimensionally regularized version and applying differential equations for the set of the corresponding 213 master integrals. To solve these linear differential equations we follow the strategy suggested by Henn and switch to a uniformly transcendental basis of master integrals. We find a solution to these equations up to weight eight in terms of multiple polylogarithms. Further, we present an analytical result for the given four-loop conformal integral considered in four-dimensional space-time in terms of single-valued harmonic polylogarithms. As a by-product, we obtain analytical results for all the other 212 master integrals within dimensional regularization, i.e. considered in D dimensions.

Localized and periodic exact solutions to the nonlinear Schroedinger equation with spatially modulated parameters: Linear and nonlinear lattices

International Nuclear Information System (INIS)

Belmonte-Beitia, Juan; Konotop, Vladimir V.; Perez-Garcia, Victor M.; Vekslerchik, Vadym E.

2009-01-01

Using similarity transformations we construct explicit solutions of the nonlinear Schroedinger equation with linear and nonlinear periodic potentials. We present explicit forms of spatially localized and periodic solutions, and study their properties. We put our results in the framework of the exploited perturbation techniques and discuss their implications on the properties of associated linear periodic potentials and on the possibilities of stabilization of gap solitons using polychromatic lattices.

A systematic method for correlating measurements of channel powers with the lattice constants in the neutron diffusion equations

International Nuclear Information System (INIS)

Buckler, A.N.

1978-10-01

The report describes the theoretical basis of the methods that have been developed for correlating measurements of spatially distributed quantities taken on the reactor with the lattice constants in the diffusion equations. The method can be used with any thermal reactor system of current interest, but the first application is to provide a replacement for the SAMSON code for Winfrith SGHW studies, where the measurements of interest are channel powers. (author)

Unconditionally stable integration of Maxwell's equations

NARCIS (Netherlands)

J.G. Verwer (Jan); M.A. Botchev

2008-01-01

htmlabstractNumerical integration of Maxwell''s equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction

Is Yang-Mills equation a totally integrable system. Lecture III

International Nuclear Information System (INIS)

Chau Wang, L.L.

1981-01-01

Topics covered include: loop-space formulation of gauge theory - loop-space chiral equation; two dimensional chiral equation - conservation laws, linear system and integrability; and parallel development for the loop-space chiral equation - subtlety

Adaptive integral equation methods in transport theory

International Nuclear Information System (INIS)

Kelley, C.T.

1992-01-01

In this paper, an adaptive multilevel algorithm for integral equations is described that has been developed with the Chandrasekhar H equation and its generalizations in mind. The algorithm maintains good performance when the Frechet derivative of the nonlinear map is singular at the solution, as happens in radiative transfer with conservative scattering and in critical neutron transport. Numerical examples that demonstrate the algorithm's effectiveness are presented

A Unified Theory of Non-Ideal Gas Lattice Boltzmann Models

Science.gov (United States)

Luo, Li-Shi

1998-01-01

A non-ideal gas lattice Boltzmann model is directly derived, in an a priori fashion, from the Enskog equation for dense gases. The model is rigorously obtained by a systematic procedure to discretize the Enskog equation (in the presence of an external force) in both phase space and time. The lattice Boltzmann model derived here is thermodynamically consistent and is free of the defects which exist in previous lattice Boltzmann models for non-ideal gases. The existing lattice Boltzmann models for non-ideal gases are analyzed and compared with the model derived here.

Integration of equations of parabolic type by the method of nets

CERN Document Server

Saul'Yev, V K; Stark, M; Ulam, S

1964-01-01

International Series of Monographs in Pure and Applied Mathematics, Volume 54: Integration of Equations of Parabolic Type by the Method of Nets deals with solving parabolic partial differential equations using the method of nets. The first part of this volume focuses on the construction of net equations, with emphasis on the stability and accuracy of the approximating net equations. The method of nets or method of finite differences (used to define the corresponding numerical method in ordinary differential equations) is one of many different approximate methods of integration of partial diff

Geometry, Heat Equation and Path Integrals on the Poincare Upper Half-Plane

OpenAIRE

Reijiro, KUBO; Research Institute for Theoretical Physics Hiroshima University

1988-01-01

Geometry, heat equation and Feynman's path integrals are studied on the Poincare upper half-plane. The fundamental solution to the heat equation ∂f/∂t=Δ_Hf is expressed in terms of a path integral defined on the upper half-plane. It is shown that Kac's statement that Feynman's path integral satisfies the Schrodinger equation is also valid for our case.

Unconditionally stable integration of Maxwell's equations

NARCIS (Netherlands)

Verwer, J.G.; Bochev, Mikhail A.

Numerical integration of Maxwell's equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction implicit finite difference

Unconditionally stable integration of Maxwell's equations

NARCIS (Netherlands)

J.G. Verwer (Jan); M.A. Botchev

2009-01-01

textabstractNumerical integration of Maxwell’s equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction implicit –

Excess Gibbs Energy for Ternary Lattice Solutions of Nonrandom Mixing

Energy Technology Data Exchange (ETDEWEB)

Jung, Hae Young [DukSung Womens University, Seoul (Korea, Republic of)

2008-12-15

It is assumed for three components lattice solution that the number of ways of arranging particles randomly on the lattice follows a normal distribution of a linear combination of N{sub 12}, N{sub 23}, N{sub 13} which are the number of the nearest neighbor interactions between different molecules. It is shown by random number simulations that this assumption is reasonable. From this distribution, an approximate equation for the excess Gibbs energy of three components lattice solution is derived. Using this equation, several liquid-vapor equilibria are calculated and compared with the results from other equations.

Numerical solution of the potential problem by integral equations without Green's functions

International Nuclear Information System (INIS)

De Mey, G.

1977-01-01

An integral equation technique will be presented to solve Laplace's equation in a two-dimensional area S. The Green's function has been replaced by a particular solution of Laplace equation in order to establish the integral equation. It is shown that accurate results can be obtained provided the pivotal elimination method is used to solve the linear algebraic set

ON THE CONSTRUCTION OF PARTIAL DIFFERENCE SCHEMES II: DISCRETE VARIABLES AND SCHWARZIAN LATTICES

Directory of Open Access Journals (Sweden)

Decio Levi

2016-06-01

Full Text Available In the process of constructing invariant difference schemes which approximate partial differential equations we write down a procedure for discretizing a partial differential equation on an arbitrary lattice. An open problem is the meaning of a lattice which does not satisfy the Clairaut–Schwarz–Young theorem. To analyze it we apply the procedure on a simple example, the potential Burgers equation with two different lattices, an orthogonal lattice which is invariant under the symmetries of the equation and satisfies the commutativity of the partial difference operators and an exponential lattice which is not invariant and does not satisfy the Clairaut–Schwarz–Young theorem. A discussion on the numerical results is presented showing the different behavior of both schemes for two different exact solutions and their numerical approximations.

Discrete breathers in honeycomb Fermi–Pasta–Ulam lattices

International Nuclear Information System (INIS)

AD Wattis, Jonathan; M James, Lauren

2014-01-01

We consider the two-dimensional Fermi–Pasta–Ulam lattice with hexagonal honeycomb symmetry, which is a Hamiltonian system describing the evolution of a scalar-valued quantity subject to nearest neighbour interactions. Using multiple-scale analysis we reduce the governing lattice equations to a nonlinear Schrödinger equation coupled to a second equation for an accompanying slow mode. Two cases in which the latter equation can be solved and so the system decoupled are considered in more detail: firstly, in the case of a symmetric potential, we derive the form of moving breathers. We find an ellipticity criterion for the wavenumbers of the carrier wave, together with asymptotic estimates for the breather energy. The minimum energy threshold depends on the wavenumber of the breather. We find that this threshold is locally maximized by stationary breathers. Secondly, for an asymmetric potential we find stationary breathers, which, even with a quadratic nonlinearity generate no second harmonic component in the breather. Plots of all our findings show clear hexagonal symmetry as we would expect from our lattice structure. Finally, we compare the properties of stationary breathers in the square, triangular and honeycomb lattices. (paper)

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)

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

Geometry, heat equation and path integrals on the Poincare upper half-plane

International Nuclear Information System (INIS)

Kubo, Reijiro.

1987-08-01

Geometry, heat equation and Feynman's path integrals are studied on the Poincare upper half-plane. The fundamental solution to the heat equation δf/δt = Δ H f is expressed in terms of a path integral defined on the upper half-plane. It is shown that Kac's proof that Feynman's path integral satisfies the Schroedinger equation is also valid for our case. (author)

Lax pair and exact solutions of a discrete coupled system related to coupled KdV and coupled mKdV equations

International Nuclear Information System (INIS)

Liu Ping; Jia Man; Lou Senyue

2007-01-01

A modified Korteweg-de Vries (mKdV) lattice is also found to be a discrete Korteweg-de Vries (KdV) equation in this paper. The Lax pair for the discrete equation is found with the help of the Lax pair for a similar discrete equation. A Lax-integrable coupled extension of the lattice is posed, which is a common discrete version of both the coupled KdV and coupled mKdV systems. Some rational expansions of the Jacobian elliptic, trigonometric and hyperbolic functions are used to construct cnoidal waves, negaton and positon solutions of the discrete coupled system

Exponential Convergence for Numerical Solution of Integral Equations Using Radial Basis Functions

Directory of Open Access Journals (Sweden)

Zakieh Avazzadeh

2014-01-01

Full Text Available We solve some different type of Urysohn integral equations by using the radial basis functions. These types include the linear and nonlinear Fredholm, Volterra, and mixed Volterra-Fredholm integral equations. Our main aim is to investigate the rate of convergence to solve these equations using the radial basis functions which have normic structure that utilize approximation in higher dimensions. Of course, the use of this method often leads to ill-posed systems. Thus we propose an algorithm to improve the results. Numerical results show that this method leads to the exponential convergence for solving integral equations as it was already confirmed for partial and ordinary differential equations.

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.''

Monograph - The Numerical Integration of Ordinary Differential Equations.

Science.gov (United States)

Hull, T. E.

The materials presented in this monograph are intended to be included in a course on ordinary differential equations at the upper division level in a college mathematics program. These materials provide an introduction to the numerical integration of ordinary differential equations, and they can be used to supplement a regular text on this…

A novel lattice energy calculation technique for simple inorganic crystals

Energy Technology Data Exchange (ETDEWEB)

Kaya, Cemal [Department of Chemistry, Faculty of Science, Cumhuriyet University, 58140 Sivas (Turkey); Kaya, Savaş, E-mail: savaskaya@cumhuriyet.edu.tr [Department of Chemistry, Faculty of Science, Cumhuriyet University, 58140 Sivas (Turkey); Banerjee, Priyabrata [Surface Engineering and Tribology Group, CSIR-Central Mechanical Engineering Research Institute, Mahatma Gandhi Avenue, Durgapur 713209 (India)

2017-01-01

In this pure theoretical study, a hitherto unexplored equation based on Shannon radii of the ions forming that crystal and chemical hardness of any crystal to calculate the lattice energies of simple inorganic ionic crystals has been presented. To prove the credibility of this equation, the results of the equation have been compared with experimental outcome obtained from Born-Fajans-Haber- cycle which is fundamentally enthalpy-based thermochemical cycle and prevalent theoretical approaches proposed for the calculation of lattice energies of ionic compounds. The results obtained and the comparisons made have demonstrated that the new equation is more useful compared to other theoretical approaches and allows to exceptionally accurate calculation of lattice energies of inorganic ionic crystals without doing any complex calculations.

Abecedarian School on Symmetries and Integrability of Difference Equations (ASIDE) & SIDE 12 International Conference Symmetries and Integrability of Difference Equations

CERN Document Server

Rebelo, Raphaël; Winternitz, Pavel

2017-01-01

This book shows how Lie group and integrability techniques, originally developed for differential equations, have been adapted to the case of difference equations. Difference equations are playing an increasingly important role in the natural sciences. Indeed, many phenomena are inherently discrete and thus naturally described by difference equations. More fundamentally, in subatomic physics, space-time may actually be discrete. Differential equations would then just be approximations of more basic discrete ones. Moreover, when using differential equations to analyze continuous processes, it is often necessary to resort to numerical methods. This always involves a discretization of the differential equations involved, thus replacing them by difference ones. Each of the nine peer-reviewed chapters in this volume serves as a self-contained treatment of a topic, containing introductory material as well as the latest research results and exercises. Each chapter is presented by one or more early career researchers...

LAMP-B: a Fortran program set for the lattice cell analysis by collision probability method

International Nuclear Information System (INIS)

Tsuchihashi, Keiichiro

1979-02-01

Nature of physical problem solved: LAMB-B solves an integral transport equation by the collision probability method for many variety of lattice cell geometries: spherical, plane and cylindrical lattice cell; square and hexagonal arrays of pin rods; annular clusters and square clusters. LAMP-B produces homogenized constants for multi and/or few group diffusion theory programs. Method of solution: LAMP-B performs an exact numerical integration to obtain the collision probabilities. Restrictions on the complexity of the problem: Not more than 68 group in the fast group calculation, and not more than 20 regions in the resonance integral calculation. Typical running time: It varies with the number of energy groups and the selection of the geometry. Unusual features of the program: Any or any combination of constituent subprograms can be used so that the partial use of this program is available. (author)

On monotonic solutions of an integral equation of Abel type

International Nuclear Information System (INIS)

Darwish, Mohamed Abdalla

2007-08-01

We present an existence theorem of monotonic solutions for a quadratic integral equation of Abel type in C[0, 1]. The famous Chandrasekhar's integral equation is considered as a special case. The concept of measure of noncompactness and a fi xed point theorem due to Darbo are the main tools in carrying out our proof. (author)

Alternative integral equations and perturbation expansions for self-coupled scalar fields

International Nuclear Information System (INIS)

Ford, L.H.

1985-01-01

It is shown that the theory of a self-coupled scalar field may be expressed in terms of a class of integral equations which include the Yang-Feldman equation as a particular case. Other integral equations in this class could be used to generate alternative perturbation expansions which contain a nonanalytic dependence upon the coupling constant and are less ultraviolet divergent than the conventional perturbation expansion. (orig.)

Geometrical-integrability constraints and equations of motion in four plus extended super spaces

International Nuclear Information System (INIS)

Chau, L.L.

1987-01-01

It is pointed out that many equations of motion in physics, including gravitational and Yang-Mills equations, have a common origin: i.e. they are the results of certain geometrical integrability conditions. These integrability conditions lead to linear systems and conservation laws that are important in integrating these equations of motion

Differential equations for loop integrals in Baikov representation

Science.gov (United States)

Bosma, Jorrit; Larsen, Kasper J.; Zhang, Yang

2018-05-01

We present a proof that differential equations for Feynman loop integrals can always be derived in Baikov representation without involving dimension-shift identities. We moreover show that in a large class of two- and three-loop diagrams it is possible to avoid squared propagators in the intermediate steps of setting up the differential equations.

TBA equations for excited states in the sine-Gordon model

International Nuclear Information System (INIS)

Balog, Janos; Hegedus, Arpad

2004-01-01

We propose thermodynamic Bethe ansatz (TBA) integral equations for multi-particle soliton (fermion) states in the sine-Gordon (massive Thirring) model. This is based on T-system and Y-system equations, which follow from the Bethe ansatz solution in the light-cone lattice formulation of the model. Even and odd charge sectors are treated on an equal footing, corresponding to periodic and twisted boundary conditions, respectively. The analytic properties of the Y-system functions are conjectured on the basis of the large volume solution of the system, which we find explicitly. A simple relation between the TBA Y-functions and the counting function variable of the alternative non-linear integral equation (Destri-de Vega equation) description of the model is given. At the special value β 2 = 6π of the sine-Gordon coupling, exact expressions for energy and momentum eigenvalues of one-particle states are found

Integral equations for free-molecule ow in MEMS: recent advancements

Directory of Open Access Journals (Sweden)

Fedeli Patrick

2017-03-01

Full Text Available We address a Boundary Integral Equation (BIE approach for the analysis of gas dissipation in near-vacuum for Micro Electro Mechanical Systems (MEMS. Inspired by an analogy with the radiosity equation in computer graphics, we discuss an efficient way to compute the visible domain of integration. Moreover, we tackle the issue of near singular integrals by developing a set of analytical formulas for planar polyhedral domains. Finally a validation with experimental results taken from the literature is presented.

Nonlinear integral equations for the sausage model

Science.gov (United States)

Ahn, Changrim; Balog, Janos; Ravanini, Francesco

2017-08-01

The sausage model, first proposed by Fateev, Onofri, and Zamolodchikov, is a deformation of the O(3) sigma model preserving integrability. The target space is deformed from the sphere to ‘sausage’ shape by a deformation parameter ν. This model is defined by a factorizable S-matrix which is obtained by deforming that of the O(3) sigma model by a parameter λ. Clues for the deformed sigma model are provided by various UV and IR information through the thermodynamic Bethe ansatz (TBA) analysis based on the S-matrix. Application of TBA to the sausage model is, however, limited to the case of 1/λ integer where the coupled integral equations can be truncated to a finite number. In this paper, we propose a finite set of nonlinear integral equations (NLIEs), which are applicable to generic value of λ. Our derivation is based on T-Q relations extracted from the truncated TBA equations. For a consistency check, we compute next-leading order corrections of the vacuum energy and extract the S-matrix information in the IR limit. We also solved the NLIE both analytically and numerically in the UV limit to get the effective central charge and compared with that of the zero-mode dynamics to obtain exact relation between ν and λ. Dedicated to the memory of Petr Petrovich Kulish.

Feynman path integral related to stochastic schroedinger equation

International Nuclear Information System (INIS)

Belavkin, V.P.; Smolyanov, O.G.

1998-01-01

The derivation of the Schroedinger equation describing the continuous measurement process is presented. The representation of the solution of the stochastic Schroedinger equation for continuous measurements is obtained by means of the Feynman path integral. The connection with the heuristic approach to the description of continuous measurements is considered. The connection with the Senon paradox is established [ru

On the structure of the commutative Z2 graded algebra valued integrable equations

International Nuclear Information System (INIS)

Konopelchenko, B.G.

1980-01-01

Partial differential equations integrable by the linear matrix spectral problem of arbitrary order are considered for the case that the 'potentials' take their values in the commutative infinte-dimensional Z 2 graded algebra (superalgebra). The general form of the integrable equations and their Baecklund transformations are found. The infinite sets of the integrals of the motion are constructed. The hamiltonian character of the integrable equations is proved. (orig.)

Representations of classical groups on the lattice and its application to the field theory on discrete space-time

OpenAIRE

Lorente, M.

2003-01-01

We explore the mathematical consequences of the assumption of a discrete space-time. The fundamental laws of physics have to be translated into the language of discrete mathematics. We find integral transformations that leave the lattice of any dimension invariant and apply these transformations to field equations.

Quantum Lattice-Gas Model for the Diffusion Equation

National Research Council Canada - National Science Library

Yepez, J

2001-01-01

.... It is a minimal model with two qubits per node of a one-dimensional lattice and it is suitable for implementation on a large array of small quantum computers interconnected by nearest-neighbor...

Set-valued and fuzzy stochastic integral equations driven by semimartingales under Osgood condition

Directory of Open Access Journals (Sweden)

Malinowski Marek T.

2015-01-01

Full Text Available We analyze the set-valued stochastic integral equations driven by continuous semimartingales and prove the existence and uniqueness of solutions to such equations in the framework of the hyperspace of nonempty, bounded, convex and closed subsets of the Hilbert space L2 (consisting of square integrable random vectors. The coefficients of the equations are assumed to satisfy the Osgood type condition that is a generalization of the Lipschitz condition. Continuous dependence of solutions with respect to data of the equation is also presented. We consider equations driven by semimartingale Z and equations driven by processes A;M from decomposition of Z, where A is a process of finite variation and M is a local martingale. These equations are not equivalent. Finally, we show that the analysis of the set-valued stochastic integral equations can be extended to a case of fuzzy stochastic integral equations driven by semimartingales under Osgood type condition. To obtain our results we use the set-valued and fuzzy Maruyama type approximations and Bihari’s inequality.

Alternate Solution to Generalized Bernoulli Equations via an Integrating Factor: An Exact Differential Equation Approach

Science.gov (United States)

Tisdell, C. C.

2017-01-01

Solution methods to exact differential equations via integrating factors have a rich history dating back to Euler (1740) and the ideas enjoy applications to thermodynamics and electromagnetism. Recently, Azevedo and Valentino presented an analysis of the generalized Bernoulli equation, constructing a general solution by linearizing the problem…

Lattice gas methods for predicting intrinsic permeability of porous media

Energy Technology Data Exchange (ETDEWEB)

Santos, L.O.E.; Philippi, P.C. [Santa Catarina Univ., Florianopolis, SC (Brazil). Dept. de Engenharia Mecanica. Lab. de Propriedades Termofisicas e Meios Porosos)]. E-mail: emerich@lmpt.ufsc.br; philippi@lmpt.ufsc.br; Damiani, M.C. [Engineering Simulation and Scientific Software (ESSS), Florianopolis, SC (Brazil). Parque Tecnologico]. E-mail: damiani@lmpt.ufsc.br

2000-07-01

This paper presents a method for predicting intrinsic permeability of porous media based on Lattice Gas Cellular Automata methods. Two methods are presented. The first is based on a Boolean model (LGA). The second is Boltzmann method (LB) based on Boltzmann relaxation equation. LGA is a relatively recent method developed to perform hydrodynamic calculations. The method, in its simplest form, consists of a regular lattice populated with particles that hop from site to site in discrete time steps in a process, called propagation. After propagation, the particles in each site interact with each other in a process called collision, in which the number of particles and momentum are conserved. An exclusion principle is imposed in order to achieve better computational efficiency. In despite of its simplicity, this model evolves in agreement with Navier-Stokes equation for low Mach numbers. LB methods were recently developed for the numerical integration of the Navier-Stokes equation based on discrete Boltzmann transport equation. Derived from LGA, LB is a powerful alternative to the standard methods in computational fluid dynamics. In recent years, it has received much attention and has been used in several applications like simulations of flows through porous media, turbulent flows and multiphase flows. It is important to emphasize some aspects that make Lattice Gas Cellular Automata methods very attractive for simulating flows through porous media. In fact, boundary conditions in flows through complex geometry structures are very easy to describe in simulations using these methods. In LGA methods simulations are performed with integers needing less resident memory capability and boolean arithmetic reduces running time. The two methods are used to simulate flows through several Brazilian reservoir petroleum rocks leading to intrinsic permeability prediction. Simulation is compared with experimental results. (author)

ICM: an Integrated Compartment Method for numerically solving partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Yeh, G.T.

1981-05-01

An integrated compartment method (ICM) is proposed to construct a set of algebraic equations from a system of partial differential equations. The ICM combines the utility of integral formulation of finite element approach, the simplicity of interpolation of finite difference approximation, and the flexibility of compartment analyses. The integral formulation eases the treatment of boundary conditions, in particular, the Neumann-type boundary conditions. The simplicity of interpolation provides great economy in computation. The flexibility of discretization with irregular compartments of various shapes and sizes offers advantages in resolving complex boundaries enclosing compound regions of interest. The basic procedures of ICM are first to discretize the region of interest into compartments, then to apply three integral theorems of vectors to transform the volume integral to the surface integral, and finally to use interpolation to relate the interfacial values in terms of compartment values to close the system. The Navier-Stokes equations are used as an example of how to derive the corresponding ICM alogrithm for a given set of partial differential equations. Because of the structure of the algorithm, the basic computer program remains the same for cases in one-, two-, or three-dimensional problems.

A Concept Lattice for Semantic Integration of Geo-Ontologies Based on Weight of Inclusion Degree Importance and Information Entropy

Directory of Open Access Journals (Sweden)

Jia Xiao

2016-11-01

Full Text Available Constructing a merged concept lattice with formal concept analysis (FCA is an important research direction in the field of integrating multi-source geo-ontologies. Extracting essential geographical properties and reducing the concept lattice are two key points of previous research. A formal integration method is proposed to address the challenges in these two areas. We first extract essential properties from multi-source geo-ontologies and use FCA to build a merged formal context. Second, the combined importance weight of each single attribute of the formal context is calculated by introducing the inclusion degree importance from rough set theory and information entropy; then a weighted formal context is built from the merged formal context. Third, a combined weighted concept lattice is established from the weighted formal context with FCA and the importance weight value of every concept is defined as the sum of weight of attributes belonging to the concept’s intent. Finally, semantic granularity of concept is defined by its importance weight; we, then gradually reduce the weighted concept lattice by setting up diminishing threshold of semantic granularity. Additionally, all of those reduced lattices are organized into a regular hierarchy structure based on the threshold of semantic granularity. A workflow is designed to demonstrate this procedure. A case study is conducted to show feasibility and validity of this method and the procedure to integrate multi-source geo-ontologies.

Multi-symplectic variational integrators for nonlinear Schrödinger equations with variable coefficients

International Nuclear Information System (INIS)

Liao Cui-Cui; Cui Jin-Chao; Liang Jiu-Zhen; Ding Xiao-Hua

2016-01-01

In this paper, we propose a variational integrator for nonlinear Schrödinger equations with variable coefficients. It is shown that our variational integrator is naturally multi-symplectic. The discrete multi-symplectic structure of the integrator is presented by a multi-symplectic form formula that can be derived from the discrete Lagrangian boundary function. As two examples of nonlinear Schrödinger equations with variable coefficients, cubic nonlinear Schrödinger equations and Gross–Pitaevskii equations are extensively studied by the proposed integrator. Our numerical simulations demonstrate that the integrator is capable of preserving the mass, momentum, and energy conservation during time evolutions. Convergence tests are presented to verify that our integrator has second-order accuracy both in time and space. (paper)

On the Asymptotic Behavior of the Kernel Function in the Generalized Langevin Equation: A One-Dimensional Lattice Model

Science.gov (United States)

Chu, Weiqi; Li, Xiantao

2018-01-01

We present some estimates for the memory kernel function in the generalized Langevin equation, derived using the Mori-Zwanzig formalism from a one-dimensional lattice model, in which the particles interactions are through nearest and second nearest neighbors. The kernel function can be explicitly expressed in a matrix form. The analysis focuses on the decay properties, both spatially and temporally, revealing a power-law behavior in both cases. The dependence on the level of coarse-graining is also studied.

Toward lattice fractional vector calculus

International Nuclear Information System (INIS)

Tarasov, Vasily E

2014-01-01

An analog of fractional vector calculus for physical lattice models is suggested. We use an approach based on the models of three-dimensional lattices with long-range inter-particle interactions. The lattice analogs of fractional partial derivatives are represented by kernels of lattice long-range interactions, where the Fourier series transformations of these kernels have a power-law form with respect to wave vector components. In the continuum limit, these lattice partial derivatives give derivatives of non-integer order with respect to coordinates. In the three-dimensional description of the non-local continuum, the fractional differential operators have the form of fractional partial derivatives of the Riesz type. As examples of the applications of the suggested lattice fractional vector calculus, we give lattice models with long-range interactions for the fractional Maxwell equations of non-local continuous media and for the fractional generalization of the Mindlin and Aifantis continuum models of gradient elasticity. (papers)

Toward lattice fractional vector calculus

Science.gov (United States)

Tarasov, Vasily E.

2014-09-01

An analog of fractional vector calculus for physical lattice models is suggested. We use an approach based on the models of three-dimensional lattices with long-range inter-particle interactions. The lattice analogs of fractional partial derivatives are represented by kernels of lattice long-range interactions, where the Fourier series transformations of these kernels have a power-law form with respect to wave vector components. In the continuum limit, these lattice partial derivatives give derivatives of non-integer order with respect to coordinates. In the three-dimensional description of the non-local continuum, the fractional differential operators have the form of fractional partial derivatives of the Riesz type. As examples of the applications of the suggested lattice fractional vector calculus, we give lattice models with long-range interactions for the fractional Maxwell equations of non-local continuous media and for the fractional generalization of the Mindlin and Aifantis continuum models of gradient elasticity.

SPLET - A program for calculating the space-lethargy distribution of epithermal neutrons in a reactor lattice cell

International Nuclear Information System (INIS)

Matausek, M.V.; Zmijatevic, I.

1981-01-01

A procedure to solve the space-single-lethargy dependent transport equation for epithermal neutrons in a cylindricised multi-region reactor lattice cell has been developed and proposed in the earlier papers. Here, the computational algorithm is comprised and the computing program SPLET, which calculates the space-lethargy distribution of the spherical harmonics neutron flux moments, as well as the related integral quantities as reaction rates and resonance integrals, is described. (author)

Localization of the eigenvalues of linear integral equations with applications to linear ordinary differential equations.

Science.gov (United States)

Sloss, J. M.; Kranzler, S. K.

1972-01-01

The equivalence of a considered integral equation form with an infinite system of linear equations is proved, and the localization of the eigenvalues of the infinite system is expressed. Error estimates are derived, and the problems of finding upper bounds and lower bounds for the eigenvalues are solved simultaneously.

Coupling-governed metamorphoses of the integrable nonlinear Schrödinger system on a triangular-lattice ribbon

Energy Technology Data Exchange (ETDEWEB)

Vakhnenko, Oleksiy O., E-mail: vakhnenko@bitp.kiev.ua

2016-05-27

Highlights: • The integrable nonlinear Schrödinger system on a triangular-lattice ribbon is inclined to metamorphoses. • The system under study is capable to incorporate the effect of external linear potential. • The system criticality against the background parameter reduces the number of independent field variables. • At critical point the system Poisson structure becomes degenerate. • The effect of criticality is elucidated by the system one-soliton solution. - Abstract: The variativity of governing coupling parameters in the integrable nonlinear Schrödinger system on a triangular-lattice ribbon is shown to ensure the important qualitative rearrangements in the system dynamics. There are at least the two types of system crucial modifications stipulated by the two types of governing parameters. Thus the longitudinal coupling parameters regulated mainly by the background values of concomitant field variables are responsible for the bifurcation of primary integrable nonlinear system into the integrable nonlinear system of Ablowitz–Ladik type. As a consequence in a critical point the number of independent field variables is reduced by a half and the system Poisson structure turns out to be degenerate. On the other hand the transverse coupling parameters regulated basically by the choice of their a priori arbitrary dependencies on time are capable to incorporate the effect of external linear potential. As a consequence the primary integrable nonlinear system with appropriately adjusted parametrical driving becomes isomorphic to the system modeling the Bloch oscillations of charged nonlinear carriers in an electrically biased ribbon of triangular lattice. The multi-component structure of basic integrable system alongside with the attractive character of system nonlinearities has predetermined the logic of whole consideration.

Integral equation approach to time-dependent kinematic dynamos in finite domains

International Nuclear Information System (INIS)

Xu Mingtian; Stefani, Frank; Gerbeth, Gunter

2004-01-01

The homogeneous dynamo effect is at the root of cosmic magnetic field generation. With only a very few exceptions, the numerical treatment of homogeneous dynamos is carried out in the framework of the differential equation approach. The present paper tries to facilitate the use of integral equations in dynamo research. Apart from the pedagogical value to illustrate dynamo action within the well-known picture of the Biot-Savart law, the integral equation approach has a number of practical advantages. The first advantage is its proven numerical robustness and stability. The second and perhaps most important advantage is its applicability to dynamos in arbitrary geometries. The third advantage is its intimate connection to inverse problems relevant not only for dynamos but also for technical applications of magnetohydrodynamics. The paper provides the first general formulation and application of the integral equation approach to time-dependent kinematic dynamos, with stationary dynamo sources, in finite domains. The time dependence is restricted to the magnetic field, whereas the velocity or corresponding mean-field sources of dynamo action are supposed to be stationary. For the spherically symmetric α 2 dynamo model it is shown how the general formulation is reduced to a coupled system of two radial integral equations for the defining scalars of the poloidal and toroidal field components. The integral equation formulation for spherical dynamos with general stationary velocity fields is also derived. Two numerical examples - the α 2 dynamo model with radially varying α and the Bullard-Gellman model - illustrate the equivalence of the approach with the usual differential equation method. The main advantage of the method is exemplified by the treatment of an α 2 dynamo in rectangular domains

Green function of the double-fractional Fokker-Planck equation: Path integral and stochastic differential equations

Science.gov (United States)

Kleinert, H.; Zatloukal, V.

2013-11-01

The statistics of rare events, the so-called black-swan events, is governed by non-Gaussian distributions with heavy power-like tails. We calculate the Green functions of the associated Fokker-Planck equations and solve the related stochastic differential equations. We also discuss the subject in the framework of path integration.

On the hierarchical lattices approximation of Bravais lattices: Specific heat and correlation length

International Nuclear Information System (INIS)

Tsallis, C.

1984-01-01

Certain types of real-space renormalization groups (which essentially approximate Bravais lattices through hierarchical ones) do not preserve standard thermodynamic convexity properties. It is pointed out that this serious defect is not intrinsic to any real-space renormalization. It can be avoided if form-invariance (under uniform translation of the energy scale) of the equation connecting the Bravais lattice (which is intended to study) to the hierarchical one (which approximates it) is demanded. In addition to that expressions for the critical exponentes ν and α corresponding to hierarchical lattices are analysed; these are consistent with Melrose recent analysis of the fractal intrinsic dimensionality. (Author) [pt

Integral transport multiregion geometrical shadowing factor for the approximate collision probability matrix calculation of infinite closely packed lattices

International Nuclear Information System (INIS)

Jowzani-Moghaddam, A.

1981-01-01

An integral transport method of calculating the geometrical shadowing factor in multiregion annular cells for infinite closely packed lattices in cylindrical geometry is developed. This analytical method has been programmed in the TPGS code. This method is based upon a consideration of the properties of the integral transport method for a nonuniform body, which together with Bonalumi's approximations allows the determination of the approximate multiregion collision probability matrix for infinite closely packed lattices with sufficient accuracy. The multiregion geometrical shadowing factors have been calculated for variations in fuel pin annular segment rings in a geometry of annular cells. These shadowing factors can then be used in the calculation of neutron transport from one annulus to another in an infinite lattice. The result of this new geometrical shadowing and collision probability matrix are compared with the Dancoff-Ginsburg correction and the probability matrix using constant shadowing on Yankee fuel elements in an infinite lattice. In these cases the Dancoff-Ginsburg correction factor and collision probability matrix using constant shadowing are in difference by at most 6.2% and 6%, respectively

Integrability and structural stability of solutions to the Ginzburg-Landau equation

Science.gov (United States)

Keefe, Laurence R.

1986-01-01

The integrability of the Ginzburg-Landau equation is studied to investigate if the existence of chaotic solutions found numerically could have been predicted a priori. The equation is shown not to possess the Painleveproperty, except for a special case of the coefficients that corresponds to the integrable, nonlinear Schroedinger (NLS) equation. Regarding the Ginzburg-Landau equation as a dissipative perturbation of the NLS, numerical experiments show all but one of a family of two-tori solutions, possessed by the NLS under particular conditions, to disappear under real perturbations to the NLS coefficients of O(10 to the -6th).

The sl (2 n|2 n) (1) super-Toda lattices and the heavenly equations as continuum limit

International Nuclear Information System (INIS)

Kuznetsova, Zhanna; Popowicz, Ziemowit; Toppan, Francesco.

2005-04-01

The n → ∞ continuum limit of super-Toda models associated with the affine sl (2 n|2 n) (1) (super)algebra series produces (2 + 1)-dimensional integrable equations in the S 1 versus R 2 spacetimes. The equations of motion of the (super)Toda hierarchies depend not only on the chosen (super)algebras but also on the specific presentation of their Cartan matrices. Four distinct series of integrable hierarchies in relation with symmetric-versus-antisymmetric, null-versus-non null presentations of the corresponding Cartan matrices are investigated. In the continuum limit we derive four classes of integrable equations of heavenly type, generalizing the results previously obtained in the literature. The systems are manifestly N = 1 supersymmetric and, for specific choices of the Cartan matrix preserving the complex structure, admit a hidden N = 2 supersymmetry. The coset reduction of the (super)-heavenly equation to the I versus R (2) (S 1 /Z 2 ) versus R 2 spacetime (with I a line segment) is illustrated. Finally, integrable N = 2,4 super symmetrically extended models in (1 + 1) dimensions are constructed through dimensional reduction of the previous systems. (author)

Nonlinear Kinetics on Lattices Based on the Kinetic Interaction Principle

Directory of Open Access Journals (Sweden)

Giorgio Kaniadakis

2018-06-01

Full Text Available Master equations define the dynamics that govern the time evolution of various physical processes on lattices. In the continuum limit, master equations lead to Fokker–Planck partial differential equations that represent the dynamics of physical systems in continuous spaces. Over the last few decades, nonlinear Fokker–Planck equations have become very popular in condensed matter physics and in statistical physics. Numerical solutions of these equations require the use of discretization schemes. However, the discrete evolution equation obtained by the discretization of a Fokker–Planck partial differential equation depends on the specific discretization scheme. In general, the discretized form is different from the master equation that has generated the respective Fokker–Planck equation in the continuum limit. Therefore, the knowledge of the master equation associated with a given Fokker–Planck equation is extremely important for the correct numerical integration of the latter, since it provides a unique, physically motivated discretization scheme. This paper shows that the Kinetic Interaction Principle (KIP that governs the particle kinetics of many body systems, introduced in G. Kaniadakis, Physica A 296, 405 (2001, univocally defines a very simple master equation that in the continuum limit yields the nonlinear Fokker–Planck equation in its most general form.

One-way spatial integration of hyperbolic equations

Science.gov (United States)

Towne, Aaron; Colonius, Tim

2015-11-01

In this paper, we develop and demonstrate a method for constructing well-posed one-way approximations of linear hyperbolic systems. We use a semi-discrete approach that allows the method to be applied to a wider class of problems than existing methods based on analytical factorization of idealized dispersion relations. After establishing the existence of an exact one-way equation for systems whose coefficients do not vary along the axis of integration, efficient approximations of the one-way operator are constructed by generalizing techniques previously used to create nonreflecting boundary conditions. When physically justified, the method can be applied to systems with slowly varying coefficients in the direction of integration. To demonstrate the accuracy and computational efficiency of the approach, the method is applied to model problems in acoustics and fluid dynamics via the linearized Euler equations; in particular we consider the scattering of sound waves from a vortex and the evolution of hydrodynamic wavepackets in a spatially evolving jet. The latter problem shows the potential of the method to offer a systematic, convergent alternative to ad hoc regularizations such as the parabolized stability equations.

Integrable discretizations of the (2+1)-dimensional sinh-Gordon equation

International Nuclear Information System (INIS)

Hu, Xing-Biao; Yu, Guo-Fu

2007-01-01

In this paper, we propose two semi-discrete equations and one fully discrete equation and study them by Hirota's bilinear method. These equations have continuum limits into a system which admits the (2+1)-dimensional generalization of the sinh-Gordon equation. As a result, two integrable semi-discrete versions and one fully discrete version for the sinh-Gordon equation are found. Baecklund transformations, nonlinear superposition formulae, determinant solution and Lax pairs for these discrete versions are presented

Lattice formulation of a two-dimensional topological field theory

International Nuclear Information System (INIS)

Ohta, Kazutoshi; Takimi, Tomohisa

2007-01-01

We investigate an integrable property and the observables of 2-dimensional N=(4,4) topological field theory defined on a discrete lattice by using the 'orbifolding' and 'deconstruction' methods. We show that our lattice model is integrable and, for this reason, the partition function reduces to matrix integrals of scalar fields on the lattice sites. We elucidate meaningful differences between a discrete lattice and a differentiable manifold. This is important for studying topological quantities on a lattice. We also propose a new construction of N=(2,2) supersymmetric lattice theory, which is realized through a suitable truncation of scalar fields from the N=(4,4) theory. (author)

A New Algorithm for System of Integral Equations

Directory of Open Access Journals (Sweden)

Abdujabar Rasulov

2014-01-01

Full Text Available We develop a new algorithm to solve the system of integral equations. In this new method no need to use matrix weights. Beacause of it, we reduce computational complexity considerable. Using the new algorithm it is also possible to solve an initial boundary value problem for system of parabolic equations. To verify the efficiency, the results of computational experiments are given.

On integration of the first order differential equations in a finite terms

International Nuclear Information System (INIS)

Malykh, M D

2017-01-01

There are several approaches to the description of the concept called briefly as integration of the first order differential equations in a finite terms or symbolical integration. In the report three of them are considered: 1.) finding of a rational integral (Beaune or Poincaré problem), 2.) integration by quadratures and 3.) integration when the general solution of given differential equation is an algebraical function of a constant (Painlevé problem). Their realizations in Sage are presented. (paper)

Reduction of lattice equations to the Painlevé equations: PIV and PV

Science.gov (United States)

Nakazono, Nobutaka

2018-02-01

In this paper, we construct a new relation between Adler-Bobenko-Suris equations and Painlevé equations. Moreover, using this connection we construct the difference-differential Lax representations of the fourth and fifth Painlevé equations.

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.

Lattice constant measurement from electron backscatter diffraction patterns

DEFF Research Database (Denmark)

Saowadee, Nath; Agersted, Karsten; Bowen, Jacob R.

2017-01-01

Kikuchi bands in election backscattered diffraction patterns (EBSP) contain information about lattice constants of crystallographic samples that can be extracted via the Bragg equation. An advantage of lattice constant measurement from EBSPs over diffraction (XRD) is the ability to perform local ...

Spherical Harmonics Treatment of Epithermal Neutron Spectra in Reactor lattices

International Nuclear Information System (INIS)

Matausek, M.V.

1972-04-01

A procedure has been developed to solve the slowing down transport equation for neutrons in a cylindrized reactor lattice cell. Treating the anisotropy of the epithermal neutron flux by the spherical harmonics formalism, which reduces the space-angle-lethargy-dependent transport equation to the matrix integrodifferential equation in space and lethargy, and replacing the lethargy transfer integrals by finite-difference forms, a set of matrix ordinary differential equations, with lethargy and space dependent coefficients, is obtained. In the resonance region this set takes a lower block triangular form and can be directly solved by forward block substitution; in the lethargy range, where the fast fission effects have to be considered, the iterative procedure is introduced. A simple and efficient approximation is then proposed, making possible the analytical solution for the spatial dependence of the spherical harmonics flux moments. The proposed procedure has been numerically examined and approved. Some typical results are presented and discussed. (author)

Thermodynamically self-consistent integral equations and the structure of liquid metals

International Nuclear Information System (INIS)

Pastore, G.; Kahl, G.

1987-01-01

We discuss the application of the new thermodynamically self-consistent integral equations for the determination of the structural properties of liquid metals. We present a detailed comparison of the structure (S(q) and g(r)) for models of liquid alkali metals as obtained from two thermodynamically self-consistent integral equations and some published exact computer simulation results; the range of states extends from the triple point to the expanded metal. The theories which only impose thermodynamic self-consistency without any fitting of external data show an excellent agreement with the simulation results, thus demonstrating that this new type of integral equation is definitely superior to the conventional ones (hypernetted chain, Percus-Yevick, mean spherical approximation, etc). (author)

Constructing soliton solutions and super-bilinear form of lattice supersymmetric KdV equation

International Nuclear Information System (INIS)

Carstea, A S

2015-01-01

The Hirota bilinear form and multisoliton solution for semidiscrete and fully discrete (difference–difference) versions of the supersymmetric Korteweg–de Vries (KdV) equation found by Xue et al (2013 J. Phys. A: Math. Theor 46 502001) are presented. The solitonic interaction term displays a fermionic dressing factor as in the continuous supersymmetric case. Using bilinear equations it is also shown that a new integrable semidiscrete (and fully discrete) version of supersymmetric KdV can be constructed with a simpler bilinear form but a more complicated interaction dressing. Its continuum limit is also computed. (paper)

Distribution theory for Schrödinger’s integral equation

NARCIS (Netherlands)

Lange, R.J.

2015-01-01

Much of the literature on point interactions in quantum mechanics has focused on the differential form of Schrödinger's equation. This paper, in contrast, investigates the integral form of Schrödinger's equation. While both forms are known to be equivalent for smooth potentials, this is not true for

Unconditionally stable integration of Maxwell’s equations

NARCIS (Netherlands)

Verwer, J.G.; Botchev, M.A.

2009-01-01

Numerical integration of Maxwell’s equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction implicit - finite

A new integral method for solving the point reactor neutron kinetics equations

International Nuclear Information System (INIS)

Li Haofeng; Chen Wenzhen; Luo Lei; Zhu Qian

2009-01-01

A numerical integral method that efficiently provides the solution of the point kinetics equations by using the better basis function (BBF) for the approximation of the neutron density in one time step integrations is described and investigated. The approach is based on an exact analytic integration of the neutron density equation, where the stiffness of the equations is overcome by the fully implicit formulation. The procedure is tested by using a variety of reactivity functions, including step reactivity insertion, ramp input and oscillatory reactivity changes. The solution of the better basis function method is compared to other analytical and numerical solutions of the point reactor kinetics equations. The results show that selecting a better basis function can improve the efficiency and accuracy of this integral method. The better basis function method can be used in real time forecasting for power reactors in order to prevent reactivity accidents.

Lagrangian structures, integrability and chaos for 3D dynamical equations

International Nuclear Information System (INIS)

Bustamante, Miguel D; Hojman, Sergio A

2003-01-01

In this paper, we consider the general setting for constructing action principles for three-dimensional first-order autonomous equations. We present the results for some integrable and non-integrable cases of the Lotka-Volterra equation, and show Lagrangian descriptions which are valid for systems satisfying Shil'nikov criteria on the existence of strange attractors, though chaotic behaviour has not been verified up to now. The Euler-Lagrange equations we get for these systems usually present 'time reparametrization' invariance, though other kinds of invariance may be found according to the kernel of the associated symplectic 2-form. The formulation of a Hamiltonian structure (Poisson brackets and Hamiltonians) for these systems from the Lagrangian viewpoint leads to a method of finding new constants of the motion starting from known ones, which is applied to some systems found in the literature known to possess a constant of the motion, to find the other and thus showing their integrability. In particular, we show that the so-called ABC system is completely integrable if it possesses one constant of the motion

Universal and integrable nonlinear evolution systems of equations in 2+1 dimensions

International Nuclear Information System (INIS)

Maccari, A.

1997-01-01

Integrable systems of nonlinear partial differential equations (PDEs) are obtained from integrable equations in 2+1 dimensions, by means of a reduction method of broad applicability based on Fourier expansion and spatio endash temporal rescalings, which is asymptotically exact in the limit of weak nonlinearity. The integrability by the spectral transform is explicitly demonstrated, because the corresponding Lax pairs have been derived, applying the same reduction method to the Lax pair of the initial equation. These systems of nonlinear PDEs are likely to be of applicative relevance and have a open-quotes universalclose quotes character, inasmuch as they may be derived from a very large class of nonlinear evolution equations with a linear dispersive part. copyright 1997 American Institute of Physics

Evolution equations for connected and disconnected sea parton distributions

Science.gov (United States)

Liu, Keh-Fei

2017-08-01

It has been revealed from the path-integral formulation of the hadronic tensor that there are connected sea and disconnected sea partons. The former is responsible for the Gottfried sum rule violation primarily and evolves the same way as the valence. Therefore, the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equations can be extended to accommodate them separately. We discuss its consequences and implications vis-á-vis lattice calculations.

Statistical mechanics view of quantum chromodynamics: Lattice gauge theory

International Nuclear Information System (INIS)

Kogut, J.B.

1984-01-01

Recent developments in lattice gauge theory are discussed from a statistial mechanics viewpoint. The basic physics problems of quantum chromodynamics (QCD) are reviewed for an audience of critical phenomena theorists. The idea of local gauge symmetry and color, the connection between statistical mechanics and field theory, asymptotic freedom and the continuum limit of lattice gauge theories, and the order parameters (confinement and chiral symmetry) of QCD are reviewed. Then recent developments in the field are discussed. These include the proof of confinement in the lattice theory, numerical evidence for confinement in the continuum limit of lattice gauge theory, and perturbative improvement programs for lattice actions. Next, we turn to the new challenges facing the subject. These include the need for a better understanding of the lattice Dirac equation and recent progress in the development of numerical methods for fermions (the pseudofermion stochastic algorithm and the microcanonical, molecular dynamics equation of motion approach). Finally, some of the applications of lattice gauge theory to QCD spectrum calculations and the thermodynamics of QCD will be discussed and a few remarks concerning future directions of the field will be made

Integration Processes of Delay Differential Equation Based on Modified Laguerre Functions

Directory of Open Access Journals (Sweden)

Yeguo Sun

2012-01-01

Full Text Available We propose long-time convergent numerical integration processes for delay differential equations. We first construct an integration process based on modified Laguerre functions. Then we establish its global convergence in certain weighted Sobolev space. The proposed numerical integration processes can also be used for systems of delay differential equations. We also developed a technique for refinement of modified Laguerre-Radau interpolations. Lastly, numerical results demonstrate the spectral accuracy of the proposed method and coincide well with analysis.

Breatherlike impurity modes in discrete nonlinear lattices

DEFF Research Database (Denmark)

Hennig, D.; Rasmussen, Kim; Tsironis, G. P.

1995-01-01

We investigate the properties of a disordered generalized discrete nonlinear Schrodinger equation, containing both diagonal and nondiagonal nonlinear terms. The equation models a Linear host lattice doped with nonlinear impurities. We find different types of impurity states that form itinerant...

Stability of negative solitary waves for an integrable modified Camassa-Holm equation

International Nuclear Information System (INIS)

Yin Jiuli; Tian Lixin; Fan Xinghua

2010-01-01

In this paper, we prove that the modified Camassa-Holm equation is Painleve integrable. We also study the orbital stability problem of negative solitary waves for this integrable equation. It is shown that the negative solitary waves are stable for arbitrary wave speed of propagation.

Fibonacci-regularization method for solving Cauchy integral equations of the first kind

Directory of Open Access Journals (Sweden)

Mohammad Ali Fariborzi Araghi

2017-09-01

Full Text Available In this paper, a novel scheme is proposed to solve the first kind Cauchy integral equation over a finite interval. For this purpose, the regularization method is considered. Then, the collocation method with Fibonacci base function is applied to solve the obtained second kind singular integral equation. Also, the error estimate of the proposed scheme is discussed. Finally, some sample Cauchy integral equations stem from the theory of airfoils in fluid mechanics are presented and solved to illustrate the importance and applicability of the given algorithm. The tables in the examples show the efficiency of the method.

Field Method for Integrating the First Order Differential Equation

Institute of Scientific and Technical Information of China (English)

JIA Li-qun; ZHENG Shi-wang; ZHANG Yao-yu

2007-01-01

An important modern method in analytical mechanics for finding the integral, which is called the field-method, is used to research the solution of a differential equation of the first order. First, by introducing an intermediate variable, a more complicated differential equation of the first order can be expressed by two simple differential equations of the first order, then the field-method in analytical mechanics is introduced for solving the two differential equations of the first order. The conclusion shows that the field-method in analytical mechanics can be fully used to find the solutions of a differential equation of the first order, thus a new method for finding the solutions of the first order is provided.

Lattice solitons in Bose-Einstein condensates

International Nuclear Information System (INIS)

Efremidis, Nikolaos K.; Christodoulides, Demetrios N.

2003-01-01

We systematically study the properties of lattice solitons in Bose-Einstein condensates with either attractive or repulsive atom interactions. This is done, by exactly solving the mean-field Gross-Pitaevskii equation in the presence of a periodic potential. We find new families of lattice soliton solutions that are characterized by the position of the energy eigenvalue within the associated band structure. These include lattice solitons in condensates with either attractive or repulsive atom interactions that exist in finite or semi-infinite gaps, as well as nonlinear modes that exhibit atomic population cutoffs

Electronic properties of graphene antidot lattices

DEFF Research Database (Denmark)

Fürst, Joachim Alexander; Pedersen, Jesper Goor; Flindt, C.

2009-01-01

Graphene antidot lattices constitute a novel class of nano-engineered graphene devices with controllable electronic and optical properties. An antidot lattice consists of a periodic array of holes that causes a band gap to open up around the Fermi level, turning graphene from a semimetal...... into a semiconductor. We calculate the electronic band structure of graphene antidot lattices using three numerical approaches with different levels of computational complexity, efficiency and accuracy. Fast finite-element solutions of the Dirac equation capture qualitative features of the band structure, while full...

On an integrable deformed affinsphären equation. A reciprocal gasdynamic connection

International Nuclear Information System (INIS)

Rogers, C.; Huang, Yehui

2012-01-01

The integrable affinsphären equation originally arose in a geometric context but has an interesting gasdynamic connection. Here, an integrable deformed version of the affinsphären equation is derived in a novel manner via the action of reciprocal transformations on a related anisentropic gasdynamics system. A linear representation for the deformed affinsphären equation is constructed by means of the reciprocal transformations. The latter are then employed to derive a class of exact solutions in parametric form. -- Highlights: ► A deformed affinsphären equation is derived via a reciprocal transformation. ► A linear representation for the deformed affinsphären equation is constructed. ► A class of exact solutions of the deformed affinsphären equation is presented.

Analysis and regularization of the thin-wire integral equation with reduced kernel

NARCIS (Netherlands)

Beurden, van M.C.; Tijhuis, A.G.

2007-01-01

For the straight wire, modeled as a hollow tube, we establish a conditional equivalence relation between the integral equations with exact and reduced kernel. This relation allows us to examine the existence and uniqueness conditions for the integral equation with reduced kernel, based on a local

Derivation of the Schrodinger Equation from the Hamilton-Jacobi Equation in Feynman's Path Integral Formulation of Quantum Mechanics

Science.gov (United States)

Field, J. H.

2011-01-01

It is shown how the time-dependent Schrodinger equation may be simply derived from the dynamical postulate of Feynman's path integral formulation of quantum mechanics and the Hamilton-Jacobi equation of classical mechanics. Schrodinger's own published derivations of quantum wave equations, the first of which was also based on the Hamilton-Jacobi…

Application of perturbation theory to lattice calculations based on method of cyclic characteristics

Science.gov (United States)

Assawaroongruengchot, Monchai

Perturbation theory is a technique used for the estimation of changes in performance functionals, such as linear reaction rate ratio and eigenvalue affected by small variations in reactor core compositions. Here the algorithm of perturbation theory is developed for the multigroup integral neutron transport problems in 2D fuel assemblies with isotropic scattering. The integral transport equation is used in the perturbative formulation because it represents the interconnecting neutronic systems of the lattice assemblies via the tracking lines. When the integral neutron transport equation is used in the formulation, one needs to solve the resulting integral transport equations for the flux importance and generalized flux importance functions. The relationship between the generalized flux importance and generalized source importance functions is defined in order to transform the generalized flux importance transport equations into the integro-differential equations for the generalized adjoints. Next we develop the adjoint and generalized adjoint transport solution algorithms based on the method of cyclic characteristics (MOCC) in DRAGON code. In the MOCC method, the adjoint characteristics equations associated with a cyclic tracking line are formulated in such a way that a closed form for the adjoint angular function can be obtained. The MOCC method then requires only one cycle of scanning over the cyclic tracking lines in each spatial iteration. We also show that the source importance function by CP method is mathematically equivalent to the adjoint function by MOCC method. In order to speed up the MOCC solution algorithm, a group-reduction and group-splitting techniques based on the structure of the adjoint scattering matrix are implemented. A combined forward flux/adjoint function iteration scheme, based on the group-splitting technique and the common use of a large number of variables storing tracking-line data and exponential values, is proposed to reduce the

Application of perturbation theory to lattice calculations based on method of cyclic characteristics

Energy Technology Data Exchange (ETDEWEB)

Assawaroongruengchot, M

2007-07-01

Perturbation theory is a technique used for the estimation of changes in performance functionals, such as linear reaction rate ratio and eigenvalue affected by small variations in reactor core compositions. Here the algorithm of perturbation theory is developed for the multigroup integral neutron transport problems in 2D fuel assemblies with isotropic scattering. The integral transport equation is used in the perturbative formulation because it represents the interconnecting neutronic systems of the lattice assemblies via the tracking lines. When the integral neutron transport equation is used in the formulation, one needs to solve the resulting integral transport equations for the flux importance and generalized flux importance functions. The relationship between the generalized flux importance and generalized source importance functions is defined in order to transform the generalized flux importance transport equations into the integro-differential equations for the generalized adjoints. Next we develop the adjoint and generalized adjoint transport solution algorithms based on the method of cyclic characteristics (MOCC) in DRAGON code. In the MOCC method, the adjoint characteristics equations associated with a cyclic tracking line are formulated in such a way that a closed form for the adjoint angular function can be obtained. The MOCC method then requires only one cycle of scanning over the cyclic tracking lines in each spatial iteration. We also show that the source importance function by CP method is mathematically equivalent to the adjoint function by MOCC method. In order to speed up the MOCC solution algorithm, a group-reduction and group-splitting techniques based on the structure of the adjoint scattering matrix are implemented. A combined forward flux/adjoint function iteration scheme, based on the group-splitting technique and the common use of a large number of variables storing tracking-line data and exponential values, is proposed to reduce the

Application of perturbation theory to lattice calculations based on method of cyclic characteristics

International Nuclear Information System (INIS)

Assawaroongruengchot, M.

2007-01-01

Perturbation theory is a technique used for the estimation of changes in performance functionals, such as linear reaction rate ratio and eigenvalue affected by small variations in reactor core compositions. Here the algorithm of perturbation theory is developed for the multigroup integral neutron transport problems in 2D fuel assemblies with isotropic scattering. The integral transport equation is used in the perturbative formulation because it represents the interconnecting neutronic systems of the lattice assemblies via the tracking lines. When the integral neutron transport equation is used in the formulation, one needs to solve the resulting integral transport equations for the flux importance and generalized flux importance functions. The relationship between the generalized flux importance and generalized source importance functions is defined in order to transform the generalized flux importance transport equations into the integro-differential equations for the generalized adjoints. Next we develop the adjoint and generalized adjoint transport solution algorithms based on the method of cyclic characteristics (MOCC) in DRAGON code. In the MOCC method, the adjoint characteristics equations associated with a cyclic tracking line are formulated in such a way that a closed form for the adjoint angular function can be obtained. The MOCC method then requires only one cycle of scanning over the cyclic tracking lines in each spatial iteration. We also show that the source importance function by CP method is mathematically equivalent to the adjoint function by MOCC method. In order to speed up the MOCC solution algorithm, a group-reduction and group-splitting techniques based on the structure of the adjoint scattering matrix are implemented. A combined forward flux/adjoint function iteration scheme, based on the group-splitting technique and the common use of a large number of variables storing tracking-line data and exponential values, is proposed to reduce the

Numerical method for solving linear Fredholm fuzzy integral equations of the second kind

Energy Technology Data Exchange (ETDEWEB)

Abbasbandy, S. [Department of Mathematics, Imam Khomeini International University, P.O. Box 288, Ghazvin 34194 (Iran, Islamic Republic of)]. E-mail: saeid@abbasbandy.com; Babolian, E. [Faculty of Mathematical Sciences and Computer Engineering, Teacher Training University, Tehran 15618 (Iran, Islamic Republic of); Alavi, M. [Department of Mathematics, Arak Branch, Islamic Azad University, Arak 38135 (Iran, Islamic Republic of)

2007-01-15

In this paper we use parametric form of fuzzy number and convert a linear fuzzy Fredholm integral equation to two linear system of integral equation of the second kind in crisp case. We can use one of the numerical method such as Nystrom and find the approximation solution of the system and hence obtain an approximation for fuzzy solution of the linear fuzzy Fredholm integral equations of the second kind. The proposed method is illustrated by solving some numerical examples.

Simulation electromagnetic scattering on bodies through integral equation and neural networks methods

Science.gov (United States)

Lvovich, I. Ya; Preobrazhenskiy, A. P.; Choporov, O. N.

2018-05-01

The paper deals with the issue of electromagnetic scattering on a perfectly conducting diffractive body of a complex shape. Performance calculation of the body scattering is carried out through the integral equation method. Fredholm equation of the second time was used for calculating electric current density. While solving the integral equation through the moments method, the authors have properly described the core singularity. The authors determined piecewise constant functions as basic functions. The chosen equation was solved through the moments method. Within the Kirchhoff integral approach it is possible to define the scattered electromagnetic field, in some way related to obtained electrical currents. The observation angles sector belongs to the area of the front hemisphere of the diffractive body. To improve characteristics of the diffractive body, the authors used a neural network. All the neurons contained a logsigmoid activation function and weighted sums as discriminant functions. The paper presents the matrix of weighting factors of the connectionist model, as well as the results of the optimized dimensions of the diffractive body. The paper also presents some basic steps in calculation technique of the diffractive bodies, based on the combination of integral equation and neural networks methods.

Behaviour of the extended Toda lattice

Science.gov (United States)

Wattis, Jonathan A. D.; Gordoa, Pilar R.; Pickering, Andrew

2015-11-01

We consider the first member of an extended Toda lattice hierarchy. This system of equations is differential with respect to one independent variable and differential-delay with respect to a second independent variable. We use asymptotic analysis to consider the long wavelength limits of the system. By considering various magnitudes for the parameters involved, we derive reduced equations related to the Korteweg-de Vries and potential Boussinesq equations.

Numerical Simulation of Antennas with Improved Integral Equation Method

International Nuclear Information System (INIS)

Ma Ji; Fang Guang-You; Lu Wei

2015-01-01

Simulating antennas around a conducting object is a challenge task in computational electromagnetism, which is concerned with the behaviour of electromagnetic fields. To analyze this model efficiently, an improved integral equation-fast Fourier transform (IE-FFT) algorithm is presented in this paper. The proposed scheme employs two Cartesian grids with different size and location to enclose the antenna and the other object, respectively. On the one hand, IE-FFT technique is used to store matrix in a sparse form and accelerate the matrix-vector multiplication for each sub-domain independently. On the other hand, the mutual interaction between sub-domains is taken as the additional exciting voltage in each matrix equation. By updating integral equations several times, the whole electromagnetic system can achieve a stable status. Finally, the validity of the presented method is verified through the analysis of typical antennas in the presence of a conducting object. (paper)

The matrix realization of affine Jacobi varieties and the extended Lotka-Volterra lattice

International Nuclear Information System (INIS)

Inoue, Rei

2004-01-01

We study completely integrable Hamiltonian systems whose monodromy matrices are related to the representatives for the set of gauge equivalence classes M F of polynomial matrices. Let X be the algebraic curve given by the common characteristic equation for M F . We construct the isomorphism from the set of representatives to an affine part of the Jacobi variety of X. This variety corresponds to the invariant manifold of the system, where the Hamiltonian flow is linearized. As an application, we discuss the algebraic complete integrability of the extended Lotka-Volterra lattice with a periodic boundary condition

Dhage Iteration Method for Generalized Quadratic Functional Integral Equations

Directory of Open Access Journals (Sweden)

Bapurao C. Dhage

2015-01-01

Full Text Available In this paper we prove the existence as well as approximations of the solutions for a certain nonlinear generalized quadratic functional integral equation. An algorithm for the solutions is developed and it is shown that the sequence of successive approximations starting at a lower or upper solution converges monotonically to the solutions of related quadratic functional integral equation under some suitable mixed hybrid conditions. We rely our main result on Dhage iteration method embodied in a recent hybrid fixed point theorem of Dhage (2014 in partially ordered normed linear spaces. An example is also provided to illustrate the abstract theory developed in the paper.

TBA-like integral equations from quantized mirror curves

Energy Technology Data Exchange (ETDEWEB)

Okuyama, Kazumi [Department of Physics, Shinshu University,Matsumoto 390-8621 (Japan); Zakany, Szabolcs [Département de Physique Théorique, Université de Genève,Genève, CH-1211 (Switzerland)

2016-03-15

Quantizing the mirror curve of certain toric Calabi-Yau (CY) three-folds leads to a family of trace class operators. The resolvent function of these operators is known to encode topological data of the CY. In this paper, we show that in certain cases, this resolvent function satisfies a system of non-linear integral equations whose structure is very similar to the Thermodynamic Bethe Ansatz (TBA) systems. This can be used to compute spectral traces, both exactly and as a semiclassical expansion. As a main example, we consider the system related to the quantized mirror curve of local ℙ{sup 2}. According to a recent proposal, the traces of this operator are determined by the refined BPS indices of the underlying CY. We use our non-linear integral equations to test that proposal.

TBA-like integral equations from quantized mirror curves

Science.gov (United States)

Okuyama, Kazumi; Zakany, Szabolcs

2016-03-01

Quantizing the mirror curve of certain toric Calabi-Yau (CY) three-folds leads to a family of trace class operators. The resolvent function of these operators is known to encode topological data of the CY. In this paper, we show that in certain cases, this resolvent function satisfies a system of non-linear integral equations whose structure is very similar to the Thermodynamic Bethe Ansatz (TBA) systems. This can be used to compute spectral traces, both exactly and as a semiclassical expansion. As a main example, we consider the system related to the quantized mirror curve of local P2. According to a recent proposal, the traces of this operator are determined by the refined BPS indices of the underlying CY. We use our non-linear integral equations to test that proposal.

Solving Abel’s Type Integral Equation with Mikusinski’s Operator of Fractional Order

Directory of Open Access Journals (Sweden)

Ming Li

2013-01-01

Full Text Available This paper gives a novel explanation of the integral equation of Abel’s type from the point of view of Mikusinski’s operational calculus. The concept of the inverse of Mikusinski’s operator of fractional order is introduced for constructing a representation of the solution to the integral equation of Abel’s type. The proof of the existence of the inverse of the fractional Mikusinski operator is presented, providing an alternative method of treating the integral equation of Abel’s type.

Integrable couplings of relativistic Toda lattice systems in polynomial form and rational form, their hierarchies and bi-Hamiltonian structures

Energy Technology Data Exchange (ETDEWEB)

Xu Xixiang [College of Science, Shandong University of Science and Technology, Qingdao 266510 (China)], E-mail: xixiang_xu@yahoo.com.cn

2009-10-02

Integrable couplings of relativistic Toda lattice systems in polynomial form and rational form, and their hierarchies, are derived from a four-by-four discrete matrix eigenvalue problem. The bi-Hamiltonian structure for every integrable coupling in the two hierarchies obtained is established by means of the discrete variational identity. Ultimately, Liouvolle integrability of the obtained integrable couplings is demonstrated.

Integrability of a system of two nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Zhukhunashvili, V.Z.

1989-01-01

In recent years the inverse scattering method has achieved significant successes in the integration of nonlinear models that arise in different branches of physics. However, its region of applicability is still restricted, i.e., not all nonlinear models can be integrated. In view of the great mathematical difficulties that arise in integration, it is clearly worth testing a model for integrability before turning to integration. Such a possibility is provided by the Zakharov-Schulman method. The question of the integrability of a system of two nonlinear Schroedinger equations is resolved. It is shown that the previously known cases exhaust all integrable variants

On integrability conditions of the equations of nonsymmetrical chiral field on SO(4)

International Nuclear Information System (INIS)

Tskhakaya, D.D.

1990-01-01

Possibility of integrating the equations of nonsymmetrical chiral field on SO(4) by means of the inverse scattering method is investigated. Maximal number of the motion integrals is found for the corresponding system of ordinary differential equations

Thermal conduction in classical low-dimensional lattices

International Nuclear Information System (INIS)

Lepri, Stefano; Livi, Roberto; Politi, Antonio

2003-01-01

Deriving macroscopic phenomenological laws of irreversible thermodynamics from simple microscopic models is one of the tasks of non-equilibrium statistical mechanics. We consider stationary energy transport in crystals with reference to simple mathematical models consisting of coupled oscillators on a lattice. The role of lattice dimensionality on the breakdown of the Fourier's law is discussed and some universal quantitative aspects are emphasized: the divergence of the finite-size thermal conductivity is characterized by universal laws in one and two dimensions. Equilibrium and non-equilibrium molecular dynamics methods are presented along with a critical survey of previous numerical results. Analytical results for the non-equilibrium dynamics can be obtained in the harmonic chain where the role of disorder and localization can be also understood. The traditional kinetic approach, based on the Boltzmann-Peierls equation is also briefly sketched with reference to one-dimensional chains. Simple toy models can be defined in which the conductivity is finite. Anomalous transport in integrable non-linear systems is briefly discussed. Finally, possible future research themes are outlined

Lattice QCD for cosmology

International Nuclear Information System (INIS)

Borsanyi, Sz.; Kampert, K.H.; Fodor, Z.; Forschungszentrum Juelich; Eoetvoes Univ., Budapest

2016-06-01

We present a full result for the equation of state (EoS) in 2+1+1 (up/down, strange and charm quarks are present) flavour lattice QCD. We extend this analysis and give the equation of state in 2+1+1+1 flavour QCD. In order to describe the evolution of the universe from temperatures several hundreds of GeV to the MeV scale we also include the known effects of the electroweak theory and give the effective degree of freedoms. As another application of lattice QCD we calculate the topological susceptibility (χ) up to the few GeV temperature region. These two results, EoS and χ, can be used to predict the dark matter axion's mass in the post-inflation scenario and/or give the relationship between the axion's mass and the universal axionic angle, which acts as a initial condition of our universe.

Lattices of dielectric resonators

CERN Document Server

Trubin, Alexander

2016-01-01

This book provides the analytical theory of complex systems composed of a large number of high-Q dielectric resonators. Spherical and cylindrical dielectric resonators with inferior and also whispering gallery oscillations allocated in various lattices are considered. A new approach to S-matrix parameter calculations based on perturbation theory of Maxwell equations, developed for a number of high-Q dielectric bodies, is introduced. All physical relationships are obtained in analytical form and are suitable for further computations. Essential attention is given to a new unified formalism of the description of scattering processes. The general scattering task for coupled eigen oscillations of the whole system of dielectric resonators is described. The equations for the  expansion coefficients are explained in an applicable way. The temporal Green functions for the dielectric resonator are presented. The scattering process of short pulses in dielectric filter structures, dielectric antennas  and lattices of d...

Lattice fluid dynamics from perfect discretizations of continuum flows

International Nuclear Information System (INIS)

Katz, E.; Wiese, U.

1998-01-01

We use renormalization group methods to derive equations of motion for large scale variables in fluid dynamics. The large scale variables are averages of the underlying continuum variables over cubic volumes and naturally exist on a lattice. The resulting lattice dynamics represents a perfect discretization of continuum physics, i.e., grid artifacts are completely eliminated. Perfect equations of motion are derived for static, slow flows of incompressible, viscous fluids. For Hagen-Poiseuille flow in a channel with a square cross section the equations reduce to a perfect discretization of the Poisson equation for the velocity field with Dirichlet boundary conditions. The perfect large scale Poisson equation is used in a numerical simulation and is shown to represent the continuum flow exactly. For nonsquare cross sections one can use a numerical iterative procedure to derive flow equations that are approximately perfect. copyright 1998 The American Physical Society

Solving differential equations for Feynman integrals by expansions near singular points

Science.gov (United States)

Lee, Roman N.; Smirnov, Alexander V.; Smirnov, Vladimir A.

2018-03-01

We describe a strategy to solve differential equations for Feynman integrals by powers series expansions near singular points and to obtain high precision results for the corresponding master integrals. We consider Feynman integrals with two scales, i.e. non-trivially depending on one variable. The corresponding algorithm is oriented at situations where canonical form of the differential equations is impossible. We provide a computer code constructed with the help of our algorithm for a simple example of four-loop generalized sunset integrals with three equal non-zero masses and two zero masses. Our code gives values of the master integrals at any given point on the real axis with a required accuracy and a given order of expansion in the regularization parameter ɛ.

CALL FOR PAPERS: Special issue on Symmetries and Integrability of Difference Equations

Science.gov (United States)

Doliwa, Adam; Korhonen, Risto; Lafortune, Stephane

2006-10-01

This is a call for contributions to a special issue of Journal of Physics A: Mathematical and General entitled `Special issue on Symmetries and Integrability of Difference Equations' as featured at the SIDE VII meeting held during July 2006 in Melbourne (http://web.maths.unsw.edu.au/%7Eschief/side/side.html). Participants at that meeting, as well as other researchers working in the field of difference equations and discrete systems, are invited to submit a research paper to this issue. This meeting was the seventh of a series of biennial meetings devoted to the study of integrable difference equations and related topics. The notion of integrability was first introduced in the 19th century in the context of classical mechanics with the definition of Liouville integrability for Hamiltonian flows. Since then, several notions of integrability have been introduced for partial and ordinary differential equations. Closely related to integrability theory is the symmetry analysis of nonlinear evolution equations. Symmetry analysis takes advantage of the Lie group structure of a given equation to study its properties. Together, integrability theory and symmetry analysis provide the main method by which nonlinear evolution equations can be solved explicitly. Difference equations, just as differential equations, are important in numerous fields of science and have a wide variety of applications in such areas as: mathematical physics, computer visualization, numerical analysis, mathematical biology, economics, combinatorics, quantum field theory, etc. It is thus crucial to develop tools to study and solve difference equations. While the theory of symmetry and integrability for differential equations is now well-established, this is not yet the case for discrete equations. The situation has undergone impressive development in recent years and has affected a broad range of fields, including the theory of special functions, quantum integrable systems, numerical analysis, cellular

Discrete breathers in a two-dimensional Fermi-Pasta-Ulam lattice

International Nuclear Information System (INIS)

Butt, Imran A; Wattis, Jonathan A D

2006-01-01

Using asymptotic methods, we investigate whether discrete breathers are supported by a two-dimensional Fermi-Pasta-Ulam lattice. A scalar (one-component) two-dimensional Fermi-Pasta-Ulam lattice is shown to model the charge stored within an electrical transmission lattice. A third-order multiple-scale analysis in the semi-discrete limit fails, since at this order, the lattice equations reduce to the (2 + 1)-dimensional cubic nonlinear Schroedinger (NLS) equation which does not support stable soliton solutions for the breather envelope. We therefore extend the analysis to higher order and find a generalized (2 + 1)-dimensional NLS equation which incorporates higher order dispersive and nonlinear terms as perturbations. We find an ellipticity criterion for the wave numbers of the carrier wave. Numerical simulations suggest that both stationary and moving breathers are supported by the system. Calculations of the energy show the expected threshold behaviour whereby the energy of breathers does not go to zero with the amplitude; we find that the energy threshold is maximized by stationary breathers, and becomes arbitrarily small as the boundary of the domain of ellipticity is approached

Efficient method for the solution of the energy dependent integral Boltzmann transport equation in the resolved resonance energy region

International Nuclear Information System (INIS)

Schwenk, G.A. Jr.

1980-01-01

The calculation of neutron-nuclei reaction rates in the lower resolved resonance region (167 eV - 1.855 eV) is considered in this dissertation. Particular emphasis is placed on the calculation of these reaction rates for tight lattices where their accuracy is most important. The results of the continuous energy Monte Carlo code, VIM, are chosen as reference values for this study. The primary objective of this work is to develop a method for calculating resonance reaction rates which agree well with the reference solution, yet is efficient enough to be used by nuclear reactor fuel cycle designers on a production basis. A very efficient multigroup solution of the two spatial region energy dependent integral transport equation is developed. This solution, denoted the Broad Group Integral Method (BGIM), uses escape probabilities to obtain the spatial coupling between regions and uses an analytical flux shape within a multigroup to obtain weighted cross sections which account for the rapidly varying resonance cross sections. The multigroup lethargy widths chosen for the numerical integration of the two region energy-dependent neutron continuity equations can be chosen much wider (a factor of 30 larger) than in the direct numerical integration methods since the analytical flux shape is used to account for fine structure effects. The BGIM solution is made highly efficient through the use of these broad groups. It is estimated that for a 10 step unit cell fuel cycle depletion calculation, the computer running time for a production code such as EPRI-LEOPARD would be increased by only 6% through the use of the more accurate and intricate BGIM method in the lower resonance energy region

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)

Lattice models and conformal field theories

International Nuclear Information System (INIS)

Saleur, H.

1988-01-01

Theoretical studies concerning the connection between critical physical systems and the conformal theories are reviewed. The conformal theory associated to a critical (integrable) lattice model is derived. The obtention of the central charge, critical exponents and torus partition function, using renormalization group arguments, is shown. The quantum group structure, in the integrable lattice models, and the theory of Visaro algebra representations are discussed. The relations between off-critical integrable models and conformal theories, in finite geometries, are studied

New block matrix spectral problem and Hamiltonian structure of the discrete integrable coupling system

Energy Technology Data Exchange (ETDEWEB)

Yu Fajun [College of Maths and Systematic Science, Shenyang Normal University, Shenyang 110034 (China)], E-mail: yufajun888@163.com

2008-06-09

In [W.X. Ma, J. Phys. A: Math. Theor. 40 (2007) 15055], Prof. Ma gave a beautiful result (a discrete variational identity). In this Letter, based on a discrete block matrix spectral problem, a new hierarchy of Lax integrable lattice equations with four potentials is derived. By using of the discrete variational identity, we obtain Hamiltonian structure of the discrete soliton equation hierarchy. Finally, an integrable coupling system of the soliton equation hierarchy and its Hamiltonian structure are obtained through the discrete variational identity.

New block matrix spectral problem and Hamiltonian structure of the discrete integrable coupling system

International Nuclear Information System (INIS)

Yu Fajun

2008-01-01

In [W.X. Ma, J. Phys. A: Math. Theor. 40 (2007) 15055], Prof. Ma gave a beautiful result (a discrete variational identity). In this Letter, based on a discrete block matrix spectral problem, a new hierarchy of Lax integrable lattice equations with four potentials is derived. By using of the discrete variational identity, we obtain Hamiltonian structure of the discrete soliton equation hierarchy. Finally, an integrable coupling system of the soliton equation hierarchy and its Hamiltonian structure are obtained through the discrete variational identity

Transverse centroid oscillations in solenoidially focused beam transport lattices

International Nuclear Information System (INIS)

Lund, Steven M.; Wootton, Christopher J.; Lee, Edward P.

2009-01-01

Transverse centroid oscillations are analyzed for a beam in a solenoid transport lattice. Linear equations of motion are derived that describe small-amplitude centroid oscillations induced by displacement and rotational misalignments of the focusing solenoids in the transport lattice, dipole steering elements, and initial centroid offset errors. These equations are analyzed in a local rotating Larmor frame to derive complex-variable 'alignment functions' and 'bending functions' that efficiently describe the characteristics of the centroid oscillations induced by both mechanical misalignments of the solenoids and dipole steering elements. The alignment and bending functions depend only on the properties of the ideal lattice in the absence of errors and steering, and have associated expansion amplitudes set by the misalignments and steering fields, respectively. Applications of this formulation are presented for statistical analysis of centroid oscillations, calculation of actual lattice misalignments from centroid measurements, and optimal beam steering.

Lax Pairs for Discrete Integrable Equations via Darboux Transformations

International Nuclear Information System (INIS)

Cao Ce-Wen; Zhang Guang-Yao

2012-01-01

A method is developed to construct discrete Lax pairs using Darboux transformations. More kinds of Lax pairs are found for some newly appeared discrete integrable equations, including the H1, the special H3 and the Q1 models in the Adler—Bobenko—Suris list and the closely related discrete and semi-discrete pKdV, pMKdV, SG and Liouville equations. (general)

A lattice Boltzmann model for solute transport in open channel flow

Science.gov (United States)

Wang, Hongda; Cater, John; Liu, Haifei; Ding, Xiangyi; Huang, Wei

2018-01-01

A lattice Boltzmann model of advection-dispersion problems in one-dimensional (1D) open channel flows is developed for simulation of solute transport and pollutant concentration. The hydrodynamics are calculated based on a previous lattice Boltzmann approach to solving the 1D Saint-Venant equations (LABSVE). The advection-dispersion model is coupled with the LABSVE using the lattice Boltzmann method. Our research recovers the advection-dispersion equations through the Chapman-Enskog expansion of the lattice Boltzmann equation. The model differs from the existing schemes in two points: (1) the lattice Boltzmann numerical method is adopted to solve the advection-dispersion problem by meso-scopic particle distribution; (2) and the model describes the relation between discharge, cross section area and solute concentration, which increases the applicability of the water quality model in practical engineering. The model is verified using three benchmark tests: (1) instantaneous solute transport within a short distance; (2) 1D point source pollution with constant velocity; (3) 1D point source pollution in a dam break flow. The model is then applied to a 50-year flood point source pollution accident on the Yongding River, which showed good agreement with a MIKE 11 solution and gauging data.

Non-integrability of time-dependent spherically symmetric Yang-Mills equations

Energy Technology Data Exchange (ETDEWEB)

Matinyan, S G; Prokhorenko, E B; Savvidy, G K

1988-03-07

The integrability of time-dependent spherically symmetric Yang-Mills equations is studied using the Fermi-Pasta-Ulam method. It is shown that the motion of this system is ergodic, while the system itself is non-integrable, i.e. manifests dynamical chaos.

Exact finite volume expectation values of \\overline{Ψ}Ψ in the massive Thirring model from light-cone lattice correlators

Science.gov (United States)

Hegedűs, Árpád

2018-03-01

In this paper, using the light-cone lattice regularization, we compute the finite volume expectation values of the composite operator \\overline{Ψ}Ψ between pure fermion states in the Massive Thirring Model. In the light-cone regularized picture, this expectation value is related to 2-point functions of lattice spin operators being located at neighboring sites of the lattice. The operator \\overline{Ψ}Ψ is proportional to the trace of the stress-energy tensor. This is why the continuum finite volume expectation values can be computed also from the set of non-linear integral equations (NLIE) governing the finite volume spectrum of the theory. Our results for the expectation values coming from the computation of lattice correlators agree with those of the NLIE computations. Previous conjectures for the LeClair-Mussardo-type series representation of the expectation values are also checked.

A New time Integration Scheme for Cahn-hilliard Equations

KAUST Repository

Schaefer, R.

2015-06-01

In this paper we present a new integration scheme that can be applied to solving difficult non-stationary non-linear problems. It is obtained by a successive linearization of the Crank- Nicolson scheme, that is unconditionally stable, but requires solving non-linear equation at each time step. We applied our linearized scheme for the time integration of the challenging Cahn-Hilliard equation, modeling the phase separation in fluids. At each time step the resulting variational equation is solved using higher-order isogeometric finite element method, with B- spline basis functions. The method was implemented in the PETIGA framework interfaced via the PETSc toolkit. The GMRES iterative solver was utilized for the solution of a resulting linear system at every time step. We also apply a simple adaptivity rule, which increases the time step size when the number of GMRES iterations is lower than 30. We compared our method with a non-linear, two stage predictor-multicorrector scheme, utilizing a sophisticated step length adaptivity. We controlled the stability of our simulations by monitoring the Ginzburg-Landau free energy functional. The proposed integration scheme outperforms the two-stage competitor in terms of the execution time, at the same time having a similar evolution of the free energy functional.

A New time Integration Scheme for Cahn-hilliard Equations

KAUST Repository

Schaefer, R.; Smol-ka, M.; Dalcin, L; Paszyn'ski, M.

2015-01-01

In this paper we present a new integration scheme that can be applied to solving difficult non-stationary non-linear problems. It is obtained by a successive linearization of the Crank- Nicolson scheme, that is unconditionally stable, but requires solving non-linear equation at each time step. We applied our linearized scheme for the time integration of the challenging Cahn-Hilliard equation, modeling the phase separation in fluids. At each time step the resulting variational equation is solved using higher-order isogeometric finite element method, with B- spline basis functions. The method was implemented in the PETIGA framework interfaced via the PETSc toolkit. The GMRES iterative solver was utilized for the solution of a resulting linear system at every time step. We also apply a simple adaptivity rule, which increases the time step size when the number of GMRES iterations is lower than 30. We compared our method with a non-linear, two stage predictor-multicorrector scheme, utilizing a sophisticated step length adaptivity. We controlled the stability of our simulations by monitoring the Ginzburg-Landau free energy functional. The proposed integration scheme outperforms the two-stage competitor in terms of the execution time, at the same time having a similar evolution of the free energy functional.

Any order approximate analytical solution of the nonlinear Volterra's integral equation for accelerator dynamic systems

International Nuclear Information System (INIS)

Liu Chunliang; Xie Xi; Chen Yinbao

1991-01-01

The universal nonlinear dynamic system equation is equivalent to its nonlinear Volterra's integral equation, and any order approximate analytical solution of the nonlinear Volterra's integral equation is obtained by exact analytical method, thus giving another derivation procedure as well as another computation algorithm for the solution of the universal nonlinear dynamic system equation

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.

Integrability of the one dimensional Schrödinger equation

Science.gov (United States)

Combot, Thierry

2018-02-01

We present a definition of integrability for the one-dimensional Schrödinger equation, which encompasses all known integrable systems, i.e., systems for which the spectrum can be explicitly computed. For this, we introduce the class of rigid functions, built as Liouvillian functions, but containing all solutions of rigid differential operators in the sense of Katz, and a notion of natural of boundary conditions. We then make a complete classification of rational integrable potentials. Many new integrable cases are found, some of them physically interesting.

On the use of the Lie group technique for differential equations with a small parameter: Approximate solutions and integrable equations

International Nuclear Information System (INIS)

Burde, G.I.

2002-01-01

A new approach to the use of the Lie group technique for partial and ordinary differential equations dependent on a small parameter is developed. In addition to determining approximate solutions to the perturbed equation, the approach allows constructing integrable equations that have solutions with (partially) prescribed features. Examples of application of the approach to partial differential equations are given

Group theory and lattice gauge fields

International Nuclear Information System (INIS)

Creutz, M.

1988-09-01

Lattice gauge theory, formulated in terms of invariant integrals over group elements on lattice bonds, benefits from many group theoretical notions. Gauge invariance provides an enormous symmetry and powerful constraints on expectation values. Strong coupling expansions require invariant integrals over polynomials in group elements, all of which can be evaluated by symmetry considerations. Numerical simulations involve random walks over the group. These walks automatically generate the invariant group measure, avoiding explicit parameterization. A recently proposed overrelaxation algorithm is particularly efficient at exploring the group manifold. These and other applications of group theory to lattice gauge fields are reviewed in this talk. 17 refs

Existence Results for Some Nonlinear Functional-Integral Equations in Banach Algebra with Applications

Directory of Open Access Journals (Sweden)

Lakshmi Narayan Mishra

2016-04-01

Full Text Available In the present manuscript, we prove some results concerning the existence of solutions for some nonlinear functional-integral equations which contains various integral and functional equations that considered in nonlinear analysis and its applications. By utilizing the techniques of noncompactness measures, we operate the fixed point theorems such as Darbo's theorem in Banach algebra concerning the estimate on the solutions. The results obtained in this paper extend and improve essentially some known results in the recent literature. We also provide an example of nonlinear functional-integral equation to show the ability of our main result.

On the absence of pentaquark states from dynamics in strongly coupled lattice QCD

Energy Technology Data Exchange (ETDEWEB)

Anjos, Petrus Henrique Ribeiro dos [Universidade Federal de Goias (UFG), Goiania, GO (Brazil); Veiga, Paulo Afonso Faria da; O' Carroll, Michael [Universidade de Sao Paulo (USP), SP (Brazil); Francisco Neto, Antonio [Universidade Federal de Ouro Preto (UFOP), MG (Brazil)

2011-07-01

Full text: We consider an imaginary time functional integral formulation of a two-flavor, 3 + 1 lattice QCD model with Wilson's action and in the strong coupling regime (with a small hopping parameter, {kappa}0, and a much smaller plaquette coupling, {beta} = 1/g{sub 0}{sup 2}, so that the quarks and glueballs are heavy). The model has local SU(3){sub c} gauge and global SU(2){sub f} flavor symmetries, and incorporates the corresponding part of the eightfold way particles: baryons (mesons) of asymptotic mass -3ln{kappa}(-2 ln {kappa}). We search for pentaquark states as meson-baryon bound states in the energy-momentum spectrum of the model, using a lattice Bethe-Salpeter equation. This equation is solved within a ladder approximation, given by the lowest nonvanishing order in {kappa} and {beta} of the Bethe-Salpeter kernel. It includes order 2 contributions with a q-barq exchange potential together with a contribution that is a local-in-space, energy-dependent potential. The attractive or repulsive nature of the exchange interaction depends on the spin of the meson-baryon states. The Bethe-Salpeter equation presents integrable singularities, forcing the couplings to be above a threshold value for the meson and the baryon to bind in a pentaquark. We analyzed all the total isospin sectors, I = 1/2/3/2/ 5/2, for the system. For all I, the net attraction resulting from the two sources of interaction is not strong enough for the meson and the baryon to bind. Thus, within our approximation, these pentaquark states are not present up to near the free meson-baryon energy threshold of - 5 ln{kappa}. This result is to be contrasted with the spinless case for which our method detects meson-baryon bound states, as well as for Yukawa effective baryon and meson field models. A physical interpretation of our results emerges from an approximate correspondence between meson-baryon bound states and negative energy states of a one-particle lattice Schroedinger Hamiltonian

Lattice Boltzmann equation calculation of internal, pressure-driven turbulent flow

International Nuclear Information System (INIS)

Hammond, L A; Halliday, I; Care, C M; Stevens, A

2002-01-01

We describe a mixing-length extension of the lattice Boltzmann approach to the simulation of an incompressible liquid in turbulent flow. The method uses a simple, adaptable, closure algorithm to bound the lattice Boltzmann fluid incorporating a law-of-the-wall. The test application, of an internal, pressure-driven and smooth duct flow, recovers correct velocity profiles for Reynolds number to 1.25 x 10 5 . In addition, the Reynolds number dependence of the friction factor in the smooth-wall branch of the Moody chart is correctly recovered. The method promises a straightforward extension to other curves of the Moody chart and to cylindrical pipe flow

Monte Carlo numerical study of lattice field theories

International Nuclear Information System (INIS)

Gan Cheekwan; Kim Seyong; Ohta, Shigemi

1997-01-01

The authors are interested in the exact first-principle calculations of quantum field theories which are indeed exact ones. For quantum chromodynamics (QCD) at low energy scale, a nonperturbation method is needed, and the only known such method is the lattice method. The path integral can be evaluated by putting a system on a finite 4-dimensional volume and discretizing space time continuum into finite points, lattice. The continuum limit is taken by making the lattice infinitely fine. For evaluating such a finite-dimensional integral, the Monte Carlo numerical estimation of the path integral can be obtained. The calculation of light hadron mass in quenched lattice QCD with staggered quarks, 3-dimensional Thirring model calculation and the development of self-test Monte Carlo method have been carried out by using the RIKEN supercomputer. The motivation of this study, lattice QCD formulation, continuum limit, Monte Carlo update, hadron propagator, light hadron mass, auto-correlation and source size dependence are described on lattice QCD. The phase structure of the 3-dimensional Thirring model for a small 8 3 lattice has been mapped. The discussion on self-test Monte Carlo method is described again. (K.I.)

The multidensity integral equation approach in the theory of complex liquids

International Nuclear Information System (INIS)

Holovko, M.F.

2001-01-01

Recent development of the multi-density integral equation approach and its application to the statistical mechanical modelling of a different type of association and clusterization in liquids and solutions are reviewed. The effects of dimerization, polymerization and network formation are discussed. The numerical and analytical solutions of the integral equations in the multi-density formalism for pair correlation functions are used for the description of structural and thermodynamical properties of ionic solutions, polymers and network forming fluids

The ATOMFT integrator - Using Taylor series to solve ordinary differential equations

Science.gov (United States)

Berryman, Kenneth W.; Stanford, Richard H.; Breckheimer, Peter J.

1988-01-01

This paper discusses the application of ATOMFT, an integration package based on Taylor series solution with a sophisticated user interface. ATOMFT has the capabilities to allow the implementation of user defined functions and the solution of stiff and algebraic equations. Detailed examples, including the solutions to several astrodynamics problems, are presented. Comparisons with its predecessor ATOMCC and other modern integrators indicate that ATOMFT is a fast, accurate, and easy method to use to solve many differential equation problems.

The Neumann Type Systems and Algebro-Geometric Solutions of a System of Coupled Integrable Equations

International Nuclear Information System (INIS)

Chen Jinbing; Qiao Zhijun

2011-01-01

A system of (1+1)-dimensional coupled integrable equations is decomposed into a pair of new Neumann type systems that separate the spatial and temporal variables for this system over a symplectic submanifold. Then, the Neumann type flows associated with the coupled integrable equations are integrated on the complex tour of a Riemann surface. Finally, the algebro-geometric solutions expressed by Riemann theta functions of the system of coupled integrable equations are obtained by means of the Jacobi inversion.

Three-dimensional lattice Boltzmann model for compressible flows.

Science.gov (United States)

Sun, Chenghai; Hsu, Andrew T

2003-07-01

A three-dimensional compressible lattice Boltzmann model is formulated on a cubic lattice. A very large particle-velocity set is incorporated in order to enable a greater variation in the mean velocity. Meanwhile, the support set of the equilibrium distribution has only six directions. Therefore, this model can efficiently handle flows over a wide range of Mach numbers and capture shock waves. Due to the simple form of the equilibrium distribution, the fourth-order velocity tensors are not involved in the formulation. Unlike the standard lattice Boltzmann model, no special treatment is required for the homogeneity of fourth-order velocity tensors on square lattices. The Navier-Stokes equations were recovered, using the Chapman-Enskog method from the Bhatnagar-Gross-Krook (BGK) lattice Boltzmann equation. The second-order discretization error of the fluctuation velocity in the macroscopic conservation equation was eliminated by means of a modified collision invariant. The model is suitable for both viscous and inviscid compressible flows with or without shocks. Since the present scheme deals only with the equilibrium distribution that depends only on fluid density, velocity, and internal energy, boundary conditions on curved wall are easily implemented by an extrapolation of macroscopic variables. To verify the scheme for inviscid flows, we have successfully simulated a three-dimensional shock-wave propagation in a box and a normal shock of Mach number 10 over a wedge. As an application to viscous flows, we have simulated a flat plate boundary layer flow, flow over a cylinder, and a transonic flow over a NACA0012 airfoil cascade.

High-precision numerical integration of equations in dynamics

Science.gov (United States)

Alesova, I. M.; Babadzanjanz, L. K.; Pototskaya, I. Yu.; Pupysheva, Yu. Yu.; Saakyan, A. T.

2018-05-01

An important requirement for the process of solving differential equations in Dynamics, such as the equations of the motion of celestial bodies and, in particular, the motion of cosmic robotic systems is high accuracy at large time intervals. One of effective tools for obtaining such solutions is the Taylor series method. In this connection, we note that it is very advantageous to reduce the given equations of Dynamics to systems with polynomial (in unknowns) right-hand sides. This allows us to obtain effective algorithms for finding the Taylor coefficients, a priori error estimates at each step of integration, and an optimal choice of the order of the approximation used. In the paper, these questions are discussed and appropriate algorithms are considered.

A predictor-corrector scheme for solving the Volterra integral equation

KAUST Repository

Al Jarro, Ahmed

2011-08-01

The occurrence of late time instabilities is a common problem of almost all time marching methods developed for solving time domain integral equations. Implicit marching algorithms are now considered stable with various efforts that have been developed for removing low and high frequency instabilities. On the other hand, literature on stabilizing explicit schemes, which might be considered more efficient since they do not require a matrix inversion at each time step, is practically non-existent. In this work, a stable but still explicit predictor-corrector scheme is proposed for solving the Volterra integral equation and its efficacy is verified numerically. © 2011 IEEE.

Path integral solution of linear second order partial differential equations I: the general construction

International Nuclear Information System (INIS)

LaChapelle, J.

2004-01-01

A path integral is presented that solves a general class of linear second order partial differential equations with Dirichlet/Neumann boundary conditions. Elementary kernels are constructed for both Dirichlet and Neumann boundary conditions. The general solution can be specialized to solve elliptic, parabolic, and hyperbolic partial differential equations with boundary conditions. This extends the well-known path integral solution of the Schroedinger/diffusion equation in unbounded space. The construction is based on a framework for functional integration introduced by Cartier/DeWitt-Morette

Irreversible stochastic processes on lattices

International Nuclear Information System (INIS)

Nord, R.S.

1986-01-01

Models for irreversible random or cooperative filling of lattices are required to describe many processes in chemistry and physics. Since the filling is assumed to be irreversible, even the stationary, saturation state is not in equilibrium. The kinetics and statistics of these processes are described by recasting the master equations in infinite hierarchical form. Solutions can be obtained by implementing various techniques: refinements in these solution techniques are presented. Programs considered include random dimer, trimer, and tetramer filling of 2D lattices, random dimer filling of a cubic lattice, competitive filling of two or more species, and the effect of a random distribution of inactive sites on the filling. Also considered is monomer filling of a linear lattice with nearest neighbor cooperative effects and solve for the exact cluster-size distribution for cluster sizes up to the asymptotic regime. Additionally, a technique is developed to directly determine the asymptotic properties of the cluster size distribution. Finally cluster growth is considered via irreversible aggregation involving random walkers. In particular, explicit results are provided for the large-lattice-size asymptotic behavior of trapping probabilities and average walk lengths for a single walker on a lattice with multiple traps. Procedures for exact calculation of these quantities on finite lattices are also developed

Analytic properties for the honeycomb lattice Green function at the origin

Science.gov (United States)

Joyce, G. S.

2018-05-01

The analytic properties of the honeycomb lattice Green function are investigated, where is a complex variable which lies in a plane. This double integral defines a single-valued analytic function provided that a cut is made along the real axis from w  =  ‑3 to . In order to analyse the behaviour of along the edges of the cut it is convenient to define the limit function where . It is shown that and can be evaluated exactly for all in terms of various hypergeometric functions, where the argument function is always real-valued and rational. The second-order linear Fuchsian differential equation satisfied by is also used to derive series expansions for and which are valid in the neighbourhood of the regular singular points and . Integral representations are established for and , where with . In particular, it is proved that where J 0(z) and Y 0(z) denote Bessel functions of the first and second kind, respectively. The results derived in the paper are utilized to evaluate the associated logarithmic integral where w lies in the cut plane. A new set of orthogonal polynomials which are connected with the honeycomb lattice Green function are also briefly discussed. Finally, a link between and the theory of Pearson random walks in a plane is established.

Numerical Study of Two-Dimensional Volterra Integral Equations by RDTM and Comparison with DTM

Directory of Open Access Journals (Sweden)

Reza Abazari

2013-01-01

Full Text Available The two-dimensional Volterra integral equations are solved using more recent semianalytic method, the reduced differential transform method (the so-called RDTM, and compared with the differential transform method (DTM. The concepts of DTM and RDTM are briefly explained, and their application to the two-dimensional Volterra integral equations is studied. The results obtained by DTM and RDTM together are compared with exact solution. As an important result, it is depicted that the RDTM results are more accurate in comparison with those obtained by DTM applied to the same Volterra integral equations. The numerical results reveal that the RDTM is very effective, convenient, and quite accurate compared to the other kind of nonlinear integral equations. It is predicted that the RDTM can be found widely applicable in engineering sciences.

An algorithm of computing inhomogeneous differential equations for definite integrals

OpenAIRE

Nakayama, Hiromasa; Nishiyama, Kenta

2010-01-01

We give an algorithm to compute inhomogeneous differential equations for definite integrals with parameters. The algorithm is based on the integration algorithm for $D$-modules by Oaku. Main tool in the algorithm is the Gr\\"obner basis method in the ring of differential operators.

On Generating Discrete Integrable Systems via Lie Algebras and Commutator Equations

International Nuclear Information System (INIS)

Zhang Yu-Feng; Tam, Honwah

2016-01-01

In the paper, we introduce the Lie algebras and the commutator equations to rewrite the Tu-d scheme for generating discrete integrable systems regularly. By the approach the various loop algebras of the Lie algebra A_1 are defined so that the well-known Toda hierarchy and a novel discrete integrable system are obtained, respectively. A reduction of the later hierarchy is just right the famous Ablowitz–Ladik hierarchy. Finally, via two different enlarging Lie algebras of the Lie algebra A_1, we derive two resulting differential-difference integrable couplings of the Toda hierarchy, of course, they are all various discrete expanding integrable models of the Toda hierarchy. When the introduced spectral matrices are higher degrees, the way presented in the paper is more convenient to generate discrete integrable equations than the Tu-d scheme by using the software Maple. (paper)

Stochastic quantization of field theories on the lattice and supersymmetrical models

International Nuclear Information System (INIS)

Aldazabal, Gerardo.

1984-01-01

Several aspects of the stochastic quantization method are considered. Specifically, field theories on the lattice and supersymmetrical models are studied. A non-linear sigma model is studied firstly, and it is shown that it is possible to obtain evolution equations written directly for invariant quantities. These ideas are generalized to obtain Langevin equations for the Wilson loops of non-abelian lattice gauge theories U (N) and SU (N). In order to write these equations, some different ways of introducing the constraints which the fields must satisfy are discussed. It is natural to have a strong coupling expansion in these equations. The correspondence with quantum field theory is established, and it is noticed that at all orders in the perturbation theory, Langevin equations reduce to Schwinger-Dyson equations. From another point of view, stochastic quantization is applied to large N matrix models on the lattice. As a result, a simple and systematic way of building reduced models is found. Referring to stochastic quantization in supersymmetric theories, a simple supersymmetric model is studied. It is shown that it is possible to write an evolution equation for the superfield wich leads to quantum field theory results in equilibrium. As the Langevin equation preserves supersymmetry, the property of dimensional reduction known for the quantum model is shown to be valid at all times. (M.E.L.) [es

Higher-Order Integral Equation Methods in Computational Electromagnetics

DEFF Research Database (Denmark)

Jørgensen, Erik; Meincke, Peter

Higher-order integral equation methods have been investigated. The study has focused on improving the accuracy and efficiency of the Method of Moments (MoM) applied to electromagnetic problems. A new set of hierarchical Legendre basis functions of arbitrary order is developed. The new basis...

Approximate solutions for the two-dimensional integral transport equation. The critically mixed methods of resolution

International Nuclear Information System (INIS)

Sanchez, Richard.

1980-11-01

This work is divided into two part the first part (note CEA-N-2165) deals with the solution of complex two-dimensional transport problems, the second one treats the critically mixed methods of resolution. These methods are applied for one-dimensional geometries with highly anisotropic scattering. In order to simplify the set of integral equation provided by the integral transport equation, the integro-differential equation is used to obtain relations that allow to lower the number of integral equation to solve; a general mathematical and numerical study is presented [fr

Stochastic integration and differential equations

CERN Document Server

Protter, Philip E

2003-01-01

It has been 15 years since the first edition of Stochastic Integration and Differential Equations, A New Approach appeared, and in those years many other texts on the same subject have been published, often with connections to applications, especially mathematical finance. Yet in spite of the apparent simplicity of approach, none of these books has used the functional analytic method of presenting semimartingales and stochastic integration. Thus a 2nd edition seems worthwhile and timely, though it is no longer appropriate to call it "a new approach". The new edition has several significant changes, most prominently the addition of exercises for solution. These are intended to supplement the text, but lemmas needed in a proof are never relegated to the exercises. Many of the exercises have been tested by graduate students at Purdue and Cornell Universities. Chapter 3 has been completely redone, with a new, more intuitive and simultaneously elementary proof of the fundamental Doob-Meyer decomposition theorem, t...

Retarded potentials and time domain boundary integral equations a road map

CERN Document Server

Sayas, Francisco-Javier

2016-01-01

This book offers a thorough and self-contained exposition of the mathematics of time-domain boundary integral equations associated to the wave equation, including applications to scattering of acoustic and elastic waves. The book offers two different approaches for the analysis of these integral equations, including a systematic treatment of their numerical discretization using Galerkin (Boundary Element) methods in the space variables and Convolution Quadrature in the time variable. The first approach follows classical work started in the late eighties, based on Laplace transforms estimates. This approach has been refined and made more accessible by tailoring the necessary mathematical tools, avoiding an excess of generality. A second approach contains a novel point of view that the author and some of his collaborators have been developing in recent years, using the semigroup theory of evolution equations to obtain improved results. The extension to electromagnetic waves is explained in one of the appendices...

Lattice Paths and the Constant Term

International Nuclear Information System (INIS)

Brak, R; Essam, J; Osborn, J; Owczarek, A L; Rechnitzer, A

2006-01-01

We firstly review the constant term method (CTM), illustrating its combinatorial connections and show how it can be used to solve a certain class of lattice path problems. We show the connection between the CTM, the transfer matrix method (eigenvectors and eigenvalues), partial difference equations, the Bethe Ansatz and orthogonal polynomials. Secondly, we solve a lattice path problem first posed in 1971. The model stated in 1971 was only solved for a special case - we solve the full model

Lattice Boltzmann simulation of antiplane shear loading of a stationary crack

Science.gov (United States)

Schlüter, Alexander; Kuhn, Charlotte; Müller, Ralf

2018-01-01

In this work, the lattice Boltzmann method is applied to study the dynamic behaviour of linear elastic solids under antiplane shear deformation. In this case, the governing set of partial differential equations reduces to a scalar wave equation for the out of plane displacement in a two dimensional domain. The lattice Boltzmann approach developed by Guangwu (J Comput Phys 161(1):61-69, 2000) in 2006 is used to solve the problem numerically. Some aspects of the scheme are highlighted, including the treatment of the boundary conditions. Subsequently, the performance of the lattice Boltzmann scheme is tested for a stationary crack problem for which an analytic solution exists. The treatment of cracks is new compared to the examples that are discussed in Guangwu's work. Furthermore, the lattice Boltzmann simulations are compared to finite element computations. Finally, the influence of the lattice Boltzmann relaxation parameter on the stability of the scheme is illustrated.

APPLICATION OF BOUNDARY INTEGRAL EQUATION METHOD FOR THERMOELASTICITY PROBLEMS

Directory of Open Access Journals (Sweden)

Vorona Yu.V.

2015-12-01

Full Text Available Boundary Integral Equation Method is used for solving analytically the problems of coupled thermoelastic spherical wave propagation. The resulting mathematical expressions coincide with the solutions obtained in a conventional manner.

Deflation of Eigenvalues for GMRES in Lattice QCD

International Nuclear Information System (INIS)

Morgan, Ronald B.; Wilcox, Walter

2002-01-01

Versions of GMRES with deflation of eigenvalues are applied to lattice QCD problems. Approximate eigenvectors corresponding to the smallest eigenvalues are generated at the same time that linear equations are solved. The eigenvectors improve convergence for the linear equations, and they help solve other right-hand sides

On Fredholm-Stieltjes quadratic integral equation with supremum

International Nuclear Information System (INIS)

Darwish, M.A.

2007-08-01

We prove an existence theorem of monotonic solutions for a quadratic integral equation of Fredholm-Stieltjes type in C[0,1]. The concept of measure of non-compactness and a fixed point theorem due to Darbo are the main tools in carrying out our proof. (author)

Fields on a random lattice

International Nuclear Information System (INIS)

Itzykson, C.

1983-10-01

We review the formulation of field theory and statistical mechanics on a Poissonian random lattice. Topics discussed include random geometry, the construction of field equations for arbitrary spin, the free field spectrum and the question of localization illustrated in the one dimensional case

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.

Canonical algorithms for numerical integration of charged particle motion equations

Science.gov (United States)

Efimov, I. N.; Morozov, E. A.; Morozova, A. R.

2017-02-01

A technique for numerically integrating the equation of charged particle motion in a magnetic field is considered. It is based on the canonical transformations of the phase space in Hamiltonian mechanics. The canonical transformations make the integration process stable against counting error accumulation. The integration algorithms contain a minimum possible amount of arithmetics and can be used to design accelerators and devices of electron and ion optics.

Integral equations with difference kernels on finite intervals

CERN Document Server

Sakhnovich, Lev A

2015-01-01

This book focuses on solving integral equations with difference kernels on finite intervals. The corresponding problem on the semiaxis was previously solved by N. Wiener–E. Hopf and by M.G. Krein. The problem on finite intervals, though significantly more difficult, may be solved using our method of operator identities. This method is also actively employed in inverse spectral problems, operator factorization and nonlinear integral equations. Applications of the obtained results to optimal synthesis, light scattering, diffraction, and hydrodynamics problems are discussed in this book, which also describes how the theory of operators with difference kernels is applied to stable processes and used to solve the famous M. Kac problems on stable processes. In this second edition these results are extensively generalized and include the case of all Levy processes. We present the convolution expression for the well-known Ito formula of the generator operator, a convolution expression that has proven to be fruitful...

Integrable discretizations and self-adaptive moving mesh method for a coupled short pulse equation

International Nuclear Information System (INIS)

Feng, Bao-Feng; Chen, Junchao; Chen, Yong; Maruno, Ken-ichi; Ohta, Yasuhiro

2015-01-01

In the present paper, integrable semi-discrete and fully discrete analogues of a coupled short pulse (CSP) equation are constructed. The key to the construction are the bilinear forms and determinant structure of the solutions of the CSP equation. We also construct N-soliton solutions for the semi-discrete and fully discrete analogues of the CSP equations in the form of Casorati determinants. In the continuous limit, we show that the fully discrete CSP equation converges to the semi-discrete CSP equation, then further to the continuous CSP equation. Moreover, the integrable semi-discretization of the CSP equation is used as a self-adaptive moving mesh method for numerical simulations. The numerical results agree with the analytical results very well. (paper)

Elliptic Euler–Poisson–Darboux equation, critical points and integrable systems

International Nuclear Information System (INIS)

Konopelchenko, B G; Ortenzi, G

2013-01-01

The structure and properties of families of critical points for classes of functions W(z, z-bar ) obeying the elliptic Euler–Poisson–Darboux equation E(1/2, 1/2) are studied. General variational and differential equations governing the dependence of critical points in variational (deformation) parameters are found. Explicit examples of the corresponding integrable quasi-linear differential systems and hierarchies are presented. There are the extended dispersionless Toda/nonlinear Schrödinger hierarchies, the ‘inverse’ hierarchy and equations associated with the real-analytic Eisenstein series E(β, β-bar ;1/2) among them. The specific bi-Hamiltonian structure of these equations is also discussed. (paper)

Solution of four-nucleon integral equations using the effective UPA

International Nuclear Information System (INIS)

Perne, R.; Sandhas, W.

1978-01-01

In the three-body case it is standard to either solve the (two-dimensional) Faddeev equations directly, or to reduce them first to one-dimensional equations by means of separable approximation (expansion) of the underlying two-body interactions. The basic four-body operator identities are reduced by the latter treatment to effective three-body equations only. These may be handled like their genuine three-body analoga, i.e., by directly solving them, or by expanding the effective interactions ocurring into separable terms. Such a procedure provides us in a second step with one-dimensional integral equations for the four-body problem, too. (orig./WL) [de

Finding a nonlinear lattice with improved integrability using Lie transform perturbation theory

International Nuclear Information System (INIS)

Sonnad, Kiran G.; Cary, John R.

2004-01-01

A condition for improved dynamic aperture for nonlinear, alternating gradient transport systems is derived using Lie transform perturbation theory. The Lie transform perturbation method is used here to perform averaging over fast oscillations by canonically transforming to slowly oscillating variables. This is first demonstrated for a linear sinusoidal focusing system. This method is then employed to average the dynamics over a lattice period for a nonlinear focusing system, provided by the use of higher order poles such as sextupoles and octupoles along with alternate gradient quadrupoles. Unlike the traditional approach, the higher order focusing is not treated as a perturbation. The Lie transform method is particularly advantageous for such a system where the form of the Hamiltonian is complex. This is because the method exploits the property of canonical invariance of Poisson brackets so that the change of variables is accomplished by just replacing the old ones with the new. The analysis shows the existence of a condition in which the system is azimuthally symmetric in the transformed, slowly oscillating frame. Such a symmetry in the time averaged frame renders the system nearly integrable in the laboratory frame. This condition leads to reduced chaos and improved confinement when compared to a system that is not close to integrability. Numerical calculations of single-particle trajectories and phase space projections of the dynamic aperture performed for a lattice with quadrupoles and sextupoles confirm that this is indeed the case

Computational Aeroacoustics Using the Generalized Lattice Boltzmann Equation, Phase I

Data.gov (United States)

National Aeronautics and Space Administration — The overall objective of the proposed project is to develop a generalized lattice Boltzmann (GLB) approach as a potential computational aeroacoustics (CAA) tool for...

Integral equations of the first kind, inverse problems and regularization: a crash course

International Nuclear Information System (INIS)

Groetsch, C W

2007-01-01

This paper is an expository survey of the basic theory of regularization for Fredholm integral equations of the first kind and related background material on inverse problems. We begin with an historical introduction to the field of integral equations of the first kind, with special emphasis on model inverse problems that lead to such equations. The basic theory of linear Fredholm equations of the first kind, paying particular attention to E. Schmidt's singular function analysis, Picard's existence criterion, and the Moore-Penrose theory of generalized inverses is outlined. The fundamentals of the theory of Tikhonov regularization are then treated and a collection of exercises and a bibliography are provided

A calderón multiplicative preconditioner for the combined field integral equation

KAUST Repository

Bagci, Hakan; Andriulli, Francesco P.; Cools, Kristof; Olyslager, Femke; Michielssen, Eric

2009-01-01

A Calderón multiplicative preconditioner (CMP) for the combined field integral equation (CFIE) is developed. Just like with previously proposed Caldern-preconditioned CFIEs, a localization procedure is employed to ensure that the equation

An efficient explicit marching on in time solver for magnetic field volume integral equation

KAUST Repository

Sayed, Sadeed Bin; Ulku, H. Arda; Bagci, Hakan

2015-01-01

An efficient explicit marching on in time (MOT) scheme for solving the magnetic field volume integral equation is proposed. The MOT system is cast in the form of an ordinary differential equation and is integrated in time using a PE(CE)m multistep

An integral equation-based numerical solver for Taylor states in toroidal geometries

Science.gov (United States)

O'Neil, Michael; Cerfon, Antoine J.

2018-04-01

We present an algorithm for the numerical calculation of Taylor states in toroidal and toroidal-shell geometries using an analytical framework developed for the solution to the time-harmonic Maxwell equations. Taylor states are a special case of what are known as Beltrami fields, or linear force-free fields. The scheme of this work relies on the generalized Debye source representation of Maxwell fields and an integral representation of Beltrami fields which immediately yields a well-conditioned second-kind integral equation. This integral equation has a unique solution whenever the Beltrami parameter λ is not a member of a discrete, countable set of resonances which physically correspond to spontaneous symmetry breaking. Several numerical examples relevant to magnetohydrodynamic equilibria calculations are provided. Lastly, our approach easily generalizes to arbitrary geometries, both bounded and unbounded, and of varying genus.

Integral equation models for image restoration: high accuracy methods and fast algorithms