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

International Nuclear Information System (INIS)

Grigoriu, Mircea; Samorodnitsky, Gennady

2004-01-01

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

A high order solver for the unbounded Poisson equation

DEFF Research Database (Denmark)

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

2013-01-01

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

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

Science.gov (United States)

Krylov, N. V.; Priola, E.

2017-09-01

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

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

Science.gov (United States)

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

1986-01-01

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

Two dimensional generalizations of the Newcomb equation

International Nuclear Information System (INIS)

Dewar, R.L.; Pletzer, A.

1989-11-01

The Bineau reduction to scalar form of the equation governing ideal, zero frequency linearized displacements from a hydromagnetic equilibrium possessing a continuous symmetry is performed in 'universal coordinates', applicable to both the toroidal and helical cases. The resulting generalized Newcomb equation (GNE) has in general a more complicated form than the corresponding one dimensional equation obtained by Newcomb in the case of circular cylindrical symmetry, but in this cylindrical case , the equation can be transformed to that of Newcomb. In the two dimensional case there is a transformation which leaves the form of the GNE invariant and simplifies the Frobenius expansion about a rational surface, especially in the limit of zero pressure gradient. The Frobenius expansions about a mode rational surface is developed and the connection with Hamiltonian transformation theory is shown. 17 refs

Classical r-matrices and Poisson bracket structures on infinite-dimensional groups

International Nuclear Information System (INIS)

Aratyn, H.; Nissimov, E.; Pacheva, S.

1992-01-01

Starting with a canonical symplectic structure defined on the contangent bundle T * G we derive, via Dirac hamiltonian reduction, Poisson brackets (PBs) on an arbitrary infinite-dimensional group G (admitting central extension). The PB structures are given in terms of an r-operator kernel related to the two-cocycle of the underlying Lie algebra and satisfying a differential classical Yang-Baxter equation. The explicit expressions of the PBs among the group variables for the (N, 0) for N=0, 1, ..., 4 (super-) Virasoro groups and the group of area-preserving diffeomorphisms on the torus are presented. (orig.)

Macroscopic Modeling of a One-Dimensional Electrochemical Cell using the Poisson-Nernst-Planck Equations

Science.gov (United States)

Yan, David

This thesis presents the one-dimensional equations, numerical method and simulations of a model to characterize the dynamical operation of an electrochemical cell. This model extends the current state-of-the art in that it accounts, in a primitive way, for the physics of the electrolyte/electrode interface and incorporates diffuse-charge dynamics, temperature coupling, surface coverage, and polarization phenomena. The one-dimensional equations account for a system with one or two mobile ions of opposite charge, and the electrode reaction we consider (when one is needed) is a one-electron electrodeposition reaction. Though the modeled system is far from representing a realistic electrochemical device, our results show a range of dynamics and behaviors which have not been observed previously, and explore the numerical challenges required when adding more complexity to a model. Furthermore, the basic transport equations (which are developed in three spatial dimensions) can in future accomodate the inclusion of additional physics, and coupling to more complex boundary conditions that incorporate two-dimensional surface phenomena and multi-rate reactions. In the model, the Poisson-Nernst-Planck equations are used to model diffusion and electromigration in an electrolyte, and the generalized Frumkin-Butler-Volmer equation is used to model reaction kinetics at electrodes. An energy balance equation is derived and coupled to the diffusion-migration equation. The model also includes dielectric polarization effects by introducing different values of the dielectric permittivity in different regions of the bulk, as well as accounting for surface coverage effects due to adsorption, and finite size "crowding", or steric effects. Advection effects are not modeled but could in future be incorporated. In order to solve the coupled PDE's, we use a variable step size second order scheme in time and finite differencing in space. Numerical tests are performed on a simplified system and

Two dimensional analytical model for a reconfigurable field effect transistor

Science.gov (United States)

Ranjith, R.; Jayachandran, Remya; Suja, K. J.; Komaragiri, Rama S.

2018-02-01

This paper presents two-dimensional potential and current models for a reconfigurable field effect transistor (RFET). Two potential models which describe subthreshold and above-threshold channel potentials are developed by solving two-dimensional (2D) Poisson's equation. In the first potential model, 2D Poisson's equation is solved by considering constant/zero charge density in the channel region of the device to get the subthreshold potential characteristics. In the second model, accumulation charge density is considered to get above-threshold potential characteristics of the device. The proposed models are applicable for the device having lightly doped or intrinsic channel. While obtaining the mathematical model, whole body area is divided into two regions: gated region and un-gated region. The analytical models are compared with technology computer-aided design (TCAD) simulation results and are in complete agreement for different lengths of the gated regions as well as at various supply voltage levels.

Poisson equation in the Kohn-Sham Coulomb problem

OpenAIRE

Manby, F. R.; Knowles, Peter James

2001-01-01

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

Poisson hierarchy of discrete strings

International Nuclear Information System (INIS)

Ioannidou, Theodora; Niemi, Antti J.

2016-01-01

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

Poisson hierarchy of discrete strings

Energy Technology Data Exchange (ETDEWEB)

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

2016-01-28

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

Control Operator for the Two-Dimensional Energized Wave Equation

Directory of Open Access Journals (Sweden)

Sunday Augustus REJU

2006-07-01

Full Text Available This paper studies the analytical model for the construction of the two-dimensional Energized wave equation. The control operator is given in term of space and time t independent variables. The integral quadratic objective cost functional is subject to the constraint of two-dimensional Energized diffusion, Heat and a source. The operator that shall be obtained extends the Conjugate Gradient method (ECGM as developed by Hestenes et al (1952, [1]. The new operator enables the computation of the penalty cost, optimal controls and state trajectories of the two-dimensional energized wave equation when apply to the Conjugate Gradient methods in (Waziri & Reju, LEJPT & LJS, Issues 9, 2006, [2-4] to appear in this series.

Coefficient Inverse Problem for Poisson's Equation in a Cylinder

NARCIS (Netherlands)

Solov'ev, V. V.

2011-01-01

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

The Poisson equation on Klein surfaces

Directory of Open Access Journals (Sweden)

Monica Rosiu

2016-04-01

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

Poisson structure of the equations of ideal multispecies fluid electrodynamics

International Nuclear Information System (INIS)

Spencer, R.G.

1984-01-01

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

Numerical simulation of transient, adiabatic, two-dimensional two-phase flow using the two-fluid model

International Nuclear Information System (INIS)

Neves Conti, T. das.

1983-01-01

A numerical method is developed to simulate adiabatic, transient, two-dimensional two-phase flow. The two-fluid model is used to obtain the mass and momentum conservation equations. These are solved by an iterative algorithm emphoying a time-marching scheme. Based on the corrective procedure of Hirt and Harlow a poisson equation is derived for the pressure field. This equation is finite-differenced and solved by a suitable matrix inversion technique. In the absence of experiment results several numerical tests were made in order to chec accuracy, convergence and stability of the proposed method. Several tests were also performed to check whether the behavior of void fraction and phasic velocities conforms with previous observations. (Author) [pt

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

Science.gov (United States)

Dick, Christian; Rogowsky, Marcus; Westermann, Rudiger

2016-11-01

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

Geometrical aspects of solvable two dimensional models

International Nuclear Information System (INIS)

Tanaka, K.

1989-01-01

It was noted that there is a connection between the non-linear two-dimensional (2D) models and the scalar curvature r, i.e., when r = -2 the equations of motion of the Liouville and sine-Gordon models were obtained. Further, solutions of various classical nonlinear 2D models can be obtained from the condition that the appropriate curvature two form Ω = 0, which suggests that these models are closely related. This relation is explored further in the classical version by obtaining the equations of motion from the evolution equations, the infinite number of conserved quantities, and the common central charge. The Poisson brackets of the solvable 2D models are specified by the Virasoro algebra. 21 refs

Poisson-Lie T-plurality

International Nuclear Information System (INIS)

Unge, Rikard von

2002-01-01

We extend the path-integral formalism for Poisson-Lie T-duality to include the case of Drinfeld doubles which can be decomposed into bi-algebras in more than one way. We give the correct shift of the dilaton, correcting a mistake in the literature. We then use the fact that the six dimensional Drinfeld doubles have been classified to write down all possible conformal Poisson-Lie T-duals of three dimensional space times and we explicitly work out two duals to the constant dilaton and zero anti-symmetric tensor Bianchi type V space time and show that they satisfy the string equations of motion. This space-time was previously thought to have no duals because of the tracefulness of the structure constants. (author)

A high order solver for the unbounded Poisson equation

DEFF Research Database (Denmark)

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

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

Poisson equation for weak gravitational lensing

International Nuclear Information System (INIS)

Kling, Thomas P.; Campbell, Bryan

2008-01-01

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

High order Poisson Solver for unbounded flows

DEFF Research Database (Denmark)

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

2015-01-01

This paper presents a high order method for solving the unbounded Poisson equation on a regular mesh using a Green’s function solution. The high order convergence was achieved by formulating mollified integration kernels, that were derived from a filter regularisation of the solution field....... The method was implemented on a rectangular domain using fast Fourier transforms (FFT) to increase computational efficiency. The Poisson solver was extended to directly solve the derivatives of the solution. This is achieved either by including the differential operator in the integration kernel...... the equations of fluid mechanics as an example, but can be used in many physical problems to solve the Poisson equation on a rectangular unbounded domain. For the two-dimensional case we propose an infinitely smooth test function which allows for arbitrary high order convergence. Using Gaussian smoothing...

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

International Nuclear Information System (INIS)

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

2017-01-01

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

A node-centered local refinement algorithm for poisson's equation in complex geometries

International Nuclear Information System (INIS)

McCorquodale, Peter; Colella, Phillip; Grote, David P.; Vay, Jean-Luc

2004-01-01

This paper presents a method for solving Poisson's equation with Dirichlet boundary conditions on an irregular bounded three-dimensional region. The method uses a nodal-point discretization and adaptive mesh refinement (AMR) on Cartesian grids, and the AMR multigrid solver of Almgren. The discrete Laplacian operator at internal boundaries comes from either linear or quadratic (Shortley-Weller) extrapolation, and the two methods are compared. It is shown that either way, solution error is second order in the mesh spacing. Error in the gradient of the solution is first order with linear extrapolation, but second order with Shortley-Weller. Examples are given with comparison with the exact solution. The method is also applied to a heavy-ion fusion accelerator problem, showing the advantage of adaptivity

Hamiltonian field description of two-dimensional vortex fluids and guiding center plasmas

International Nuclear Information System (INIS)

Morrison, P.J.

1981-03-01

The equations that describe the motion of two-dimensional vortex fluids and guiding center plasmas are shown to possess underlying field Hamiltonian structure. A Poisson bracket which is given in terms of the vorticity, the physical although noncanonical dynamical variable, casts these equations into Heisenberg form. The Hamiltonian density is the kinetic energy density of the fluid. The well-known conserved quantities are seen to be in involution with respect to this Poisson bracket. Expanding the vorticity in terms of a Fourier-Dirac series transforms the field description given here into the usual canonical equations for discrete vortex motion. A Clebsch potential representation of the vorticity transforms the noncanonical field description into a canonical description

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

Energy Technology Data Exchange (ETDEWEB)

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

2016-06-15

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

Recognition of Equations Using a Two-Dimensional Stochastic Context-Free Grammar

Science.gov (United States)

Chou, Philip A.

1989-11-01

We propose using two-dimensional stochastic context-free grammars for image recognition, in a manner analogous to using hidden Markov models for speech recognition. The value of the approach is demonstrated in a system that recognizes printed, noisy equations. The system uses a two-dimensional probabilistic version of the Cocke-Younger-Kasami parsing algorithm to find the most likely parse of the observed image, and then traverses the corresponding parse tree in accordance with translation formats associated with each production rule, to produce eqn I troff commands for the imaged equation. In addition, it uses two-dimensional versions of the Inside/Outside and Baum re-estimation algorithms for learning the parameters of the grammar from a training set of examples. Parsing the image of a simple noisy equation currently takes about one second of cpu time on an Alliant FX/80.

Cellular solutions for the Poisson equation in extended systems

International Nuclear Information System (INIS)

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

1994-01-01

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

Exact solution for the Poisson field in a semi-infinite strip.

Science.gov (United States)

Cohen, Yossi; Rothman, Daniel H

2017-04-01

The Poisson equation is associated with many physical processes. Yet exact analytic solutions for the two-dimensional Poisson field are scarce. Here we derive an analytic solution for the Poisson equation with constant forcing in a semi-infinite strip. We provide a method that can be used to solve the field in other intricate geometries. We show that the Poisson flux reveals an inverse square-root singularity at a tip of a slit, and identify a characteristic length scale in which a small perturbation, in a form of a new slit, is screened by the field. We suggest that this length scale expresses itself as a characteristic spacing between tips in real Poisson networks that grow in response to fluxes at tips.

Stabilizing local boundary conditions for two-dimensional shallow water equations

KAUST Repository

Dia, Ben Mansour; Oppelstrup, Jesper

2018-01-01

In this article, we present a sub-critical two-dimensional shallow water flow regulation. From the energy estimate of a set of one-dimensional boundary stabilization problems, we obtain a set of polynomial equations with respect to the boundary

Global Existence and Large Time Behavior of Solutions to the Bipolar Nonisentropic Euler-Poisson Equations

Directory of Open Access Journals (Sweden)

Min Chen

2014-01-01

Full Text Available We study the one-dimensional bipolar nonisentropic Euler-Poisson equations which can model various physical phenomena, such as the propagation of electron and hole in submicron semiconductor devices, the propagation of positive ion and negative ion in plasmas, and the biological transport of ions for channel proteins. We show the existence and large time behavior of global smooth solutions for the initial value problem, when the difference of two particles’ initial mass is nonzero, and the far field of two particles’ initial temperatures is not the ambient device temperature. This result improves that of Y.-P. Li, for the case that the difference of two particles’ initial mass is zero, and the far field of the initial temperature is the ambient device temperature.

Nambu-Poisson reformulation of the finite dimensional dynamical systems

International Nuclear Information System (INIS)

Baleanu, D.; Makhaldiani, N.

1998-01-01

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

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.

Multivariate fractional Poisson processes and compound sums

OpenAIRE

Beghin, Luisa; Macci, Claudio

2015-01-01

In this paper we present multivariate space-time fractional Poisson processes by considering common random time-changes of a (finite-dimensional) vector of independent classical (non-fractional) Poisson processes. In some cases we also consider compound processes. We obtain some equations in terms of some suitable fractional derivatives and fractional difference operators, which provides the extension of known equations for the univariate processes.

Zero-dimensional limit of the two-dimensional Lugiato-Lefever equation

Science.gov (United States)

Cardoso, Wesley B.; Salasnich, Luca; Malomed, Boris A.

2017-05-01

We study effects of tight harmonic-oscillator confinement on the electromagnetic field in a laser cavity by solving the two-dimensional Lugiato-Lefever (2D LL) equation, taking into account self-focusing or defocusing nonlinearity, losses, pump, and the trapping potential. Tightly confined (quasi-zero-dimensional) optical modes (pixels), produced by this model, are analyzed by means of the variational approximation, which provides a qualitative picture of the ensuing phenomena. This is followed by systematic simulations of the time-dependent 2D LL equation, which reveal the shape, stability, and dynamical behavior of the resulting localized patterns. In this way, we produce stability diagrams for the expected pixels. Then, we consider the LL model with the vortical pump, showing that it can produce stable pixels with embedded vorticity (vortex solitons) in remarkably broad stability areas. Alongside confined vortices with the simple single-ring structure, in the latter case the LL model gives rise to stable multi-ring states, with a spiral phase field. In addition to the numerical results, a qualitatively correct description of the vortex solitons is provided by the Thomas-Fermi approximation. Contribution to the Topical Issue: "Theory and Applications of the Lugiato-Lefever Equation", edited by Yanne K. Chembo, Damia Gomila, Mustapha Tlidi, Curtis R. Menyuk.

A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains

Science.gov (United States)

Johansen, Hans; Colella, Phillip

1998-11-01

We present a numerical method for solving Poisson's equation, with variable coefficients and Dirichlet boundary conditions, on two-dimensional regions. The approach uses a finite-volume discretization, which embeds the domain in a regular Cartesian grid. We treat the solution as a cell-centered quantity, even when those centers are outside the domain. Cells that contain a portion of the domain boundary use conservative differencing of second-order accurate fluxes on each cell volume. The calculation of the boundary flux ensures that the conditioning of the matrix is relatively unaffected by small cell volumes. This allows us to use multigrid iterations with a simple point relaxation strategy. We have combined this with an adaptive mesh refinement (AMR) procedure. We provide evidence that the algorithm is second-order accurate on various exact solutions and compare the adaptive and nonadaptive calculations.

Two-dimensional computer simulation of high intensity proton beams

CERN Document Server

Lapostolle, Pierre M

1972-01-01

A computer program has been developed which simulates the two- dimensional transverse behaviour of a proton beam in a focusing channel. The model is represented by an assembly of a few thousand 'superparticles' acted upon by their own self-consistent electric field and an external focusing force. The evolution of the system is computed stepwise in time by successively solving Poisson's equation and Newton's law of motion. Fast Fourier transform techniques are used for speed in the solution of Poisson's equation, while extensive area weighting is utilized for the accurate evaluation of electric field components. A computer experiment has been performed on the CERN CDC 6600 computer to study the nonlinear behaviour of an intense beam in phase space, showing under certain circumstances a filamentation due to space charge and an apparent emittance growth. (14 refs).

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

International Nuclear Information System (INIS)

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

2006-01-01

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

Robust iterative observer for source localization for Poisson equation

KAUST Repository

Majeed, Muhammad Usman

2017-01-05

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

Robust iterative observer for source localization for Poisson equation

KAUST Repository

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

2017-01-01

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

An energy principle for two-dimensional collisionless relativistic plasmas

International Nuclear Information System (INIS)

Otto, A.; Schindler, K.

1984-01-01

Using relativistic Vlasov theory an energy principle for two-dimensional plasmas is derived, which provides a sufficient and necessary criterion for the stability of relativistic plasma equilibria. This energy principle includes charge separating effects since the exact Poisson equation was taken into consideration. Applying the variational principle to the case of the relativistic plane plasma sheet, the same marginal wave length is found as in the non-relativistic case. (author)

Solution and Study of the Two-Dimensional Nodal Neutron Transport Equation

International Nuclear Information System (INIS)

Panta Pazos, Ruben; Biasotto Hauser, Eliete; Tullio de Vilhena, Marco

2002-01-01

In the last decade Vilhena and coworkers reported an analytical solution to the two-dimensional nodal discrete-ordinates approximations of the neutron transport equation in a convex domain. The key feature of these works was the application of the combined collocation method of the angular variable and nodal approach in the spatial variables. By nodal approach we mean the transverse integration of the SN equations. This procedure leads to a set of one-dimensional S N equations for the average angular fluxes in the variables x and y. These equations were solved by the old version of the LTS N method, which consists in the application of the Laplace transform to the set of nodal S N equations and solution of the resulting linear system by symbolic computation. It is important to recall that this procedure allow us to increase N the order of S N up to 16. To overcome this drawback we step forward performing a spectral painstaking analysis of the nodal S N equations for N up to 16 and we begin the convergence of the S N nodal equations defining an error for the angular flux and estimating the error in terms of the truncation error of the quadrature approximations of the integral term. Furthermore, we compare numerical results of this approach with those of other techniques used to solve the two-dimensional discrete approximations of the neutron transport equation. (authors)

Reference manual for the POISSON/SUPERFISH Group of Codes

Energy Technology Data Exchange (ETDEWEB)

1987-01-01

The POISSON/SUPERFISH Group codes were set up to solve two separate problems: the design of magnets and the design of rf cavities in a two-dimensional geometry. The first stage of either problem is to describe the layout of the magnet or cavity in a way that can be used as input to solve the generalized Poisson equation for magnets or the Helmholtz equations for cavities. The computer codes require that the problems be discretized by replacing the differentials (dx,dy) by finite differences ({delta}X,{delta}Y). Instead of defining the function everywhere in a plane, the function is defined only at a finite number of points on a mesh in the plane.

Exact solutions and conservation laws of the system of two-dimensional viscous Burgers equations

Science.gov (United States)

Abdulwahhab, Muhammad Alim

2016-10-01

Fluid turbulence is one of the phenomena that has been studied extensively for many decades. Due to its huge practical importance in fluid dynamics, various models have been developed to capture both the indispensable physical quality and the mathematical structure of turbulent fluid flow. Among the prominent equations used for gaining in-depth insight of fluid turbulence is the two-dimensional Burgers equations. Its solutions have been studied by researchers through various methods, most of which are numerical. Being a simplified form of the two-dimensional Navier-Stokes equations and its wide range of applicability in various fields of science and engineering, development of computationally efficient methods for the solution of the two-dimensional Burgers equations is still an active field of research. In this study, Lie symmetry method is used to perform detailed analysis on the system of two-dimensional Burgers equations. Optimal system of one-dimensional subalgebras up to conjugacy is derived and used to obtain distinct exact solutions. These solutions not only help in understanding the physical effects of the model problem but also, can serve as benchmarks for constructing algorithms and validation of numerical solutions of the system of Burgers equations under consideration at finite Reynolds numbers. Independent and nontrivial conserved vectors are also constructed.

Self-energy-modified Poisson-Nernst-Planck equations: WKB approximation and finite-difference approaches.

Science.gov (United States)

Xu, Zhenli; Ma, Manman; Liu, Pei

2014-07-01

We propose a modified Poisson-Nernst-Planck (PNP) model to investigate charge transport in electrolytes of inhomogeneous dielectric environment. The model includes the ionic polarization due to the dielectric inhomogeneity and the ion-ion correlation. This is achieved by the self energy of test ions through solving a generalized Debye-Hückel (DH) equation. We develop numerical methods for the system composed of the PNP and DH equations. Particularly, toward the numerical challenge of solving the high-dimensional DH equation, we developed an analytical WKB approximation and a numerical approach based on the selective inversion of sparse matrices. The model and numerical methods are validated by simulating the charge diffusion in electrolytes between two electrodes, for which effects of dielectrics and correlation are investigated by comparing the results with the prediction by the classical PNP theory. We find that, at the length scale of the interface separation comparable to the Bjerrum length, the results of the modified equations are significantly different from the classical PNP predictions mostly due to the dielectric effect. It is also shown that when the ion self energy is in weak or mediate strength, the WKB approximation presents a high accuracy, compared to precise finite-difference results.

Whitham modulation theory for the two-dimensional Benjamin-Ono equation.

Science.gov (United States)

Ablowitz, Mark; Biondini, Gino; Wang, Qiao

2017-09-01

Whitham modulation theory for the two-dimensional Benjamin-Ono (2DBO) equation is presented. A system of five quasilinear first-order partial differential equations is derived. The system describes modulations of the traveling wave solutions of the 2DBO equation. These equations are transformed to a singularity-free hydrodynamic-like system referred to here as the 2DBO-Whitham system. Exact reductions of this system are discussed, the formulation of initial value problems is considered, and the system is used to study the transverse stability of traveling wave solutions of the 2DBO equation.

Electromagnetic wave propagation over an inhomogeneous flat earth (two-dimensional integral equation formulation)

International Nuclear Information System (INIS)

de Jong, G.

1975-01-01

With the aid of a two-dimensional integral equation formulation, the ground wave propagation of electromagnetic waves transmitted by a vertical electric dipole over an inhomogeneous flat earth is investigated. For the configuration in which a ground wave is propagating across an ''island'' on a flat earth, the modulus and argument of the attenuation function have been computed. The results for the two-dimensional treatment are significantly more accurate in detail than the calculations using a one-dimensional integral equation

The Poisson equation in axisymmetric domains with conical points

International Nuclear Information System (INIS)

Nkemzi, B.

2003-01-01

This paper analyzes the application of the Fourier-finite-element method (FFEM) for the resolution of the Derichlet problem for the Poisson equation -Δu-circumflex = f-circumflex in axisymmetric domains Ω-circumflex subset of R 3 with conical points on the rotation axis. The FFEM combines the approximate Fourier method with respect to one space direction with the finite element method for the approximate calculation of the Fourier coefficients of the solution. Here, the influence of the conical points on the regularity of the Fourier coefficients of the solution is analyzed and the asymptotic behaviour of the coefficients near the conical points is described by some singularity functions and treated numerically by mesh grading in the two-dimensional meridian of Ω-circumflex. It is proved that for f-circumflex in L 2 (Ω-circumflex), the rate of convergence of the combined approximations in the Sobolev space W 2 1 (Ω-circumflex) is of the order O(h + N -1 ), where h and N represent, respectively, the parameters of the finite-element- and the Fourier-approximation, with h → 0 and n → ∞. (author)

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

International Nuclear Information System (INIS)

Yang Xiao; Du Dianlou

2010-01-01

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

Poisson's theorem and integrals of KdV equation

International Nuclear Information System (INIS)

Tasso, H.

1978-01-01

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

Poisson's equation in de Sitter space-time

Energy Technology Data Exchange (ETDEWEB)

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

1980-11-01

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

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

International Nuclear Information System (INIS)

Sharifi, M. J.; Adibi, A.

2000-01-01

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

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

Science.gov (United States)

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

2017-06-01

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

A fast Poisson solver for unsteady incompressible Navier-Stokes equations on the half-staggered grid

Science.gov (United States)

Golub, G. H.; Huang, L. C.; Simon, H.; Tang, W. -P.

1995-01-01

In this paper, a fast Poisson solver for unsteady, incompressible Navier-Stokes equations with finite difference methods on the non-uniform, half-staggered grid is presented. To achieve this, new algorithms for diagonalizing a semi-definite pair are developed. Our fast solver can also be extended to the three dimensional case. The motivation and related issues in using this second kind of staggered grid are also discussed. Numerical testing has indicated the effectiveness of this algorithm.

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

Generalized master equations for non-Poisson dynamics on networks.

Science.gov (United States)

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

2012-10-01

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

Potential formation in a one-dimensional bounded plasma system containing a two-electron temperature plasma: Kinetic model and PIC simulation

International Nuclear Information System (INIS)

Gyergyek, T.; Jurcic-Zlobec, B.; Cercek, M.

2008-01-01

Potential formation in a bounded plasma system that contains electrons with a two-temperature velocity distribution and is terminated by a floating, electron emitting electrode (collector) is studied by a one-dimensional kinetic model. A method on how to determine the boundary conditions at the collector for the numerical solution of the Poisson equation is presented. The difference between the regular and the irregular numerical solutions of the Poisson equation is explained. The regular numerical solution of the Poisson equation fulfills the boundary conditions at the source and can be computed for any distance from the collector. The irregular solution does not fulfill the source boundary conditions and the computation breaks down at some distance from the collector. An excellent agreement of the values of the potential at the inflection point found from the numerical solution of the Poisson equation with the values predicted by the analytical model is obtained. Potential, electric field, and particle density profiles found by the numerical solution of the Poisson equation are compared to the profiles obtained with the particle in cell computer simulation. A very good quantitative agreement of the potential and electric field profiles is obtained. For certain values of the parameters the analytical model predicts three possible values of the potential at the inflection point. In such cases always only one of the corresponding numerical solutions of the Poisson equation is regular, while the other two are irregular. The regular numerical solution of the Poisson equation always corresponds to the solution of the model that predicts the largest ion flux to the collector

A mixed method Poisson solver for three-dimensional self-gravitating astrophysical fluid dynamical systems

Science.gov (United States)

Duncan, Comer; Jones, Jim

1993-01-01

A key ingredient in the simulation of self-gravitating astrophysical fluid dynamical systems is the gravitational potential and its gradient. This paper focuses on the development of a mixed method multigrid solver of the Poisson equation formulated so that both the potential and the Cartesian components of its gradient are self-consistently and accurately generated. The method achieves this goal by formulating the problem as a system of four equations for the gravitational potential and the three Cartesian components of the gradient and solves them using a distributed relaxation technique combined with conventional full multigrid V-cycles. The method is described, some tests are presented, and the accuracy of the method is assessed. We also describe how the method has been incorporated into our three-dimensional hydrodynamics code and give an example of an application to the collision of two stars. We end with some remarks about the future developments of the method and some of the applications in which it will be used in astrophysics.

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

Science.gov (United States)

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

2011-11-01

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

Methods for the solution of the two-dimensional radiation-transfer equation

International Nuclear Information System (INIS)

Weaver, R.; Mihalas, D.; Olson, G.

1982-01-01

We use the variable Eddington factor (VEF) approximation to solve the time-dependent two-dimensional radiation transfer equation. The transfer equation and its moments are derived for an inertial frame of reference in cylindrical geometry. Using the VEF tensor to close the moment equations, we manipulate them into a combined moment equation that results in an energy equation, which is automatically flux limited. There are two separable facets in this method of solution. First, given the variable Eddington tensor, we discuss the efficient solution of the combined moment matrix equation. The second facet of the problem is the calculation of the variable Eddington tensor. Several options for this calculation, as well as physical limitations on the use of locally-calculated Eddington factors, are discussed

On integrability of a noncommutative q-difference two-dimensional Toda lattice equation

Energy Technology Data Exchange (ETDEWEB)

Li, C.X., E-mail: trisha_li2001@163.com [School of Mathematical Sciences, Capital Normal University, Beijing 100048 (China); Department of Mathematics, College of Charleston, Charleston, SC 29401 (United States); Nimmo, J.J.C., E-mail: jonathan.nimmo@glasgow.ac.uk [School of Mathematics and Statistics, University of Glasgow, Glasgow G12 8QW (United Kingdom); Shen, Shoufeng, E-mail: mathssf@zjut.edu.cn [Department of Applied Mathematics, Zhejiang University of Technology, Hangzhou 310023 (China)

2015-12-18

In our previous work (C.X. Li and J.J.C. Nimmo, 2009 [18]), we presented a generalized type of Darboux transformations in terms of a twisted derivation in a unified form. The twisted derivation includes ordinary derivatives, forward difference operators, super derivatives and q-difference operators as its special cases. This result not only enables one to recover the known Darboux transformations and quasideterminant solutions to the noncommutative KP equation, the non-Abelian two-dimensional Toda lattice equation, the non-Abelian Hirota–Miwa equation and the super KdV equation, but also inspires us to investigate quasideterminant solutions to q-difference soliton equations. In this paper, we first construct the bilinear Bäcklund transformations for the known bilinear q-difference two-dimensional Toda lattice equation (q-2DTL) and then derive a Lax pair whose compatibility gives a formally different nonlinear q-2DTL equation and finally obtain its quasideterminant solutions by iterating its Darboux transformations. - Highlights: • Examples are given to illustrate the extensive applications of twisted derivations. • Bilinear Bäcklund transformation is constructed for the known q-2DTL equation. • Lax pair is obtained for an equivalent q-2DTL equation. • Quasideterminant solutions are found for the nc q-2DTL equation.

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

International Nuclear Information System (INIS)

Cavazos-Cadena, Rolando; Hernandez-Hernandez, Daniel

2006-01-01

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

Two-dimensional nonlinear string-type equations and their exact integration

International Nuclear Information System (INIS)

Leznov, A.N.; Saveliev, M.V.

1982-01-01

On the base of group-theoretical formulation for exactly integrable two-dimensional non-linear dynamical systems associated with a local part of an arbitrary graded Lie algebra we study a string-type subclass of the equations. Explicit expressions have been obtained for their general solutions

Discrete formulation for two-dimensional multigroup neutron diffusion equations

Energy Technology Data Exchange (ETDEWEB)

Vosoughi, Naser E-mail: vosoughi@mehr.sharif.edu; Salehi, Ali A.; Shahriari, Majid

2003-02-01

The objective of this paper is to introduce a new numerical method for neutronic calculation in a reactor core. This method can produce the final finite form of the neutron diffusion equation by classifying the neutronic variables and using two kinds of cell complexes without starting from the conventional differential form of the neutron diffusion equation. The method with linear interpolation produces the same convergence as the linear continuous finite element method. The quadratic interpolation is proven; the convergence order depends on the shape of the dual cell. The maximum convergence order is achieved by choosing the dual cell based on two Gauss' points. The accuracy of the method was examined with a well-known IAEA two-dimensional benchmark problem. The numerical results demonstrate the effectiveness of the new method.

Discrete formulation for two-dimensional multigroup neutron diffusion equations

International Nuclear Information System (INIS)

Vosoughi, Naser; Salehi, Ali A.; Shahriari, Majid

2003-01-01

The objective of this paper is to introduce a new numerical method for neutronic calculation in a reactor core. This method can produce the final finite form of the neutron diffusion equation by classifying the neutronic variables and using two kinds of cell complexes without starting from the conventional differential form of the neutron diffusion equation. The method with linear interpolation produces the same convergence as the linear continuous finite element method. The quadratic interpolation is proven; the convergence order depends on the shape of the dual cell. The maximum convergence order is achieved by choosing the dual cell based on two Gauss' points. The accuracy of the method was examined with a well-known IAEA two-dimensional benchmark problem. The numerical results demonstrate the effectiveness of the new method

Solution of the two-dimensional space-time reactor kinetics equation by a locally one-dimensional method

International Nuclear Information System (INIS)

Chen, G.S.; Christenson, J.M.

1985-01-01

In this paper, the authors present some initial results from an investigation of the application of a locally one-dimensional (LOD) finite difference method to the solution of the two-dimensional, two-group reactor kinetics equations. Although the LOD method is relatively well known, it apparently has not been previously applied to the space-time kinetics equations. In this investigation, the LOD results were benchmarked against similar computational results (using the same computing environment, the same programming structure, and the same sample problems) obtained by the TWIGL program. For all of the problems considered, the LOD method provided accurate results in one-half to one-eight of the time required by the TWIGL program

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

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

Science.gov (United States)

Yang, Xiao; Du, Dianlou

2010-08-01

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

A high order solver for the unbounded Poisson equation

DEFF Research Database (Denmark)

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

2012-01-01

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

Solution of the two- dimensional heat equation for a rectangular plate

Directory of Open Access Journals (Sweden)

Nurcan BAYKUŞ SAVAŞANERİL

2015-11-01

Full Text Available Laplace equation is a fundamental equation of applied mathematics. Important phenomena in engineering and physics, such as steady-state temperature distribution, electrostatic potential and fluid flow, are modeled by means of this equation. The Laplace equation which satisfies boundary values is known as the Dirichlet problem. The solutions to the Dirichlet problem form one of the most celebrated topics in the area of applied mathematics. In this study, a novel method is presented for the solution of two-dimensional heat equation for a rectangular plate. In this alternative method, the solution function of the problem is based on the Green function, and therefore on elliptic functions.

Travelling wave solutions and proper solutions to the two-dimensional Burgers-Korteweg-de Vries equation

International Nuclear Information System (INIS)

Feng Zhaosheng

2003-01-01

In this paper, we study the two-dimensional Burgers-Korteweg-de Vries (2D-BKdV) equation by analysing an equivalent two-dimensional autonomous system, which indicates that under some particular conditions, the 2D-BKdV equation has a unique bounded travelling wave solution. Then by using a direct method, a travelling solitary wave solution to the 2D-BKdV equation is expressed explicitly, which appears to be more efficient than the existing methods proposed in the literature. At the end of the paper, the asymptotic behaviour of the proper solutions of the 2D-BKdV equation is established by applying the qualitative theory of differential equations

Two routes to the one-dimensional discrete nonpolynomial Schroedinger equation

International Nuclear Information System (INIS)

Gligoric, G.; Hadzievski, Lj.; Maluckov, A.; Salasnich, L.; Malomed, B. A.

2009-01-01

The Bose-Einstein condensate (BEC), confined in a combination of the cigar-shaped trap and axial optical lattice, is studied in the framework of two models described by two versions of the one-dimensional (1D) discrete nonpolynomial Schroedinger equation (NPSE). Both models are derived from the three-dimensional Gross-Pitaevskii equation (3D GPE). To produce 'model 1' (which was derived in recent works), the 3D GPE is first reduced to the 1D continual NPSE, which is subsequently discretized. 'Model 2,' which was not considered before, is derived by first discretizing the 3D GPE, which is followed by the reduction in the dimension. The two models seem very different; in particular, model 1 is represented by a single discrete equation for the 1D wave function, while model 2 includes an additional equation for the transverse width. Nevertheless, numerical analyses show similar behaviors of fundamental unstaggered solitons in both systems, as concerns their existence region and stability limits. Both models admit the collapse of the localized modes, reproducing the fundamental property of the self-attractive BEC confined in tight traps. Thus, we conclude that the fundamental properties of discrete solitons predicted for the strongly trapped self-attracting BEC are reliable, as the two distinct models produce them in a nearly identical form. However, a difference between the models is found too, as strongly pinned (very narrow) discrete solitons, which were previously found in model 1, are not generated by model 2--in fact, in agreement with the continual 1D NPSE, which does not have such solutions either. In that respect, the newly derived model provides for a more accurate approximation for the trapped BEC.

Analytic energies and wave functions of the two-dimensional Schrodinger equation: ground state of two-dimensional quartic potential and classification of solutions

Czech Academy of Sciences Publication Activity Database

Tichý, V.; Kuběna, Aleš Antonín; Skála, L.

2012-01-01

Roč. 90, č. 6 (2012), s. 503-513 ISSN 0008-4204 Institutional support: RVO:67985556 Keywords : Schroninger equation * partial differential equation * analytic solution * anharmonic oscilator * double-well Subject RIV: BE - Theoretical Physics Impact factor: 0.902, year: 2012 http://library.utia.cas.cz/separaty/2012/E/kubena-analytic energies and wave functions of the two-dimensional schrodinger equation.pdf

Stabilizing local boundary conditions for two-dimensional shallow water equations

KAUST Repository

Dia, Ben Mansour

2018-03-27

In this article, we present a sub-critical two-dimensional shallow water flow regulation. From the energy estimate of a set of one-dimensional boundary stabilization problems, we obtain a set of polynomial equations with respect to the boundary values as a requirement for the energy decrease. Using the Riemann invariant analysis, we build stabilizing local boundary conditions that guarantee the stability of the hydrodynamical state around a given steady state. Numerical results for the controller applied to the nonlinear problem demonstrate the performance of the method.

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

Science.gov (United States)

Gao, Yi; Bouix, Sylvain

2016-05-01

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

Contravariant gravity on Poisson manifolds and Einstein gravity

International Nuclear Information System (INIS)

Kaneko, Yukio; Watamura, Satoshi; Muraki, Hisayoshi

2017-01-01

A relation between gravity on Poisson manifolds proposed in Asakawa et al (2015 Fortschr. Phys . 63 683–704) and Einstein gravity is investigated. The compatibility of the Poisson and Riemann structures defines a unique connection, the contravariant Levi-Civita connection, and leads to the idea of the contravariant gravity. The Einstein–Hilbert-type action yields an equation of motion which is written in terms of the analog of the Einstein tensor, and it includes couplings between the metric and the Poisson tensor. The study of the Weyl transformation reveals properties of those interactions. It is argued that this theory can have an equivalent description as a system of Einstein gravity coupled to matter. As an example, it is shown that the contravariant gravity on a two-dimensional Poisson manifold can be described by a real scalar field coupled to the metric in a specific manner. (paper)

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

Science.gov (United States)

Koehl, Patrice; Delarue, Marc

2010-02-14

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

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

Science.gov (United States)

Gavish, Nir

2018-04-01

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

A New Auto-Baecklund Transformation and Two-Soliton Solution for (3+1)-Dimensional Jimbo-Miwa Equation

International Nuclear Information System (INIS)

Liu Chunping; Zhou Ling

2011-01-01

By improving the extended homogeneous balance method, a general method is suggested to derive a new auto-Baecklund transformation (BT) for (3+1)-Dimensional Jimbo-Miwa (JM) equation. The auto-BT obtained by using our method only involves one quadratic homogeneity equation written as a bilinear equation. Based on the auto-BT, two-soliton solution of the (3+1)-Dimensional JM equation is obtained. (general)

Approximate analytical solution of two-dimensional multigroup P-3 equations

International Nuclear Information System (INIS)

Matausek, M.V.; Milosevic, M.

1981-01-01

Iterative solution of multigroup spherical harmonics equations reduces, in the P-3 approximation and in two-dimensional geometry, to a problem of solving an inhomogeneous system of eight ordinary first order differential equations. With appropriate boundary conditions, these equations have to be solved for each energy group and in each iteration step. The general solution of the corresponding homogeneous system of equations is known in analytical form. The present paper shows how the right-hand side of the system can be approximated in order to derive a particular solution and thus an approximate analytical expression for the general solution of the inhomogeneous system. This combined analytical-numerical approach was shown to have certain advantages compared to the finite-difference method or the Lie-series expansion method, which have been used to solve similar problems. (author)

Comparative study of the two-fluid momentum equations for multi-dimensional bubbly flows: Modification of Reynolds stress

Energy Technology Data Exchange (ETDEWEB)

Lee, Seung Jun; Park, Ik Kyu; Yoon, Han Young [Thermal-Hydraulic Safety Research Division, Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of); Jae, Byoung [School of Mechanical Engineering, Chungnam National University, Daejeon (Korea, Republic of)

2017-01-15

Two-fluid equations are widely used to obtain averaged behaviors of two-phase flows. This study addresses a problem that may arise when the two-fluid equations are used for multi-dimensional bubbly flows. If steady drag is the only accounted force for the interfacial momentum transfer, the disperse-phase velocity would be the same as the continuous-phase velocity when the flow is fully developed without gravity. However, existing momentum equations may show unphysical results in estimating the relative velocity of the disperse phase against the continuous-phase. First, we examine two types of existing momentum equations. One is the standard two-fluid momentum equation in which the disperse-phase is treated as a continuum. The other is the averaged momentum equation derived from a solid/ fluid particle motion. We show that the existing equations are not proper for multi-dimensional bubbly flows. To resolve the problem mentioned above, we modify the form of the Reynolds stress terms in the averaged momentum equation based on the solid/fluid particle motion. The proposed equation shows physically correct results for both multi-dimensional laminar and turbulent flows.

Two-dimensional electron flow in pulsed power transmission lines and plasma opening switches

International Nuclear Information System (INIS)

Church, B.W.; Longcope, D.W.; Ng, C.K.; Sudan, R.N.

1991-01-01

The operation of magnetically insulated transmission lines (MITL) and the interruption of current in a plasma opening switch (POS) are determined by the physics of the electrons emitted by the cathode surface. A mathematical model describes the self-consistent two-dimensional flow of an electron fluid. A finite element code, FERUS, has been developed to solve the two equations which describe Poisson's and Ampere's law in two dimensions. The solutions from this code are obtained for parameters where the electron orbits are considerably modified by the self-magnetic field of the current. Next, the self-insulated electron flow in a MITL with a step change in cross-section is studied using a conventional two-dimensional fully electromagnetic particle-in-cell code, MASK. The equations governing two-dimensional quasi-static electron flow are solved numerically by a third technique which is suitable for predicting current interruption in a POS. The object of the study is to determine the critical load impedance, Z CL , required for current interruption for a given applied voltage, cathode voltage and plasma length. (author). 9 refs, 5 figs

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.

Appearance of eigen modes for the linearized Vlasov-Poisson equation

International Nuclear Information System (INIS)

Degond, P.

1983-01-01

In order to determine the asymptotic behaviour, when the time goes to infinity, of the solution of the linearized Vlasov-Poisson equation, we use eigen modes, associated to continuous linear functionals on a Banach space of analytic functions [fr

Family of two-dimensional Born-Infeld equations and a system of conservation laws

International Nuclear Information System (INIS)

Koiv, M.; Rosenhaus, V.

1979-01-01

Lower-order conserved quantities, the ''currents'', for two-dimensional Lorentz-invariant Born-Infeld equation are considered. The currents are divided into pairs, which contain a class (basic currents) leading to the family equations. The basic currents determine the transformations between the solutions of the Born-Infeld eqution family. The equations belonging to the family are fully hodograph-invariant, partly hodograph-invariant, and effectively linear, i.e. non-linear equations with linear image of hodograph transformation

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

DEFF Research Database (Denmark)

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

2008-01-01

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

A modified Poisson-Boltzmann equation applied to protein adsorption.

Science.gov (United States)

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

2018-01-05

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

Approximate analytical solution of two-dimensional multigroup P-3 equations

International Nuclear Information System (INIS)

Matausek, M.V.; Milosevic, M.

1981-01-01

Iterative solution of multigroup spherical harmonics equations reduces, in the P-3 approximation and in two-dimensional geometry, to a problem of solving an inhomogeneous system of eight ordinary first order differential equations. With appropriate boundary conditions, these equations have to be solved for each energy group and in each iteration step. The general solution of the corresponding homogeneous system of equations is known in analytical form. The present paper shows how the right-hand side of the system can be approximated in order to derive a particular solution and thus an approximate analytical expression for the general solution of the inhomogeneous system. This combined analytical-numerical approach was shown to have certain advantages compared to the finite-difference method or the Lie-series expansion method, which have been used to solve similar problems. (orig./RW) [de

The Poisson equation at second order in relativistic cosmology

International Nuclear Information System (INIS)

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

2013-01-01

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

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

Science.gov (United States)

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

2010-09-20

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

A fast semi-discrete Kansa method to solve the two-dimensional spatiotemporal fractional diffusion equation

Science.gov (United States)

Sun, HongGuang; Liu, Xiaoting; Zhang, Yong; Pang, Guofei; Garrard, Rhiannon

2017-09-01

Fractional-order diffusion equations (FDEs) extend classical diffusion equations by quantifying anomalous diffusion frequently observed in heterogeneous media. Real-world diffusion can be multi-dimensional, requiring efficient numerical solvers that can handle long-term memory embedded in mass transport. To address this challenge, a semi-discrete Kansa method is developed to approximate the two-dimensional spatiotemporal FDE, where the Kansa approach first discretizes the FDE, then the Gauss-Jacobi quadrature rule solves the corresponding matrix, and finally the Mittag-Leffler function provides an analytical solution for the resultant time-fractional ordinary differential equation. Numerical experiments are then conducted to check how the accuracy and convergence rate of the numerical solution are affected by the distribution mode and number of spatial discretization nodes. Applications further show that the numerical method can efficiently solve two-dimensional spatiotemporal FDE models with either a continuous or discrete mixing measure. Hence this study provides an efficient and fast computational method for modeling super-diffusive, sub-diffusive, and mixed diffusive processes in large, two-dimensional domains with irregular shapes.

Steady state solution of the Poisson-Nernst-Planck equations

International Nuclear Information System (INIS)

Golovnev, A.; Trimper, S.

2010-01-01

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

New lumps of Veselov-Novikov integrable nonlinear equation and new exact rational potentials of two-dimensional stationary Schroedinger equation via ∂-macron-dressing method

International Nuclear Information System (INIS)

Dubrovsky, V.G.; Formusatik, I.B.

2003-01-01

The scheme for calculating via Zakharov-Manakov ∂-macron-dressing method of new rational solutions with constant asymptotic values at infinity of the famous two-dimensional Veselov-Novikov (VN) integrable nonlinear evolution equation and new exact rational potentials of two-dimensional stationary Schroedinger (2DSchr) equation with multiple pole wave functions is developed. As examples new lumps of VN nonlinear equation and new exact rational potentials of 2DSchr equation with multiple pole of order two wave functions are calculated. Among the constructed rational solutions are as nonsingular and also singular

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

DEFF Research Database (Denmark)

Hejlesen, Mads Mølholm; Walther, Jens Honore

2016-01-01

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

Vapour-liquid equilibrium properties for two- and three-dimensional Lennard-Jones fluids from equations of state

International Nuclear Information System (INIS)

Mulero, A.; Cuadros, F; Faundez, C.A.

1999-01-01

Vapour-liquid equilibrium properties for both three- and two-dimensional Lennard-Jones fluids were obtained using simple cubic-in-density equations of state proposed by the authors. Results were compared with those obtained by other workers from computer simulations and also with results given by other more complex semi-theoretical or semi-empirical equations of state. In the three-dimensional case good agreement is found for all properties and all temperatures. In the two-dimensional case only the coexistence densities were compared, producing good agreement for low temperatures only. The present work is the first to give numerical data for the vapour-liquid equilibrium properties of Lennard-Jones fluids calculated from equations of state. Copyright (1999) CSIRO Australia

A Kronecker product splitting preconditioner for two-dimensional space-fractional diffusion equations

Science.gov (United States)

Chen, Hao; Lv, Wen; Zhang, Tongtong

2018-05-01

We study preconditioned iterative methods for the linear system arising in the numerical discretization of a two-dimensional space-fractional diffusion equation. Our approach is based on a formulation of the discrete problem that is shown to be the sum of two Kronecker products. By making use of an alternating Kronecker product splitting iteration technique we establish a class of fixed-point iteration methods. Theoretical analysis shows that the new method converges to the unique solution of the linear system. Moreover, the optimal choice of the involved iteration parameters and the corresponding asymptotic convergence rate are computed exactly when the eigenvalues of the system matrix are all real. The basic iteration is accelerated by a Krylov subspace method like GMRES. The corresponding preconditioner is in a form of a Kronecker product structure and requires at each iteration the solution of a set of discrete one-dimensional fractional diffusion equations. We use structure preserving approximations to the discrete one-dimensional fractional diffusion operators in the action of the preconditioning matrix. Numerical examples are presented to illustrate the effectiveness of this approach.

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

Directory of Open Access Journals (Sweden)

Diem Dang Huan

2015-12-01

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

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

Science.gov (United States)

Chadha, Alka; Bora, Swaroop Nandan

2017-11-01

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

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

Science.gov (United States)

Mullenmeister, Paul

1988-01-01

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

Modified Splitting FDTD Methods for Two-Dimensional Maxwell’s Equations

Directory of Open Access Journals (Sweden)

Liping Gao

2017-01-01

Full Text Available In this paper, we develop a new method to reduce the error in the splitting finite-difference method of Maxwell’s equations. By this method two modified splitting FDTD methods (MS-FDTDI, MS-FDTDII for the two-dimensional Maxwell equations are proposed. It is shown that the two methods are second-order accurate in time and space and unconditionally stable by Fourier methods. By energy method, it is proved that MS-FDTDI is second-order convergent. By deriving the numerical dispersion (ND relations, we prove rigorously that MS-FDTDI has less ND errors than the ADI-FDTD method and the ND errors of ADI-FDTD are less than those of MS-FDTDII. Numerical experiments for computing ND errors and simulating a wave guide problem and a scattering problem are carried out and the efficiency of the MS-FDTDI and MS-FDTDII methods is confirmed.

Comparison of preconditioned generalized conjugate gradient methods to two-dimensional neutron and photon transport equation

International Nuclear Information System (INIS)

Chen, G.S.

1997-01-01

We apply and compare the preconditioned generalized conjugate gradient methods to solve the linear system equation that arises in the two-dimensional neutron and photon transport equation in this paper. Several subroutines are developed on the basis of preconditioned generalized conjugate gradient methods for time-independent, two-dimensional neutron and photon transport equation in the transport theory. These generalized conjugate gradient methods are used. TFQMR (transpose free quasi-minimal residual algorithm), CGS (conjuage gradient square algorithm), Bi-CGSTAB (bi-conjugate gradient stabilized algorithm) and QMRCGSTAB (quasi-minimal residual variant of bi-conjugate gradient stabilized algorithm). These sub-routines are connected to computer program DORT. Several problems are tested on a personal computer with Intel Pentium CPU. (author)

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

KAUST Repository

Majeed, Muhammad Usman

2017-07-10

A state-observer based method is developed to solve point source localization problem for Poisson equation in a 3D rectangular prism with available boundary data. The technique requires a weighted sum of solutions of multiple boundary data estimation problems for Laplace equation over the 3D domain. The solution of each of these boundary estimation problems involves writing down the mathematical problem in state-space-like representation using one of the space variables as time-like. First, system observability result for 3D boundary estimation problem is recalled in an infinite dimensional setting. Then, based on the observability result, the boundary estimation problem is decomposed into a set of independent 2D sub-problems. These 2D problems are then solved using an iterative observer to obtain the solution. Theoretical results are provided. The method is implemented numerically using finite difference discretization schemes. Numerical illustrations along with simulation results are provided.

One- and Two-dimensional Solitary Wave States in the Nonlinear Kramers Equation with Movement Direction as a Variable

Science.gov (United States)

Sakaguchi, Hidetsugu; Ishibashi, Kazuya

2018-06-01

We study self-propelled particles by direct numerical simulation of the nonlinear Kramers equation for self-propelled particles. In our previous paper, we studied self-propelled particles with velocity variables in one dimension. In this paper, we consider another model in which each particle exhibits directional motion. The movement direction is expressed with a variable ϕ. We show that one-dimensional solitary wave states appear in direct numerical simulations of the nonlinear Kramers equation in one- and two-dimensional systems, which is a generalization of our previous result. Furthermore, we find two-dimensionally localized states in the case that each self-propelled particle exhibits rotational motion. The center of mass of the two-dimensionally localized state exhibits circular motion, which implies collective rotating motion. Finally, we consider a simple one-dimensional model equation to qualitatively understand the formation of the solitary wave state.

Numerical solution of multi group-Two dimensional- Adjoint equation with finite element method

International Nuclear Information System (INIS)

Poursalehi, N.; Khalafi, H.; Shahriari, M.; Minoochehr

2008-01-01

Adjoint equation is used for perturbation theory in nuclear reactor design. For numerical solution of adjoint equation, usually two methods are applied. These are Finite Element and Finite Difference procedures. Usually Finite Element Procedure is chosen for solving of adjoint equation, because it is more use able in variety of geometries. In this article, Galerkin Finite Element method is discussed. This method is applied for numerical solving multi group, multi region and two dimensional (X, Y) adjoint equation. Typical reactor geometry is partitioned with triangular meshes and boundary condition for adjoint flux is considered zero. Finally, for a case of defined parameters, Finite Element Code was applied and results were compared with Citation Code

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

Science.gov (United States)

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

2001-09-01

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

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

International Nuclear Information System (INIS)

Balian, R.; Veneroni, M.

1989-01-01

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

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

International Nuclear Information System (INIS)

Veneroni, M.; Balian, R.

1989-01-01

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

ColDICE: A parallel Vlasov–Poisson solver using moving adaptive simplicial tessellation

International Nuclear Information System (INIS)

Sousbie, Thierry; Colombi, Stéphane

2016-01-01

Resolving numerically Vlasov–Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm consisting in representing the phase-space sheet with a conforming, self-adaptive simplicial tessellation of which the vertices follow the Lagrangian equations of motion. The algorithm is implemented both in six- and four-dimensional phase-space. Refinement of the tessellation mesh is performed using the bisection method and a local representation of the phase-space sheet at second order relying on additional tracers created when needed at runtime. In order to preserve in the best way the Hamiltonian nature of the system, refinement is anisotropic and constrained by measurements of local Poincaré invariants. Resolution of Poisson equation is performed using the fast Fourier method on a regular rectangular grid, similarly to particle in cells codes. To compute the density projected onto this grid, the intersection of the tessellation and the grid is calculated using the method of Franklin and Kankanhalli [65–67] generalised to linear order. As preliminary tests of the code, we study in four dimensional phase-space the evolution of an initially small patch in a chaotic potential and the cosmological collapse of a fluctuation composed of two sinusoidal waves. We also perform a “warm” dark matter simulation in six-dimensional phase-space that we use to check the parallel scaling of the code.

ColDICE: A parallel Vlasov–Poisson solver using moving adaptive simplicial tessellation

Energy Technology Data Exchange (ETDEWEB)

Sousbie, Thierry, E-mail: tsousbie@gmail.com [Institut d' Astrophysique de Paris, CNRS UMR 7095 and UPMC, 98bis, bd Arago, F-75014 Paris (France); Department of Physics, The University of Tokyo, Tokyo 113-0033 (Japan); Research Center for the Early Universe, School of Science, The University of Tokyo, Tokyo 113-0033 (Japan); Colombi, Stéphane, E-mail: colombi@iap.fr [Institut d' Astrophysique de Paris, CNRS UMR 7095 and UPMC, 98bis, bd Arago, F-75014 Paris (France); Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502 (Japan)

2016-09-15

Resolving numerically Vlasov–Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm consisting in representing the phase-space sheet with a conforming, self-adaptive simplicial tessellation of which the vertices follow the Lagrangian equations of motion. The algorithm is implemented both in six- and four-dimensional phase-space. Refinement of the tessellation mesh is performed using the bisection method and a local representation of the phase-space sheet at second order relying on additional tracers created when needed at runtime. In order to preserve in the best way the Hamiltonian nature of the system, refinement is anisotropic and constrained by measurements of local Poincaré invariants. Resolution of Poisson equation is performed using the fast Fourier method on a regular rectangular grid, similarly to particle in cells codes. To compute the density projected onto this grid, the intersection of the tessellation and the grid is calculated using the method of Franklin and Kankanhalli [65–67] generalised to linear order. As preliminary tests of the code, we study in four dimensional phase-space the evolution of an initially small patch in a chaotic potential and the cosmological collapse of a fluctuation composed of two sinusoidal waves. We also perform a “warm” dark matter simulation in six-dimensional phase-space that we use to check the parallel scaling of the code.

The linearized pressure Poisson equation for global instability analysis of incompressible flows

Science.gov (United States)

Theofilis, Vassilis

2017-12-01

The linearized pressure Poisson equation (LPPE) is used in two and three spatial dimensions in the respective matrix-forming solution of the BiGlobal and TriGlobal eigenvalue problem in primitive variables on collocated grids. It provides a disturbance pressure boundary condition which is compatible with the recovery of perturbation velocity components that satisfy exactly the linearized continuity equation. The LPPE is employed to analyze instability in wall-bounded flows and in the prototype open Blasius boundary layer flow. In the closed flows, excellent agreement is shown between results of the LPPE and those of global linear instability analyses based on the time-stepping nektar++, Semtex and nek5000 codes, as well as with those obtained from the FreeFEM++ matrix-forming code. In the flat plate boundary layer, solutions extracted from the two-dimensional LPPE eigenvector at constant streamwise locations are found to be in very good agreement with profiles delivered by the NOLOT/PSE space marching code. Benchmark eigenvalue data are provided in all flows analyzed. The performance of the LPPE is seen to be superior to that of the commonly used pressure compatibility (PC) boundary condition: at any given resolution, the discrete part of the LPPE eigenspectrum contains converged and not converged, but physically correct, eigenvalues. By contrast, the PC boundary closure delivers some of the LPPE eigenvalues and, in addition, physically wrong eigenmodes. It is concluded that the LPPE should be used in place of the PC pressure boundary closure, when BiGlobal or TriGlobal eigenvalue problems are solved in primitive variables by the matrix-forming approach on collocated grids.

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

Science.gov (United States)

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

2013-01-01

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

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

International Nuclear Information System (INIS)

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

1974-01-01

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

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

International Nuclear Information System (INIS)

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

1997-01-01

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

Backlund transformations and three-dimensional lattice equations

NARCIS (Netherlands)

Nijhoff, F.W.; Capel, H.W.; Wiersma, G.L.; Quispel, G.R.W.

1984-01-01

A (nonlocal) linear integral equation is studied, which allows for Bäcklund transformations in the measure. The compatibility of three of these transformations leads to an integrable nonlinear three-dimensional lattice equation. In appropriate continuum limits the two-dimensional Toda-lattice

Generalized results on the role of new-time transformations in finite-dimensional Poisson systems

Energy Technology Data Exchange (ETDEWEB)

Hernandez-Bermejo, Benito, E-mail: benito.hernandez@urjc.e [Departamento de Fisica, Escuela Superior de Ciencias Experimentales y Tecnologia, Universidad Rey Juan Carlos, Calle Tulipan S/N, 28933 Mostoles, Madrid (Spain)

2010-01-25

The problem of characterizing all new-time transformations preserving the Poisson structure of a finite-dimensional Poisson system is completely solved in a constructive way. As a corollary, this leads to a broad generalization of previously known results. Examples are given.

Poisson-Box Sampling algorithms for three-dimensional Markov binary mixtures

Science.gov (United States)

Larmier, Coline; Zoia, Andrea; Malvagi, Fausto; Dumonteil, Eric; Mazzolo, Alain

2018-02-01

Particle transport in Markov mixtures can be addressed by the so-called Chord Length Sampling (CLS) methods, a family of Monte Carlo algorithms taking into account the effects of stochastic media on particle propagation by generating on-the-fly the material interfaces crossed by the random walkers during their trajectories. Such methods enable a significant reduction of computational resources as opposed to reference solutions obtained by solving the Boltzmann equation for a large number of realizations of random media. CLS solutions, which neglect correlations induced by the spatial disorder, are faster albeit approximate, and might thus show discrepancies with respect to reference solutions. In this work we propose a new family of algorithms (called 'Poisson Box Sampling', PBS) aimed at improving the accuracy of the CLS approach for transport in d-dimensional binary Markov mixtures. In order to probe the features of PBS methods, we will focus on three-dimensional Markov media and revisit the benchmark problem originally proposed by Adams, Larsen and Pomraning [1] and extended by Brantley [2]: for these configurations we will compare reference solutions, standard CLS solutions and the new PBS solutions for scalar particle flux, transmission and reflection coefficients. PBS will be shown to perform better than CLS at the expense of a reasonable increase in computational time.

Comparison of three-dimensional poisson solution methods for particle-based simulation and inhomogeneous dielectrics.

Science.gov (United States)

Berti, Claudio; Gillespie, Dirk; Bardhan, Jaydeep P; Eisenberg, Robert S; Fiegna, Claudio

2012-07-01

Particle-based simulation represents a powerful approach to modeling physical systems in electronics, molecular biology, and chemical physics. Accounting for the interactions occurring among charged particles requires an accurate and efficient solution of Poisson's equation. For a system of discrete charges with inhomogeneous dielectrics, i.e., a system with discontinuities in the permittivity, the boundary element method (BEM) is frequently adopted. It provides the solution of Poisson's equation, accounting for polarization effects due to the discontinuity in the permittivity by computing the induced charges at the dielectric boundaries. In this framework, the total electrostatic potential is then found by superimposing the elemental contributions from both source and induced charges. In this paper, we present a comparison between two BEMs to solve a boundary-integral formulation of Poisson's equation, with emphasis on the BEMs' suitability for particle-based simulations in terms of solution accuracy and computation speed. The two approaches are the collocation and qualocation methods. Collocation is implemented following the induced-charge computation method of D. Boda et al. [J. Chem. Phys. 125, 034901 (2006)]. The qualocation method is described by J. Tausch et al. [IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems 20, 1398 (2001)]. These approaches are studied using both flat and curved surface elements to discretize the dielectric boundary, using two challenging test cases: a dielectric sphere embedded in a different dielectric medium and a toy model of an ion channel. Earlier comparisons of the two BEM approaches did not address curved surface elements or semiatomistic models of ion channels. Our results support the earlier findings that for flat-element calculations, qualocation is always significantly more accurate than collocation. On the other hand, when the dielectric boundary is discretized with curved surface elements, the

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

International Nuclear Information System (INIS)

Milgrom, M.

1997-01-01

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

On the biophysics of cathodal galvanotaxis in rat prostate cancer cells: Poisson-Nernst-Planck equation approach.

Science.gov (United States)

Borys, Przemysław

2012-06-01

Rat prostate cancer cells have been previously investigated using two cell lines: a highly metastatic one (Mat-Ly-Lu) and a nonmetastatic one (AT-2). It turns out that the highly metastatic Mat-Ly-Lu cells exhibit a phenomenon of cathodal galvanotaxis in an electric field which can be blocked by interrupting the voltage-gated sodium channel (VGSC) activity. The VGSC activity is postulated to be characteristic for metastatic cells and seems to be a reasonable driving force for motile behavior. However, the classical theory of cellular motion depends on calcium ions rather than sodium ions. The current research provides a theoretical connection between cellular sodium inflow and cathodal galvanotaxis of Mat-Ly-Lu cells. Electrical repulsion of intracellular calcium ions by entering sodium ions is proposed after depolarization starting from the cathodal side. The disturbance in the calcium distribution may then drive actin polymerization and myosin contraction. The presented modeling is done within a continuous one-dimensional Poisson-Nernst-Planck equation framework.

Matrix-type multiple reciprocity boundary element method for solving three-dimensional two-group neutron diffusion equations

International Nuclear Information System (INIS)

Itagaki, Masafumi; Sahashi, Naoki.

1997-01-01

The multiple reciprocity boundary element method has been applied to three-dimensional two-group neutron diffusion problems. A matrix-type boundary integral equation has been derived to solve the first and the second group neutron diffusion equations simultaneously. The matrix-type fundamental solutions used here satisfy the equation which has a point source term and is adjoint to the neutron diffusion equations. A multiple reciprocity method has been employed to transform the matrix-type domain integral related to the fission source into an equivalent boundary one. The higher order fundamental solutions required for this formulation are composed of a series of two types of analytic functions. The eigenvalue itself is also calculated using only boundary integrals. Three-dimensional test calculations indicate that the present method provides stable and accurate solutions for criticality problems. (author)

The Fractional Poisson Process and the Inverse Stable Subordinator

OpenAIRE

Meerschaert, Mark; Nane, Erkan; Vellaisamy, P.

2011-01-01

The fractional Poisson process is a renewal process with Mittag-Leffler waiting times. Its distributions solve a time-fractional analogue of the Kolmogorov forward equation for a Poisson process. This paper shows that a traditional Poisson process, with the time variable replaced by an independent inverse stable subordinator, is also a fractional Poisson process. This result unifies the two main approaches in the stochastic theory of time-fractional diffusion equations. The equivalence extend...

Comparison of Poisson structures and Poisson-Lie dynamical r-matrices

OpenAIRE

Enriquez, B.; Etingof, P.; Marshall, I.

2004-01-01

We construct a Poisson isomorphism between the formal Poisson manifolds g^* and G^*, where g is a finite dimensional quasitriangular Lie bialgebra. Here g^* is equipped with its Lie-Poisson (or Kostant-Kirillov-Souriau) structure, and G^* with its Poisson-Lie structure. We also quantize Poisson-Lie dynamical r-matrices of Balog-Feher-Palla.

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

Czech Academy of Sciences Publication Activity Database

Medková, Dagmar

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

Soliton solutions of the two-dimensional KdV-Burgers equation by homotopy perturbation method

International Nuclear Information System (INIS)

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

2007-01-01

In this Letter, the He's homotopy perturbation method (HPM) to finding the soliton solutions of the two-dimensional Korteweg-de Vries Burgers' equation (tdKdVB) for the initial conditions was applied. Numerical solutions of the equation were obtained. The obtained solutions, in comparison with the exact solutions admit a remarkable accuracy. The results reveal that the HPM is very effective and simple

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

Science.gov (United States)

Fellner, Klemens; Kovtunenko, Victor A

2016-01-01

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

Dynamics in discrete two-dimensional nonlinear Schrödinger equations in the presence of point defects

DEFF Research Database (Denmark)

Christiansen, Peter Leth; Gaididei, Yuri Borisovich; Rasmussen, Kim

1996-01-01

The dynamics of two-dimensional discrete structures is studied in the framework of the generalized two-dimensional discrete nonlinear Schrodinger equation. The nonlinear coupling in the form of the Ablowitz-Ladik nonlinearity and point impurities is taken into account. The stability properties...... of the stationary solutions are examined. The essential importance of the existence of stable immobile solitons in the two-dimensional dynamics of the traveling pulses is demonstrated. The typical scenario of the two-dimensional quasicollapse of a moving intense pulse represents the formation of standing trapped...... narrow spikes. The influence of the point impurities on this dynamics is also investigated....

Parallel SOR methods with a parabolic-diffusion acceleration technique for solving an unstructured-grid Poisson equation on 3D arbitrary geometries

Science.gov (United States)

Zapata, M. A. Uh; Van Bang, D. Pham; Nguyen, K. D.

2016-05-01

This paper presents a parallel algorithm for the finite-volume discretisation of the Poisson equation on three-dimensional arbitrary geometries. The proposed method is formulated by using a 2D horizontal block domain decomposition and interprocessor data communication techniques with message passing interface. The horizontal unstructured-grid cells are reordered according to the neighbouring relations and decomposed into blocks using a load-balanced distribution to give all processors an equal amount of elements. In this algorithm, two parallel successive over-relaxation methods are presented: a multi-colour ordering technique for unstructured grids based on distributed memory and a block method using reordering index following similar ideas of the partitioning for structured grids. In all cases, the parallel algorithms are implemented with a combination of an acceleration iterative solver. This solver is based on a parabolic-diffusion equation introduced to obtain faster solutions of the linear systems arising from the discretisation. Numerical results are given to evaluate the performances of the methods showing speedups better than linear.

Effective electrodiffusion equation for non-uniform nanochannels.

Science.gov (United States)

Marini Bettolo Marconi, Umberto; Melchionna, Simone; Pagonabarraga, Ignacio

2013-06-28

We derive a one-dimensional formulation of the Planck-Nernst-Poisson equation to describe the dynamics of a symmetric binary electrolyte in channels whose section is nanometric and varies along the axial direction. The approach is in the spirit of the Fick-Jacobs diffusion equation and leads to a system of coupled equations for the partial densities which depends on the charge sitting at the walls in a non-trivial fashion. We consider two kinds of non-uniformities, those due to the spatial variation of charge distribution and those due to the shape variation of the pore and report one- and three-dimensional solutions of the electrokinetic equations.

Travelling wave solutions of two-dimensional Korteweg-de Vries-Burgers and Kadomtsev-Petviashvili equations

International Nuclear Information System (INIS)

Estevez, P G; Kuru, S; Negro, J; Nieto, L M

2006-01-01

The travelling wave solutions of the two-dimensional Korteweg-de Vries-Burgers and Kadomtsev-Petviashvili equations are studied from two complementary points of view. The first one is an adaptation of the factorization technique that provides particular as well as general solutions. The second one applies the Painleve analysis to both equations, throwing light on some aspects of the first method and giving an explanation to some restriction on the coefficients, as well as the relation between factorizations and integrals of motion

Symplectic and Hamiltonian structures of nonlinear evolution equations

International Nuclear Information System (INIS)

Dorfman, I.Y.

1993-01-01

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

A two-dimensional analytical subthreshold behavior model for junctionless dual-material cylindrical surrounding-gate MOSFETs

International Nuclear Information System (INIS)

Li Cong; Zhuang Yi-Qi; Zhang Li; Jin Gang

2014-01-01

A two-dimensional analytical subthreshold behavior model for junctionless dual-material cylindrical surrounding-gate (JLDMCSG) metal-oxide-semiconductor field-effect transistors (MOSFETs) is proposed. It is derived by solving the two-dimensional Poisson's equation in two continuous cylindrical regions with any simplifying assumption. Using this analytical model, the subthreshold characteristics of JLDMCSG MOSFETs are investigated in terms of channel electrostatic potential, horizontal electric field, and subthreshold current. Compared to junctionless single-material cylindrical surrounding-gate MOSFETs, JLDMCSG MOSFETs can effectively suppress short-channel effects and simultaneously improve carrier transport efficiency. It is found that the subthreshold current of JLDMCSG MOSFETs can be significantly reduced by adopting both a thin oxide and thin silicon channel. The accuracy of the analytical model is verified by its good agreement with the three-dimensional numerical simulator ISE TCAD

Fourth-order poisson solver for the simulation of bounded plasmas

International Nuclear Information System (INIS)

Knorr, G.; Joyce, G.; Marcus, A.J.

1980-01-01

The solution of the two-dimensional Poisson equation in a rectangle with periodic boundaries in one direction and Dirichlet or Neumann boundaries in the other can be handled by a Fast Fourier Transform in one dimension and a fast nonperiodic procedure such as splines in the other. Such a solution is necessary for the simulation of semiperiodic plasma systems. A method is presented which is direct and of fourth order in both the electric potential and the electric fields

Two-dimensional hydrodynamics of uniform ion plasma in electrostatic field

International Nuclear Information System (INIS)

Mahdieh, M. H.; Gavili, A.

2005-01-01

Two-dimensional hydrodynamics of ion extraction from uniform quasi-neutral plasma, in electrostatic field has been simulated numerically. Experimentally, tunable pulsed lasers produce non-uniform plasma through stepwise photo-excitation and photo-ionization or multi-photo-ionization processes. Poisson's equation was solved simultaneously with the equations of mass, and momentum, assuming the Maxwell-Boltzmann distribution for electrons. In the calculation, the initial density profile at the boundaries has been assumed to be very steep for the ion plasma. In these calculations dynamics of electric potential and the ions density were assessed. The ion extraction time was also estimated from the calculation. The knowledge of spatial distribution of the ions across the cathode is very important for the practical purposes. In this simulation, the spatial distribution of the ion current density across the cathode as well as its temporal distribution was calculated

Solution of Schroedinger Equation for Two-Dimensional Complex Quartic Potentials

International Nuclear Information System (INIS)

Singh, Ram Mehar; Chand, Fakir; Mishra, S. C.

2009-01-01

We investigate the quasi-exact solutions of the Schroedinger wave equation for two-dimensional non-hermitian complex Hamiltonian systems within the frame work of an extended complex phase space characterized by x = x 1 + ip 3 , y = x 2 + ip 4 , p x = p 1 + ix 3 , p y = p 2 + ix 4 . Explicit expressions of the energy eigenvalues and the eigenfunctions for ground and first excited states for a complex quartic potential are obtained. Eigenvalue spectra of some variants of the complex quartic potential, including PT-symmetric one, are also worked out. (general)

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

International Nuclear Information System (INIS)

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

2016-01-01

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

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

Science.gov (United States)

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

2016-01-07

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

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

Energy Technology Data Exchange (ETDEWEB)

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

2016-01-07

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

A Seemingly Unrelated Poisson Regression Model

OpenAIRE

King, Gary

1989-01-01

This article introduces a new estimator for the analysis of two contemporaneously correlated endogenous event count variables. This seemingly unrelated Poisson regression model (SUPREME) estimator combines the efficiencies created by single equation Poisson regression model estimators and insights from "seemingly unrelated" linear regression models.

Lagrangian analysis of invariant third-order equations of motion in relativistic classical particle mechanics

International Nuclear Information System (INIS)

Matsyuk, R.Ya.

1985-01-01

The problem on the existence of the invariant third-order Euler-Poisson equations in the pseudo-Euclidean space is investigated. The locally variational problem is determined by the Lagrangian density over the space of the second-order jets. The one - parameter family of the invariant third-order Euler-Poisson equations is groved to be the only one in the three-dimensional pseudo-Euclidean space. No invariant third-order Euler-Poisson equations exist in the four-dimensional pseudo-Euclidean space. It is shown that the Mathisson equation and the equation of geodesic circles in particular cases may be considered in the context of the Ostrogradiskij mechanics and the Kavaguchi geometry

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

Science.gov (United States)

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

1990-01-01

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

Non-chiral, molecular model of negative Poisson ratio in two dimensions

International Nuclear Information System (INIS)

Wojciechowski, K W

2003-01-01

A two-dimensional model of tri-atomic molecules (in which 'atoms' are distributed on vertices of equilateral triangles, and which are further referred to as cyclic trimers) is solved exactly in the static (zero-temperature) limit for the nearest-neighbour site-site interactions. It is shown that the cyclic trimers form a mechanically stable and elastically isotropic non-chiral phase of negative Poisson ratio. The properties of the system are illustrated by three examples of atom-atom interaction potentials: (i) the purely repulsive (n-inverse-power) potential, (ii) the purely attractive (n-power) potential and (iii) the Lennard-Jones potential which shows both the repulsive and the attractive part. The analytic form of the dependence of the Poisson ratio on the interatomic potential is obtained. It is shown that the Poisson ratio depends, in a universal way, only on the trimer anisotropy parameter both (1) in the limit of n → ∞ for cases (i) and (ii), as well as (2) at the zero external pressure for any potential with a doubly differentiable minimum, case (iii) is an example

Solutions of diffusion equations in two-dimensional cylindrical geometry by series expansions

International Nuclear Information System (INIS)

Ohtani, Nobuo

1976-01-01

A solution of the multi-group multi-regional diffusion equation in two-dimensional cylindrical (rho-z) geometry is obtained in the form of a regionwise double series composed of Bessel and trigonometrical functions. The diffusion equation is multiplied by weighting functions, which satisfy the homogeneous part of the diffusion equation, and the products are integrated over the region for obtaining the equations to determine the fluxes and their normal derivatives at the region boundaries. Multiplying the diffusion equation by each function of the set used for the flux expansion, then integrating the products, the coefficients of the double series of the flux inside each region are calculated using the boundary values obtained above. Since the convergence of the series thus obtained is slow especially near the region boundaries, a method for improving the convergence has been developed. The double series of the flux is separated into two parts. The normal derivative at the region boundary of the first part is zero, and that of the second part takes the value which is obtained in the first stage of this method. The second part is replaced by a continuous function, and the flux is represented by the sum of the continuous function and the double series. A sample critical problem of a two-group two-region system is numerically studied. The results show that the present method yields very accurately the flux integrals in each region with only a small number of expansion terms. (auth.)

An Auxiliary Equation for the Bellman Equation in a One-Dimensional Ergodic Control

International Nuclear Information System (INIS)

Fujita, Y.

2001-01-01

In this paper we consider the Bellman equation in a one-dimensional ergodic control. Our aim is to show the existence and the uniqueness of its solution under general assumptions. For this purpose we introduce an auxiliary equation whose solution gives the invariant measure of the diffusion corresponding to an optimal control. Using this solution, we construct a solution to the Bellman equation. Our method of using this auxiliary equation has two advantages in the one-dimensional case. First, we can solve the Bellman equation under general assumptions. Second, this auxiliary equation gives an optimal Markov control explicitly in many examples

Solution of two-dimensional equations of neutron transport in 4P0-approximation of spherical harmonics method

International Nuclear Information System (INIS)

Polivanskij, V.P.

1989-01-01

The method to solve two-dimensional equations of neutron transport using 4P 0 -approximation is presented. Previously such approach was efficiently used for the solution of one-dimensional problems. New an attempt is made to apply the approach to solution of two-dimensional problems. Algorithm of the solution is given, as well as results of test neutron-physical calculations. A considerable as compared with diffusion approximation is shown. 11 refs

Asymptotic solution of the Vlasov and Poisson equations for an inhomogeneous plasma

International Nuclear Information System (INIS)

Croci, R.

1991-01-01

The asymptotic solutions to a class of inhomogeneous integral equations that reduce to algebraic equations when a parameter η goes to zero (the kernel becoming proportional to a Dirac δ function) are derived. This class includes the integral equations obtained from the system of Vlasov and Poisson equations for the Fourier transform in space and the Laplace transform in time of the electrostatic potential, when the equilibrium magnetic field is uniform and the equilibrium plasma density depends on ηx, with the co-ordinate z being the direction of the magnetic field. In this case the inhomogeneous term is given by the initial conditions and possibly by sources, and the Laplace-transform variable ω is the eigenvalue parameter. (Author)

Particular solutions of generalized Euler-Poisson-Darboux equation

Directory of Open Access Journals (Sweden)

Rakhila B. Seilkhanova

2015-01-01

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

On a method of construction of exact solutions for equations of two-dimensional hydrodynamics of incompressible liquids

International Nuclear Information System (INIS)

Yurov, A.V.; Yurova, A.A.

2006-01-01

The simple algebraic method for construction of exact solutions of two-dimensional hydrodynamic equations of incompressible flow is proposed. This method can be applied both to nonviscous flow (Euler equations) and to viscous flow (Navier-Stokes equations). In the case of nonviscous flow, the problem is reduced to sequential solving of three linear partial differential equations. In the case of viscous flow, the Navier-Stokes equations are reduced to three linear partial differential equations and one differential equation of the first order [ru

Solution of the multigroup diffusion equation for two-dimensional triangular regions by finite Fourier transformation

International Nuclear Information System (INIS)

Takeshi, Y.; Keisuke, K.

1983-01-01

The multigroup neutron diffusion equation for two-dimensional triangular geometry is solved by the finite Fourier transformation method. Using the zero-th-order equation of the integral equation derived by this method, simple algebraic expressions for the flux are derived and solved by the alternating direction implicit method. In sample calculations for a benchmark problem of a fast breeder reactor, it is shown that the present method gives good results with fewer mesh points than the usual finite difference method

Team behaviour analysis in sports using the poisson equation

OpenAIRE

Direkoglu, Cem; O'Connor, Noel E.

2012-01-01

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

Two-Dimensional Space-Time Dependent Multi-group Diffusion Equation with SLOR Method

International Nuclear Information System (INIS)

Yulianti, Y.; Su'ud, Z.; Waris, A.; Khotimah, S. N.

2010-01-01

The research of two-dimensional space-time diffusion equations with SLOR (Successive-Line Over Relaxation) has been done. SLOR method is chosen because this method is one of iterative methods that does not required to defined whole element matrix. The research is divided in two cases, homogeneous case and heterogeneous case. Homogeneous case has been inserted by step reactivity. Heterogeneous case has been inserted by step reactivity and ramp reactivity. In general, the results of simulations are agreement, even in some points there are differences.

Coupling Navier-stokes and Cahn-hilliard Equations in a Two-dimensional Annular flow Configuration

KAUST Repository

Vignal, Philippe

2015-06-01

In this work, we present a novel isogeometric analysis discretization for the Navier-Stokes- Cahn-Hilliard equation, which uses divergence-conforming spaces. Basis functions generated with this method can have higher-order continuity, and allow to directly discretize the higher- order operators present in the equation. The discretization is implemented in PetIGA-MF, a high-performance framework for discrete differential forms. We present solutions in a two- dimensional annulus, and model spinodal decomposition under shear flow.

Solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

International Nuclear Information System (INIS)

Rosenfeld, M.; Kwak, D.; Vinokur, M.

1988-01-01

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method. 23 references

Trefftz method in solving Fourier-Kirchhoff equation for two-phase flow boiling in a vertical rectangular minichannel

Directory of Open Access Journals (Sweden)

Hożejowska Sylwia

2017-01-01

Full Text Available This paper presents the results of investigations into flow boiling heat transfer in an asymmetrically heated vertical minichannel of 1.7 mm depth. The heated element for FC-72 flowing in the minichannel was an alloy plate 0.45 mm thick, microstructured on one side, in direct contact with the flowing fluid. The computational part of the study contains approximate steady state solutions of the heat transfer problems described by Poisson.s equation and the energy equation for the heated plate and the fluid, respectively. For both equations, the boundary conditions were specified on the basis of experimental data. Temperature of the outer plate surface, measured by infrared thermography, and heat losses to ambient air were included in the calculations. For the energy equation we assumed parabolic profile of fluid velocity and the equality of temperatures and heat fluxes at the interface between the heated surface and the fluid. The void fraction was taken from a single-phase flow model. Two-dimensional temperature distributions were obtained by the Trefftz method and, due to the Robin condition at the interface between them, it was possible to calculate the heat transfer coefficient. Its values were compared to those obtained by other correlations known from literature.

Solution of two-dimensional neutron diffusion equation for triangular region by finite Fourier transformation

International Nuclear Information System (INIS)

Kobayashi, Keisuke; Ishibashi, Hideo

1978-01-01

A two-dimensional neutron diffusion equation for a triangular region is shown to be solved by the finite Fourier transformation. An application of the Fourier transformation to the diffusion equation for triangular region yields equations whose unknowns are the expansion coefficients of the neutron flux and current in Fourier series or Legendre polynomials expansions only at the region boundary. Some numerical calculations have revealed that the present technique gives accurate results. It is shown also that the solution using the expansion in Legendre polynomials converges with relatively few terms even if the solution in Fourier series exhibits the Gibbs' phenomenon. (auth.)

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

International Nuclear Information System (INIS)

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

1982-01-01

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

Iterative solution of the semiconductor device equations

Energy Technology Data Exchange (ETDEWEB)

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

1996-12-31

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

2D sigma models and differential Poisson algebras

International Nuclear Information System (INIS)

Arias, Cesar; Boulanger, Nicolas; Sundell, Per; Torres-Gomez, Alexander

2015-01-01

We construct a two-dimensional topological sigma model whose target space is endowed with a Poisson algebra for differential forms. The model consists of an equal number of bosonic and fermionic fields of worldsheet form degrees zero and one. The action is built using exterior products and derivatives, without any reference to a worldsheet metric, and is of the covariant Hamiltonian form. The equations of motion define a universally Cartan integrable system. In addition to gauge symmetries, the model has one rigid nilpotent supersymmetry corresponding to the target space de Rham operator. The rigid and local symmetries of the action, respectively, are equivalent to the Poisson bracket being compatible with the de Rham operator and obeying graded Jacobi identities. We propose that perturbative quantization of the model yields a covariantized differential star product algebra of Kontsevich type. We comment on the resemblance to the topological A model.

Decomposition of almost-Poisson structure of generalised Chaplygin's nonholonomic systems

International Nuclear Information System (INIS)

Chang, Liu; Peng, Chang; Shi-Xing, Liu; Yong-Xin, Guo

2010-01-01

This paper constructs an almost-Poisson structure for the non-self-adjoint dynamical systems, which can be decomposed into a sum of a Poisson bracket and the other almost-Poisson bracket. The necessary and sufficient condition for the decomposition of the almost-Poisson bracket to be two Poisson ones is obtained. As an application, the almost-Poisson structure for generalised Chaplygin's systems is discussed in the framework of the decomposition theory. It proves that the almost-Poisson bracket for the systems can be decomposed into the sum of a canonical Poisson bracket and another two noncanonical Poisson brackets in some special cases, which is useful for integrating the equations of motion

A Local Approximation of Fundamental Measure Theory Incorporated into Three Dimensional Poisson-Nernst-Planck Equations to Account for Hard Sphere Repulsion Among Ions

Science.gov (United States)

Qiao, Yu; Liu, Xuejiao; Chen, Minxin; Lu, Benzhuo

2016-04-01

The hard sphere repulsion among ions can be considered in the Poisson-Nernst-Planck (PNP) equations by combining the fundamental measure theory (FMT). To reduce the nonlocal computational complexity in 3D simulation of biological systems, a local approximation of FMT is derived, which forms a local hard sphere PNP (LHSPNP) model. In the derivation, the excess chemical potential from hard sphere repulsion is obtained with the FMT and has six integration components. For the integrands and weighted densities in each component, Taylor expansions are performed and the lowest order approximations are taken, which result in the final local hard sphere (LHS) excess chemical potential with four components. By plugging the LHS excess chemical potential into the ionic flux expression in the Nernst-Planck equation, the three dimensional LHSPNP is obtained. It is interestingly found that the essential part of free energy term of the previous size modified model (Borukhov et al. in Phys Rev Lett 79:435-438, 1997; Kilic et al. in Phys Rev E 75:021502, 2007; Lu and Zhou in Biophys J 100:2475-2485, 2011; Liu and Eisenberg in J Chem Phys 141:22D532, 2014) has a very similar form to one term of the LHS model, but LHSPNP has more additional terms accounting for size effects. Equation of state for one component homogeneous fluid is studied for the local hard sphere approximation of FMT and is proved to be exact for the first two virial coefficients, while the previous size modified model only presents the first virial coefficient accurately. To investigate the effects of LHS model and the competitions among different counterion species, numerical experiments are performed for the traditional PNP model, the LHSPNP model, the previous size modified PNP (SMPNP) model and the Monte Carlo simulation. It's observed that in steady state the LHSPNP results are quite different from the PNP results, but are close to the SMPNP results under a wide range of boundary conditions. Besides, in both

Characterization and global analysis of a family of Poisson structures

International Nuclear Information System (INIS)

Hernandez-Bermejo, Benito

2006-01-01

A three-dimensional family of solutions of the Jacobi equations for Poisson systems is characterized. In spite of its general form it is possible the explicit and global determination of its main features, such as the symplectic structure and the construction of the Darboux canonical form. Examples are given

Characterization and global analysis of a family of Poisson structures

Energy Technology Data Exchange (ETDEWEB)

Hernandez-Bermejo, Benito [Escuela Superior de Ciencias Experimentales y Tecnologia, Edificio Departamental II, Universidad Rey Juan Carlos, Calle Tulipan S/N, 28933 (Mostoles), Madrid (Spain)]. E-mail: benito.hernandez@urjc.es

2006-06-26

A three-dimensional family of solutions of the Jacobi equations for Poisson systems is characterized. In spite of its general form it is possible the explicit and global determination of its main features, such as the symplectic structure and the construction of the Darboux canonical form. Examples are given.

Investigation of the use of Prandtl/Navier--Stokes equation procedures for two-dimensional incompressible flows

International Nuclear Information System (INIS)

Anderson, C.R.; Reider, M.B.

1994-01-01

The technique of combining solutions of the Prandtl equations with solutions of the Navier--Stokes equations to compute incompressible flow around two-dimensional bodies is investigated herein. Computational evidence is presented which shows that if the ''obvious'' coupling is used to combine the solutions, then the resulting solution is not accurate. An alternate coupling procedure is described which greatly improves the accuracy of the solutions obtained with the combined equation approach. An alternate coupling that can be used to create a more accurate vortex sheet/vortex blob method is then shown

First error bounds for the porous media approximation of the Poisson-Nernst-Planck equations

Energy Technology Data Exchange (ETDEWEB)

Schmuck, Markus [Imperial College, London (United Kingdom). Depts. of Chemical Engineering and Mathematics

2012-04-15

We study the well-accepted Poisson-Nernst-Planck equations modeling transport of charged particles. By formal multiscale expansions we rederive the porous media formulation obtained by two-scale convergence in [42, 43]. The main result is the derivation of the error which occurs after replacing a highly heterogeneous solid-electrolyte composite by a homogeneous one. The derived estimates show that the approximation errors for both, the ion densities quantified in L{sup 2}-norm and the electric potential measured in H{sup 1}-norm, are of order O(s{sup 1/2}). (orig.)

Solving the two-dimensional stationary transport equation with the aid of the nodal method

International Nuclear Information System (INIS)

Mesina, M.

1976-07-01

In this document the two-dimensional stationary transport equation for the geometry of a fuel assembly or for a system of square boxes has been formulated as an algebraic eigenvalue problem, and the solution was achieved with the computer code NODE 2 which was developed for this purpose. (orig.) [de

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

Science.gov (United States)

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

2018-03-13

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

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

Science.gov (United States)

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

2010-11-01

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

An integrable (2+1)-dimensional Toda equation with two discrete variables

International Nuclear Information System (INIS)

Cao Cewen; Cao Jianli

2007-01-01

An integrable (2+1)-dimensional Toda equation with two discrete variables is presented from the compatible condition of a Lax triad composed of the ZS-AKNS (Zakharov, Shabat; Ablowitz, Kaup, Newell, Segur) eigenvalue problem and two discrete spectral problems. Through the nonlinearization technique, the Lax triad is transformed into a Hamiltonian system and two symplectic maps, respectively, which are integrable in the Liouville sense, sharing the same set of integrals, functionally independent and involutive with each other. In the Jacobi variety of the associated algebraic curve, both the continuous and the discrete flows are straightened out by the Abel-Jacobi coordinates, and are integrated by quadratures. An explicit algebraic-geometric solution in the original variable is obtained by the Riemann-Jacobi inversion

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

KAUST Repository

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

2017-01-01

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

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

DEFF Research Database (Denmark)

Johannesson, Björn

2010-01-01

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

Optimal Control Strategies in a Two Dimensional Differential Game Using Linear Equation under a Perturbed System

Directory of Open Access Journals (Sweden)

Musa Danjuma SHEHU

2008-06-01

Full Text Available This paper lays emphasis on formulation of two dimensional differential games via optimal control theory and consideration of control systems whose dynamics is described by a system of Ordinary Differential equation in the form of linear equation under the influence of two controls U(. and V(.. Base on this, strategies were constructed. Hence we determine the optimal strategy for a control say U(. under a perturbation generated by the second control V(. within a given manifold M.

Large-Time Behavior of Solutions to Vlasov-Poisson-Fokker-Planck Equations: From Evanescent Collisions to Diffusive Limit

Science.gov (United States)

Herda, Maxime; Rodrigues, L. Miguel

2018-03-01

The present contribution investigates the dynamics generated by the two-dimensional Vlasov-Poisson-Fokker-Planck equation for charged particles in a steady inhomogeneous background of opposite charges. We provide global in time estimates that are uniform with respect to initial data taken in a bounded set of a weighted L^2 space, and where dependencies on the mean-free path τ and the Debye length δ are made explicit. In our analysis the mean free path covers the full range of possible values: from the regime of evanescent collisions τ → ∞ to the strongly collisional regime τ → 0. As a counterpart, the largeness of the Debye length, that enforces a weakly nonlinear regime, is used to close our nonlinear estimates. Accordingly we pay a special attention to relax as much as possible the τ -dependent constraint on δ ensuring exponential decay with explicit τ -dependent rates towards the stationary solution. In the strongly collisional limit τ → 0, we also examine all possible asymptotic regimes selected by a choice of observation time scale. Here also, our emphasis is on strong convergence, uniformity with respect to time and to initial data in bounded sets of a L^2 space. Our proofs rely on a detailed study of the nonlinear elliptic equation defining stationary solutions and a careful tracking and optimization of parameter dependencies of hypocoercive/hypoelliptic estimates.

Finite element method with quadratic quadrilateral unit for solving two dimensional incompressible N-S equation

International Nuclear Information System (INIS)

Tao Ganqiang; Yu Qing; Xiao Xiao

2011-01-01

Viscous and incompressible fluid flow is important for numerous engineering mechanics problems. Because of high non linear and incompressibility for Navier-Stokes equation, it is very difficult to solve Navier-Stokes equation by numerical method. According to its characters of Navier-Stokes equation, quartic derivation controlling equation of the two dimensional incompressible Navier-Stokes equation is set up firstly. The method solves the problem for dealing with vorticity boundary and automatically meets incompressibility condition. Then Finite Element equation for Navier-Stokes equation is proposed by using quadratic quadrilateral unit with 8 nodes in which the unit function is quadratic and non linear.-Based on it, the Finite Element program of quadratic quadrilateral unit with 8 nodes is developed. Lastly, numerical experiment proves the accuracy and dependability of the method and also shows the method has good application prospect in computational fluid mechanics. (authors)

A GPU accelerated and error-controlled solver for the unbounded Poisson equation in three dimensions

Science.gov (United States)

Exl, Lukas

2017-12-01

An efficient solver for the three dimensional free-space Poisson equation is presented. The underlying numerical method is based on finite Fourier series approximation. While the error of all involved approximations can be fully controlled, the overall computation error is driven by the convergence of the finite Fourier series of the density. For smooth and fast-decaying densities the proposed method will be spectrally accurate. The method scales with O(N log N) operations, where N is the total number of discretization points in the Cartesian grid. The majority of the computational costs come from fast Fourier transforms (FFT), which makes it ideal for GPU computation. Several numerical computations on CPU and GPU validate the method and show efficiency and convergence behavior. Tests are performed using the Vienna Scientific Cluster 3 (VSC3). A free MATLAB implementation for CPU and GPU is provided to the interested community.

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

Generalization of the separation of variables in the Jacobi identities for finite-dimensional Poisson systems

International Nuclear Information System (INIS)

Hernandez-Bermejo, Benito

2011-01-01

A new n-dimensional family of Poisson structures is globally characterized and analyzed, including the construction of its main features: the symplectic structure and the reduction to the Darboux canonical form. Examples are given that include the generalization of previously known solution families such as the separable Poisson structures. - Highlights: → A new family of Poisson structures is globally characterized and analyzed. → Such family is globally defined for arbitrary values of the dimension and the rank. → Global construction of Casimir invariants and Darboux canonical form is provided. → Very diverse and previously known solutions of physical interest are generalized.

Generalization of the separation of variables in the Jacobi identities for finite-dimensional Poisson systems

Energy Technology Data Exchange (ETDEWEB)

Hernandez-Bermejo, Benito, E-mail: benito.hernandez@urjc.e [Departamento de Fisica, Escuela Superior de Ciencias Experimentales y Tecnologia, Universidad Rey Juan Carlos, Calle Tulipan S/N, 28933 Mostoles, Madrid (Spain)

2011-05-09

A new n-dimensional family of Poisson structures is globally characterized and analyzed, including the construction of its main features: the symplectic structure and the reduction to the Darboux canonical form. Examples are given that include the generalization of previously known solution families such as the separable Poisson structures. - Highlights: A new family of Poisson structures is globally characterized and analyzed. Such family is globally defined for arbitrary values of the dimension and the rank. Global construction of Casimir invariants and Darboux canonical form is provided. Very diverse and previously known solutions of physical interest are generalized.

Localization of Point Sources for Poisson Equation using State Observers

KAUST Repository

Majeed, Muhammad Usman

2016-08-09

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

Localization of Point Sources for Poisson Equation using State Observers

KAUST Repository

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

2016-01-01

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

A Fibonacci collocation method for solving a class of Fredholm–Volterra integral equations in two-dimensional spaces

Directory of Open Access Journals (Sweden)

Farshid Mirzaee

2014-06-01

Full Text Available In this paper, we present a numerical method for solving two-dimensional Fredholm–Volterra integral equations (F-VIE. The method reduces the solution of these integral equations to the solution of a linear system of algebraic equations. The existence and uniqueness of the solution and error analysis of proposed method are discussed. The method is computationally very simple and attractive. Finally, numerical examples illustrate the efficiency and accuracy of the method.

A Fokker-Planck-Landau collision equation solver on two-dimensional velocity grid and its application to particle-in-cell simulation

Energy Technology Data Exchange (ETDEWEB)

Yoon, E. S.; Chang, C. S., E-mail: cschang@pppl.gov [Princeton Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08543 (United States); Korea Advanced Institute of Science and Technology, Yuseong-gu, DaeJeon 305-701 (Korea, Republic of)

2014-03-15

An approximate two-dimensional solver of the nonlinear Fokker-Planck-Landau collision operator has been developed using the assumption that the particle probability distribution function is independent of gyroangle in the limit of strong magnetic field. The isotropic one-dimensional scheme developed for nonlinear Fokker-Planck-Landau equation by Buet and Cordier [J. Comput. Phys. 179, 43 (2002)] and for linear Fokker-Planck-Landau equation by Chang and Cooper [J. Comput. Phys. 6, 1 (1970)] have been modified and extended to two-dimensional nonlinear equation. In addition, a method is suggested to apply the new velocity-grid based collision solver to Lagrangian particle-in-cell simulation by adjusting the weights of marker particles and is applied to a five dimensional particle-in-cell code to calculate the neoclassical ion thermal conductivity in a tokamak plasma. Error verifications show practical aspects of the present scheme for both grid-based and particle-based kinetic codes.

Reduction of infinite dimensional equations

Directory of Open Access Journals (Sweden)

Zhongding Li

2006-02-01

Full Text Available In this paper, we use the general Legendre transformation to show the infinite dimensional integrable equations can be reduced to a finite dimensional integrable Hamiltonian system on an invariant set under the flow of the integrable equations. Then we obtain the periodic or quasi-periodic solution of the equation. This generalizes the results of Lax and Novikov regarding the periodic or quasi-periodic solution of the KdV equation to the general case of isospectral Hamiltonian integrable equation. And finally, we discuss the AKNS hierarchy as a special example.

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

Directory of Open Access Journals (Sweden)

Yibin Zhang

2013-01-01

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

RTk/SN Solutions of the Two-Dimensional Multigroup Transport Equations in Hexagonal Geometry

International Nuclear Information System (INIS)

Valle, Edmundo del; Mund, Ernest H.

2004-01-01

This paper describes an extension to the hexagonal geometry of some weakly discontinuous nodal finite element schemes developed by Hennart and del Valle for the two-dimensional discrete ordinates transport equation in quadrangular geometry. The extension is carried out in a way similar to the extension to the hexagonal geometry of nodal element schemes for the diffusion equation using a composite mapping technique suggested by Hennart, Mund, and del Valle. The combination of the weakly discontinuous nodal transport scheme and the composite mapping is new and is detailed in the main section of the paper. The algorithm efficiency is shown numerically through some benchmark calculations on classical problems widely referred to in the literature

Homogeneous Poisson structures

International Nuclear Information System (INIS)

Shafei Deh Abad, A.; Malek, F.

1993-09-01

We provide an algebraic definition for Schouten product and give a decomposition for any homogenenous Poisson structure in any n-dimensional vector space. A large class of n-homogeneous Poisson structures in R k is also characterized. (author). 4 refs

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

International Nuclear Information System (INIS)

Parra, Felix I.; Catto, Peter J.

2009-01-01

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

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

NARCIS (Netherlands)

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

1985-01-01

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

Poisson-Boltzmann versus Size-Modified Poisson-Boltzmann Electrostatics Applied to Lipid Bilayers.

Science.gov (United States)

Wang, Nuo; Zhou, Shenggao; Kekenes-Huskey, Peter M; Li, Bo; McCammon, J Andrew

2014-12-26

Mean-field methods, such as the Poisson-Boltzmann equation (PBE), are often used to calculate the electrostatic properties of molecular systems. In the past two decades, an enhancement of the PBE, the size-modified Poisson-Boltzmann equation (SMPBE), has been reported. Here, the PBE and the SMPBE are reevaluated for realistic molecular systems, namely, lipid bilayers, under eight different sets of input parameters. The SMPBE appears to reproduce the molecular dynamics simulation results better than the PBE only under specific parameter sets, but in general, it performs no better than the Stern layer correction of the PBE. These results emphasize the need for careful discussions of the accuracy of mean-field calculations on realistic systems with respect to the choice of parameters and call for reconsideration of the cost-efficiency and the significance of the current SMPBE formulation.

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

Directory of Open Access Journals (Sweden)

Tsugio Fukuchi

2014-06-01

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

Quaternionic and Poisson-Lie structures in three-dimensional gravity: The cosmological constant as deformation parameter

International Nuclear Information System (INIS)

Meusburger, C.; Schroers, B. J.

2008-01-01

Each of the local isometry groups arising in three-dimensional (3d) gravity can be viewed as a group of unit (split) quaternions over a ring which depends on the cosmological constant. In this paper we explain and prove this statement and use it as a unifying framework for studying Poisson structures associated with the local isometry groups. We show that, in all cases except for the case of Euclidean signature with positive cosmological constant, the local isometry groups are equipped with the Poisson-Lie structure of a classical double. We calculate the dressing action of the factor groups on each other and find, among others, a simple and unified description of the symplectic leaves of SU(2) and SL(2,R). We also compute the Poisson structure on the dual Poisson-Lie groups of the local isometry groups and on their Heisenberg doubles; together, they determine the Poisson structure of the phase space of 3d gravity in the so-called combinatorial description

An inverse source problem of the Poisson equation with Cauchy data

Directory of Open Access Journals (Sweden)

Ji-Chuan Liu

2017-05-01

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

Two dimensional heat transfer problem in flow boiling in a rectangular minichannel

Directory of Open Access Journals (Sweden)

Hożejowska Sylwia

2015-01-01

Full Text Available The paper presents mathematical modelling of flow boiling heat transfer in a rectangular minichannel asymmetrically heated by a thin and one-sided enhanced foil. Both surfaces are available for observations due to the openings covered with glass sheets. Thus, changes in the colour of the plain foil surface can be registered and then processed. Plain side of the heating foil is covered with a base coat and liquid crystal paint. Observation of the opposite, enhanced surface of the minichannel allows for identification of the gas-liquid two-phase flow patterns and vapour quality. A two-dimensional mathematical model of heat transfer in three subsequent layers (sheet glass, heating foil, liquid was proposed. Heat transfer in all these layers was described with the respective equations: Laplace equation, Poisson equation and energy equation, subject to boundary conditions corresponding to the observed physical process. The solutions (temperature distributions in all three layers were obtained by Trefftz method. Additionally, the temperature of the boiling liquid was obtained by homotopy perturbation method (HPM combined with Trefftz method. The heat transfer coefficient, derived from Robin boundary condition, was estimated in both approaches. In comparison, the results by both methods show very good agreement especially when restricted to the thermal sublayer.

Parametric study of nonlinear electrostatic waves in two-dimensional quantum dusty plasmas

International Nuclear Information System (INIS)

Ali, S; Moslem, W M; Kourakis, I; Shukla, P K

2008-01-01

The nonlinear properties of two-dimensional cylindrical quantum dust-ion-acoustic (QDIA) and quantum dust-acoustic (QDA) waves are studied in a collisionless, unmagnetized and dense (quantum) dusty plasma. For this purpose, the reductive perturbation technique is employed to the quantum hydrodynamical equations and the Poisson equation, obtaining the cylindrical Kadomtsev-Petviashvili (CKP) equations. The effects of quantum diffraction, as well as quantum statistical and geometric effects on the profiles of QDIA and QDA solitary waves are examined. It is found that the amplitudes and widths of the nonplanar QDIA and QDA waves are significantly affected by the quantum electron tunneling effect. The addition of a dust component to a quantum plasma is seen to affect the propagation characteristics of localized QDIA excitations. In the case of low-frequency QDA waves, this effect is even stronger, since the actual form of the potential solitary waves, in fact, depends on the dust charge polarity (positive/negative) itself (allowing for positive/negative potential forms, respectively). The relevance of the present investigation to metallic nanostructures is highlighted

Stability of plane wave solutions of the two-space-dimensional nonlinear Schroedinger equation

International Nuclear Information System (INIS)

Martin, D.U.; Yuen, H.C.; Saffman, P.G.

1980-01-01

The stability of plane, periodic solutions of the two-dimensional nonlinear Schroedinger equation to infinitesimal, two-dimensional perturbation has been calculated and verified numerically. For standing wave disturbances, instability is found for both odd and even modes; as the period of the unperturbed solution increases, the instability associated with the odd modes remains but that associated with the even mode disappears, which is consistent with the results of Zakharov and Rubenchik, Saffman and Yuen and Ablowitz and Segur on the stability of solitons. In addition, we have identified travelling wave instabilities for the even mode perturbations which are absent in the long-wave limit. Extrapolation to the case of an unperturbed solution with infinite period suggests that these instabilities may also be present for the soliton. In other words, the soliton is unstable to odd, standing-wave perturbations, and very likely also to even, travelling-wave perturbations. (orig.)

Inverse Problem for Two-Dimensional Discrete Schr`dinger Equation

CERN Document Server

Serdyukova, S I

2000-01-01

For two-dimensional discrete Schroedinger equation the boundary-value problem in rectangle M times N with zero boundary conditions is solved. It's stated in this work, that inverse problem reduces to reconstruction of C symmetric five-diagonal matrix with given spectrum and given first k(M,N), 1<-k

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

Science.gov (United States)

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

2007-02-01

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

Birth–death process of local structures in defect turbulence described by the one-dimensional complex Ginzburg–Landau equation

Energy Technology Data Exchange (ETDEWEB)

Uchiyama, Yusuke, E-mail: r1230160@risk.tsukuba.ac.jp; Konno, Hidetoshi

2014-04-01

Defect turbulence described by the one-dimensional complex Ginzburg–Landau equation is investigated and analyzed via a birth–death process of the local structures composed of defects, holes, and modulated amplitude waves (MAWs). All the number statistics of each local structure, in its stationary state, are subjected to Poisson statistics. In addition, the probability density functions of interarrival times of defects, lifetimes of holes, and MAWs show the existence of long-memory and some characteristic time scales caused by zigzag motions of oscillating traveling holes. The corresponding stochastic process for these observations is fully described by a non-Markovian master equation.

TWO-DIMENSIONAL STELLAR EVOLUTION CODE INCLUDING ARBITRARY MAGNETIC FIELDS. II. PRECISION IMPROVEMENT AND INCLUSION OF TURBULENCE AND ROTATION

International Nuclear Information System (INIS)

Li Linghuai; Sofia, Sabatino; Basu, Sarbani; Demarque, Pierre; Ventura, Paolo; Penza, Valentina; Bi Shaolan

2009-01-01

In the second paper of this series we pursue two objectives. First, in order to make the code more sensitive to small effects, we remove many approximations made in Paper I. Second, we include turbulence and rotation in the two-dimensional framework. The stellar equilibrium is described by means of a set of five differential equations, with the introduction of a new dependent variable, namely the perturbation to the radial gravity, that is found when the nonradial effects are considered in the solution of the Poisson equation. Following the scheme of the first paper, we write the equations in such a way that the two-dimensional effects can be easily disentangled. The key concept introduced in this series is the equipotential surface. We use the underlying cause-effect relation to develop a recurrence relation to calculate the equipotential surface functions for uniform rotation, differential rotation, rotation-like toroidal magnetic fields, and turbulence. We also develop a more precise code to numerically solve the two-dimensional stellar structure and evolution equations based on the equipotential surface calculations. We have shown that with this formulation we can achieve the precision required by observations by appropriately selecting the convergence criterion. Several examples are presented to show that the method works well. Since we are interested in modeling the effects of a dynamo-type field on the detailed envelope structure and global properties of the Sun, the code has been optimized for short timescales phenomena (down to 1 yr). The time dependence of the code has so far been tested exclusively to address such problems.

dimensional Jaulent–Miodek equations

Indian Academy of Sciences (India)

(2+1)-dimensional Jaulent–Miodek equation; the first integral method; kinks; ... and effective method for solving nonlinear partial differential equations which can ... of the method employed and exact kink and soliton solutions are constructed ...

One-dimensional reduction of the three-dimenstional Gross-Pitaevskii equation with two- and three-body interactions

International Nuclear Information System (INIS)

Cardoso, W. B.; Avelar, A. T.; Bazeia, D.

2011-01-01

We deal with the three-dimensional Gross-Pitaevskii equation which is used to describe a cloud of dilute bosonic atoms that interact under competing two- and three-body scattering potentials. We study the case where the cloud of atoms is strongly confined in two spatial dimensions, allowing us to build an unidimensional nonlinear equation,controlled by the nonlinearities and the confining potentials that trap the system along the longitudinal coordinate. We focus attention on specific limits dictated by the cubic and quintic coefficients, and we implement numerical simulations to help us to quantify the validity of the procedure.

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

Science.gov (United States)

Prentice, J. S. C.

2012-01-01

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

Decoupling of the nernst-planck and poisson equations. Application to a membrane system at overlimiting currents.

Science.gov (United States)

Urtenov, Mahamet A-Kh; Kirillova, Evgeniya V; Seidova, Natalia M; Nikonenko, Victor V

2007-12-27

This paper deals with one-dimensional stationary Nernst-Planck and Poisson (NPP) equations describing ion electrodiffusion in multicomponent solution/electrode or ion-conductive membrane systems. A general method for resolving ordinary and singularly perturbed problems with these equations is developed. This method is based on the decoupling of NPP equations that results in deduction of an equation containing only the terms with different powers of the electrical field and its derivatives. Then, the solution of this equation, analytical in several cases or numerical, is substituted into the Nernst-Planck equations for calculating the concentration profile for each ion present in the system. Different ionic species are grouped in valency classes that allows one to reduce the dimension of the original set of equations and leads to a relatively easy treatment of multi-ion systems. When applying the method developed, the main attention is paid to ion transfer at limiting and overlimiting currents, where a significant deviation from local electroneutrality occurs. The boundary conditions and different approximations are examined: the local electroneutrality (LEN) assumption and the original assumption of quasi-uniform distribution of the space charge density (QCD). The relations between the ion fluxes at limiting and overlimiting currents are discussed. In particular, attention is paid to the "exaltation" of counterion transfer toward an ion-exchange membrane by co-ion flux leaking through the membrane or generated at the membrane/solution interface. The structure of the multi-ion concentration field in a depleted diffusion boundary layer (DBL) near an ion-exchange membrane at overlimiting currents is analyzed. The presence of salt ions and hydrogen and hydroxyl ions generated in the course of the water "splitting" reaction is considered. The thickness of the DBL and its different zones, as functions of applied current density, are found by fitting experimental current

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

Science.gov (United States)

Kopp, Michael; Vattis, Kyriakos; Skordis, Constantinos

2017-12-01

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

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

DEFF Research Database (Denmark)

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

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

Solving the Helmholtz equation in conformal mapped ARROWstructures using homotopy perturbation method

DEFF Research Database (Denmark)

Reck, Kasper; Thomsen, Erik Vilain; Hansen, Ole

2011-01-01

. The solution of the mapped Helmholtz equation is found by solving an infinite series of Poisson equations using two dimensional Fourier series. The solution is entirely based on analytical expressions and is not mesh dependent. The analytical results are compared to a numerical (finite element method) solution......The scalar wave equation, or Helmholtz equation, describes within a certain approximation the electromagnetic field distribution in a given system. In this paper we show how to solve the Helmholtz equation in complex geometries using conformal mapping and the homotopy perturbation method...

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

Directory of Open Access Journals (Sweden)

J. Izadian

2014-01-01

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

On a Poisson homogeneous space of bilinear forms with a Poisson-Lie action

Science.gov (United States)

Chekhov, L. O.; Mazzocco, M.

2017-12-01

Let \\mathscr A be the space of bilinear forms on C^N with defining matrices A endowed with a quadratic Poisson structure of reflection equation type. The paper begins with a short description of previous studies of the structure, and then this structure is extended to systems of bilinear forms whose dynamics is governed by the natural action A\\mapsto B ABT} of the {GL}_N Poisson-Lie group on \\mathscr A. A classification is given of all possible quadratic brackets on (B, A)\\in {GL}_N× \\mathscr A preserving the Poisson property of the action, thus endowing \\mathscr A with the structure of a Poisson homogeneous space. Besides the product Poisson structure on {GL}_N× \\mathscr A, there are two other (mutually dual) structures, which (unlike the product Poisson structure) admit reductions by the Dirac procedure to a space of bilinear forms with block upper triangular defining matrices. Further generalisations of this construction are considered, to triples (B,C, A)\\in {GL}_N× {GL}_N× \\mathscr A with the Poisson action A\\mapsto B ACT}, and it is shown that \\mathscr A then acquires the structure of a Poisson symmetric space. Generalisations to chains of transformations and to the quantum and quantum affine algebras are investigated, as well as the relations between constructions of Poisson symmetric spaces and the Poisson groupoid. Bibliography: 30 titles.

Computing stationary solutions of the two-dimensional Gross-Pitaevskii equation with deflated continuation

Science.gov (United States)

Charalampidis, E. G.; Kevrekidis, P. G.; Farrell, P. E.

2018-01-01

In this work we employ a recently proposed bifurcation analysis technique, the deflated continuation algorithm, to compute steady-state solitary waveforms in a one-component, two-dimensional nonlinear Schrödinger equation with a parabolic trap and repulsive interactions. Despite the fact that this system has been studied extensively, we discover a wide variety of previously unknown branches of solutions. We analyze the stability of the newly discovered branches and discuss the bifurcations that relate them to known solutions both in the near linear (Cartesian, as well as polar) and in the highly nonlinear regimes. While deflated continuation is not guaranteed to compute the full bifurcation diagram, this analysis is a potent demonstration that the algorithm can discover new nonlinear states and provide insights into the energy landscape of complex high-dimensional Hamiltonian dynamical systems.

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

Science.gov (United States)

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

2017-01-01

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

Two-dimensional differential transform method for solving linear and non-linear Schroedinger equations

International Nuclear Information System (INIS)

Ravi Kanth, A.S.V.; Aruna, K.

2009-01-01

In this paper, we propose a reliable algorithm to develop exact and approximate solutions for the linear and nonlinear Schroedinger equations. The approach rest mainly on two-dimensional differential transform method which is one of the approximate methods. The method can easily be applied to many linear and nonlinear problems and is capable of reducing the size of computational work. Exact solutions can also be achieved by the known forms of the series solutions. Several illustrative examples are given to demonstrate the effectiveness of the present method.

Equivalence of two-dimensional gravities

International Nuclear Information System (INIS)

Mohammedi, N.

1990-01-01

The authors find the relationship between the Jackiw-Teitelboim model of two-dimensional gravity and the SL(2,R) induced gravity. These are shown to be related to a two-dimensional gauge theory obtained by dimensionally reducing the Chern-Simons action of the 2 + 1 dimensional gravity. The authors present an explicit solution to the equations of motion of the auxiliary field of the Jackiw-Teitelboim model in the light-cone gauge. A renormalization of the cosmological constant is also given

Equilibrium: two-dimensional configurations

International Nuclear Information System (INIS)

Anon.

1987-01-01

In Chapter 6, the problem of toroidal force balance is addressed in the simplest, nontrivial two-dimensional geometry, that of an axisymmetric torus. A derivation is presented of the Grad-Shafranov equation, the basic equation describing axisymmetric toroidal equilibrium. The solutions to equations provide a complete description of ideal MHD equilibria: radial pressure balance, toroidal force balance, equilibrium Beta limits, rotational transform, shear, magnetic wall, etc. A wide number of configurations are accurately modeled by the Grad-Shafranov equation. Among them are all types of tokamaks, the spheromak, the reversed field pinch, and toroidal multipoles. An important aspect of the analysis is the use of asymptotic expansions, with an inverse aspect ratio serving as the expansion parameter. In addition, an equation similar to the Grad-Shafranov equation, but for helically symmetric equilibria, is presented. This equation represents the leading-order description low-Beta and high-Beta stellarators, heliacs, and the Elmo bumpy torus. The solutions all correspond to infinitely long straight helices. Bending such a configuration into a torus requires a full three-dimensional calculation and is discussed in Chapter 7

Analytic solution of the two-dimensional Fokker-Planck equation governing stochastic ion heating by a lower hybrid wave

International Nuclear Information System (INIS)

Malescio, G.

1981-04-01

The two-dimensional Fokker-Planck equation describing the ion motion in a coherent lower hybrid wave above the stochasticity threshold is analytically solved. An expression is given for the steady state power dissipation

Critical behavior in two-dimensional quantum gravity and equations of motion of the string

International Nuclear Information System (INIS)

Das, S.R.; Dhar, A.; Wadia, S.R.

1990-01-01

The authors show how consistent quantization determines the renormalization of couplings in a quantum field theory coupled to gravity in two dimensions. The special status of couplings corresponding to conformally invariant matter is discussed. In string theory, where the dynamical degree of freedom of the two-dimensional metric plays the role of time in target space, these renormalization group equations are themselves the classical equations of motion. Time independent solutions, like classical vacuua, correspond to the situation in which matter is conformally invariant. Time dependent solutions, like tunnelling configurations between vacuua, correspond to special trajectories in theory space. The authors discuss an example of such a trajectory in the space containing the c ≤ 1 minimal models. The authors also discuss the connection between this work and the recent attempts to construct non-pertubative string theories based on matrix models

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

Directory of Open Access Journals (Sweden)

José Colmenares

2014-01-01

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

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

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

International Nuclear Information System (INIS)

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

2009-01-01

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

Bifurcations of solitary wave solutions for (two and three)-dimensional nonlinear partial differential equation in quantum and magnetized plasma by using two different methods

Science.gov (United States)

Khater, Mostafa M. A.; Seadawy, Aly R.; Lu, Dianchen

2018-06-01

In this research, we study new two techniques that called the extended simple equation method and the novel (G‧/G) -expansion method. The extended simple equation method depend on the auxiliary equation (dϕ/dξ = α + λϕ + μϕ2) which has three ways for solving depends on the specific condition on the parameters as follow: When (λ = 0) this auxiliary equation reduces to Riccati equation, when (α = 0) this auxiliary equation reduces to Bernoulli equation and when (α ≠ 0, λ ≠ 0, μ ≠ 0) we the general solutions of this auxiliary equation while the novel (G‧/G) -expansion method depends also on similar auxiliary equation (G‧/G)‧ = μ + λ(G‧/G) + (v - 1)(G‧/G) 2 which depend also on the value of (λ2 - 4 (v - 1) μ) and the specific condition on the parameters as follow: When (λ = 0) this auxiliary equation reduces to Riccati equation, when (μ = 0) this auxiliary equation reduces to Bernoulli equation and when (λ2 ≠ 4 (v - 1) μ) we the general solutions of this auxiliary equation. This show how both of these auxiliary equation are special cases of Riccati equation. We apply these methods on two dimensional nonlinear Kadomtsev-Petviashvili Burgers equation in quantum plasma and three-dimensional nonlinear modified Zakharov-Kuznetsov equation of ion-acoustic waves in a magnetized plasma. We obtain the exact traveling wave solutions of these important models and under special condition on the parameters, we get solitary traveling wave solutions. All calculations in this study have been established and verified back with the aid of the Maple package program. The executed method is powerful, effective and straightforward for solving nonlinear partial differential equations to obtain more and new solutions.

Baicklund transformation and multiple soliton solutions for the (3+1)-dimensional Jimbo-Miwa equation

Institute of Scientific and Technical Information of China (English)

张解放; 吴锋民

2002-01-01

We study an approach to constructing multiple soliton solutions of the (3+1)-dimensional nonlinear evolution equation. We take the (3+1)-dimensional Jimbo-Miwa (JM) equation as an example. Using the extended homogeneous balance method, one can find a Backlund transformation to decompose the (3+1)-dimensional JM equation into a linear partial differential equation and two bilinear partial differential equations. Starting from these linear and bilinear partial differential equations, some multiple soliton solutions for the (3+1)-dimensional JM equation are obtained by introducing a class of formal solutions.

New method for solving three-dimensional Schroedinger equation

International Nuclear Information System (INIS)

Melezhik, V.S.

1992-01-01

A new method is developed for solving the multidimensional Schroedinger equation without the variable separation. To solve the Schroedinger equation in a multidimensional coordinate space X, a difference grid Ω i (i=1,2,...,N) for some of variables, Ω, from X={R,Ω} is introduced and the initial partial-differential equation is reduced to a system of N differential-difference equations in terms of one of the variables R. The arising multi-channel scattering (or eigenvalue) problem is solved by the algorithm based on a continuous analog of the Newton method. The approach has been successfully tested for several two-dimensional problems (scattering on a nonspherical potential well and 'dipole' scatterer, a hydrogen atom in a homogenous magnetic field) and for a three-dimensional problem of the helium-atom bound states. (author)

A meshless method for solving two-dimensional variable-order time fractional advection-diffusion equation

Science.gov (United States)

Tayebi, A.; Shekari, Y.; Heydari, M. H.

2017-07-01

Several physical phenomena such as transformation of pollutants, energy, particles and many others can be described by the well-known convection-diffusion equation which is a combination of the diffusion and advection equations. In this paper, this equation is generalized with the concept of variable-order fractional derivatives. The generalized equation is called variable-order time fractional advection-diffusion equation (V-OTFA-DE). An accurate and robust meshless method based on the moving least squares (MLS) approximation and the finite difference scheme is proposed for its numerical solution on two-dimensional (2-D) arbitrary domains. In the time domain, the finite difference technique with a θ-weighted scheme and in the space domain, the MLS approximation are employed to obtain appropriate semi-discrete solutions. Since the newly developed method is a meshless approach, it does not require any background mesh structure to obtain semi-discrete solutions of the problem under consideration, and the numerical solutions are constructed entirely based on a set of scattered nodes. The proposed method is validated in solving three different examples including two benchmark problems and an applied problem of pollutant distribution in the atmosphere. In all such cases, the obtained results show that the proposed method is very accurate and robust. Moreover, a remarkable property so-called positive scheme for the proposed method is observed in solving concentration transport phenomena.

Density of states, Poisson's formula of summation and Walfisz's formula

International Nuclear Information System (INIS)

Fucho, P.

1980-06-01

Using Poisson's formula for summation, we obtain an expression for density of states of d-dimensional scalar Helmoholtz's equation under various boundary conditions. Likewise, we also obtain formulas of Walfisz's type. It becomes evident that the formulas obtained by Pathria et al. in connection with ideal bosons in a finite system are exactly the same as those obtained by utilizing the formulas for density of states. (author)

Poisson Coordinates.

Science.gov (United States)

Li, Xian-Ying; Hu, Shi-Min

2013-02-01

Harmonic functions are the critical points of a Dirichlet energy functional, the linear projections of conformal maps. They play an important role in computer graphics, particularly for gradient-domain image processing and shape-preserving geometric computation. We propose Poisson coordinates, a novel transfinite interpolation scheme based on the Poisson integral formula, as a rapid way to estimate a harmonic function on a certain domain with desired boundary values. Poisson coordinates are an extension of the Mean Value coordinates (MVCs) which inherit their linear precision, smoothness, and kernel positivity. We give explicit formulas for Poisson coordinates in both continuous and 2D discrete forms. Superior to MVCs, Poisson coordinates are proved to be pseudoharmonic (i.e., they reproduce harmonic functions on n-dimensional balls). Our experimental results show that Poisson coordinates have lower Dirichlet energies than MVCs on a number of typical 2D domains (particularly convex domains). As well as presenting a formula, our approach provides useful insights for further studies on coordinates-based interpolation and fast estimation of harmonic functions.

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

International Nuclear Information System (INIS)

Pabst, M.

2014-01-01

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

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

Science.gov (United States)

Egan, Raphael; Gibou, Frédéric

2017-10-01

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

On Riemann boundary value problems for null solutions of the two dimensional Helmholtz equation

Science.gov (United States)

Bory Reyes, Juan; Abreu Blaya, Ricardo; Rodríguez Dagnino, Ramón Martin; Kats, Boris Aleksandrovich

2018-01-01

The Riemann boundary value problem (RBVP to shorten notation) in the complex plane, for different classes of functions and curves, is still widely used in mathematical physics and engineering. For instance, in elasticity theory, hydro and aerodynamics, shell theory, quantum mechanics, theory of orthogonal polynomials, and so on. In this paper, we present an appropriate hyperholomorphic approach to the RBVP associated to the two dimensional Helmholtz equation in R^2 . Our analysis is based on a suitable operator calculus.

Approximate solutions for the two-dimensional integral transport equation. Solution of complex two-dimensional transport problems

International Nuclear Information System (INIS)

Sanchez, Richard.

1980-11-01

This work is divided into two parts: the first part deals with the solution of complex two-dimensional transport problems, the second one (note CEA-N-2166) treats the critically mixed methods of resolution. A set of approximate solutions for the isotropic two-dimensional neutron transport problem has been developed using the interface current formalism. The method has been applied to regular lattices of rectangular cells containing a fuel pin, cladding, and water, or homogenized structural material. The cells are divided into zones that are homogeneous. A zone-wise flux expansion is used to formulate a direct collision probability problem within a cell. The coupling of the cells is effected by making extra assumptions on the currents entering and leaving the interfaces. Two codes have been written: CALLIOPE uses a cylindrical cell model and one or three terms for the flux expansion, and NAUSICAA uses a two-dimensional flux representation and does a truly two-dimensional calculation inside each cell. In both codes, one or three terms can be used to make a space-independent expansion of the angular fluxes entering and leaving each side of the cell. The accuracies and computing times achieved with the different approximations are illustrated by numerical studies on two benchmark problems and by calculations performed in the APOLLO multigroup code [fr

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

The Analysis of Corporate Bond Valuation under an Infinite Dimensional Compound Poisson Framework

Directory of Open Access Journals (Sweden)

Sheng Fan

2014-01-01

Full Text Available This paper analyzes the firm bond valuation and credit spread with an endogenous model for the pure default and callable default corporate bond. Regarding the stochastic instantaneous forward rates and the firm value as an infinite dimensional Poisson process, we provide some analytical results for the embedded American options and firm bond valuations.

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

Czech Academy of Sciences Publication Activity Database

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

2009-01-01

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

Analytic study of solutions for a (3 + 1) -dimensional generalized KP equation

Science.gov (United States)

Gao, Hui; Cheng, Wenguang; Xu, Tianzhou; Wang, Gangwei

2018-03-01

The (3 + 1) -dimensional generalized KP (gKP) equation is an important nonlinear partial differential equation in theoretical and mathematical physics which can be used to describe nonlinear wave motion. Through the Hirota bilinear method, one-solition, two-solition and N-solition solutions are derived via symbolic computation. Two classes of lump solutions, rationally localized in all directions in space, to the dimensionally reduced cases in (2 + 1)-dimensions, are constructed by using a direct method based on the Hirota bilinear form of the equation. It implies that we can derive the lump solutions of the reduced gKP equation from positive quadratic function solutions to the aforementioned bilinear equation. Meanwhile, we get interaction solutions between a lump and a kink of the gKP equation. The lump appears from a kink and is swallowed by it with the change of time. This work offers a possibility which can enrich the variety of the dynamical features of solutions for higher-dimensional nonlinear evolution equations.

Darboux transformations for (1+2)-dimensional Fokker-Planck equations with constant diffusion matrix

International Nuclear Information System (INIS)

Schulze-Halberg, Axel

2012-01-01

We construct a Darboux transformation for (1+2)-dimensional Fokker-Planck equations with constant diffusion matrix. Our transformation is based on the two-dimensional supersymmetry formalism for the Schrödinger equation. The transformed Fokker-Planck equation and its solutions are obtained in explicit form.

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

Science.gov (United States)

Xie, Yang; Ying, Jinyong; Xie, Dexuan

2017-03-30

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

Pentadiagonal alternating-direction-implicit finite-difference time-domain method for two-dimensional Schrödinger equation

Science.gov (United States)

Tay, Wei Choon; Tan, Eng Leong

2014-07-01

In this paper, we have proposed a pentadiagonal alternating-direction-implicit (Penta-ADI) finite-difference time-domain (FDTD) method for the two-dimensional Schrödinger equation. Through the separation of complex wave function into real and imaginary parts, a pentadiagonal system of equations for the ADI method is obtained, which results in our Penta-ADI method. The Penta-ADI method is further simplified into pentadiagonal fundamental ADI (Penta-FADI) method, which has matrix-operator-free right-hand-sides (RHS), leading to the simplest and most concise update equations. As the Penta-FADI method involves five stencils in the left-hand-sides (LHS) of the pentadiagonal update equations, special treatments that are required for the implementation of the Dirichlet's boundary conditions will be discussed. Using the Penta-FADI method, a significantly higher efficiency gain can be achieved over the conventional Tri-ADI method, which involves a tridiagonal system of equations.

Ring-shaped quasi-soliton solutions to the two-and three-dimensional Sine-Gordon equation

International Nuclear Information System (INIS)

Christiansen, P.L.; Olsen, O.H.

1979-01-01

Ring-shaped solitary wave solutions to the Sine-Gordon equation in two and three spatial dimensions are investigated by numerical computation. Each expanding wave exhibits a return effect. The reflection of the shrinking wave at the singularity at the center of the wave is investigated in a particular case. Collision experiments in numero for expanding and shrinking concentric ring waves show that the solutions possess quasisoliton properties. A Baecklund transformation for the non-symmetric three-dimensional case is given. (Auth.)

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

Directory of Open Access Journals (Sweden)

Elina L. Shishkina

2017-07-01

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

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

Directory of Open Access Journals (Sweden)

Hua Yang

2012-01-01

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

A meshless local radial basis function method for two-dimensional incompressible Navier-Stokes equations

KAUST Repository

Wang, Zhiheng

2014-12-10

A meshless local radial basis function method is developed for two-dimensional incompressible Navier-Stokes equations. The distributed nodes used to store the variables are obtained by the philosophy of an unstructured mesh, which results in two main advantages of the method. One is that the unstructured nodes generation in the computational domain is quite simple, without much concern about the mesh quality; the other is that the localization of the obtained collocations for the discretization of equations is performed conveniently with the supporting nodes. The algebraic system is solved by a semi-implicit pseudo-time method, in which the convective and source terms are explicitly marched by the Runge-Kutta method, and the diffusive terms are implicitly solved. The proposed method is validated by several benchmark problems, including natural convection in a square cavity, the lid-driven cavity flow, and the natural convection in a square cavity containing a circular cylinder, and very good agreement with the existing results are obtained.

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

Science.gov (United States)

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

2007-01-01

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

Numerical Solution of the Fractional Partial Differential Equations by the Two-Dimensional Fractional-Order Legendre Functions

Directory of Open Access Journals (Sweden)

Fukang Yin

2013-01-01

Full Text Available A numerical method is presented to obtain the approximate solutions of the fractional partial differential equations (FPDEs. The basic idea of this method is to achieve the approximate solutions in a generalized expansion form of two-dimensional fractional-order Legendre functions (2D-FLFs. The operational matrices of integration and derivative for 2D-FLFs are first derived. Then, by these matrices, a system of algebraic equations is obtained from FPDEs. Hence, by solving this system, the unknown 2D-FLFs coefficients can be computed. Three examples are discussed to demonstrate the validity and applicability of the proposed method.

Two-dimensional Haar wavelet Collocation Method for the solution of Stationary Neutron Transport Equation in a homogeneous isotropic medium

International Nuclear Information System (INIS)

Patra, A.; Saha Ray, S.

2014-01-01

Highlights: • A stationary transport equation has been solved using the technique of Haar wavelet Collocation Method. • This paper intends to provide the great utility of Haar wavelets to nuclear science problem. • In the present paper, two-dimensional Haar wavelets are applied. • The proposed method is mathematically very simple, easy and fast. - Abstract: This paper emphasizes on finding the solution for a stationary transport equation using the technique of Haar wavelet Collocation Method (HWCM). Haar wavelet Collocation Method is efficient and powerful in solving wide class of linear and nonlinear differential equations. Recently Haar wavelet transform has gained the reputation of being a very effective tool for many practical applications. This paper intends to provide the great utility of Haar wavelets to nuclear science problem. In the present paper, two-dimensional Haar wavelets are applied for solution of the stationary Neutron Transport Equation in homogeneous isotropic medium. The proposed method is mathematically very simple, easy and fast. To demonstrate about the efficiency of the method, one test problem is discussed. It can be observed from the computational simulation that the numerical approximate solution is much closer to the exact solution

Comparison of preconditioned generalized conjugate gradient methods to two-dimensional neutron and photon transport equation

International Nuclear Information System (INIS)

Chen, G.S.; Yang, D.Y.

1998-01-01

We apply and compare the preconditioned generalized conjugate gradient methods to solve the linear system equation that arises in the two-dimensional neutron and photon transport equation in this paper. Several subroutines are developed on the basis of preconditioned generalized conjugate gradient methods for time-independent, two-dimensional neutron and photon transport equation in the transport theory. These generalized conjugate gradient methods are used: TFQMR (transpose free quasi-minimal residual algorithm) CGS (conjugate gradient square algorithm), Bi-CGSTAB (bi-conjugate gradient stabilized algorithm) and QMRCGSTAB (quasi-minimal residual variant of bi-conjugate gradient stabilized algorithm). These subroutines are connected to computer program DORT. Several problems are tested on a personal computer with Intel Pentium CPU. The reasons to choose the generalized conjugate gradient methods are that the methods have better residual (equivalent to error) control procedures in the computation and have better convergent rate. The pointwise incomplete LU factorization ILU, modified pointwise incomplete LU factorization MILU, block incomplete factorization BILU and modified blockwise incomplete LU factorization MBILU are the preconditioning techniques used in the several testing problems. In Bi-CGSTAB, CGS, TFQMR and QMRCGSTAB method, we find that either CGS or Bi-CGSTAB method combined with preconditioner MBILU is the most efficient algorithm in these methods in the several testing problems. The numerical solution of flux by preconditioned CGS and Bi-CGSTAB methods has the same result as those from Cray computer, obtained by either the point successive relaxation method or the line successive relaxation method combined with Gaussian elimination

A closed-form solution for the two-dimensional transport equation by the LTSN nodal method in the range of Compton Effect

International Nuclear Information System (INIS)

Rodriguez, Barbara D.A.; Tullio de Vilhena, Marco; Hoff, Gabriela

2008-01-01

In this paper we report a two-dimensional LTS N nodal solution for homogeneous and heterogeneous rectangular domains, assuming the Klein-Nishina scattering kernel and multigroup model. The main idea relies on the solution of the two one-dimensional S N equations resulting from transverse integration of the S N equations in the rectangular domain by the LTS N nodal method, considering the leakage angular fluxes approximated by exponential, which allow us to determine a closed-form solution for the photons transport equation. The angular flux and the parameters of the medium are used for the calculation of the absorbed energy in rectangular domains with different dimensions and compositions. The incoming photons will be tracked until their whole energy is deposited and/or they leave the domain of interest. In this study, the absorbed energy by Compton Effect will be considered. The remaining effects will not be taken into account. We present numerical simulations and comparisons with results obtained by using Geant4 (version 9.1) program which applies the Monte Carlo's technique to low energy libraries for a two-dimensional problem assuming the Klein-Nishina scattering kernel. (authors)

Dynamics of one- and two-dimensional fronts in a bistable equation with time-delayed global feedback: Propagation failure and control mechanisms

International Nuclear Information System (INIS)

Boubendir, Yassine; Mendez, Vicenc; Rotstein, Horacio G.

2010-01-01

We study the evolution of fronts in a bistable equation with time-delayed global feedback in the fast reaction and slow diffusion regime. This equation generalizes the Hodgkin-Grafstein and Allen-Cahn equations. We derive a nonlinear equation governing the motion of fronts, which includes a term with delay. In the one-dimensional case this equation is linear. We study the motion of one- and two-dimensional fronts, finding a much richer dynamics than for the previously studied cases (without time-delayed global feedback). We explain the mechanism by which localized fronts created by inhibitory global coupling loose stability in a Hopf bifurcation as the delay time increases. We show that for certain delay times, the prevailing phase is different from that corresponding to the system in the absence of global coupling. Numerical simulations of the partial differential equation are in agreement with the analytical predictions.

Mathematical modeling and the two-phase constitutive equations

International Nuclear Information System (INIS)

Boure, J.A.

1975-01-01

The problems raised by the mathematical modeling of two-phase flows are summarized. The models include several kinds of equations, which cannot be discussed independently, such as the balance equations and the constitutive equations. A review of the various two-phase one-dimensional models proposed to date, and of the constitutive equations they imply, is made. These models are either mixture models or two-fluid models. Due to their potentialities, the two-fluid models are discussed in more detail. To avoid contradictions, the form of the constitutive equations involved in two-fluid models must be sufficiently general. A special form of the two-fluid models, which has particular advantages, is proposed. It involves three mixture balance equations, three balance equations for slip and thermal non-equilibriums, and the necessary constitutive equations [fr

Singular Poisson tensors

International Nuclear Information System (INIS)

Littlejohn, R.G.

1982-01-01

The Hamiltonian structures discovered by Morrison and Greene for various fluid equations were obtained by guessing a Hamiltonian and a suitable Poisson bracket formula, expressed in terms of noncanonical (but physical) coordinates. In general, such a procedure for obtaining a Hamiltonian system does not produce a Hamiltonian phase space in the usual sense (a symplectic manifold), but rather a family of symplectic manifolds. To state the matter in terms of a system with a finite number of degrees of freedom, the family of symplectic manifolds is parametrized by a set of Casimir functions, which are characterized by having vanishing Poisson brackets with all other functions. The number of independent Casimir functions is the corank of the Poisson tensor J/sup ij/, the components of which are the Poisson brackets of the coordinates among themselves. Thus, these Casimir functions exist only when the Poisson tensor is singular

Solving Two -Dimensional Diffusion Equations with Nonlocal Boundary Conditions by a Special Class of Padé Approximants

Directory of Open Access Journals (Sweden)

Mohammad Siddique

2010-08-01

Full Text Available Parabolic partial differential equations with nonlocal boundary conditions arise in modeling of a wide range of important application areas such as chemical diffusion, thermoelasticity, heat conduction process, control theory and medicine science. In this paper, we present the implementation of positivity- preserving Padé numerical schemes to the two-dimensional diffusion equation with nonlocal time dependent boundary condition. We successfully implemented these numerical schemes for both Homogeneous and Inhomogeneous cases. The numerical results show that these Padé approximation based numerical schemes are quite accurate and easily implemented.

Two-dimensional transport of tokamak plasmas

International Nuclear Information System (INIS)

Hirshman, S.P.; Jardin, S.C.

1979-01-01

A reduced set of two-fluid transport equations is obtained from the conservation equations describing the time evolution of the differential particle number, entropy, and magnetic fluxes in an axisymmetric toroidal plasma with nested magnetic surfaces. Expanding in the small ratio of perpendicular to parallel mobilities and thermal conductivities yields as solubility constraints one-dimensional equations for the surface-averaged thermodynamic variables and magnetic fluxes. Since Ohm's law E +u x B =R', where R' accounts for any nonideal effects, only determines the particle flow relative to the diffusing magnetic surfaces, it is necessary to solve a single two-dimensional generalized differential equation, (partial/partialt) delpsi. (delp - J x B) =0, to find the absolute velocity of a magnetic surface enclosing a fixed toroidal flux. This equation is linear but nonstandard in that it involves flux surface averages of the unknown velocity. Specification of R' and the cross-field ion and electron heat fluxes provides a closed system of equations. A time-dependent coordinate transformation is used to describe the diffusion of plasma quantities through magnetic surfaces of changing shape

Exactly solvable model of the two-dimensional electrical double layer.

Science.gov (United States)

Samaj, L; Bajnok, Z

2005-12-01

We consider equilibrium statistical mechanics of a simplified model for the ideal conductor electrode in an interface contact with a classical semi-infinite electrolyte, modeled by the two-dimensional Coulomb gas of pointlike unit charges in the stability-against-collapse regime of reduced inverse temperatures 0layer) carries some nonzero surface charge density. The model is mappable onto an integrable semi-infinite sine-Gordon theory with Dirichlet boundary conditions. The exact form-factor and boundary state information gained from the mapping provide asymptotic forms of the charge and number density profiles of electrolyte particles at large distances from the interface. The result for the asymptotic behavior of the induced electric potential, related to the charge density via the Poisson equation, confirms the validity of the concept of renormalized charge and the corresponding saturation hypothesis. It is documented on the nonperturbative result for the asymptotic density profile at a strictly nonzero beta that the Debye-Hückel beta-->0 limit is a delicate issue.

Numerical solutions of the Vlasov equation

International Nuclear Information System (INIS)

Satofuka, Nobuyuki; Morinishi, Koji; Nishida, Hidetoshi

1985-01-01

A numerical procedure is derived for the solutions of the one- and two-dimensional Vlasov-Poisson system equations. This numerical procedure consists of the phase space discretization and the integration of the resulting set of ordinary differential equations. In the phase space discretization, derivatives with respect to the phase space variable are approximated by a weighted sum of the values of the distribution function at properly chosen neighboring points. Then, the resulting set of ordinary differential equations is solved by using an appropriate time integration scheme. The results for linear Landau damping, nonlinear Landau damping and counter-streaming plasmas are investigated and compared with those of the splitting scheme. The proposed method is found to be very accurate and efficient. (author)

Two-dimensional topological field theories coupled to four-dimensional BF theory

International Nuclear Information System (INIS)

Montesinos, Merced; Perez, Alejandro

2008-01-01

Four-dimensional BF theory admits a natural coupling to extended sources supported on two-dimensional surfaces or string world sheets. Solutions of the theory are in one to one correspondence with solutions of Einstein equations with distributional matter (cosmic strings). We study new (topological field) theories that can be constructed by adding extra degrees of freedom to the two-dimensional world sheet. We show how two-dimensional Yang-Mills degrees of freedom can be added on the world sheet, producing in this way, an interactive (topological) theory of Yang-Mills fields with BF fields in four dimensions. We also show how a world sheet tetrad can be naturally added. As in the previous case the set of solutions of these theories are contained in the set of solutions of Einstein's equations if one allows distributional matter supported on two-dimensional surfaces. These theories are argued to be exactly quantizable. In the context of quantum gravity, one important motivation to study these models is to explore the possibility of constructing a background-independent quantum field theory where local degrees of freedom at low energies arise from global topological (world sheet) degrees of freedom at the fundamental level

Stochastic resonance induced by the novel random transitions of two-dimensional weak damping bistable duffing oscillator and bifurcation of moment equation

International Nuclear Information System (INIS)

Zhang Guangjun; Xu Jianxue; Wang Jue; Yue Zhifeng; Zou Hailin

2009-01-01

In this paper stochastic resonance induced by the novel random transitions of two-dimensional weak damping bistable Duffing oscillator is analyzed by moment method. This kind of novel transition refers to the one among three potential well on two sides of bifurcation point of original system at the presence of internal noise. Several conclusions are drawn. First, the semi-analytical result of stochastic resonance induced by the novel random transitions of two-dimensional weak damping bistable Duffing oscillator can be obtained, and the semi-analytical result is qualitatively compatible with the one of Monte Carlo simulation. Second, a bifurcation of double-branch fixed point curves occurs in the moment equations with noise intensity as their bifurcation parameter. Third, the bifurcation of moment equations corresponds to stochastic resonance of original system. Finally, the mechanism of stochastic resonance is presented from another viewpoint through analyzing the energy transfer induced by the bifurcation of moment equation.

2D Poisson sigma models with gauged vectorial supersymmetry

Energy Technology Data Exchange (ETDEWEB)

Bonezzi, Roberto [Dipartimento di Fisica ed Astronomia, Università di Bologna and INFN, Sezione di Bologna,via Irnerio 46, I-40126 Bologna (Italy); Departamento de Ciencias Físicas, Universidad Andres Bello,Republica 220, Santiago (Chile); Sundell, Per [Departamento de Ciencias Físicas, Universidad Andres Bello,Republica 220, Santiago (Chile); Torres-Gomez, Alexander [Departamento de Ciencias Físicas, Universidad Andres Bello,Republica 220, Santiago (Chile); Instituto de Ciencias Físicas y Matemáticas, Universidad Austral de Chile-UACh,Valdivia (Chile)

2015-08-12

In this note, we gauge the rigid vectorial supersymmetry of the two-dimensional Poisson sigma model presented in arXiv:1503.05625. We show that the consistency of the construction does not impose any further constraints on the differential Poisson algebra geometry than those required for the ungauged model. We conclude by proposing that the gauged model provides a first-quantized framework for higher spin gravity.

Macroscopic balance equations for two-phase flow models

International Nuclear Information System (INIS)

Hughes, E.D.

1979-01-01

The macroscopic, or overall, balance equations of mass, momentum, and energy are derived for a two-fluid model of two-phase flows in complex geometries. These equations provide a base for investigating methods of incorporating improved analysis methods into computer programs, such as RETRAN, which are used for transient and steady-state thermal-hydraulic analyses of nuclear steam supply systems. The equations are derived in a very general manner so that three-dimensional, compressible flows can be analysed. The equations obtained supplement the various partial differential equation two-fluid models of two-phase flow which have recently appeared in the literature. The primary objective of the investigation is the macroscopic balance equations. (Auth.)

General solution of the Dirac equation for quasi-two-dimensional electrons

Energy Technology Data Exchange (ETDEWEB)

Eremko, Alexander, E-mail: eremko@bitp.kiev.ua [Bogolyubov Institute for Theoretical Physics, Metrologichna Str., 14-b, Kyiv, 03680 (Ukraine); Brizhik, Larissa, E-mail: brizhik@bitp.kiev.ua [Bogolyubov Institute for Theoretical Physics, Metrologichna Str., 14-b, Kyiv, 03680 (Ukraine); Loktev, Vadim, E-mail: vloktev@bitp.kiev.ua [Bogolyubov Institute for Theoretical Physics, Metrologichna Str., 14-b, Kyiv, 03680 (Ukraine); National Technical University of Ukraine “KPI”, Peremohy av., 37, Kyiv, 03056 (Ukraine)

2016-06-15

The general solution of the Dirac equation for quasi-two-dimensional electrons confined in an asymmetric quantum well, is found. The energy spectrum of such a system is exactly calculated using special unitary operator and is shown to depend on the electron spin polarization. This solution contains free parameters, whose variation continuously transforms one known particular solution into another. As an example, two different cases are considered in detail: electron in a deep and in a strongly asymmetric shallow quantum well. The effective mass renormalized by relativistic corrections and Bychkov–Rashba coefficients are analytically obtained for both cases. It is demonstrated that the general solution transforms to the particular solutions, found previously (Eremko et al., 2015) with the use of spin invariants. The general solution allows to establish conditions at which a specific (accompanied or non-accompanied by Rashba splitting) spin state can be realized. These results can prompt the ways to control the spin degree of freedom via the synthesis of spintronic heterostructures with the required properties.

Poisson traces, D-modules, and symplectic resolutions.

Science.gov (United States)

Etingof, Pavel; Schedler, Travis

2018-01-01

We survey the theory of Poisson traces (or zeroth Poisson homology) developed by the authors in a series of recent papers. The goal is to understand this subtle invariant of (singular) Poisson varieties, conditions for it to be finite-dimensional, its relationship to the geometry and topology of symplectic resolutions, and its applications to quantizations. The main technique is the study of a canonical D-module on the variety. In the case the variety has finitely many symplectic leaves (such as for symplectic singularities and Hamiltonian reductions of symplectic vector spaces by reductive groups), the D-module is holonomic, and hence, the space of Poisson traces is finite-dimensional. As an application, there are finitely many irreducible finite-dimensional representations of every quantization of the variety. Conjecturally, the D-module is the pushforward of the canonical D-module under every symplectic resolution of singularities, which implies that the space of Poisson traces is dual to the top cohomology of the resolution. We explain many examples where the conjecture is proved, such as symmetric powers of du Val singularities and symplectic surfaces and Slodowy slices in the nilpotent cone of a semisimple Lie algebra. We compute the D-module in the case of surfaces with isolated singularities and show it is not always semisimple. We also explain generalizations to arbitrary Lie algebras of vector fields, connections to the Bernstein-Sato polynomial, relations to two-variable special polynomials such as Kostka polynomials and Tutte polynomials, and a conjectural relationship with deformations of symplectic resolutions. In the appendix we give a brief recollection of the theory of D-modules on singular varieties that we require.

Poisson traces, D-modules, and symplectic resolutions

Science.gov (United States)

Etingof, Pavel; Schedler, Travis

2018-03-01

We survey the theory of Poisson traces (or zeroth Poisson homology) developed by the authors in a series of recent papers. The goal is to understand this subtle invariant of (singular) Poisson varieties, conditions for it to be finite-dimensional, its relationship to the geometry and topology of symplectic resolutions, and its applications to quantizations. The main technique is the study of a canonical D-module on the variety. In the case the variety has finitely many symplectic leaves (such as for symplectic singularities and Hamiltonian reductions of symplectic vector spaces by reductive groups), the D-module is holonomic, and hence, the space of Poisson traces is finite-dimensional. As an application, there are finitely many irreducible finite-dimensional representations of every quantization of the variety. Conjecturally, the D-module is the pushforward of the canonical D-module under every symplectic resolution of singularities, which implies that the space of Poisson traces is dual to the top cohomology of the resolution. We explain many examples where the conjecture is proved, such as symmetric powers of du Val singularities and symplectic surfaces and Slodowy slices in the nilpotent cone of a semisimple Lie algebra. We compute the D-module in the case of surfaces with isolated singularities and show it is not always semisimple. We also explain generalizations to arbitrary Lie algebras of vector fields, connections to the Bernstein-Sato polynomial, relations to two-variable special polynomials such as Kostka polynomials and Tutte polynomials, and a conjectural relationship with deformations of symplectic resolutions. In the appendix we give a brief recollection of the theory of D-modules on singular varieties that we require.

Exact Solutions for Certain Nonlinear Autonomous Ordinary Differential Equations of the Second Order and Families of Two-Dimensional Autonomous Systems

Directory of Open Access Journals (Sweden)

M. P. Markakis

2010-01-01

Full Text Available Certain nonlinear autonomous ordinary differential equations of the second order are reduced to Abel equations of the first kind ((Ab-1 equations. Based on the results of a previous work, concerning a closed-form solution of a general (Ab-1 equation, and introducing an arbitrary function, exact one-parameter families of solutions are derived for the original autonomous equations, for the most of which only first integrals (in closed or parametric form have been obtained so far. Two-dimensional autonomous systems of differential equations of the first order, equivalent to the considered herein autonomous forms, are constructed and solved by means of the developed analysis.

On nonlinear equations associated with Lie algebras of diffeomorphism groups of two-dimensional manifolds

International Nuclear Information System (INIS)

Kashaev, R.M.; Savel'ev, M.V.; Savel'eva, S.A.

1990-01-01

Nonlinear equations associated through a zero curvature type representation with Lie algebras S 0 Diff T 2 and of infinitesimal diffeomorphisms of (S 1 ) 2 , and also with a new infinite-dimensional Lie algebras. In particular, the general solution (in the sense of the Goursat problem) of the heavently equation which describes self-dual Einstein spaces with one rotational Killing symmetry is discussed, as well as the solutions to a generalized equation. The paper is supplied with Appendix containing the definition of the continuum graded Lie algebras and the general construction of the nonlinear equations associated with them. 11 refs

Poisson solvers for self-consistent multi-particle simulations

International Nuclear Information System (INIS)

Qiang, J; Paret, S

2014-01-01

Self-consistent multi-particle simulation plays an important role in studying beam-beam effects and space charge effects in high-intensity beams. The Poisson equation has to be solved at each time-step based on the particle density distribution in the multi-particle simulation. In this paper, we review a number of numerical methods that can be used to solve the Poisson equation efficiently. The computational complexity of those numerical methods will be O(N log(N)) or O(N) instead of O(N2), where N is the total number of grid points used to solve the Poisson equation

On the confinement of a Dirac particle to a two-dimensional ring

International Nuclear Information System (INIS)

Bakke, K.; Furtado, C.

2012-01-01

In this contribution, we propose a new model for studying the confinement of a spin-half particle to a two-dimensional quantum ring for systems described by the Dirac equation by introducing a new coupling into the Dirac equation. We show that the introduction of this new coupling into the Dirac equation yields a generalization of the two-dimensional quantum ring model proposed by Tan and Inkson [W.-C. Tan, J.C. Inkson, Semicond. Sci. Technol. 11 (1996) 1635] for relativistic spin-half quantum particles. -- Highlights: ► Two-dimensional ring model for condensed matter systems described by the Dirac equation. ► Exact solutions of the Dirac equation. ► Persistent currents for Dirac-like systems confined to a two-dimensional quantum ring.

Quasistationary model of high-current relativistic electron beam. 1. Exact solution of Poisson equations

International Nuclear Information System (INIS)

Brenner, S.E.; Gandyl', E.M.; Podkopaev, A.P.

1995-01-01

The dynamics of high-current relativistic electron beam moving trough the cylindrical drift space has been modelled by the large particles, the shape of which allows to solve the Poisson equations exactly, and in such a way to avoid the linearization being usually used in those problems. The expressions for the components of own electric field of electron beam passing through the cylindrical drift space have been obtained. (author). 11 refs., 1 fig

Characteristics and stability analyses of transient one-dimensional two-phase flow equations and their finite difference approximations

International Nuclear Information System (INIS)

Lyczkowski, R.W.; Gidaspow, D.; Solbrig, C.W.; Hughes, E.D.

1975-01-01

Equation systems describing one-dimensional, transient, two-phase flow with separate continuity, momentum, and energy equations for each phase are classified by use of the method of characteristics. Little attempt is made to justify the physics of these equations. Many of the equation systems possess complex-valued characteristics and hence, according to well-known mathematical theorems, are not well-posed as initial-value problems (IVPs). Real-valued characteristics are necessary but not sufficient to insure well-posedness. In the absence of lower order source or sink terms (potential type flows), which can affect the well-posedness of IVPs, the complex characteristics associated with these two-phase flow equations imply unbounded exponential growth for disturbances of all wavelengths. Analytical and numerical examples show that the ill-posedness of IVPs for the two-phase flow partial differential equations which possess complex characteristics produce unstable numerical schemes. These unstable numerical schemes can produce apparently stable and even accurate results if the growth rate resulting from the complex characteristics remains small throughout the time span of the numerical experiment or if sufficient numerical damping is present for the increment size used. Other examples show that clearly nonphysical numerical instabilities resulting from the complex characteristics can be produced. These latter types of numerical instabilities are shown to be removed by the addition of physically motivated differential terms which eliminate the complex characteristics. (auth)

Lie algebra contractions on two-dimensional hyperboloid

International Nuclear Information System (INIS)

Pogosyan, G. S.; Yakhno, A.

2010-01-01

The Inoenue-Wigner contraction from the SO(2, 1) group to the Euclidean E(2) and E(1, 1) group is used to relate the separation of variables in Laplace-Beltrami (Helmholtz) equations for the four corresponding two-dimensional homogeneous spaces: two-dimensional hyperboloids and two-dimensional Euclidean and pseudo-Euclidean spaces. We show how the nine systems of coordinates on the two-dimensional hyperboloids contracted to the four systems of coordinates on E 2 and eight on E 1,1 . The text was submitted by the authors in English.

Poisson equations of rotational motion for a rigid triaxial body with application to a tumbling artificial satellite

Science.gov (United States)

Liu, J. J. F.; Fitzpatrick, P. M.

1975-01-01

A mathematical model is developed for studying the effects of gravity gradient torque on the attitude stability of a tumbling triaxial rigid satellite. Poisson equations are used to investigate the rotation of the satellite (which is in elliptical orbit about an attracting point mass) about its center of mass. An averaging method is employed to obtain an intermediate set of differential equations for the nonresonant, secular behavior of the osculating elements which describe the rotational motions of the satellite, and the averaged equations are then integrated to obtain long-term secular solutions for the osculating elements.

Dirac equation in 5- and 6-dimensional curved space-time manifolds

International Nuclear Information System (INIS)

Vladimirov, Yu.S.; Popov, A.D.

1984-01-01

The program of plotting unified multidimensional theory of gravitation, electromagnetism and electrically charged matter with transition from 5-dimensional variants to 6-dimensional theory possessing signature (+----+) is developed. For recording the Dirac equation in 5- and 6-dimensional curved space-time manifolds the tetrad formalism and γ-matrix formulation of the General Relativity Theory are used. It is shown that the 6-dimensional theory case unifies the two private cases of 5-dimensional theory and corresponds to two possibilities of the theory developed by Kadyshevski

Comparison between two bivariate Poisson distributions through the ...

African Journals Online (AJOL)

These two models express themselves by their probability mass function. ... To remedy this problem, Berkhout and Plug proposed a bivariate Poisson distribution accepting the correlation as well negative, equal to zero, that positive.

Darboux transformation and explicit solutions for some (2+1)-dimensional equations

International Nuclear Information System (INIS)

Wang Yan; Shen Lijuan; Du Dianlou

2007-01-01

Three systems of (2+1)-dimensional soliton equations and their decompositions into the (1+1)-dimensional soliton equations are proposed. These equations include KPI, CKP, MKPI. With the help of Darboux transformation of (1+1)-dimensional equations, we get the explicit solutions of the (2+1)-dimensional equations

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

Science.gov (United States)

Lu, Benzhuo; Zhou, Y.C.

2011-01-01

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

Solution of two-dimensional diffusion equation for hexagonal cells by the finite Fourier transformation

International Nuclear Information System (INIS)

Kobayashi, Keisuke

1975-01-01

A method of solution is presented for a monoenergetic diffusion equation in two-dimensional hexagonal cells by a finite Fourier transformation. Up to the present, the solution by the finite Fourier transformation has been developed for x-y, r-z and x-y-z geometries, and the flux and current at the boundary are obtained in terms of Fourier series. It is shown here that the method can be applied to hexagonal cells and the expansion of boundary values in a Legendre polynomials gives numerically a higher accuracy than is obtained by a Fourier series. (orig.) [de

General method and exact solutions to a generalized variable-coefficient two-dimensional KdV equation

International Nuclear Information System (INIS)

Chen, Yong; Shanghai Jiao-Tong Univ., Shangai; Chinese Academy of sciences, Beijing

2005-01-01

A general method to uniformly construct exact solutions in terms of special function of nonlinear partial differential equations is presented by means of a more general ansatz and symbolic computation. Making use of the general method, we can successfully obtain the solutions found by the method proposed by Fan (J. Phys. A., 36 (2003) 7009) and find other new and more general solutions, which include polynomial solutions, exponential solutions, rational solutions, triangular periodic wave solution, soliton solutions, soliton-like solutions and Jacobi, Weierstrass doubly periodic wave solutions. A general variable-coefficient two-dimensional KdV equation is chosen to illustrate the method. As a result, some new exact soliton-like solutions are obtained. planets. The numerical results are given in tables. The results are discussed in the conclusion

Combined analytical-numerical procedure to solve multigroup spherical harmonics equations in two-dimensional r-z geometry

International Nuclear Information System (INIS)

Matausek, M.V.; Milosevic, M.

1986-01-01

In the present paper a generalization is performed of a procedure to solve multigroup spherical harmonics equations, which has originally been proposed and developed for one-dimensional systems in cylindrical or spherical geometry, and later extended for a special case of a two-dimensional system in r-z geometry. The expressions are derived for the axial and the radial dependence of the group values of the neutron flux moments, in the P-3 approximation of the spherical harmonics method, in a cylindrically symmetrical system with an arbitrary number of material regions in both r- and z-directions. In the special case of an axially homogeneous system, these expressions reduce to the relations derived previously. (author)

Renormalized charge in a two-dimensional model of colloidal suspension from hypernetted chain approach.

Science.gov (United States)

Camargo, Manuel; Téllez, Gabriel

2008-04-07

The renormalized charge of a simple two-dimensional model of colloidal suspension was determined by solving the hypernetted chain approximation and Ornstein-Zernike equations. At the infinite dilution limit, the asymptotic behavior of the correlation functions is used to define the effective interactions between the components of the system and these effective interactions were compared to those derived from the Poisson-Boltzmann theory. The results we obtained show that, in contrast to the mean-field theory, the renormalized charge does not saturate, but exhibits a maximum value and then decays monotonically as the bare charge increases. The results also suggest that beyond the counterion layer near to the macroion surface, the ionic cloud is not a diffuse layer which can be handled by means of the linearized theory, as the two-state model claims, but a more complex structure is settled by the correlations between microions.

Whitham modulation theory for (2  +  1)-dimensional equations of Kadomtsev–Petviashvili type

Science.gov (United States)

Ablowitz, Mark J.; Biondini, Gino; Rumanov, Igor

2018-05-01

Whitham modulation theory for certain two-dimensional evolution equations of Kadomtsev–Petviashvili (KP) type is presented. Three specific examples are considered in detail: the KP equation, the two-dimensional Benjamin–Ono (2DBO) equation and a modified KP (m2KP) equation. A unified derivation is also provided. In the case of the m2KP equation, the corresponding Whitham modulation system exhibits features different from the other two. The approach presented here does not require integrability of the original evolution equation. Indeed, while the KP equation is known to be a completely integrable equation, the 2DBO equation and the m2KP equation are not known to be integrable. In each of the cases considered, the Whitham modulation system obtained consists of five first-order quasilinear partial differential equations. The Riemann problem (i.e. the analogue of the Gurevich–Pitaevskii problem) for the one-dimensional reduction of the m2KP equation is studied. For the m2KP equation, the system of modulation equations is used to analyze the linear stability of traveling wave solutions.

Two-dimensional time dependent Riemann solvers for neutron transport

International Nuclear Information System (INIS)

Brunner, Thomas A.; Holloway, James Paul

2005-01-01

A two-dimensional Riemann solver is developed for the spherical harmonics approximation to the time dependent neutron transport equation. The eigenstructure of the resulting equations is explored, giving insight into both the spherical harmonics approximation and the Riemann solver. The classic Roe-type Riemann solver used here was developed for one-dimensional problems, but can be used in multidimensional problems by treating each face of a two-dimensional computation cell in a locally one-dimensional way. Several test problems are used to explore the capabilities of both the Riemann solver and the spherical harmonics approximation. The numerical solution for a simple line source problem is compared to the analytic solution to both the P 1 equation and the full transport solution. A lattice problem is used to test the method on a more challenging problem

Non-isothermal Smoluchowski-Poisson equation as a singular limit of the Navier-Stokes-Fourier-Poisson system

Czech Academy of Sciences Publication Activity Database

Feireisl, Eduard; Laurençot, P.

2007-01-01

Roč. 88, - (2007), s. 325-349 ISSN 0021-7824 R&D Projects: GA ČR GA201/05/0164 Institutional research plan: CEZ:AV0Z10190503 Keywords : Navier-Stokes-Fourier- Poisson system * Smoluchowski- Poisson system * singular limit Subject RIV: BA - General Mathematics Impact factor: 1.118, year: 2007

An accessible four-dimensional treatment of Maxwell's equations in terms of differential forms

Science.gov (United States)

Sá, Lucas

2017-03-01

Maxwell’s equations are derived in terms of differential forms in the four-dimensional Minkowski representation, starting from the three-dimensional vector calculus differential version of these equations. Introducing all the mathematical and physical concepts needed (including the tool of differential forms), using only knowledge of elementary vector calculus and the local vector version of Maxwell’s equations, the equations are reduced to a simple and elegant set of two equations for a unified quantity, the electromagnetic field. The treatment should be accessible for students taking a first course on electromagnetism.

Universal equations of unsteady two-dimensional MHD boundary layer whose temperature varies with time

Directory of Open Access Journals (Sweden)

Boričić Zoran

2009-01-01

Full Text Available This paper concerns with unsteady two-dimensional temperature laminar magnetohydrodynamic (MHD boundary layer of incompressible fluid. It is assumed that induction of outer magnetic field is function of longitudinal coordinate with force lines perpendicular to the body surface on which boundary layer forms. Outer electric filed is neglected and magnetic Reynolds number is significantly lower then one i.e. considered problem is in inductionless approximation. Characteristic properties of fluid are constant because velocity of flow is much lower than speed of light and temperature difference is small enough (under 50ºC . Introduced assumptions simplify considered problem in sake of mathematical solving, but adopted physical model is interesting from practical point of view, because its relation with large number of technically significant MHD flows. Obtained partial differential equations can be solved with modern numerical methods for every particular problem. Conclusions based on these solutions are related only with specific temperature MHD boundary layer problem. In this paper, quite different approach is used. First new variables are introduced and then sets of similarity parameters which transform equations on the form which don't contain inside and in corresponding boundary conditions characteristics of particular problems and in that sense equations are considered as universal. Obtained universal equations in appropriate approximation can be solved numerically once for all. So-called universal solutions of equations can be used to carry out general conclusions about temperature MHD boundary layer and for calculation of arbitrary particular problems. To calculate any particular problem it is necessary also to solve corresponding momentum integral equation.

Closed-form solution of a two-dimensional fuel temperature model for TRIGA-type reactors

Energy Technology Data Exchange (ETDEWEB)

Rivard, J B [Sandia Laboratories (United States)

1974-07-01

If azimuthal power density variations are ignored, the steady-state temperature distribution within a TRIGA-type fuel element is given by the solution of the Poisson equation in two dimensions (r and z) . This paper presents a closed-form solution of this equation as a function of the axial and radial power density profiles, the conductivity of the U-ZrH, the inlet temperature, specific heat and flow rate of the coolant, and the overall heat transfer coefficient. The method begins with the development of a system of linear ordinary differential equations describing mass and energy balances in the fuel and coolant. From the solution of this system, an expression for the second derivative of the fuel temperature distribution in the axial (z) direction is found. Substitution of this expression into the Poisson equation for T(r,z) reduces it from a partial differential equation to an ordinary differential equation in r, which is subsequently solved in closed-form. The results of typical calculations using the model are presented. (author)

Dimensional reduction of a general advection–diffusion equation in 2D channels

Science.gov (United States)

Kalinay, Pavol; Slanina, František

2018-06-01

Diffusion of point-like particles in a two-dimensional channel of varying width is studied. The particles are driven by an arbitrary space dependent force. We construct a general recurrence procedure mapping the corresponding two-dimensional advection-diffusion equation onto the longitudinal coordinate x. Unlike the previous specific cases, the presented procedure enables us to find the one-dimensional description of the confined diffusion even for non-conservative (vortex) forces, e.g. caused by flowing solvent dragging the particles. We show that the result is again the generalized Fick–Jacobs equation. Despite of non existing scalar potential in the case of vortex forces, the effective one-dimensional scalar potential, as well as the corresponding quasi-equilibrium and the effective diffusion coefficient can be always found.

Conservation laws for two (2 + 1)-dimensional differential-difference systems

International Nuclear Information System (INIS)

Yu Guofu; Tam, H.-W.

2006-01-01

Two integrable differential-difference equations are considered. One is derived from the discrete BKP equation and the other is a symmetric (2 + 1)-dimensional Lotka-Volterra equation. An infinite number of conservation laws for the two differential-difference equations are deduced

Study of the nonlinear three-dimensional Debye screening in plasmas

International Nuclear Information System (INIS)

Lin Chang; Zhao Jinbao; Zhang Xiulian

2000-01-01

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

Heat transport in two-dimensional materials by directly solving the phonon Boltzmann equation under Callaway's dual relaxation model

Science.gov (United States)

Guo, Yangyu; Wang, Moran

2017-10-01

The single mode relaxation time approximation has been demonstrated to greatly underestimate the lattice thermal conductivity of two-dimensional materials due to the collective effect of phonon normal scattering. Callaway's dual relaxation model represents a good approximation to the otherwise ab initio solution of the phonon Boltzmann equation. In this work we develop a discrete-ordinate-method (DOM) scheme for the numerical solution of the phonon Boltzmann equation under Callaway's model. Heat transport in a graphene ribbon with different geometries is modeled by our scheme, which produces results quite consistent with the available molecular dynamics, Monte Carlo simulations, and experimental measurements. Callaway's lattice thermal conductivity model with empirical boundary scattering rates is examined and shown to overestimate or underestimate the direct DOM solution. The length convergence of the lattice thermal conductivity of a rectangular graphene ribbon is explored and found to depend appreciably on the ribbon width, with a semiquantitative correlation provided between the convergence length and the width. Finally, we predict the existence of a phonon Knudsen minimum in a graphene ribbon only at a low system temperature and isotope concentration so that the average normal scattering rate is two orders of magnitude stronger than the intrinsic resistive one. The present work will promote not only the methodology for the solution of the phonon Boltzmann equation but also the theoretical modeling and experimental detection of hydrodynamic phonon transport in two-dimensional materials.

Polynomial Poisson algebras: Gel'fand-Kirillov problem and Poisson spectra

OpenAIRE

Lecoutre, César

2014-01-01

We study the fields of fractions and the Poisson spectra of polynomial Poisson algebras.\\ud \\ud First we investigate a Poisson birational equivalence problem for polynomial Poisson algebras over a field of arbitrary characteristic. Namely, the quadratic Poisson Gel'fand-Kirillov problem asks whether the field of fractions of a Poisson algebra is isomorphic to the field of fractions of a Poisson affine space, i.e. a polynomial algebra such that the Poisson bracket of two generators is equal to...

On the Fedosov deformation quantization beyond the regular Poisson manifolds

International Nuclear Information System (INIS)

Dolgushev, V.A.; Isaev, A.P.; Lyakhovich, S.L.; Sharapov, A.A.

2002-01-01

A simple iterative procedure is suggested for the deformation quantization of (irregular) Poisson brackets associated to the classical Yang-Baxter equation. The construction is shown to admit a pure algebraic reformulation giving the Universal Deformation Formula (UDF) for any triangular Lie bialgebra. A simple proof of classification theorem for inequivalent UDF's is given. As an example the explicit quantization formula is presented for the quasi-homogeneous Poisson brackets on two-plane

A six-mode truncation of the Navier-Stokes equations on a two-dimensional torus: a numerical study

International Nuclear Information System (INIS)

Angelo, P.M.; Riela, G.

1981-01-01

We study a model obtained from a six-mode truncation of the Navier-Stokes equations for a two-dimensional incompressible fluid on a torus. We find that at low values of the Reynolds number R the dynamics is characterized by fixed points and, at large values of R, by two stable periodic orbits; at intermediate values of R two infinite sequences of bifurcations of periodic orbits into periodic orbits of doubled period lead to two regions of ''turbulent'' or ''chaotic'' behaviour. The turbulent regions end up for values of R for which stable periodic orbits appear. (author)

Penyelesaian Persamaan Poisson 2D dengan Menggunakan Metode Gauss-Seidel dan Conjugate Gradien

OpenAIRE

Mahmudah, Dewi Erla; Naf'an, Muhammad Zidny

2017-01-01

In this paper we focus on solution of 2D Poisson equation numerically. 2D Poisson equation is a partial differential equation of second order elliptical type. This equation is a particular form or non-homogeneous form of the Laplace equation. The solution of 2D Poisson equation is performed numerically using Gauss Seidel method and Conjugate Gradient method. The result is the value using Gauss Seidel method and Conjugate Gradient method is same. But, consider the iteration process, the conver...

Five-dimensional Monopole Equation with Hedge-Hog Ansatz and Abel's Differential Equation

OpenAIRE

Kihara, Hironobu

2008-01-01

We review the generalized monopole in the five-dimensional Euclidean space. A numerical solution with the Hedge-Hog ansatz is studied. The Bogomol'nyi equation becomes a second order autonomous non-linear differential equation. The equation can be translated into the Abel's differential equation of the second kind and is an algebraic differential equation.

Poisson structure and symmetry in the Chern-Simons formulation of (2 + 1)-dimensional gravity

International Nuclear Information System (INIS)

Meusburger, C; Schroers, B J

2003-01-01

In the formulation of (2 + 1)-dimensional gravity as a Chern-Simons gauge theory, the phase space is the moduli space of flat Poincare group connections. Using the combinatorial approach developed by Fock and Rosly, we give an explicit description of the phase space and its Poisson structure for the general case of a genus g oriented surface with punctures representing particles and a boundary playing the role of spatial infinity. We give a physical interpretation and explain how the degrees of freedom associated with each handle and each particle can be decoupled. The symmetry group of the theory combines an action of the mapping class group with asymptotic Poincare transformations in a nontrivial fashion. We derive the conserved quantities associated with the latter and show that the mapping class group of the surface acts on the phase space via Poisson isomorphisms

Solution of 3-dimensional diffusion equation by finite Fourier transformation

International Nuclear Information System (INIS)

Krishnani, P.D.

1978-01-01

Three dimensional diffusion equation in Cartesian co-ordinates is solved by using the finite Fourier transformation. This method is different from the usual Fourier transformation method in the sense that the solutions are obtained without performing the inverse Fourier transformation. The advantage has been taken of the fact that the flux is finite and integrable in the finite region. By applying this condition, a two-dimensional integral equation, involving flux and its normal derivative at the boundary, is obtained. By solving this equation with given boundary conditions, all of the boundary values are determined. In order to calculate the flux inside the region, flux is expanded into three-dimensional Fourier series. The Fourier coefficients of the flux in the region are calculated from the boundary values. The advantage of this method is that the integrated flux is obtained without knowing the fluxes inside the region as in the case of finite difference method. (author)

Computation of solar perturbations with Poisson series

Science.gov (United States)

Broucke, R.

1974-01-01

Description of a project for computing first-order perturbations of natural or artificial satellites by integrating the equations of motion on a computer with automatic Poisson series expansions. A basic feature of the method of solution is that the classical variation-of-parameters formulation is used rather than rectangular coordinates. However, the variation-of-parameters formulation uses the three rectangular components of the disturbing force rather than the classical disturbing function, so that there is no problem in expanding the disturbing function in series. Another characteristic of the variation-of-parameters formulation employed is that six rather unusual variables are used in order to avoid singularities at the zero eccentricity and zero (or 90 deg) inclination. The integration process starts by assuming that all the orbit elements present on the right-hand sides of the equations of motion are constants. These right-hand sides are then simple Poisson series which can be obtained with the use of the Bessel expansions of the two-body problem in conjunction with certain interation methods. These Poisson series can then be integrated term by term, and a first-order solution is obtained.

Symmetries, integrals, and three-dimensional reductions of Plebanski's second heavenly equation

International Nuclear Information System (INIS)

Neyzi, F.; Sheftel, M. B.; Yazici, D.

2007-01-01

We study symmetries and conservation laws for Plebanski's second heavenly equation written as a first-order nonlinear evolutionary system which admits a multi-Hamiltonian structure. We construct an optimal system of one-dimensional subalgebras and all inequivalent three-dimensional symmetry reductions of the original four-dimensional system. We consider these two-component evolutionary systems in three dimensions as natural candidates for integrable systems

Zakharov-Shabat-Mikhailov scheme of construction of two-dimensional completely integrable field theories

International Nuclear Information System (INIS)

Chudnovsky, D.V.; Columbia Univ., New York; Chudnovsky, G.V.; Columbia Univ., New York

1980-01-01

General algebraic and analytic formalism for derivation and solution of general two dimensional field theory equations of Zakharov-Shabat-Mikhailov type is presented. The examples presented show that this class of equations covers most of the known two-dimensional completely integrable equations. Possible generalizations for four dimensional systems require detailed analysis of Baecklund transformation of these equations. Baecklund transformation is presented in the form of Riemann problem and one special case of dual symmetry is worked out. (orig.)

Exact Solutions to (2+1)-Dimensional Kaup-Kupershmidt Equation

International Nuclear Information System (INIS)

Lu Hailing; Liu Xiqiang

2009-01-01

In this paper, by using the symmetry method, the relationships between new explicit solutions and old ones of the (2+1)-dimensional Kaup-Kupershmidt (KK) equation are presented. We successfully obtain more general exact travelling wave solutions for (2+1)-dimensional KK equation by the symmetry method and the (G'/G)-expansion method. Consequently, we find some new solutions of (2+1)-dimensional KK equation, including similarity solutions, solitary wave solutions, and periodic solutions. (general)

dimensional nonlinear evolution equations

Indian Academy of Sciences (India)

in real-life situations, it is important to find their exact solutions. Further, in ... But only little work is done on the high-dimensional equations. .... Similarly, to determine the values of d and q, we balance the linear term of the lowest order in eq.

Two-dimensional model for subthreshold current and subthreshold swing of graded-channel dual-material double-gate (GCDMDG) MOSFETs

Science.gov (United States)

Goel, Ekta; Kumar, Sanjay; Singh, Balraj; Singh, Kunal; Jit, Satyabrata

2017-06-01

The subthreshold performance of graded-channel dual-material double-gate (GCDMDG) MOSFETs is examined through two-dimensional (2D) analytical modeling of subthreshold-current (SC) and subthreshold-swing (SS). The potential function obtained by using the parabolic approach to solve the 2D Poisson's equation, has been used to formulate SC and SS characteristics of the device. The variations of SS against different device parameters have been obtained with the help of effective conduction path parameter. The SC and SS characteristics of the GCDMDG MOS transistor have been compared with those of the dual-material double-gate (DMDG) and simple graded-channel double-gate (GCDG) MOS structures to show its better subthreshold characteristics over the latter two devices. The results of the developed model are well-agreed with the commercially available SILVACO ATLAS™ simulator data.

Prolongation Loop Algebras for a Solitonic System of Equations

Directory of Open Access Journals (Sweden)

Maria A. Agrotis

2006-11-01

Full Text Available We consider an integrable system of reduced Maxwell-Bloch equations that describes the evolution of an electromagnetic field in a two-level medium that is inhomogeneously broadened. We prove that the relevant Bäcklund transformation preserves the reality of the n-soliton potentials and establish their pole structure with respect to the broadening parameter. The natural phase space of the model is embedded in an infinite dimensional loop algebra. The dynamical equations of the model are associated to an infinite family of higher order Hamiltonian systems that are in involution. We present the Hamiltonian functions and the Poisson brackets between the extended potentials.

Generalized Runge-Kutta method for two- and three-dimensional space-time diffusion equations with a variable time step

International Nuclear Information System (INIS)

Aboanber, A.E.; Hamada, Y.M.

2008-01-01

An extensive knowledge of the spatial power distribution is required for the design and analysis of different types of current-generation reactors, and that requires the development of more sophisticated theoretical methods. Therefore, the need to develop new methods for multidimensional transient reactor analysis still exists. The objective of this paper is to develop a computationally efficient numerical method for solving the multigroup, multidimensional, static and transient neutron diffusion kinetics equations. A generalized Runge-Kutta method has been developed for the numerical integration of the stiff space-time diffusion equations. The method is fourth-order accurate, using an embedded third-order solution to arrive at an estimate of the truncation error for automatic time step control. In addition, the A(α)-stability properties of the method are investigated. The analyses of two- and three-dimensional benchmark problems as well as static and transient problems, demonstrate that very accurate solutions can be obtained with assembly-sized spatial meshes. Preliminary numerical evaluations using two- and three-dimensional finite difference codes showed that the presented generalized Runge-Kutta method is highly accurate and efficient when compared with other optimized iterative numerical and conventional finite difference methods

Two-dimensional wave propagation in layered periodic media

KAUST Repository

Quezada de Luna, Manuel

2014-09-16

We study two-dimensional wave propagation in materials whose properties vary periodically in one direction only. High order homogenization is carried out to derive a dispersive effective medium approximation. One-dimensional materials with constant impedance exhibit no effective dispersion. We show that a new kind of effective dispersion may arise in two dimensions, even in materials with constant impedance. This dispersion is a macroscopic effect of microscopic diffraction caused by spatial variation in the sound speed. We analyze this dispersive effect by using highorder homogenization to derive an anisotropic, dispersive effective medium. We generalize to two dimensions a homogenization approach that has been used previously for one-dimensional problems. Pseudospectral solutions of the effective medium equations agree to high accuracy with finite volume direct numerical simulations of the variable-coeffi cient equations.

Poisson's ratio of fiber-reinforced composites

Science.gov (United States)

Christiansson, Henrik; Helsing, Johan

1996-05-01

Poisson's ratio flow diagrams, that is, the Poisson's ratio versus the fiber fraction, are obtained numerically for hexagonal arrays of elastic circular fibers in an elastic matrix. High numerical accuracy is achieved through the use of an interface integral equation method. Questions concerning fixed point theorems and the validity of existing asymptotic relations are investigated and partially resolved. Our findings for the transverse effective Poisson's ratio, together with earlier results for random systems by other authors, make it possible to formulate a general statement for Poisson's ratio flow diagrams: For composites with circular fibers and where the phase Poisson's ratios are equal to 1/3, the system with the lowest stiffness ratio has the highest Poisson's ratio. For other choices of the elastic moduli for the phases, no simple statement can be made.

Spectral/hp least-squares finite element formulation for the Navier-Stokes equations

International Nuclear Information System (INIS)

Pontaza, J.P.; Reddy, J.N.

2003-01-01

We consider the application of least-squares finite element models combined with spectral/hp methods for the numerical solution of viscous flow problems. The paper presents the formulation, validation, and application of a spectral/hp algorithm to the numerical solution of the Navier-Stokes equations governing two- and three-dimensional stationary incompressible and low-speed compressible flows. The Navier-Stokes equations are expressed as an equivalent set of first-order equations by introducing vorticity or velocity gradients as additional independent variables and the least-squares method is used to develop the finite element model. High-order element expansions are used to construct the discrete model. The discrete model thus obtained is linearized by Newton's method, resulting in a linear system of equations with a symmetric positive definite coefficient matrix that is solved in a fully coupled manner by a preconditioned conjugate gradient method. Spectral convergence of the L 2 least-squares functional and L 2 error norms is verified using smooth solutions to the two-dimensional stationary Poisson and incompressible Navier-Stokes equations. Numerical results for flow over a backward-facing step, steady flow past a circular cylinder, three-dimensional lid-driven cavity flow, and compressible buoyant flow inside a square enclosure are presented to demonstrate the predictive capability and robustness of the proposed formulation

New high accuracy super stable alternating direction implicit methods for two and three dimensional hyperbolic damped wave equations

Directory of Open Access Journals (Sweden)

R.K. Mohanty

2014-01-01

Full Text Available In this paper, we report new three level implicit super stable methods of order two in time and four in space for the solution of hyperbolic damped wave equations in one, two and three space dimensions subject to given appropriate initial and Dirichlet boundary conditions. We use uniform grid points both in time and space directions. Our methods behave like fourth order accurate, when grid size in time-direction is directly proportional to the square of grid size in space-direction. The proposed methods are super stable. The resulting system of algebraic equations is solved by the Gauss elimination method. We discuss new alternating direction implicit (ADI methods for two and three dimensional problems. Numerical results and the graphical representation of numerical solution are presented to illustrate the accuracy of the proposed methods.

The (2+1)-dimensional axial universes—solutions to the Einstein equations, dimensional reduction points and Klein–Fock–Gordon waves

International Nuclear Information System (INIS)

Fiziev, P P; Shirkov, D V

2012-01-01

The paper presents a generalization and further development of our recent publications, where solutions of the Klein–Fock–Gordon equation defined on a few particular D = (2 + 1)-dimensional static spacetime manifolds were considered. The latter involve toy models of two-dimensional spaces with axial symmetry, including dimensional reduction to the one-dimensional space as a singular limiting case. Here, the non-static models of space geometry with axial symmetry are under consideration. To make these models closer to physical reality, we define a set of ‘admissible’ shape functions ρ(t, z) as the (2 + 1)-dimensional Einstein equation solutions in the vacuum spacetime, in the presence of the Λ-term and for the spacetime filled with the standard ‘dust’. It is curious that in the last case the Einstein equations reduce to the well-known Monge–Ampère equation, thus enabling one to obtain the general solution of the Cauchy problem, as well as a set of other specific solutions involving one arbitrary function. A few explicit solutions of the Klein–Fock–Gordon equation in this set are given. An interesting qualitative feature of these solutions relates to the dimensional reduction points, their classification and time behavior. In particular, these new entities could provide us with novel insight into the nature of P- and T-violations and of the Big Bang. A short comparison with other attempts to utilize the dimensional reduction of the spacetime is given. (paper)

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

Science.gov (United States)

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

2017-09-01

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

Poisson structures for reduced non-holonomic systems

International Nuclear Information System (INIS)

Ramos, Arturo

2004-01-01

Borisov, Mamaev and Kilin have recently found certain Poisson structures with respect to which the reduced and rescaled systems of certain non-holonomic problems, involving rolling bodies without slipping, become Hamiltonian, the Hamiltonian function being the reduced energy. We study further the algebraic origin of these Poisson structures, showing that they are of rank 2 and therefore the mentioned rescaling is not necessary. We show that they are determined, up to a non-vanishing factor function, by the existence of a system of first-order differential equations providing two integrals of motion. We generalize the form of the Poisson structures and extend their domain of definition. We apply the theory to the rolling disc, the Routh's sphere, the ball rolling on a surface of revolution, and its special case of a ball rolling inside a cylinder

Solutions and Conservation Laws of a (2+1-Dimensional Boussinesq Equation

Directory of Open Access Journals (Sweden)

Letlhogonolo Daddy Moleleki

2013-01-01

Full Text Available We study a nonlinear evolution partial differential equation, namely, the (2+1-dimensional Boussinesq equation. For the first time Lie symmetry method together with simplest equation method is used to find the exact solutions of the (2+1-dimensional Boussinesq equation. Furthermore, the new conservation theorem due to Ibragimov will be utilized to construct the conservation laws of the (2+1-dimensional Boussinesq equation.

Exact Solutions for Two Equation Hierarchies

International Nuclear Information System (INIS)

Song-Lin, Zhao; Da-Jun, Zhang; Jie, Ji

2010-01-01

Bilinear forms and double-Wronskian solutions are given for two hierarchies, the (2+1)-dimensional breaking Ablowitz–Kaup–Newell–Segur (AKNS) hierarchy and the negative order AKNS hierarchy. According to some choices of the coefficient matrix in the Wronskian condition equation set, we obtain some kinds of solutions for these two hierarchies, such as solitons, Jordan block solutions, rational solutions, complexitons and mixed solutions. (general)

MARG2D code. 1. Eigenvalue problem for two dimensional Newcomb equation

Energy Technology Data Exchange (ETDEWEB)

Tokuda, Shinji [Japan Atomic Energy Research Inst., Naka, Ibaraki (Japan). Naka Fusion Research Establishment; Watanabe, Tomoko

1997-10-01

A new method and a code MARG2D have been developed to solve the 2-dimensional Newcomb equation which plays an important role in the magnetohydrodynamic (MHD) stability analysis in an axisymmetric toroidal plasma such as a tokamak. In the present formulation, an eigenvalue problem is posed for the 2-D Newcomb equation, where the weight function (the kinetic energy integral) and the boundary conditions at rational surfaces are chosen so that an eigenfunction correctly behaves as the linear combination of the small solution and the analytical solutions around each of the rational surfaces. Thus, the difficulty on solving the 2-D Newcomb equation has been resolved. By using the MARG2D code, the ideal MHD marginally stable state can be identified for a 2-D toroidal plasma. The code is indispensable on computing the outer-region matching data necessary for the resistive MHD stability analysis. Benchmark with ERATOJ, an ideal MHD stability code, has been carried out and the MARG2D code demonstrates that it indeed identifies both stable and marginally stable states against ideal MHD motion. (author)

Finite element solution of two dimensional time dependent heat equation

International Nuclear Information System (INIS)

Maaz

1999-01-01

A Microsoft Windows based computer code, named FHEAT, has been developed for solving two dimensional heat problems in Cartesian and Cylindrical geometries. The programming language is Microsoft Visual Basic 3.0. The code makes use of Finite element formulation for spatial domain and Finite difference formulation for time domain. Presently the code is capable of solving two dimensional steady state and transient problems in xy- and rz-geometries. The code is capable excepting both triangular and rectangular elements. Validation and benchmarking was done against hand calculations and published results. (author)

Global solubility of the three-dimensional Navier-Stokes equations with uniformly large initial vorticity

International Nuclear Information System (INIS)

Makhalov, A S; Nikolaenko, V P

2003-01-01

This paper is a survey of results concerning the three-dimensional Navier-Stokes and Euler equations with initial data characterized by uniformly large vorticity. The existence of regular solutions of the three-dimensional Navier-Stokes equations on an unbounded time interval is proved for large initial data both in R 3 and in bounded cylindrical domains. Moreover, the existence of smooth solutions on large finite time intervals is established for the three-dimensional Euler equations. These results are obtained without additional assumptions on the behaviour of solutions for t>0. Any smooth solution is not close to any two-dimensional manifold. Our approach is based on the computation of singular limits of rapidly oscillating operators, non-linear averaging, and a consideration of the mutual absorption of non-linear oscillations of the vorticity field. The use of resonance conditions, methods from the theory of small divisors, and non-linear averaging of almost periodic functions leads to the limit resonant Navier-Stokes equations. Global solubility of these equations is proved without any conditions on the three-dimensional initial data. The global regularity of weak solutions of three-dimensional Navier-Stokes equations with uniformly large vorticity at t=0 is proved by using the regularity of weak solutions and the strong convergence

A parallel algorithm for the two-dimensional time fractional diffusion equation with implicit difference method.

Science.gov (United States)

Gong, Chunye; Bao, Weimin; Tang, Guojian; Jiang, Yuewen; Liu, Jie

2014-01-01

It is very time consuming to solve fractional differential equations. The computational complexity of two-dimensional fractional differential equation (2D-TFDE) with iterative implicit finite difference method is O(M(x)M(y)N(2)). In this paper, we present a parallel algorithm for 2D-TFDE and give an in-depth discussion about this algorithm. A task distribution model and data layout with virtual boundary are designed for this parallel algorithm. The experimental results show that the parallel algorithm compares well with the exact solution. The parallel algorithm on single Intel Xeon X5540 CPU runs 3.16-4.17 times faster than the serial algorithm on single CPU core. The parallel efficiency of 81 processes is up to 88.24% compared with 9 processes on a distributed memory cluster system. We do think that the parallel computing technology will become a very basic method for the computational intensive fractional applications in the near future.

On the Sodium Concentration Diffusion with Three-Dimensional Extracellular Stimulation

Directory of Open Access Journals (Sweden)

Luisa Consiglieri

2011-01-01

Full Text Available We deal with the transmembrane sodium diffusion in a nerve. We study a mathematical model of a nerve fibre in response to an imposed extracellular stimulus. The presented model is constituted by a diffusion-drift vectorial equation in a bidomain, that is, two parabolic equations defined in each of the intra- and extra-regions. This system of partial differential equations can be understood as a reduced three-dimensional Poisson-Nernst-Planck model of the sodium concentration. The representation of the membrane includes a jump boundary condition describing the mechanisms involved in the excitation-contraction couple. Our first novelty comes from this general dynamical boundary condition. The second one is the three-dimensional behaviour of the extracellular stimulus. An analytical solution to the mathematical model is proposed depending on the morphology of the excitation.

Spherical harmonics solutions of multi-dimensional neutron transport equation by finite Fourier transformation

International Nuclear Information System (INIS)

Kobayashi, Keisuke

1977-01-01

A method of solution of a monoenergetic neutron transport equation in P sub(L) approximation is presented for x-y and x-y-z geometries using the finite Fourier transformation. A reactor system is assumed to consist of multiregions in each of which the nuclear cross sections are spatially constant. Since the unknown functions of this method are the spherical harmonics components of the neutron angular flux at the material boundaries alone, the three- and two-dimensional equations are reduced to two- and one-dimensional equations, respectively. The present approach therefore gives fewer unknowns than in the usual series expansion method or in the finite difference method. Some numerical examples are shown for the criticality problem. (auth.)

Hodograph solutions of the dispersionless coupled KdV hierarchies, critical points and the Euler-Poisson-Darboux equation

International Nuclear Information System (INIS)

Konopelchenko, B; Alonso, L MartInez; Medina, E

2010-01-01

It is shown that the hodograph solutions of the dispersionless coupled KdV (dcKdV) hierarchies describe critical and degenerate critical points of a scalar function which obeys the Euler-Poisson-Darboux equation. Singular sectors of each dcKdV hierarchy are found to be described by solutions of higher genus dcKdV hierarchies. Concrete solutions exhibiting shock-type singularities are presented.

Resolvent approach for two-dimensional scattering problems. Application to the nonstationary Schroedinger problem and the KPI equation

International Nuclear Information System (INIS)

Boiti, M.; Pempinelli, F.; Pogrebkov, A.K.; Polivanov, M.C.

1993-01-01

The resolvent operator of the linear problem is determined as the full Green function continued in the complex domain in two variables. An analog of the known Hilbert identity is derived. The authors demonstrate the role of this identity in the study of two-dimensional scattering. Considering the nonstationary Schroedinger equation as an example, it is shown that all types of solutions of the linear problem, as well as spectral data known in the literature, are given as specific values of this unique function - the resolvent function. A new form of the inverse problem is formulated. 7 refs

Lagrangian coherent structures at the onset of hyperchaos in the two-dimensional Navier-Stokes equations.

Science.gov (United States)

Miranda, Rodrigo A; Rempel, Erico L; Chian, Abraham C-L; Seehafer, Norbert; Toledo, Benjamin A; Muñoz, Pablo R

2013-09-01

We study a transition to hyperchaos in the two-dimensional incompressible Navier-Stokes equations with periodic boundary conditions and an external forcing term. Bifurcation diagrams are constructed by varying the Reynolds number, and a transition to hyperchaos (HC) is identified. Before the onset of HC, there is coexistence of two chaotic attractors and a hyperchaotic saddle. After the transition to HC, the two chaotic attractors merge with the hyperchaotic saddle, generating random switching between chaos and hyperchaos, which is responsible for intermittent bursts in the time series of energy and enstrophy. The chaotic mixing properties of the flow are characterized by detecting Lagrangian coherent structures. After the transition to HC, the flow displays complex Lagrangian patterns and an increase in the level of Lagrangian chaoticity during the bursty periods that can be predicted statistically by the hyperchaotic saddle prior to HC transition.

One-Dimensional Fokker-Planck Equation with Quadratically Nonlinear Quasilocal Drift

Science.gov (United States)

Shapovalov, A. V.

2018-04-01

The Fokker-Planck equation in one-dimensional spacetime with quadratically nonlinear nonlocal drift in the quasilocal approximation is reduced with the help of scaling of the coordinates and time to a partial differential equation with a third derivative in the spatial variable. Determining equations for the symmetries of the reduced equation are derived and the Lie symmetries are found. A group invariant solution having the form of a traveling wave is found. Within the framework of Adomian's iterative method, the first iterations of an approximate solution of the Cauchy problem are obtained. Two illustrative examples of exact solutions are found.

Multisoliton formula for completely integrable two-dimensional systems

International Nuclear Information System (INIS)

Chudnovsky, D.V.; Chudnovsky, G.V.

1979-01-01

For general two-dimensional completely integrable systems, the exact formulae for multisoliton type solutions are given. The formulae are obtained algebrically from solutions of two linear partial differential equations

Prescription-induced jump distributions in multiplicative Poisson processes.

Science.gov (United States)

Suweis, Samir; Porporato, Amilcare; Rinaldo, Andrea; Maritan, Amos

2011-06-01

Generalized Langevin equations (GLE) with multiplicative white Poisson noise pose the usual prescription dilemma leading to different evolution equations (master equations) for the probability distribution. Contrary to the case of multiplicative Gaussian white noise, the Stratonovich prescription does not correspond to the well-known midpoint (or any other intermediate) prescription. By introducing an inertial term in the GLE, we show that the Itô and Stratonovich prescriptions naturally arise depending on two time scales, one induced by the inertial term and the other determined by the jump event. We also show that, when the multiplicative noise is linear in the random variable, one prescription can be made equivalent to the other by a suitable transformation in the jump probability distribution. We apply these results to a recently proposed stochastic model describing the dynamics of primary soil salinization, in which the salt mass balance within the soil root zone requires the analysis of different prescriptions arising from the resulting stochastic differential equation forced by multiplicative white Poisson noise, the features of which are tailored to the characters of the daily precipitation. A method is finally suggested to infer the most appropriate prescription from the data.

Prescription-induced jump distributions in multiplicative Poisson processes

Science.gov (United States)

Suweis, Samir; Porporato, Amilcare; Rinaldo, Andrea; Maritan, Amos

2011-06-01

Generalized Langevin equations (GLE) with multiplicative white Poisson noise pose the usual prescription dilemma leading to different evolution equations (master equations) for the probability distribution. Contrary to the case of multiplicative Gaussian white noise, the Stratonovich prescription does not correspond to the well-known midpoint (or any other intermediate) prescription. By introducing an inertial term in the GLE, we show that the Itô and Stratonovich prescriptions naturally arise depending on two time scales, one induced by the inertial term and the other determined by the jump event. We also show that, when the multiplicative noise is linear in the random variable, one prescription can be made equivalent to the other by a suitable transformation in the jump probability distribution. We apply these results to a recently proposed stochastic model describing the dynamics of primary soil salinization, in which the salt mass balance within the soil root zone requires the analysis of different prescriptions arising from the resulting stochastic differential equation forced by multiplicative white Poisson noise, the features of which are tailored to the characters of the daily precipitation. A method is finally suggested to infer the most appropriate prescription from the data.

A modified Poisson-Boltmann model including charge regulation for the adsorption of ionizable polyelectrolytes to charged interfaces, applied to lysozyme adsorption on silica

NARCIS (Netherlands)

Biesheuvel, P.M.; Veen, van der M.; Norde, W.

2005-01-01

The equilibrium adsorption of polyelectrolytes with multiple types of ionizable groups is described using a modified Poisson-Boltzmann equation including charge regulation of both the polymer and the interface. A one-dimensional mean-field model is used in which the electrostatic potential is

Analytical approach to (2+1)-dimensional Boussinesq equation and (3+1)-dimensional Kadomtsev-Petviashvili equation

Energy Technology Data Exchange (ETDEWEB)

Sariaydin, Selin; Yildirim, Ahmet [Ege Univ., Dept. of Mathematics, Bornova-Izmir (Turkey)

2010-05-15

In this paper, we studied the solitary wave solutions of the (2+1)-dimensional Boussinesq equation u{sub tt} - u{sub xx} - u{sub yy} - (u{sup 2}){sub xx} - u{sub xxxx} = 0 and the (3+1)-dimensional Kadomtsev-Petviashvili (KP) equation u{sub xt} - 6u{sub x}{sup 2} + 6uu{sub xx} - u{sub xxxx} - u{sub yy} - u{sub zz} = 0. By using this method, an explicit numerical solution is calculated in the form of a convergent power series with easily computable components. To illustrate the application of this method numerical results are derived by using the calculated components of the homotopy perturbation series. The numerical solutions are compared with the known analytical solutions. Results derived from our method are shown graphically. (orig.)

An accessible four-dimensional treatment of Maxwell's equations in terms of differential forms

International Nuclear Information System (INIS)

Sá, Lucas

2017-01-01

Maxwell’s equations are derived in terms of differential forms in the four-dimensional Minkowski representation, starting from the three-dimensional vector calculus differential version of these equations. Introducing all the mathematical and physical concepts needed (including the tool of differential forms), using only knowledge of elementary vector calculus and the local vector version of Maxwell’s equations, the equations are reduced to a simple and elegant set of two equations for a unified quantity, the electromagnetic field. The treatment should be accessible for students taking a first course on electromagnetism. (paper)

Advanced diffusion model in compacted bentonite based on modified Poisson-Boltzmann equations

International Nuclear Information System (INIS)

Yotsuji, K.; Tachi, Y.; Nishimaki, Y.

2012-01-01

Document available in extended abstract form only. Diffusion and sorption of radionuclides in compacted bentonite are the key processes in the safe geological disposal of radioactive waste. JAEA has developed the integrated sorption and diffusion (ISD) model for compacted bentonite by coupling the pore water chemistry, sorption and diffusion processes in consistent way. The diffusion model accounts consistently for cation excess and anion exclusion in narrow pores in compacted bentonite by the electric double layer (EDL) theory. The firstly developed ISD model could predict the diffusivity of the monovalent cation/anion in compacted bentonite as a function of dry density. This ISD model was modified by considering the visco-electric effect, and applied for diffusion data for various radionuclides measured under wide range of conditions (salinity, density, etc.). This modified ISD model can give better quantitative agreement with diffusion data for monovalent cation/anion, however, the model predictions still disagree with experimental data for multivalent cation and complex species. In this study we extract the additional key factors influencing diffusion model in narrow charged pores, and the effects of these factors were investigated to reach a better understanding of diffusion processes in compacted bentonite. We investigated here the dielectric saturation effect and the excluded volume effect into the present ISD model and numerically solved these modified Poisson-Boltzmann equations. In the vicinity of the negatively charged clay surfaces, it is necessary to evaluate concentration distribution of electrolytes considering the dielectric saturation effects. The Poisson-Boltzmann (P-B) equation coupled with the dielectric saturation effects was solved numerically by using Runge-Kutta and Shooting methods. Figure 1(a) shows the concentration distributions of Na + as numerical solutions of the modified and original P-B equations for 0.01 M pore water, 800 kg m -3

Exact numerical solutions of the Schrödinger equation for a two-dimensional exciton in a constant magnetic field of arbitrary strength

Energy Technology Data Exchange (ETDEWEB)

Hoang-Do, Ngoc-Tram [Department of Physics, Ho Chi Minh City University of Pedagogy 280, An Duong Vuong Street, District 5, Ho Chi Minh City (Viet Nam); Pham, Dang-Lan [Institute for Computational Science and Technology, Quang Trung Software Town, District 12, Ho Chi Minh City (Viet Nam); Le, Van-Hoang, E-mail: hoanglv@hcmup.edu.vn [Department of Physics, Ho Chi Minh City University of Pedagogy 280, An Duong Vuong Street, District 5, Ho Chi Minh City (Viet Nam)

2013-08-15

Exact numerical solutions of the Schrödinger equation for a two-dimensional exciton in a constant magnetic field of arbitrary strength are obtained for not only the ground state but also high excited states. Toward this goal, the operator method is developed by combining with the Levi-Civita transformation which transforms the problem under investigation into that of a two-dimensional anharmonic oscillator. This development of the non-perturbation method is significant because it can be applied to other problems of two-dimensional atomic systems. The obtained energies and wave functions set a new record for their precision of up to 20 decimal places. Analyzing the obtained data we also find an interesting result that exact analytical solutions exist at some values of magnetic field intensity.

Exact numerical solutions of the Schrödinger equation for a two-dimensional exciton in a constant magnetic field of arbitrary strength

International Nuclear Information System (INIS)

Hoang-Do, Ngoc-Tram; Pham, Dang-Lan; Le, Van-Hoang

2013-01-01

Exact numerical solutions of the Schrödinger equation for a two-dimensional exciton in a constant magnetic field of arbitrary strength are obtained for not only the ground state but also high excited states. Toward this goal, the operator method is developed by combining with the Levi-Civita transformation which transforms the problem under investigation into that of a two-dimensional anharmonic oscillator. This development of the non-perturbation method is significant because it can be applied to other problems of two-dimensional atomic systems. The obtained energies and wave functions set a new record for their precision of up to 20 decimal places. Analyzing the obtained data we also find an interesting result that exact analytical solutions exist at some values of magnetic field intensity

Numerical Simulations of Scattering of Light from Two-Dimensional Rough Surfaces Using the Reduced Rayleigh Equation

Directory of Open Access Journals (Sweden)

Tor eNordam

2013-09-01

Full Text Available A formalism is introduced for the non-perturbative, purely numerical, solution of the reduced Rayleigh equation for the scattering of light from two-dimensional penetrable rough surfaces. Implementation and performance issues of the method, and various consistency checks of it, are presented and discussed. The proposed method is found, within the validity of the Rayleigh hypothesis, to give reliable results. For a non-absorbing metal surface the conservation of energy was explicitly checked, and found to be satisfied to within 0.03%, or better, for the parameters assumed. This testifies to the accuracy of the approach and a satisfactory discretization. As an illustration, we calculate the full angular distribution of the mean differential reflection coefficient for the scattering of p- or s-polarized light incident on two-dimensional dielectric or metallic randomly rough surfaces defined by (isotropic or anisotropic Gaussian and cylindrical power spectra. Simulation results obtained by the proposed method agree well with experimentally measured scattering data taken from similar well-characterized, rough metal samples, or to results obtained by other numerical methods.

The one-dimensional Gross-Pitaevskii equation and its some excitation states

Energy Technology Data Exchange (ETDEWEB)

Prayitno, T. B., E-mail: trunk-002@yahoo.com [Physics Department, Faculty of Mathematics and Natural Science, Universitas Negeri Jakarta, Jl. Pemuda Rawamangun no. 10, Jakarta, 13220 (Indonesia)

2015-04-16

We have derived some excitation states of the one-dimensional Gross-Pitaevskii equation coupled by the gravitational potential. The methods that we have used here are taken by pursuing the recent work of Kivshar et. al. by considering the equation as a macroscopic quantum oscillator. To obtain the states, we have made the appropriate transformation to reduce the three-dimensional Gross-Pitaevskii equation into the one-dimensional Gross-Pitaevskii equation and applying the time-independent perturbation theory in the general solution of the one-dimensional Gross-Pitaevskii equation as a linear superposition of the normalized eigenfunctions of the Schrödinger equation for the harmonic oscillator potential. Moreover, we also impose the condition by assuming that some terms in the equation should be so small in order to preserve the use of the perturbation method.

TURBO: a computer program for two-dimensional incompressible fluid flow analysis using a two-equations turbulence model

International Nuclear Information System (INIS)

Botelho, D.A.; Moreira, M.L.

1991-06-01

The Reynolds turbulent transport equations for an incompressible fluid are integrated on a bi-dimensional staggered grid, for velocity and pressure, using the SIMPLER method. With the resulting algebraic relations it was developed the TURBO program, which final objectives are the thermal stratification and natural convection analysis of nuclear reactor pools. This program was tested in problems applications with analytic or experimental solutions previously known. (author)

Solution of the equations of motion for a super non-Abelian sigma model in curved background by the super Poisson-Lie T-duality

International Nuclear Information System (INIS)

Eghbali, Ali

2015-01-01

The equations of motion of a super non-Abelian T-dual sigma model on the Lie supergroup (C_1"1+A) in the curved background are explicitly solved by the super Poisson-Lie T-duality. To find the solution of the flat model we use the transformation of supercoordinates, transforming the metric into a constant one, which is shown to be a supercanonical transformation. Then, using the super Poisson-Lie T-duality transformations and the dual decomposition of elements of Drinfel’d superdouble, the solution of the equations of motion for the dual sigma model is obtained. The general form of the dilaton fields satisfying the vanishing β−function equations of the sigma models is found. In this respect, conformal invariance of the sigma models built on the Drinfel’d superdouble ((C_1"1+A) , I_(_2_|_2_)) is guaranteed up to one-loop, at least.

A closed-form solution for the two-dimensional transport equation by the LTS{sub N} nodal method in the energy range of Compton effect

Energy Technology Data Exchange (ETDEWEB)

Rodriguez, B.D.A., E-mail: barbararodriguez@furg.b [Universidade Federal do Rio Grande, Instituto de Matematica, Estatistica e Fisica, Rio Grande, RS (Brazil); Vilhena, M.T., E-mail: vilhena@mat.ufrgs.b [Universidade Federal do Rio Grande do Sul, Departamento de Matematica Pura e Aplicada, Porto Alegre, RS (Brazil); Hoff, G., E-mail: hoff@pucrs.b [Pontificia Universidade Catolica do Rio Grande do Sul, Faculdade de Fisica, Porto Alegre, RS (Brazil); Bodmann, B.E.J., E-mail: bardo.bodmann@ufrgs.b [Universidade Federal do Rio Grande do Sul, Departamento de Matematica Pura e Aplicada, Porto Alegre, RS (Brazil)

2011-01-15

In the present work we report on a closed-form solution for the two-dimensional Compton transport equation by the LTS{sub N} nodal method in the energy range of Compton effect. The solution is determined using the LTS{sub N} nodal approach for homogeneous and heterogeneous rectangular domains, assuming the Klein-Nishina scattering kernel and a multi-group model. The solution is obtained by two one-dimensional S{sub N} equation systems resulting from integrating out one of the orthogonal variables of the S{sub N} equations in the rectangular domain. The leakage angular fluxes are approximated by exponential forms, which allows to determine a closed-form solution for the photons transport equation. The angular flux and the parameters of the medium are used for the calculation of the absorbed energy in rectangular domains with different dimensions and compositions. In this study, only the absorbed energy by Compton effect is considered. We present numerical simulations and comparisons with results obtained by using the simulation platform GEANT4 (version 9.1) with its low energy libraries.

Solving the two-dimensional Schrödinger equation using basis ...

Indian Academy of Sciences (India)

Ihab H Naeim

2017-10-19

Oct 19, 2017 ... We shall study the case of a two-dimensional quantum system .... Solving (6) for ck,l is tantamount to pro- ... case, the final computational problem becomes quite ..... matrix approach fails in the case of two particles con-.

N-dimensional integrability from two-photon coalgebra symmetry

International Nuclear Information System (INIS)

Ballesteros, Angel; Blasco, Alfonso; Herranz, Francisco J

2009-01-01

A wide class of Hamiltonian systems with N degrees of freedom and endowed with, at least, (N - 2) functionally independent integrals of motion in involution is constructed by making use of the two-photon Lie-Poisson coalgebra (h 6 , Δ). The set of (N - 2) constants of the motion is shown to be a universal one for all these Hamiltonians, irrespective of the dependence of the latter on several arbitrary functions and N free parameters. Within this large class of quasi-integrable N-dimensional Hamiltonians, new families of completely integrable systems are identified by finding explicitly a new independent integral I through the analysis of the sub-coalgebra structure of h 6 . In particular, new completely integrable N-dimensional Hamiltonians describing natural systems, geodesic flows and static electromagnetic Hamiltonians are presented

Fractal electrodynamics via non-integer dimensional space approach

Science.gov (United States)

Tarasov, Vasily E.

2015-09-01

Using the recently suggested vector calculus for non-integer dimensional space, we consider electrodynamics problems in isotropic case. This calculus allows us to describe fractal media in the framework of continuum models with non-integer dimensional space. We consider electric and magnetic fields of fractal media with charges and currents in the framework of continuum models with non-integer dimensional spaces. An application of the fractal Gauss's law, the fractal Ampere's circuital law, the fractal Poisson equation for electric potential, and equation for fractal stream of charges are suggested. Lorentz invariance and speed of light in fractal electrodynamics are discussed. An expression for effective refractive index of non-integer dimensional space is suggested.

Relativistic two-and three-particle scattering equations using instant and light-front dynamics

International Nuclear Information System (INIS)

Adhikari, S.K.; Tomio, L.; Frederico, T.

1992-01-01

Starting from the Bethe-Salpeter equation for two particles in the ladder approximation and integrating over the time component of momentum we derive three dimensional scattering integral equations satisfying constraints of unitarity and relativity, both employing the light-front and instant-form variables. The equations we arrive at are those first derived by Weinberg and by Blankenbecler and Sugar, and are shown to be related by a transformation of variables. Hence we show how to perform and relate identical dynamical calculation using these two equations. We extends this procedure to the case of three particles interacting via two-particle separable potentials. Using light-front and instant form variables we suggest a couple of three dimensional three-particle scattering equations satisfying constraints of two and three-particle unitarity and relativity. The three-particle light-front equation is shown to be approximately related by a transformation of variables to one of the instant-form three-particle equations. (author)

Central-Upwind Schemes for Two-Layer Shallow Water Equations

KAUST Repository

Kurganov, Alexander; Petrova, Guergana

2009-01-01

We derive a second-order semidiscrete central-upwind scheme for one- and two-dimensional systems of two-layer shallow water equations. We prove that the presented scheme is well-balanced in the sense that stationary steady-state solutions

Computational time analysis of the numerical solution of 3D electrostatic Poisson's equation

Science.gov (United States)

Kamboh, Shakeel Ahmed; Labadin, Jane; Rigit, Andrew Ragai Henri; Ling, Tech Chaw; Amur, Khuda Bux; Chaudhary, Muhammad Tayyab

2015-05-01

3D Poisson's equation is solved numerically to simulate the electric potential in a prototype design of electrohydrodynamic (EHD) ion-drag micropump. Finite difference method (FDM) is employed to discretize the governing equation. The system of linear equations resulting from FDM is solved iteratively by using the sequential Jacobi (SJ) and sequential Gauss-Seidel (SGS) methods, simulation results are also compared to examine the difference between the results. The main objective was to analyze the computational time required by both the methods with respect to different grid sizes and parallelize the Jacobi method to reduce the computational time. In common, the SGS method is faster than the SJ method but the data parallelism of Jacobi method may produce good speedup over SGS method. In this study, the feasibility of using parallel Jacobi (PJ) method is attempted in relation to SGS method. MATLAB Parallel/Distributed computing environment is used and a parallel code for SJ method is implemented. It was found that for small grid size the SGS method remains dominant over SJ method and PJ method while for large grid size both the sequential methods may take nearly too much processing time to converge. Yet, the PJ method reduces computational time to some extent for large grid sizes.

Two-dimensional effects in nonlinear Kronig-Penney models

DEFF Research Database (Denmark)

Gaididei, Yuri Borisovich; Christiansen, Peter Leth; Rasmussen, Kim

1997-01-01

An analysis of two-dimensional (2D) effects in the nonlinear Kronig-Penney model is presented. We establish an effective one-dimensional description of the 2D effects, resulting in a set of pseudodifferential equations. The stationary states of the 2D system and their stability is studied...

Two-Dimensional Motions of Rockets

Science.gov (United States)

Kang, Yoonhwan; Bae, Saebyok

2007-01-01

We analyse the two-dimensional motions of the rockets for various types of rocket thrusts, the air friction and the gravitation by using a suitable representation of the rocket equation and the numerical calculation. The slope shapes of the rocket trajectories are discussed for the three types of rocket engines. Unlike the projectile motions, the…

Topological aspect of disclinations in two-dimensional crystals

International Nuclear Information System (INIS)

Wei-Kai, Qi; Tao, Zhu; Yong, Chen; Ji-Rong, Ren

2009-01-01

By using topological current theory, this paper studies the inner topological structure of disclinations during the melting of two-dimensional systems. From two-dimensional elasticity theory, it finds that there are topological currents for topological defects in homogeneous equation. The evolution of disclinations is studied, and the branch conditions for generating, annihilating, crossing, splitting and merging of disclinations are given. (the physics of elementary particles and fields)

Exact solutions of the vacuum Einstein's equations allowing for two noncommuting Killing vectors

International Nuclear Information System (INIS)

Aliev, V.N.; Leznov, A.N.

1990-01-01

Einstein's equations are written in the form of covariant gauge theory in two-dimensional space with binomial solvable gauge group, with respect to two noncommutative of Killing vectors. The theory is exact integrable in one-dimensional case and series of partial exact solutions are constructed in two-dimensional. 5 refs

Spectral analysis and multigrid preconditioners for two-dimensional space-fractional diffusion equations

Science.gov (United States)

Moghaderi, Hamid; Dehghan, Mehdi; Donatelli, Marco; Mazza, Mariarosa

2017-12-01

Fractional diffusion equations (FDEs) are a mathematical tool used for describing some special diffusion phenomena arising in many different applications like porous media and computational finance. In this paper, we focus on a two-dimensional space-FDE problem discretized by means of a second order finite difference scheme obtained as combination of the Crank-Nicolson scheme and the so-called weighted and shifted Grünwald formula. By fully exploiting the Toeplitz-like structure of the resulting linear system, we provide a detailed spectral analysis of the coefficient matrix at each time step, both in the case of constant and variable diffusion coefficients. Such a spectral analysis has a very crucial role, since it can be used for designing fast and robust iterative solvers. In particular, we employ the obtained spectral information to define a Galerkin multigrid method based on the classical linear interpolation as grid transfer operator and damped-Jacobi as smoother, and to prove the linear convergence rate of the corresponding two-grid method. The theoretical analysis suggests that the proposed grid transfer operator is strong enough for working also with the V-cycle method and the geometric multigrid. On this basis, we introduce two computationally favourable variants of the proposed multigrid method and we use them as preconditioners for Krylov methods. Several numerical results confirm that the resulting preconditioning strategies still keep a linear convergence rate.

RBSDE's with jumps and the related obstacle problems for integral-partial differential equations

Institute of Scientific and Technical Information of China (English)

FAN; Yulian

2006-01-01

The author proves, when the noise is driven by a Brownian motion and an independent Poisson random measure, the one-dimensional reflected backward stochastic differential equation with a stopping time terminal has a unique solution. And in a Markovian framework, the solution can provide a probabilistic interpretation for the obstacle problem for the integral-partial differential equation.

Stable multiple-layer stationary solutions of a semilinear parabolic equation in two-dimensional domains

Directory of Open Access Journals (Sweden)

Arnaldo Simal do Nascimento

1997-12-01

Full Text Available We use $Gamma$--convergence to prove existence of stable multiple--layer stationary solutions (stable patterns to the reaction--diffusion equation. $$ eqalign{ {partial v_varepsilon over partial t} =& varepsilon^2, hbox{div}, (k_1(xabla v_varepsilon + k_2(x(v_varepsilon -alpha(Beta-v_varepsilon (v_varepsilon -gamma_varepsilon(x,,hbox{ in }Omegaimes{Bbb R}^+ cr &v_varepsilon(x,0 = v_0 quad {partial v_varepsilon over partial widehat{n}} = 0,, quadhbox{ for } xin partialOmega,, t >0,.} $$ Given nested simple closed curves in ${Bbb R}^2$, we give sufficient conditions on their curvature so that the reaction--diffusion problem possesses a family of stable patterns. In particular, we extend to two-dimensional domains and to a spatially inhomogeneous source term, a previous result by Yanagida and Miyata.

Alternative Forms of Compound Fractional Poisson Processes

Directory of Open Access Journals (Sweden)

Luisa Beghin

2012-01-01

Full Text Available We study here different fractional versions of the compound Poisson process. The fractionality is introduced in the counting process representing the number of jumps as well as in the density of the jumps themselves. The corresponding distributions are obtained explicitly and proved to be solution of fractional equations of order less than one. Only in the final case treated in this paper, where the number of jumps is given by the fractional-difference Poisson process defined in Orsingher and Polito (2012, we have a fractional driving equation, with respect to the time argument, with order greater than one. Moreover, in this case, the compound Poisson process is Markovian and this is also true for the corresponding limiting process. All the processes considered here are proved to be compositions of continuous time random walks with stable processes (or inverse stable subordinators. These subordinating relationships hold, not only in the limit, but also in the finite domain. In some cases the densities satisfy master equations which are the fractional analogues of the well-known Kolmogorov one.

Interaction phenomenon to dimensionally reduced p-gBKP equation

Science.gov (United States)

Zhang, Runfa; Bilige, Sudao; Bai, Yuexing; Lü, Jianqing; Gao, Xiaoqing

2018-02-01

Based on searching the combining of quadratic function and exponential (or hyperbolic cosine) function from the Hirota bilinear form of the dimensionally reduced p-gBKP equation, eight class of interaction solutions are derived via symbolic computation with Mathematica. The submergence phenomenon, presented to illustrate the dynamical features concerning these obtained solutions, is observed by three-dimensional plots and density plots with particular choices of the involved parameters between the exponential (or hyperbolic cosine) function and the quadratic function. It is proved that the interference between the two solitary waves is inelastic.

The BRST complex of homological Poisson reduction

Science.gov (United States)

Müller-Lennert, Martin

2017-02-01

BRST complexes are differential graded Poisson algebras. They are associated with a coisotropic ideal J of a Poisson algebra P and provide a description of the Poisson algebra (P/J)^J as their cohomology in degree zero. Using the notion of stable equivalence introduced in Felder and Kazhdan (Contemporary Mathematics 610, Perspectives in representation theory, 2014), we prove that any two BRST complexes associated with the same coisotropic ideal are quasi-isomorphic in the case P = R[V] where V is a finite-dimensional symplectic vector space and the bracket on P is induced by the symplectic structure on V. As a corollary, the cohomology of the BRST complexes is canonically associated with the coisotropic ideal J in the symplectic case. We do not require any regularity assumptions on the constraints generating the ideal J. We finally quantize the BRST complex rigorously in the presence of infinitely many ghost variables and discuss the uniqueness of the quantization procedure.

Positioning in a flat two-dimensional space-time: The delay master equation

International Nuclear Information System (INIS)

Coll, Bartolome; Ferrando, Joan Josep; Morales-Lladosa, Juan Antonio

2010-01-01

The basic theory on relativistic positioning systems in a two-dimensional space-time has been presented in two previous papers [B. Coll, J. J. Ferrando, and J. A. Morales, Phys. Rev. D 73, 084017 (2006); ibid.74, 104003 (2006)], where the possibility of making relativistic gravimetry with these systems has been analyzed by considering specific examples. Here, generic relativistic positioning systems in the Minkowski plane are studied. The information that can be obtained from the data received by a user of the positioning system is analyzed in detail. In particular, it is shown that the accelerations of the emitters and of the user along their trajectories are determined by the sole knowledge of the emitter positioning data and of the acceleration of only one of the emitters. Moreover, as a consequence of the so-called master delay equation, the knowledge of this acceleration is only required during an echo interval, i.e., the interval between the emission time of a signal by an emitter and its reception time after being reflected by the other emitter. These results are illustrated with the obtention of the dynamics of the emitters and of the user from specific sets of data received by the user.

Classical solutions of two dimensional Stokes problems on non smooth domains. 2: Collocation method for the Radon equation

International Nuclear Information System (INIS)

Lubuma, M.S.

1991-05-01

The non uniquely solvable Radon boundary integral equation for the two-dimensional Stokes-Dirichlet problem on a non smooth domain is transformed into a well posed one by a suitable compact perturbation of the velocity double layer potential operator. The solution to the modified equation is decomposed into a regular part and a finite linear combination of intrinsic singular functions whose coefficients are computed from explicit formulae. Using these formulae, the classical collocation method, defined by continuous piecewise linear vector-valued basis functions, which converges slowly because of the lack of regularity of the solution, is improved into a collocation dual singular function method with optimal rates of convergence for the solution and for the coefficients of singularities. (author). 34 refs

Wedge-Shaped GaN Nanowalls: A Potential Candidate for Two-Dimensional Electronics and Spintronics

Science.gov (United States)

Deb, Swarup; Dhar, Subhabrata

Schrödingerand Poisson equations are solved self-consistently in order to obtain the potential and charge density distribution in n-type GaN nanowalls tapered along c-axis by different angles. The study shows two-dimensional (2D) quantum confinement of electrons in the central vertical plane of the wall for the entire range of tapering. Calculation of room temperature electron mobility in the 2D channel shows a steady decrease with the increase of the inclination angle of the side facets with respect to the base. However, it is interesting to note that the mobility remains to be much larger than that of bulk GaN even for the inclination angle of 65∘. The properties of high mobility and the vertical orientation of the 2DEG plane in this system can be exploited in fabricating highly conducting transparent interconnects and field effect transistors, which can lead to large scale integration of 2D devices in future.

Comparing two Poisson populations sequentially: an application

International Nuclear Information System (INIS)

Halteman, E.J.

1986-01-01

Rocky Flats Plant in Golden, Colorado monitors each of its employees for radiation exposure. Excess exposure is detected by comparing the means of two Poisson populations. A sequential probability ratio test (SPRT) is proposed as a replacement for the fixed sample normal approximation test. A uniformly most efficient SPRT exists, however logistics suggest using a truncated SPRT. The truncated SPRT is evaluated in detail and shown to possess large potential savings in average time spent by employees in the monitoring process

On the late phase of relaxation of two-dimensional fluids: turbulence of unitons

International Nuclear Information System (INIS)

Spineanu, F; Vlad, M

2017-01-01

The two-dimensional ideal fluid and the plasma confined by a strong magnetic field exhibit an intrinsic tendency to organization due to the inverse spectral cascade. In the asymptotic states reached at relaxation the turbulence has vanished and there are only coherent vortical structures. We are interested in the regime that precedes these ordered flow patterns, in which there still is turbulence and imperfect but robust structures have emerged. To develop an analytical description we propose to start from the stationary coherent states and (in the direction opposite to relaxation) explore the space of configurations before the extremum of the functional that defines the structures has been reached. We find necessary to assemble different but related models: point-like vortices, its field theoretical formulation as interacting matter and gauge fields, chiral model and surfaces with constant mean curvature. These models are connected by the similar ability to described randomly interacting coherent structures. They derive exactly the same equation for the asymptotic state (sinh-Poisson equation, confirmed by numerical calculation of fluid flows). The chiral model, to which one can arrive from self-duality equation of the field theoretical model for fluid and from constant mean curvature surface equations, appears to be the suitable analytical framework. Its solutions, the unitons, aquire dynamics when the system is not at the extremum of the action. In the present work we provide arguments that the underlying common nature of these models can be used to develop an approach to fluid and plasma states of turbulence interacting with structures. (paper)

Five-dimensional truncation of the plane incompressible navier-stokes equations

Energy Technology Data Exchange (ETDEWEB)

Boldrighini, C [Camerino Univ. (Italy). Istituto di Matematica; Franceschini, V [Modena Univ. (Italy). Istituto Matematico

1979-01-01

A five-modes truncation of the Navier-Stokes equations for a two dimensional incompressible fluid on a torus is considered. A computer analysis shows that for a certain range of the Reynolds number the system exhibits a stochastic behaviour, approached through an involved sequence of bifurcations.

Two-level schemes for the advection equation

Science.gov (United States)

Vabishchevich, Petr N.

2018-06-01

The advection equation is the basis for mathematical models of continuum mechanics. In the approximate solution of nonstationary problems it is necessary to inherit main properties of the conservatism and monotonicity of the solution. In this paper, the advection equation is written in the symmetric form, where the advection operator is the half-sum of advection operators in conservative (divergent) and non-conservative (characteristic) forms. The advection operator is skew-symmetric. Standard finite element approximations in space are used. The standard explicit two-level scheme for the advection equation is absolutely unstable. New conditionally stable regularized schemes are constructed, on the basis of the general theory of stability (well-posedness) of operator-difference schemes, the stability conditions of the explicit Lax-Wendroff scheme are established. Unconditionally stable and conservative schemes are implicit schemes of the second (Crank-Nicolson scheme) and fourth order. The conditionally stable implicit Lax-Wendroff scheme is constructed. The accuracy of the investigated explicit and implicit two-level schemes for an approximate solution of the advection equation is illustrated by the numerical results of a model two-dimensional problem.

Application of Green's differential equation to the analysis of ion-matrix sheaths around wedge-shaped cathodes

International Nuclear Information System (INIS)

Donolato, C

2005-01-01

A relation between the gradient of the electric field and mean curvature of equipotential surfaces (Green's differential equation) is applied to a two-dimensional free-boundary problem arising in the study of ion sheaths around wedge-shaped cathodes. With the assumption that the equipotential lines are hyperbolae, this relation leads to a nonlinear ordinary differential equation for the potential along the bisector line of the wedge. An approximate solution is found, which yields, in particular, the sheath width along this line as a function of the wedge angle. The resulting values are in good agreement with published results obtained by numerically solving Poisson's equation

Accuracy assessment of the linear Poisson-Boltzmann equation and reparametrization of the OBC generalized Born model for nucleic acids and nucleic acid-protein complexes.

Science.gov (United States)

Fogolari, Federico; Corazza, Alessandra; Esposito, Gennaro

2015-04-05

The generalized Born model in the Onufriev, Bashford, and Case (Onufriev et al., Proteins: Struct Funct Genet 2004, 55, 383) implementation has emerged as one of the best compromises between accuracy and speed of computation. For simulations of nucleic acids, however, a number of issues should be addressed: (1) the generalized Born model is based on a linear model and the linearization of the reference Poisson-Boltmann equation may be questioned for highly charged systems as nucleic acids; (2) although much attention has been given to potentials, solvation forces could be much less sensitive to linearization than the potentials; and (3) the accuracy of the Onufriev-Bashford-Case (OBC) model for nucleic acids depends on fine tuning of parameters. Here, we show that the linearization of the Poisson Boltzmann equation has mild effects on computed forces, and that with optimal choice of the OBC model parameters, solvation forces, essential for molecular dynamics simulations, agree well with those computed using the reference Poisson-Boltzmann model. © 2015 Wiley Periodicals, Inc.

Self-focusing instability of two-dimensional solitons and vortices

DEFF Research Database (Denmark)

Kuznetsov, E.A.; Juul Rasmussen, J.

1995-01-01

The instability of two-dimensional solitons and vortices is demonstrated in the framework of the three-dimensional nonlinear Schrodinger equation (NLSE). The instability can be regarded as the analog of the Kadomtsev-Petviashvili instability [B. B. Kadomtsev and V. I. Petviashvili, Sov. Phys. Dokl...

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

International Nuclear Information System (INIS)

Rasulova, M.Yu

1998-01-01

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

Numerical treatment for solving two-dimensional space-fractional advection-dispersion equation using meshless method

Science.gov (United States)

Cheng, Rongjun; Sun, Fengxin; Wei, Qi; Wang, Jufeng

2018-02-01

Space-fractional advection-dispersion equation (SFADE) can describe particle transport in a variety of fields more accurately than the classical models of integer-order derivative. Because of nonlocal property of integro-differential operator of space-fractional derivative, it is very challenging to deal with fractional model, and few have been reported in the literature. In this paper, a numerical analysis of the two-dimensional SFADE is carried out by the element-free Galerkin (EFG) method. The trial functions for the SFADE are constructed by the moving least-square (MLS) approximation. By the Galerkin weak form, the energy functional is formulated. Employing the energy functional minimization procedure, the final algebraic equations system is obtained. The Riemann-Liouville operator is discretized by the Grünwald formula. With center difference method, EFG method and Grünwald formula, the fully discrete approximation schemes for SFADE are established. Comparing with exact results and available results by other well-known methods, the computed approximate solutions are presented in the format of tables and graphs. The presented results demonstrate the validity, efficiency and accuracy of the proposed techniques. Furthermore, the error is computed and the proposed method has reasonable convergence rates in spatial and temporal discretizations.

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

New exact solutions of (2 + 1)-dimensional Gardner equation via the new sine-Gordon equation expansion method

International Nuclear Information System (INIS)

Chen Yong; Yan Zhenya

2005-01-01

In this paper (2 + 1)-dimensional Gardner equation is investigated using a sine-Gordon equation expansion method, which was presented via a generalized sine-Gordon reduction equation and a new transformation. As a consequence, it is shown that the method is more powerful to obtain many types of new doubly periodic solutions of (2 + 1)-dimensional Gardner equation. In particular, solitary wave solutions are also given as simple limits of doubly periodic solutions

Collapse in a forced three-dimensional nonlinear Schrodinger equation

DEFF Research Database (Denmark)

Lushnikov, P.M.; Saffman, M.

2000-01-01

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

A finite-difference method for the variable coefficient Poisson equation on hierarchical Cartesian meshes

Science.gov (United States)

Raeli, Alice; Bergmann, Michel; Iollo, Angelo

2018-02-01

We consider problems governed by a linear elliptic equation with varying coefficients across internal interfaces. The solution and its normal derivative can undergo significant variations through these internal boundaries. We present a compact finite-difference scheme on a tree-based adaptive grid that can be efficiently solved using a natively parallel data structure. The main idea is to optimize the truncation error of the discretization scheme as a function of the local grid configuration to achieve second-order accuracy. Numerical illustrations are presented in two and three-dimensional configurations.

Flows of non-smooth vector fields and degenerate elliptic equations with applications to the Vlasov-Poisson and semigeostrophic systems

CERN Document Server

Colombo, Maria

2017-01-01

The first part of the book is devoted to the transport equation for a given vector field, exploiting the lagrangian structure of solutions. It also treats the regularity of solutions of some degenerate elliptic equations, which appear in the eulerian counterpart of some transport models with congestion. The second part of the book deals with the lagrangian structure of solutions of the Vlasov-Poisson system, which describes the evolution of a system of particles under the self-induced gravitational/electrostatic field, and the existence of solutions of the semigeostrophic system, used in meteorology to describe the motion of large-scale oceanic/atmospheric flows.

A generalization of the simplest equation method and its application to (3+1)-dimensional KP equation and generalized Fisher equation

International Nuclear Information System (INIS)

Zhao, Zhonglong; Zhang, Yufeng; Han, Zhong; Rui, Wenjuan

2014-01-01

In this paper, the simplest equation method is used to construct exact traveling solutions of the (3+1)-dimensional KP equation and generalized Fisher equation. We summarize the main steps of the simplest equation method. The Bernoulli and Riccati equation are used as simplest equations. This method is straightforward and concise, and it can be applied to other nonlinear partial differential equations

Invariants and labels for Lie-Poisson Systems

International Nuclear Information System (INIS)

Thiffeault, J.L.; Morrison, P.J.

1998-04-01

Reduction is a process that uses symmetry to lower the order of a Hamiltonian system. The new variables in the reduced picture are often not canonical: there are no clear variables representing positions and momenta, and the Poisson bracket obtained is not of the canonical type. Specifically, we give two examples that give rise to brackets of the noncanonical Lie-Poisson form: the rigid body and the two-dimensional ideal fluid. From these simple cases, we then use the semidirect product extension of algebras to describe more complex physical systems. The Casimir invariants in these systems are examined, and some are shown to be linked to the recovery of information about the configuration of the system. We discuss a case in which the extension is not a semidirect product, namely compressible reduced MHD, and find for this case that the Casimir invariants lend partial information about the configuration of the system

On the well-posedness of the stochastic Allen–Cahn equation in two dimensions

International Nuclear Information System (INIS)

Ryser, Marc D.; Nigam, Nilima; Tupper, Paul F.

2012-01-01

White noise-driven nonlinear stochastic partial differential equations (SPDEs) of parabolic type are frequently used to model physical systems in space dimensions d = 1, 2, 3. Whereas existence and uniqueness of weak solutions to these equations are well established in one dimension, the situation is different for d ⩾ 2. Despite their popularity in the applied sciences, higher dimensional versions of these SPDE models are generally assumed to be ill-posed by the mathematics community. We study this discrepancy on the specific example of the two dimensional Allen–Cahn equation driven by additive white noise. Since it is unclear how to define the notion of a weak solution to this equation, we regularize the noise and introduce a family of approximations. Based on heuristic arguments and numerical experiments, we conjecture that these approximations exhibit divergent behavior in the continuum limit. The results strongly suggest that shrinking the mesh size in simulations of the two-dimensional white noise-driven Allen–Cahn equation does not lead to the recovery of a physically meaningful limit.

Application of fast Fourier transforms to the direct solution of a class of two-dimensional separable elliptic equations on the sphere

Science.gov (United States)

Moorthi, Shrinivas; Higgins, R. W.

1993-01-01

An efficient, direct, second-order solver for the discrete solution of a class of two-dimensional separable elliptic equations on the sphere (which generally arise in implicit and semi-implicit atmospheric models) is presented. The method involves a Fourier transformation in longitude and a direct solution of the resulting coupled second-order finite-difference equations in latitude. The solver is made efficient by vectorizing over longitudinal wave-number and by using a vectorized fast Fourier transform routine. It is evaluated using a prescribed solution method and compared with a multigrid solver and the standard direct solver from FISHPAK.

Poisson-Boltzmann-Nernst-Planck model

International Nuclear Information System (INIS)

Zheng Qiong; Wei Guowei

2011-01-01

The Poisson-Nernst-Planck (PNP) model is based on a mean-field approximation of ion interactions and continuum descriptions of concentration and electrostatic potential. It provides qualitative explanation and increasingly quantitative predictions of experimental measurements for the ion transport problems in many areas such as semiconductor devices, nanofluidic systems, and biological systems, despite many limitations. While the PNP model gives a good prediction of the ion transport phenomenon for chemical, physical, and biological systems, the number of equations to be solved and the number of diffusion coefficient profiles to be determined for the calculation directly depend on the number of ion species in the system, since each ion species corresponds to one Nernst-Planck equation and one position-dependent diffusion coefficient profile. In a complex system with multiple ion species, the PNP can be computationally expensive and parameter demanding, as experimental measurements of diffusion coefficient profiles are generally quite limited for most confined regions such as ion channels, nanostructures and nanopores. We propose an alternative model to reduce number of Nernst-Planck equations to be solved in complex chemical and biological systems with multiple ion species by substituting Nernst-Planck equations with Boltzmann distributions of ion concentrations. As such, we solve the coupled Poisson-Boltzmann and Nernst-Planck (PBNP) equations, instead of the PNP equations. The proposed PBNP equations are derived from a total energy functional by using the variational principle. We design a number of computational techniques, including the Dirichlet to Neumann mapping, the matched interface and boundary, and relaxation based iterative procedure, to ensure efficient solution of the proposed PBNP equations. Two protein molecules, cytochrome c551 and Gramicidin A, are employed to validate the proposed model under a wide range of bulk ion concentrations and external

A generalized Poisson solver for first-principles device simulations

Energy Technology Data Exchange (ETDEWEB)

Bani-Hashemian, Mohammad Hossein; VandeVondele, Joost, E-mail: joost.vandevondele@mat.ethz.ch [Nanoscale Simulations, ETH Zürich, 8093 Zürich (Switzerland); Brück, Sascha; Luisier, Mathieu [Integrated Systems Laboratory, ETH Zürich, 8092 Zürich (Switzerland)

2016-01-28

Electronic structure calculations of atomistic systems based on density functional theory involve solving the Poisson equation. In this paper, we present a plane-wave based algorithm for solving the generalized Poisson equation subject to periodic or homogeneous Neumann conditions on the boundaries of the simulation cell and Dirichlet type conditions imposed at arbitrary subdomains. In this way, source, drain, and gate voltages can be imposed across atomistic models of electronic devices. Dirichlet conditions are enforced as constraints in a variational framework giving rise to a saddle point problem. The resulting system of equations is then solved using a stationary iterative method in which the generalized Poisson operator is preconditioned with the standard Laplace operator. The solver can make use of any sufficiently smooth function modelling the dielectric constant, including density dependent dielectric continuum models. For all the boundary conditions, consistent derivatives are available and molecular dynamics simulations can be performed. The convergence behaviour of the scheme is investigated and its capabilities are demonstrated.

Estimation of a Non-homogeneous Poisson Model: An Empirical ...

African Journals Online (AJOL)

This article aims at applying the Nonhomogeneous Poisson process to trends of economic development. For this purpose, a modified Nonhomogeneous Poisson process is derived when the intensity rate is considered as a solution of stochastic differential equation which satisfies the geometric Brownian motion. The mean ...

Resolvent approach for two-dimensional scattering problems. Application to the nonstationary Schrödinger problem and the KPI equation

Science.gov (United States)

Boiti, M.; Pempinelli, F.; Pogrebkov, A. K.; Polivanov, M. C.

1992-11-01

The resolvent operator of the linear problem is determined as the full Green function continued in the complex domain in two variables. An analog of the known Hilbert identity is derived. We demonstrate the role of this identity in the study of two-dimensional scattering. Considering the nonstationary Schrödinger equation as an example, we show that all types of solutions of the linear problems, as well as spectral data known in the literature, are given as specific values of this unique function — the resolvent function. A new form of the inverse problem is formulated.

On the conservation laws and solutions of a (2+1) dimensional KdV-mKdV equation of mathematical physics

Science.gov (United States)

Motsepa, Tanki; Masood Khalique, Chaudry

2018-05-01

In this paper, we study a (2+1) dimensional KdV-mKdV equation, which models many physical phenomena of mathematical physics. This equation has two integral terms in it. By an appropriate substitution, we convert this equation into two partial differential equations, which do not have integral terms and are equivalent to the original equation. We then work with the system of two equations and obtain its exact travelling wave solutions in form of Jacobi elliptic functions. Furthermore, we employ the multiplier method to construct conservation laws for the system. Finally, we revert the results obtained into the original variables of the (2+1) dimensional KdV-mKdV equation.

Unsteady free surface flow in porous media: One-dimensional model equations including vertical effects and seepage face

Science.gov (United States)

Di Nucci, Carmine

2018-05-01

This note examines the two-dimensional unsteady isothermal free surface flow of an incompressible fluid in a non-deformable, homogeneous, isotropic, and saturated porous medium (with zero recharge and neglecting capillary effects). Coupling a Boussinesq-type model for nonlinear water waves with Darcy's law, the two-dimensional flow problem is solved using one-dimensional model equations including vertical effects and seepage face. In order to take into account the seepage face development, the system equations (given by the continuity and momentum equations) are completed by an integral relation (deduced from the Cauchy theorem). After testing the model against data sets available in the literature, some numerical simulations, concerning the unsteady flow through a rectangular dam (with an impermeable horizontal bottom), are presented and discussed.

A one-level FETI method for the drift–diffusion-Poisson system with discontinuities at an interface

KAUST Repository

Baumgartner, Stefan

2013-06-01

A 3d feti method for the drift-diffusion-Poisson system including discontinuities at a 2d interface is developed. The motivation for this work is to provide a parallel numerical algorithm for a system of PDEs that are the basic model equations for the simulation of semiconductor devices such as transistors and sensors. Moreover, discontinuities or jumps in the potential and its normal derivative at a 2d surface are included for the simulation of nanowire sensors based on a homogenized model. Using the feti method, these jump conditions can be included with the usual numerical properties and the original Farhat-Roux feti method is extended to the drift-diffusion-Poisson equations including discontinuities. We show two numerical examples. The first example verifies the correct implementation including the discontinuities on a 2d grid divided into eight subdomains. The second example is 3d and shows the application of the algorithm to the simulation of nanowire sensors with high aspect ratios. The Poisson-Boltzmann equation and the drift-diffusion-Poisson system with jump conditions are solved on a 3d grid with real-world boundary conditions. © 2013 Elsevier Inc..

The coupling of Poisson sigma models to topological backgrounds

Energy Technology Data Exchange (ETDEWEB)

Rosa, Dario [School of Physics, Korea Institute for Advanced Study,Seoul 02455 (Korea, Republic of)

2016-12-13

We extend the coupling to the topological backgrounds, recently worked out for the 2-dimensional BF-model, to the most general Poisson sigma models. The coupling involves the choice of a Casimir function on the target manifold and modifies the BRST transformations. This in turn induces a change in the BRST cohomology of the resulting theory. The observables of the coupled theory are analyzed and their geometrical interpretation is given. We finally couple the theory to 2-dimensional topological gravity: this is the first step to study a topological string theory in propagation on a Poisson manifold. As an application, we show that the gauge-fixed vectorial supersymmetry of the Poisson sigma models has a natural explanation in terms of the theory coupled to topological gravity.

Non-linear instability analysis of the two-dimensional Navier-Stokes equation: The Taylor-Green vortex problem

Science.gov (United States)

Sengupta, Tapan K.; Sharma, Nidhi; Sengupta, Aditi

2018-05-01

An enstrophy-based non-linear instability analysis of the Navier-Stokes equation for two-dimensional (2D) flows is presented here, using the Taylor-Green vortex (TGV) problem as an example. This problem admits a time-dependent analytical solution as the base flow, whose instability is traced here. The numerical study of the evolution of the Taylor-Green vortices shows that the flow becomes turbulent, but an explanation for this transition has not been advanced so far. The deviation of the numerical solution from the analytical solution is studied here using a high accuracy compact scheme on a non-uniform grid (NUC6), with the fourth-order Runge-Kutta method. The stream function-vorticity (ψ, ω) formulation of the governing equations is solved here in a periodic square domain with four vortices at t = 0. Simulations performed at different Reynolds numbers reveal that numerical errors in computations induce a breakdown of symmetry and simultaneous fragmentation of vortices. It is shown that the actual physical instability is triggered by the growth of disturbances and is explained by the evolution of disturbance mechanical energy and enstrophy. The disturbance evolution equations have been traced by looking at (a) disturbance mechanical energy of the Navier-Stokes equation, as described in the work of Sengupta et al., "Vortex-induced instability of an incompressible wall-bounded shear layer," J. Fluid Mech. 493, 277-286 (2003), and (b) the creation of rotationality via the enstrophy transport equation in the work of Sengupta et al., "Diffusion in inhomogeneous flows: Unique equilibrium state in an internal flow," Comput. Fluids 88, 440-451 (2013).

Two-dimensional quantisation of the quasi-Landau hydrogenic spectrum

International Nuclear Information System (INIS)

Gallas, J.A.C.; O'Connell, R.F.

1982-01-01

Based on the two-dimensional WKB model, an equation is derived from which the non-relativistic quasi-Landau energy spectrum of hydrogen-like atoms may be easily obtained. In addition, the solution of radial equations in the WKB approximation and its relation with models recently used to fit experimental data are discussed. (author)

An analytical approach for a nodal scheme of two-dimensional neutron transport problems

International Nuclear Information System (INIS)

Barichello, L.B.; Cabrera, L.C.; Prolo Filho, J.F.

2011-01-01

Research highlights: → Nodal equations for a two-dimensional neutron transport problem. → Analytical Discrete Ordinates Method. → Numerical results compared with the literature. - Abstract: In this work, a solution for a two-dimensional neutron transport problem, in cartesian geometry, is proposed, on the basis of nodal schemes. In this context, one-dimensional equations are generated by an integration process of the multidimensional problem. Here, the integration is performed for the whole domain such that no iterative procedure between nodes is needed. The ADO method is used to develop analytical discrete ordinates solution for the one-dimensional integrated equations, such that final solutions are analytical in terms of the spatial variables. The ADO approach along with a level symmetric quadrature scheme, lead to a significant order reduction of the associated eigenvalues problems. Relations between the averaged fluxes and the unknown fluxes at the boundary are introduced as the usually needed, in nodal schemes, auxiliary equations. Numerical results are presented and compared with test problems.

Soliton solutions of the (2 + 1)-dimensional Harry Dym equation via Darboux transformation

International Nuclear Information System (INIS)

Halim, A.A.

2008-01-01

This work introduces solitons solutions for the (2 + 1)-dimensional Harry Dym equation using Darboux transformation. The link between the (2 + 1)-dimensional Harry Dym equation and the linear system associated with the modified Kadomtzev-Patvishvili equation is used. Namely, soliton solutions for the linear system associated with the later equation are produced using Darboux transformation. These solutions are inserted in the mentioned link to produce soliton solutions for the (2 + 1)-dimensional Harry Dym equation

Iterative Two- and One-Dimensional Methods for Three-Dimensional Neutron Diffusion Calculations

International Nuclear Information System (INIS)

Lee, Hyun Chul; Lee, Deokjung; Downar, Thomas J.

2005-01-01

Two methods are proposed for solving the three-dimensional neutron diffusion equation by iterating between solutions of the two-dimensional (2-D) radial and one-dimensional (1-D) axial solutions. In the first method, the 2-D/1-D equations are coupled using a current correction factor (CCF) with the average fluxes of the lower and upper planes and the axial net currents at the plane interfaces. In the second method, an analytic expression for the axial net currents at the interface of the planes is used for planar coupling. A comparison of the new methods is made with two previously proposed methods, which use interface net currents and partial currents for planar coupling. A Fourier convergence analysis of the four methods was performed, and results indicate that the two new methods have at least three advantages over the previous methods. First, the new methods are unconditionally stable, whereas the net current method diverges for small axial mesh size. Second, the new methods provide better convergence performance than the other methods in the range of practical mesh sizes. Third, the spectral radii of the new methods asymptotically approach zero as the mesh size increases, while the spectral radius of the partial current method approaches a nonzero value as the mesh size increases. Of the two new methods proposed here, the analytic method provides a smaller spectral radius than the CCF method, but the CCF method has several advantages over the analytic method in practical applications

Analysis of the gravitational coupled collisionless Boltzmann-poisson equations and numerical simulations of the formation of self-gravitating systems

International Nuclear Information System (INIS)

Roy, Fabrice

2004-01-01

We study the formation of self-gravitating systems and their properties by means of N-body simulations of gravitational collapse. First, we summarize the major analytical results concerning the collisionless Boltzmann equation and the Poisson's equation which describe the dynamics of collisionless gravitational systems. We present a study of some analytical solutions of this coupled system of equations. We then present the software used to perform the simulations. Some of this has been parallelized and implemented with the aid of MPI. For this reason we give a brief overview of it. Finally, we present the results of the numerical simulations. Analysis of these results allows us to explain some features of self-gravitating systems and the initial conditions needed to trigger the Antonov instability and the radial orbit instability. (author) [fr

Continuum model of the two-component Becker-Döring equations

OpenAIRE

Soheili, Ali Reza

2004-01-01

The process of collision between particles is a subject of interest in many fields of physics, astronomy, polymer physics, atmospheric physics, and colloid chemistry. If two types of particles are allowed to participate in the cluster coalescence, then the time evolution of the cluster distribution has been described by an infinite system of ordinary differential equations. In this paper, we describe the model with a second-order two-dimensional partial differential equation, as a continuum m...

Continuum model of the two-component Becker-Döring equations

Directory of Open Access Journals (Sweden)

Ali Reza Soheili

2004-01-01

Full Text Available The process of collision between particles is a subject of interest in many fields of physics, astronomy, polymer physics, atmospheric physics, and colloid chemistry. If two types of particles are allowed to participate in the cluster coalescence, then the time evolution of the cluster distribution has been described by an infinite system of ordinary differential equations. In this paper, we describe the model with a second-order two-dimensional partial differential equation, as a continuum model.

Dynamics of two-dimensional solitary vortices in a low-β plasma with convective motion

International Nuclear Information System (INIS)

Makino, Mitsuhiro; Kamimura, Tetsuo; Taniuti, Tosiya.

1980-12-01

Numerical studies of the Hasegawa-Mima equation, derived in the context of drift waves but equivalent to the quasigeostrophic vortex potential equation for Rossby waves, show the stable properties of solitary vortices which are two dimensional, localized, steady and translating solutions of this same equation. A solitary vortex can propagate only in the direction (x-direction) perpendicular to the density gradient. When this solitary vortex solution is inclined at some angle with respect to the x-axis, its propagation direction oscillates in the x and y plane. In two dimensional collisions, i.e. head-on collision and overtaking, solitary vortices interact two-dimensionally and recover their initial shapes at the end of both types of collisions. (author)

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

KAUST Repository

Bourantas, Georgios

2013-07-01

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

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

KAUST Repository

Bourantas, Georgios; Burganos, Vasilis N.

2013-01-01

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

Collapse arresting in an inhomogeneous two-dimensional nonlinear Schrodinger model

DEFF Research Database (Denmark)

Schjødt-Eriksen, Jens; Gaididei, Yuri Borisovich; Christiansen, Peter Leth

2001-01-01

Collapse of (2 + 1)-dimensional beams in the inhomogeneous two-dimensional cubic nonlinear Schrodinger equation is analyzed numerically and analytically. It is shown that in the vicinity of a narrow attractive inhomogeneity, the collapse of beams that in a homogeneous medium would collapse may...

Analytical simulation of two dimensional advection dispersion ...

African Journals Online (AJOL)

The study was designed to investigate the analytical simulation of two dimensional advection dispersion equation of contaminant transport. The steady state flow condition of the contaminant transport where inorganic contaminants in aqueous waste solutions are disposed of at the land surface where it would migrate ...

Analytical Simulation of Two Dimensional Advection Dispersion ...

African Journals Online (AJOL)

ADOWIE PERE

ABSTRACT: The study was designed to investigate the analytical simulation of two dimensional advection dispersion equation of contaminant transport. The steady state flow condition of the contaminant transport where inorganic contaminants in aqueous waste solutions are disposed of at the land surface where it would ...

The Penalty Cost Functional for the Two-Dimensional

Directory of Open Access Journals (Sweden)

Victor Onomza WAZIRI

2006-07-01

Full Text Available This paper constructs the penalty cost functional for optimizing the two-dimensional control operator of the energized wave equation. In some multiplier methods such as the Lagrange multipliers and Pontrygean maximum principle, the cost of merging the constraint equation to the integral quadratic objective functional to obtain an unconstraint equation is normally guessed or obtained from the first partial derivatives of the unconstrained equation. The Extended Conjugate Gradient Method (ECGM necessitates that the penalty cost be sequentially obtained algebraically. The ECGM problem contains a functional which is completely given in terms of state and time spatial dependent variables.

Remarks for one-dimensional fractional equations

Directory of Open Access Journals (Sweden)

Massimiliano Ferrara

2014-01-01

Full Text Available In this paper we study a class of one-dimensional Dirichlet boundary value problems involving the Caputo fractional derivatives. The existence of infinitely many solutions for this equations is obtained by exploiting a recent abstract result. Concrete examples of applications are presented.

Modifications to POISSON

International Nuclear Information System (INIS)

Harwood, L.H.

1981-01-01

At MSU we have used the POISSON family of programs extensively for magnetic field calculations. In the presently super-saturated computer situation, reducing the run time for the program is imperative. Thus, a series of modifications have been made to POISSON to speed up convergence. Two of the modifications aim at having the first guess solution as close as possible to the final solution. The other two aim at increasing the convergence rate. In this discussion, a working knowledge of POISSON is assumed. The amount of new code and expected time saving for each modification is discussed

Fractional calculus phenomenology in two-dimensional plasma models

Science.gov (United States)

Gustafson, Kyle; Del Castillo Negrete, Diego; Dorland, Bill

2006-10-01

Transport processes in confined plasmas for fusion experiments, such as ITER, are not well-understood at the basic level of fully nonlinear, three-dimensional kinetic physics. Turbulent transport is invoked to describe the observed levels in tokamaks, which are orders of magnitude greater than the theoretical predictions. Recent results show the ability of a non-diffusive transport model to describe numerical observations of turbulent transport. For example, resistive MHD modeling of tracer particle transport in pressure-gradient driven turbulence for a three-dimensional plasma reveals that the superdiffusive (2̂˜t^α where α> 1) radial transport in this system is described quantitatively by a fractional diffusion equation Fractional calculus is a generalization involving integro-differential operators, which naturally describe non-local behaviors. Our previous work showed the quantitative agreement of special fractional diffusion equation solutions with numerical tracer particle flows in time-dependent linearized dynamics of the Hasegawa-Mima equation (for poloidal transport in a two-dimensional cold-ion plasma). In pursuit of a fractional diffusion model for transport in a gyrokinetic plasma, we now present numerical results from tracer particle transport in the nonlinear Hasegawa-Mima equation and a planar gyrokinetic model. Finite Larmor radius effects will be discussed. D. del Castillo Negrete, et al, Phys. Rev. Lett. 94, 065003 (2005).

Formulation of Hamiltonian mechanics with even and odd Poisson brackets

International Nuclear Information System (INIS)

Khudaverdyan, O.M.; Nersesyan, A.P.

1987-01-01

A possibility is studied as to constrict the odd Poisson bracket and odd Hamiltonian by the given dynamics in phase superspace - the even Poisson bracket and even Hamiltonian so the transition to the new structure does not change the equations of motion. 9 refs

Numerical solutions of incompressible Navier-Stokes equations using modified Bernoulli's law

Science.gov (United States)

Shatalov, A.; Hafez, M.

2003-11-01

Simulations of incompressible flows are important for many practical applications in aeronautics and beyond, particularly in the high Reynolds number regime. The present formulation is based on Helmholtz velocity decomposition where the velocity is presented as the gradient of a potential plus a rotational component. Substituting in the continuity equation yields a Poisson equation for the potential which is solved with a zero normal derivative at solid surfaces. The momentum equation is used to update the rotational component with no slip/no penetration surface boundary conditions. The pressure is related to the potential function through a special relation which is a generalization of Bernoulli's law, with a viscous term included. Results of calculations for two- and three-dimensional problems prove that the present formulation is a valid approach, with some possible benefits compared to existing methods.

K-FIX: a computer program for transient, two-dimensional, two-fluid flow

International Nuclear Information System (INIS)

Rivard, W.C.; Torrey, M.D.

1976-11-01

The transient dynamics of two-dimensional, two-phase flow with interfacial exchange are calculated at all flow speeds using the K-FIX program. Each phase is described in terms of its own density, velocity, and temperature. The six field equations for the two phases couple through mass, momentum, and energy exchange. The equations are solved using an Eulerian finite difference technique that implicitly couples the rates of phase transitions, momentum, and energy exchange to determination of the pressure, density, and velocity fields. The implicit solution is accomplished iteratively without linearizing the equations, thus eliminating the need for numerous derivative terms. K-FIX is written in a highly modular form to be easily adaptable to a variety of problems. It is applied to growth of an isolated steam bubble in a superheated water pool

Zero singularities of codimension two and three in delay differential equations

International Nuclear Information System (INIS)

Campbell, Sue Ann; Yuan Yuan

2008-01-01

We give conditions under which a general class of delay differential equations has a point of Bogdanov–Takens or a triple zero bifurcation. We show how a centre manifold projection of the delay equations reduces the dynamics to two- or three-dimensional systems of ordinary differential equations. We put these equations in normal form and determine how the coefficients of the normal forms depend on the original parameters in the model. Finally we apply our results to two neural models and compare the predictions of the theory with numerical bifurcation analysis of the full equations. One model involves a transcritical bifurcation, hence we derive and analyse the appropriate unfoldings for this case

Exact solutions of the one-dimensional generalized modified complex Ginzburg-Landau equation

International Nuclear Information System (INIS)

Yomba, Emmanuel; Kofane, Timoleon Crepin

2003-01-01

The one-dimensional (1D) generalized modified complex Ginzburg-Landau (MCGL) equation for the traveling wave systems is analytically studied. Exact solutions of this equation are obtained using a method which combines the Painleve test for integrability in the formalism of Weiss-Tabor-Carnevale and Hirota technique of bilinearization. We show that pulses, fronts, periodic unbounded waves, sources, sinks and solution as collision between two fronts are the important coherent structures that organize much of the dynamical properties of these traveling wave systems. The degeneracies of the 1D generalized MCGL equation are examined as well as several of their solutions. These degeneracies include two important equations: the 1D generalized modified Schroedinger equation and the 1D generalized real modified Ginzburg-Landau equation. We obtain that the one parameter family of traveling localized source solutions called 'Nozaki-Bekki holes' become a subfamily of the dark soliton solutions in the 1D generalized modified Schroedinger limit

Stabilization of the solution of a two-dimensional system of Navier-Stokes equations in an unbounded domain with several exits to infinity

International Nuclear Information System (INIS)

Khisamutdinova, N A

2003-01-01

The behaviour as t→∞ of the solution of the mixed problem for the system of Navier-Stokes equations with a Dirichlet condition at the boundary is studied in an unbounded two-dimensional domain with several exits to infinity. A class of domains is distinguished in which an estimate characterizing the decay of solutions in terms of the geometry of the domain is proved for exponentially decreasing initial velocities. A similar estimate of the solution of the first mixed problem for the heat equation is sharp in a broad class of domains with several exits to infinity

Two-dimensional model of coupled heat and moisture transport in frost-heaving soils

International Nuclear Information System (INIS)

Guymon, G.L.; Berg, R.L.; Hromadka, T.V.

1984-01-01

A two-dimensional model of coupled heat and moisture flow in frost-heaving soils is developed based upon well known equations of heat and moisture flow in soils. Numerical solution is by the nodal domain integration method which includes the integrated finite difference and the Galerkin finite element methods. Solution of the phase change process is approximated by an isothermal approach and phenomenological equations are assumed for processes occurring in freezing or thawing zones. The model has been verified against experimental one-dimensional freezing soil column data and experimental two-dimensional soil thawing tank data as well as two-dimensional soil seepage data. The model has been applied to several simple but useful field problems such as roadway embankment freezing and frost heaving

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.

The connection of two-particle relativistic quantum mechanics with the Bethe-Salpeter equation

International Nuclear Information System (INIS)

Sazdjian, H.

1986-02-01

We show the formal equivalence between the wave equations of two-particle relativistic quantum mechanics, based on the manifestly covariant hamiltonian formalism with constraints, and the Bethe-Salpeter equation. This is achieved by algebraically transforming the latter so as to separate it into two independent equations which match the equations of hamiltonian relativistic quantum mechanics. The first equation determines the relative time evolution of the system, while the second one yields a three-dimensional eigenvalue equation. A connection is thus established between the Bethe-Salpeter wave function and its kernel on the one hand and the quantum mechanical wave function and interaction potential on the other. For the sector of solutions of the Bethe-Salpeter equation having non-relativistic limits, this relationship can be evaluated in perturbation theory. We also device a generalized form of the instantaneous approximation which simplifies the various expressions involved in the above relations. It also permits the evaluation of the normalization condition of the quantum mechanical wave function as a three-dimensional integral

Second invariant for two-dimensional classical super systems

Indian Academy of Sciences (India)

Construction of superpotentials for two-dimensional classical super systems (for N. 2) is carried ... extensively used for the case of non-linear partial differential equation by various authors. [3,4–7,12 ..... found to be integrable just by accident.

Kinetic Theory of a Confined Quasi-Two-Dimensional Gas of Hard Spheres

Directory of Open Access Journals (Sweden)

J. Javier Brey

2017-02-01

Full Text Available The dynamics of a system of hard spheres enclosed between two parallel plates separated a distance smaller than two particle diameters is described at the level of kinetic theory. The interest focuses on the behavior of the quasi-two-dimensional fluid seen when looking at the system from above or below. In the first part, a collisional model for the effective two-dimensional dynamics is analyzed. Although it is able to describe quite well the homogeneous evolution observed in the experiments, it is shown that it fails to predict the existence of non-equilibrium phase transitions, and in particular, the bimodal regime exhibited by the real system. A critical revision analysis of the model is presented , and as a starting point to get a more accurate description, the Boltzmann equation for the quasi-two-dimensional gas has been derived. In the elastic case, the solutions of the equation verify an H-theorem implying a monotonic tendency to a non-uniform steady state. As an example of application of the kinetic equation, here the evolution equations for the vertical and horizontal temperatures of the system are derived in the homogeneous approximation, and the results compared with molecular dynamics simulation results.

Exact results in the large system size limit for the dynamics of the chemical master equation, a one dimensional chain of equations.

Science.gov (United States)

Martirosyan, A; Saakian, David B

2011-08-01

We apply the Hamilton-Jacobi equation (HJE) formalism to solve the dynamics of the chemical master equation (CME). We found exact analytical expressions (in large system-size limit) for the probability distribution, including explicit expression for the dynamics of variance of distribution. We also give the solution for some simple cases of the model with time-dependent rates. We derived the results of the Van Kampen method from the HJE approach using a special ansatz. Using the Van Kampen method, we give a system of ordinary differential equations (ODEs) to define the variance in a two-dimensional case. We performed numerics for the CME with stationary noise. We give analytical criteria for the disappearance of bistability in the case of stationary noise in one-dimensional CMEs.

A three-dimensional viscous topography mesoscale model

Energy Technology Data Exchange (ETDEWEB)

Eichhorn, J; Flender, M; Kandlbinder, T; Panhans, W G; Trautmann, T; Zdunkowski, W G [Mainz Univ. (Germany). Inst. fuer Physik der Atmosphaere; Cui, K; Ries, R; Siebert, J; Wedi, N

1997-11-01

This study describes the theoretical foundation and applications of a newly designed mesoscale model named CLIMM (climate model Mainz). In contrast to terrain following coordinates, a cartesian grid is used to keep the finite difference equations as simple as possible. The method of viscous topography is applied to the flow part of the model. Since the topography intersects the cartesian grid cells, the new concept of boundary weight factors is introduced for the solution of Poisson`s equation. A three-dimensional radiosity model was implemented to handle radiative transfer at the ground. The model is applied to study thermally induced circulations and gravity waves at an idealized mountain. Furthermore, CLIMM was used to simulate typical wind and temperature distributions for the city of Mainz and its rural surroundings. It was found that the model in all cases produced realistic results. (orig.) 38 refs.

Electromagnetic-field equations in the six-dimensional space-time R6

International Nuclear Information System (INIS)

Teli, M.T.; Palaskar, D.

1984-01-01

Maxwell's equations (without monopoles) for electromagnetic fields are obtained in six-dimensional space-time. The equations possess structural symmetry in space and time, field and source densities. Space-time-symmetric conservation laws and field solutions are obtained. The results are successfully correlated with their four-dimensional space-time counterparts

Application of Exp-function method for (2 + 1)-dimensional nonlinear evolution equations