On the Solution of the Neutron Transport Equation

Energy Technology Data Exchange (ETDEWEB)

Depken, S

1962-12-15

The neutron transport equation has occupied the attention of many authors since Placzek, Wick and others made their first attempts to solve it, Even in the simple case of energy independent cross-sections, and disregarding the motion of the scattering nucleons, it is difficult to find a solution in an analytical form which is easily surveyable and fitted for numerical calculations. In Part I of this paper some new viewpoints will be introduced which enable the solution to be presented in its simplest possible form. Part II is devoted to an investigation of some functions introduced in Part I. In Part III the results are applied to the case of large energy lethargy, and the validity of derived formulas is discussed.

The solution of the multigroup neutron transport equation using spherical harmonics

International Nuclear Information System (INIS)

Fletcher, K.

1981-01-01

A solution of the multi-group neutron transport equation in up to three space dimensions is presented. The flux is expanded in a series of unnormalised spherical harmonics. Using the various recurrence formulae a linked set of first order differential equations is obtained for the moments psisup(g)sub(lm)(r), γsup(g)sub(lm)(r). Terms with odd l are eliminated resulting in a second order system which is solved by two methods. The first is a finite difference formulation using an iterative procedure, secondly, in XYZ and XY geometry a finite element solution is given. Results for a test problem using both methods are exhibited and compared. (orig./RW) [de

A numerical solution of the coupled proton-H atom transport equations for the proton aurora

International Nuclear Information System (INIS)

Basu, B.; Jasperse, J.R.; Grossbard, N.J.

1990-01-01

A numerical code has been developed to solve the coupled proton-H atom linear transport equations for the proton aurora. The transport equations have been simplified by using plane-parallel geometry and the forward-scattering approximations only. Otherwise, the equations and their numerical solutions are exact. Results are presented for the particle fluxes and the energy deposition rates, and they are compared with the previous analytical results that were obtained by using additional simplifying approximations. It is found that although the analytical solutions for the particle fluxes differ somewhat from the numerical solutions, the energy deposition rates calculated by the two methods agree to within a few percent. The accurate particle fluxes given by the numerical code are useful for accurate calculation of the characteristic quantities of the proton aurora, such as the ionization rates and the emission rates

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

Application of synthetic diffusion method in the numerical solution of the equations of neutron transport in slab geometry

International Nuclear Information System (INIS)

Valdes Parra, J.J.

1986-01-01

One of the main problems in reactor physics is to determine the neutron distribution in reactor core, since knowing that, it is possible to calculate the rapidity of occurrence of different nuclear reaction inside the reactor core. Within different theories existing in nuclear reactor physics, is neutron transport the one in which equation who govern the exact behavior of neutronic distribution are developed even inside the proper neutron transport theory, there exist different methods of solution which are approximations to exact solution; still more, with the purpose to reach a more precise solution, the majority of methods have been approached to the obtention of solutions in numerical form with the aim of take the advantages of modern computers, and for this reason a great deal of effort is dedicated to numerical solution of the equations of neutron transport. In agreement with the above mentioned, in this work has been developed a computer program which uses a relatively new techniques known as 'acceleration of synthetic diffusion' which has been applied to solve the neutron transport equation with 'classical schemes of spatial integration' obtaining results with a smaller quantity of interactions, if they compare to done without using such equation (Author)

Solution and study of nodal neutron transport equation applying the LTSN-DiagExp method

International Nuclear Information System (INIS)

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

2003-01-01

In this paper we report advances about the three-dimensional nodal discrete-ordinates approximations of neutron transport equation for Cartesian geometry. We use the combined collocation method of the angular variables and nodal approach for the spatial variables. By nodal approach we mean the iterated transverse integration of the S N equations. This procedure leads to the set of one-dimensional averages angular fluxes in each spatial variable. The resulting system of equations is solved with the LTS N method, first applying the Laplace transform to the set of the nodal S N equations and then obtained the solution by symbolic computation. We include the LTS N method by diagonalization to solve the nodal neutron transport equation and then we outline the convergence of these nodal-LTS N approximations with the help of a norm associated to the quadrature formula used to approximate the integral term of the neutron transport equation. (author)

Application of finite element method in the solution of transport equation

International Nuclear Information System (INIS)

Maiorino, J.R.; Vieira, W.J.

1985-01-01

It is presented the application of finite element method in the solution of second order transport equation (self-adjoint) for the even parity flux. The angular component is treated by expansion in Legendre polinomials uncoupled of the spatial component, which is approached by an expansion in base functions, interpolated in each spatial element. (M.C.K.) [pt

Analytical solutions of a fractional diffusion-advection equation for solar cosmic-ray transport

International Nuclear Information System (INIS)

Litvinenko, Yuri E.; Effenberger, Frederic

2014-01-01

Motivated by recent applications of superdiffusive transport models to shock-accelerated particle distributions in the heliosphere, we analytically solve a one-dimensional fractional diffusion-advection equation for the particle density. We derive an exact Fourier transform solution, simplify it in a weak diffusion approximation, and compare the new solution with previously available analytical results and with a semi-numerical solution based on a Fourier series expansion. We apply the results to the problem of describing the transport of energetic particles, accelerated at a traveling heliospheric shock. Our analysis shows that significant errors may result from assuming an infinite initial distance between the shock and the observer. We argue that the shock travel time should be a parameter of a realistic superdiffusive transport model.

Numerical Solution of the Electron Transport Equation in the Upper Atmosphere

Energy Technology Data Exchange (ETDEWEB)

Woods, Mark Christopher [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Holmes, Mark [Rensselaer Polytechnic Inst., Troy, NY (United States); Sailor, William C [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2017-07-01

A new approach for solving the electron transport equation in the upper atmosphere is derived. The problem is a very stiff boundary value problem, and to obtain an accurate numerical solution, matrix factorizations are used to decouple the fast and slow modes. A stable finite difference method is applied to each mode. This solver is applied to a simplifieed problem for which an exact solution exists using various versions of the boundary conditions that might arise in a natural auroral display. The numerical and exact solutions are found to agree with each other to at least two significant digits.

Simulation of transport equations with Monte Carlo

International Nuclear Information System (INIS)

Matthes, W.

1975-09-01

The main purpose of the report is to explain the relation between the transport equation and the Monte Carlo game used for its solution. The introduction of artificial particles carrying a weight provides one with high flexibility in constructing many different games for the solution of the same equation. This flexibility opens a way to construct a Monte Carlo game for the solution of the adjoint transport equation. Emphasis is laid mostly on giving a clear understanding of what to do and not on the details of how to do a specific game

The use of non-dimensional representation of the solute transport equations

International Nuclear Information System (INIS)

Laurens, J.-M.

1988-07-01

This report presents the results obtained in a pilot investigation into the use of non-dimensional representations of the solute transport equations, so as to improve the efficiency of the PRA codes used by the DoE and its contractors. A reduced set of parameters was obtained for a single layer transport model. As expected, the response was shown to be highly sensitive on the new parameters. A faster convergence of the system was observed when the sampling technique used was changed to take into account the properties of the new parameters, such that uniform coverage of the reduced parameter hyperspace was achieved. (author)

A mass conservative numerical solution of vertical water flow and mass transport equations in unsaturated porous media

International Nuclear Information System (INIS)

Lim, S.C.; Lee, K.J.

1993-01-01

The Galerkin finite element method is used to solve the problem of one-dimensional, vertical flow of water and mass transport of conservative-nonconservative solutes in unsaturated porous media. Numerical approximations based on different forms of the governing equation, although they are equivalent in continuous forms, can result in remarkably different solutions in an unsaturated flow problem. Solutions given by a simple Galerkin method based on the h-based Richards equation yield a large mass balance error and an underestimation of the infiltration depth. With the employment of the ROMV (restoration of main variable) concept in the discretization step, the mass conservative numerical solution algorithm for water flow has been derived. The resulting computational schemes for water flow and mass transport are applied to sandy soil. The ROMV method shows good mass conservation in water flow analysis, whereas it seems to have a minor effect on mass transport. However, it may relax the time-step size restriction and so ensure an improved calculation output. (author)

Solution and study of nodal neutron transport equation applying the LTS{sub N}-DiagExp method

Energy Technology Data Exchange (ETDEWEB)

Hauser, Eliete Biasotto; Pazos, Ruben Panta [Pontificia Univ. Catolica do Rio Grande do Sul, Porto Alegre, RS (Brazil). Faculdade de Matematica]. E-mail: eliete@pucrs.br; rpp@mat.pucrs.br; Vilhena, Marco Tullio de [Pontificia Univ. Catolica do Rio Grande do Sul, Porto Alegre, RS (Brazil). Instituto de Matematica]. E-mail: vilhena@mat.ufrgs.br; Barros, Ricardo Carvalho de [Universidade do Estado, Nova Friburgo, RJ (Brazil). Instituto Politecnico]. E-mail: ricardo@iprj.uerj.br

2003-07-01

In this paper we report advances about the three-dimensional nodal discrete-ordinates approximations of neutron transport equation for Cartesian geometry. We use the combined collocation method of the angular variables and nodal approach for the spatial variables. By nodal approach we mean the iterated transverse integration of the S{sub N} equations. This procedure leads to the set of one-dimensional averages angular fluxes in each spatial variable. The resulting system of equations is solved with the LTS{sub N} method, first applying the Laplace transform to the set of the nodal S{sub N} equations and then obtained the solution by symbolic computation. We include the LTS{sub N} method by diagonalization to solve the nodal neutron transport equation and then we outline the convergence of these nodal-LTS{sub N} approximations with the help of a norm associated to the quadrature formula used to approximate the integral term of the neutron transport equation. (author)

An inexact Newton method for fully-coupled solution of the Navier-Stokes equations with heat and mass transport

Energy Technology Data Exchange (ETDEWEB)

Shadid, J.N.; Tuminaro, R.S. [Sandia National Labs., Albuquerque, NM (United States); Walker, H.F. [Utah State Univ., Logan, UT (United States). Dept. of Mathematics and Statistics

1997-02-01

The solution of the governing steady transport equations for momentum, heat and mass transfer in flowing fluids can be very difficult. These difficulties arise from the nonlinear, coupled, nonsymmetric nature of the system of algebraic equations that results from spatial discretization of the PDEs. In this manuscript the authors focus on evaluating a proposed nonlinear solution method based on an inexact Newton method with backtracking. In this context they use a particular spatial discretization based on a pressure stabilized Petrov-Galerkin finite element formulation of the low Mach number Navier-Stokes equations with heat and mass transport. The discussion considers computational efficiency, robustness and some implementation issues related to the proposed nonlinear solution scheme. Computational results are presented for several challenging CFD benchmark problems as well as two large scale 3D flow simulations.

Presentation of some methods for the solution of the monoenergetic neutrons transport equation

International Nuclear Information System (INIS)

Valle G, E. del.

1978-01-01

The neutrons transport theory problems whose solution has been reached were collected in order to show that the transport equation is so complicated that different techniques were developed so as to give approximative numerical solutions to problems concerning the practical application. Such a technique, which had not been investigated in the literature dealing with these problems, is described here. The results which were obtained through this technique in undimensional problems of criticity are satisfactory and speaking in a conceptual way this method is extremely simple because it times. There is no limitation to deal with problems related neutrons sources with an arbitrary distribution and in principle the application of this technique can be extended to unhomogeneous environments. (author)

Determination of a closed-form solution for the multidimensional transport equation using a fractional derivative

International Nuclear Information System (INIS)

Zabadal, J.; Vilhena, M.T.; Segatto, C.F.; Pazos, R.P.Ruben Panta.

2002-01-01

In this work we construct a closed-form solution for the multidimensional transport equation rewritten in integral form which is expressed in terms of a fractional derivative of the angular flux. We determine the unknown order of the fractional derivative comparing the kernel of the integral equation with the one of the Riemann-Liouville definition of fractional derivative. We report numerical simulations

Determination of a closed-form solution for the multidimensional transport equation using a fractional derivative

Energy Technology Data Exchange (ETDEWEB)

Zabadal, J. E-mail: jorge.zabadal@ufrgs.br; Vilhena, M.T. E-mail: vilhena@mat.ufrgs.br; Segatto, C.F. E-mail: cynthia@mat.ufrgs.br; Pazos, R.P.Ruben Panta. E-mail: rpp@mat.pucrgs.br

2002-07-01

In this work we construct a closed-form solution for the multidimensional transport equation rewritten in integral form which is expressed in terms of a fractional derivative of the angular flux. We determine the unknown order of the fractional derivative comparing the kernel of the integral equation with the one of the Riemann-Liouville definition of fractional derivative. We report numerical simulations.

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

Approximate solution to neutron transport equation with linear anisotropic scattering

International Nuclear Information System (INIS)

Coppa, G.; Ravetto, P.; Sumini, M.

1983-01-01

A method to obtain an approximate solution to the transport equation, when both sources and collisions show a linearly anisotropic behavior, is outlined and the possible implications for numerical calculations in applied neutronics as well as shielding evaluations are investigated. The form of the differential system of equations taken by the method is quite handy and looks simpler and more manageable than any other today available technique. To go deeper into the efficiency of the method, some typical calculations concerning critical dimension of multiplying systems are then performed and the results are compared with the ones coming from the classical Ssub(N) approximations. The outcome of such calculations leads us to think of interesting developments of the method which could be quite useful in alternative to other today widespread approximate procedures, for any geometry, but especially for curved ones. (author)

Solution to the monoenergetic time-dependent neutron transport equation with a time-varying source

International Nuclear Information System (INIS)

Ganapol, B.D.

1986-01-01

Even though fundamental time-dependent neutron transport problems have existed since the inception of neutron transport theory, it has only been recently that a reliable numerical solution to one of the basic problems has been obtained. Experience in generating numerical solutions to time-dependent transport equations has indicated that the multiple collision formulation is the most versatile numerical technique for model problems. The formulation coupled with a moment reconstruction of each collided flux component has led to benchmark-quality (four- to five-digit accuracy) numerical evaluation of the neutron flux in plane infinite geometry for any degree of scattering anisotropy and for both pulsed isotropic and beam sources. As will be shown in this presentation, this solution can serve as a Green's function, thus extending the previous results to more complicated source situations. Here we will be concerned with a time-varying source at the center of an infinite medium. If accurate, such solutions have both pedagogical and practical uses as benchmarks against which other more approximate solutions designed for a wider class of problems can be compared

Exact analytical solution of time-independent neutron transport equation, and its applications to systems with a point source

International Nuclear Information System (INIS)

Mikata, Y.

2014-01-01

Highlights: • An exact solution for the one-speed neutron transport equation is obtained. • This solution as well as its derivation are believed to be new. • Neutron flux for a purely absorbing material with a point neutron source off the origin is obtained. • Spherically as well as cylindrically piecewise constant cross sections are studied. • Neutron flux expressions for a point neutron source off the origin are believed to be new. - Abstract: An exact analytical solution of the time-independent monoenergetic neutron transport equation is obtained in this paper. The solution is applied to systems with a point source. Systematic analysis of the solution of the time-independent neutron transport equation, and its applications represent the primary goal of this paper. To the best of the author’s knowledge, certain key results on the scalar neutron flux as well as their derivations are new. As an application of these results, a scalar neutron flux for a purely absorbing medium with a spherically piecewise constant cross section and an isotropic point neutron source off the origin as well as that for a cylindrically piecewise constant cross section with a point neutron source off the origin are obtained. Both of these results are believed to be new

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

A solution of the monoenergetic neutral particle transport equation for adjacent half-spaces with anisotropic scattering

Energy Technology Data Exchange (ETDEWEB)

Ganapol, B.D., E-mail: ganapol@cowboy.ame.arizona.edu [Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ (United States); Mostacci, D.; Previti, A. [Montecuccolino Laboratory, University of Bologna, Via dei Colli, 16, I-40136 Bologna (Italy)

2016-07-01

We present highly accurate solutions to the neutral particle transport equation in a half-space. While our initial motivation was in response to a recently published solution based on Chandrasekhar's H-function, the presentation to follow has taken on a more comprehensive tone. The solution by H-functions certainly did achieved high accuracy but was limited to isotropic scattering and emission from spatially uniform and linear sources. Moreover, the overly complicated nature of the H-function approach strongly suggests that its extension to anisotropic scattering and general sources is not at all practical. For this reason, an all encompassing theory for the determination of highly precise benchmarks, including anisotropic scattering for a variety of spatial source distributions, is presented for particle transport in a half-space. We illustrate the approach via a collection of cases including tables of 7-place flux benchmarks to guide transport methods developers. The solution presented can be applied to a considerable number of one and two half-space transport problems with variable sources and represents a state-of-the-art benchmark solution.

Numerical solution of neutron transport equations in discrete ordinates and slab geometry

International Nuclear Information System (INIS)

Serrano Pedraza, F.

1985-01-01

An unified formalism to solve numerically, between other equation, the neutron transport in discrete ordinates, slab geometry, several energy groups and independents of time, has been developed recently. Such a formalism cover some of the conventional schemes as diamond difference, (WDD) characteristic step (SC) lineal characteristic (LC), quadratic characteristic (QC) and lineal discontinuous. Unified formation gives before hand the convergence order of the previously selected scheme. In fact it allows besides to generate a big amount of numerical schemes, with which is also possible to solve numerical equations as soon as neutron transport. The essential purpose of this work was to solve the neutron transport equations in slab geometry and discrete ordinates considering several energy groups without to take under advisement time dependence based in the above mentioned unified formalism. To reach this purpose it was necesary to design a computer code with the name TNOD1 (Neutron transport in discrete ordinates and 1 dimension) which includes each one of the schemes already pointed out. there exist two numerical schemes, also recently developed, quadratic continuous (QC) and cubic continuous (CN), although covered by unified formalism, it has been possible to include them inside this computer code without make substantial changes in its structure. In chapter I, derivative of neutron transport equation independent of time is taken, for angular flux, including boundary conditions and discontinuity. In chapter II the neutron transport equations are obtained in multigroups, independents of time, for approximation of discrete ordinates. Description of theory related with unified formalism and its relationship with mentioned discretization schemes is presented in chapter III. Chapter IV describes the computer code developed and finally, in chapter V different numerical results obtained with TNOD1 program are shown. In Appendix A theorems and mathematical arguments used

A comparison of numerical solutions of partial differential equations with probabilistic and possibilistic parameters for the quantification of uncertainty in subsurface solute transport.

Science.gov (United States)

Zhang, Kejiang; Achari, Gopal; Li, Hua

2009-11-03

Traditionally, uncertainty in parameters are represented as probabilistic distributions and incorporated into groundwater flow and contaminant transport models. With the advent of newer uncertainty theories, it is now understood that stochastic methods cannot properly represent non random uncertainties. In the groundwater flow and contaminant transport equations, uncertainty in some parameters may be random, whereas those of others may be non random. The objective of this paper is to develop a fuzzy-stochastic partial differential equation (FSPDE) model to simulate conditions where both random and non random uncertainties are involved in groundwater flow and solute transport. Three potential solution techniques namely, (a) transforming a probability distribution to a possibility distribution (Method I) then a FSPDE becomes a fuzzy partial differential equation (FPDE), (b) transforming a possibility distribution to a probability distribution (Method II) and then a FSPDE becomes a stochastic partial differential equation (SPDE), and (c) the combination of Monte Carlo methods and FPDE solution techniques (Method III) are proposed and compared. The effects of these three methods on the predictive results are investigated by using two case studies. The results show that the predictions obtained from Method II is a specific case of that got from Method I. When an exact probabilistic result is needed, Method II is suggested. As the loss or gain of information during a probability-possibility (or vice versa) transformation cannot be quantified, their influences on the predictive results is not known. Thus, Method III should probably be preferred for risk assessments.

Saturation and linear transport equation

International Nuclear Information System (INIS)

Kutak, K.

2009-03-01

We show that the GBW saturation model provides an exact solution to the one dimensional linear transport equation. We also show that it is motivated by the BK equation considered in the saturated regime when the diffusion and the splitting term in the diffusive approximation are balanced by the nonlinear term. (orig.)

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.

An Equation-Type Approach for the Numerical Solution of the Partial Differential Equations Governing Transport Phenomena in Porous Media

KAUST Repository

Sun, Shuyu

2012-06-02

A new technique for the numerical solution of the partial differential equations governing transport phenomena in porous media is introduced. In this technique, the governing equations as depicted from the physics of the problem are used without extra manipulations. In other words, there is no need to reduce the number of governing equations by some sort of mathematical manipulations. This technique enables the separation of the physics part of the problem and the solver part, which makes coding more robust and could be used in several other applications with little or no modifications (e.g., multi-phase flow in porous media). In this method, one abandons the need to construct the coefficient matrix for the pressure equation. Alternatively, the coefficients are automatically generated within the solver routine. We show examples of using this technique to solving several flow problems in porous media.

An Equation-Type Approach for the Numerical Solution of the Partial Differential Equations Governing Transport Phenomena in Porous Media

KAUST Repository

Sun, Shuyu; Salama, Amgad; El-Amin, Mohamed

2012-01-01

A new technique for the numerical solution of the partial differential equations governing transport phenomena in porous media is introduced. In this technique, the governing equations as depicted from the physics of the problem are used without extra manipulations. In other words, there is no need to reduce the number of governing equations by some sort of mathematical manipulations. This technique enables the separation of the physics part of the problem and the solver part, which makes coding more robust and could be used in several other applications with little or no modifications (e.g., multi-phase flow in porous media). In this method, one abandons the need to construct the coefficient matrix for the pressure equation. Alternatively, the coefficients are automatically generated within the solver routine. We show examples of using this technique to solving several flow problems in porous media.

Global existence of weak solutions to dissipative transport equations with nonlocal velocity

Science.gov (United States)

Bae, Hantaek; Granero-Belinchón, Rafael; Lazar, Omar

2018-04-01

We consider 1D dissipative transport equations with nonlocal velocity field: where is a nonlocal operator given by a Fourier multiplier. We especially consider two types of nonlocal operators: (1) , the Hilbert transform, (2) . In this paper, we show several global existence of weak solutions depending on the range of γ, δ and α. When , we take initial data having finite energy, while we take initial data in weighted function spaces (in the real variables or in the Fourier variables), which have infinite energy, when .

The secret to successful solute-transport modeling

Science.gov (United States)

Konikow, Leonard F.

2011-01-01

Modeling subsurface solute transport is difficult—more so than modeling heads and flows. The classical governing equation does not always adequately represent what we see at the field scale. In such cases, commonly used numerical models are solving the wrong equation. Also, the transport equation is hyperbolic where advection is dominant, and parabolic where hydrodynamic dispersion is dominant. No single numerical method works well for all conditions, and for any given complex field problem, where seepage velocity is highly variable, no one method will be optimal everywhere. Although we normally expect a numerically accurate solution to the governing groundwater-flow equation, errors in concentrations from numerical dispersion and/or oscillations may be large in some cases. The accuracy and efficiency of the numerical solution to the solute-transport equation are more sensitive to the numerical method chosen than for typical groundwater-flow problems. However, numerical errors can be kept within acceptable limits if sufficient computational effort is expended. But impractically long

Correspondence Between One- and Two-Equation Models for Solute Transport in Two-Region Heterogeneous Porous Media

KAUST Repository

Davit, Y.; Wood, B. D.; Debenest, G.; Quintard, M.

2012-01-01

In this work, we study the transient behavior of homogenized models for solute transport in two-region porous media. We focus on the following three models: (1) a time non-local, two-equation model (2eq-nlt). This model does not rely on time

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

Exact harmonic solutions to Guyer-Krumhansl-type equation and application to heat transport in thin films

Science.gov (United States)

Zhukovsky, K.; Oskolkov, D.

2018-03-01

A system of hyperbolic-type inhomogeneous differential equations (DE) is considered for non-Fourier heat transfer in thin films. Exact harmonic solutions to Guyer-Krumhansl-type heat equation and to the system of inhomogeneous DE are obtained in Cauchy- and Dirichlet-type conditions. The contribution of the ballistic-type heat transport, of the Cattaneo heat waves and of the Fourier heat diffusion is discussed and compared with each other in various conditions. The application of the study to the ballistic heat transport in thin films is performed. Rapid evolution of the ballistic quasi-temperature component in low-dimensional systems is elucidated and compared with slow evolution of its diffusive counterpart. The effect of the ballistic quasi-temperature component on the evolution of the complete quasi-temperature is explored. In this context, the influence of the Knudsen number and of Cauchy- and Dirichlet-type conditions on the evolution of the temperature distribution is explored. The comparative analysis of the obtained solutions is performed.

Hot electrons in superlattices: quantum transport versus Boltzmann equation

DEFF Research Database (Denmark)

Wacker, Andreas; Jauho, Antti-Pekka; Rott, S.

1999-01-01

A self-consistent solution of the transport equation is presented for semiconductor superlattices within different approaches: (i) a full quantum transport model based on nonequilibrium Green functions, (ii) the semiclassical Boltzmann equation for electrons in a miniband, and (iii) Boltzmann...

On the spectral analysis of iterative solutions of the discretized one-group transport equation

International Nuclear Information System (INIS)

Sanchez, Richard

2004-01-01

We analyze the Fourier-mode technique used for the spectral analysis of iterative solutions of the one-group discretized transport equation. We introduce a direct spectral analysis for the iterative solution of finite difference approximations for finite slabs composed of identical layers, providing thus a complementary analysis that is more appropriate for reactor applications. Numerical calculations for the method of characteristics and with the diamond difference approximation show the appearance of antisymmetric modes generated by the iteration on boundary data. We have also utilized the discrete Fourier transform to compute the spectrum for a periodic slab containing N identical layers and shown that at the limit N → ∞ one obtains the familiar Fourier-mode solution

dispersion equation parameters of solute transport in agricultural

African Journals Online (AJOL)

Jane

2011-08-31

Aug 31, 2011 ... fields for predicting soil quality property. Key words: ... The classical approach of modeling solute transport in porous media uses the deterministic ... concentration of the solution in the liquid phase, u0 is the mean velocity and ...

Adjoint P1 equations solution for neutron slowing down

International Nuclear Information System (INIS)

Cardoso, Carlos Eduardo Santos; Martinez, Aquilino Senra; Silva, Fernando Carvalho da

2002-01-01

In some applications of perturbation theory, it is necessary know the adjoint neutron flux, which is obtained by the solution of adjoint neutron diffusion equation. However, the multigroup constants used for this are weighted in only the direct neutron flux, from the solution of direct P1 equations. In this work, the adjoint P1 equations are derived by the neutron transport equation, the reversion operators rules and analogies between direct and adjoint parameters. The direct and adjoint neutron fluxes resulting from the solution of P 1 equations were used to three different weighting processes, to obtain the macrogroup macroscopic cross sections. It was found out noticeable differences among them. (author)

Method of solution of the neutron transport equation in multidimensional cartesian geometries using spherical harmonics and spatially orthogonal polynomials

International Nuclear Information System (INIS)

Fenstermacher, T.E.

1981-01-01

The solution of the neutron transport equation has long been a subject of intense interest to nuclear engineers. Present computer codes for the solution of this equation, however, are expensive to run for large, multidimensional problems, and also suffer from computational problems such as the ray effect. A method has been developed which eliminates many of these problems. It consists of transforming the transport equation into a set of linear partial differential equations by the use of spherical harmonics. The problem volume is divided into mesh boxes, and the flux components are approximated within each mesh box by spatially orthogonal quadratic polynomials, which need not be continuous at mesh box interfaces. A variational principle is developed, and used to solve for the unknown coefficients of these polynomials. Both one dimensional and two dimensional computer codes using this method have been written. The codes have each been tested on several test cases, and the solutions checked against solutions obtained by other methods. While the codes have some difficulty in modeling sharp transients, they produce excellent results on problems where the characteristic lengths are many mean free paths. On one test case, the two dimensional code, SHOP/2D, required only one-fourth the computer time required by the finite difference, discrete ordinates code TWOTRAN to produce a solution. In addition, SHOP/2D converged much better than TWOTRAN and produced more physical-appearing results

Solution algorithms for a PN-1 - Equivalent SN angular discretization of the transport equation in one-dimensional spherical coordinates

International Nuclear Information System (INIS)

Warsa, J. S.; Morel, J. E.

2007-01-01

Angular discretizations of the S N transport equation in curvilinear coordinate systems may result in a streaming-plus-removal operator that is dense in the angular variable or that is not lower-triangular. We investigate numerical solution algorithms for such angular discretizations using relationships given by Chandrasekhar to compute the angular derivatives in the one-dimensional S N transport equation in spherical coordinates with Gauss quadrature. This discretization makes the S N transport equation P N-1 - equivalent, but it also makes the sweep operator dense at every spatial point because the N angular derivatives are expressed in terms of the N angular fluxes. To avoid having to invert the sweep operator directly, we must work with the angular fluxes to solve the equations iteratively. We show how we can use approximations to the sweep operator to precondition the full P N-1 equivalent S N equations. We show that these pre-conditioners affect the operator enough such that convergence of a Krylov iterative method improves. (authors)

Coupling between solute transport and chemical reactions models

International Nuclear Information System (INIS)

Samper, J.; Ajora, C.

1993-01-01

During subsurface transport, reactive solutes are subject to a variety of hydrodynamic and chemical processes. The major hydrodynamic processes include advection and convection, dispersion and diffusion. The key chemical processes are complexation including hydrolysis and acid-base reactions, dissolution-precipitation, reduction-oxidation, adsorption and ion exchange. The combined effects of all these processes on solute transport must satisfy the principle of conservation of mass. The statement of conservation of mass for N mobile species leads to N partial differential equations. Traditional solute transport models often incorporate the effects of hydrodynamic processes rigorously but oversimplify chemical interactions among aqueous species. Sophisticated chemical equilibrium models, on the other hand, incorporate a variety of chemical processes but generally assume no-flow systems. In the past decade, coupled models accounting for complex hydrological and chemical processes, with varying degrees of sophistication, have been developed. The existing models of reactive transport employ two basic sets of equations. The transport of solutes is described by a set of partial differential equations, and the chemical processes, under the assumption of equilibrium, are described by a set of nonlinear algebraic equations. An important consideration in any approach is the choice of primary dependent variables. Most existing models cannot account for the complete set of chemical processes, cannot be easily extended to include mixed chemical equilibria and kinetics, and cannot handle practical two and three dimensional problems. The difficulties arise mainly from improper selection of the primary variables in the transport equations. (Author) 38 refs

The plasma transport equations derived by multiple time-scale expansions and turbulent transport. I. General theory

International Nuclear Information System (INIS)

Edenstrasser, J.W.

1995-01-01

A multiple time-scale derivative expansion scheme is applied to the dimensionless Fokker--Planck equation and to Maxwell's equations, where the parameter range of a typical fusion plasma was assumed. Within kinetic theory, the four time scales considered are those of Larmor gyration, particle transit, collisions, and classical transport. The corresponding magnetohydrodynamic (MHD) time scales are those of ion Larmor gyration, Alfven, MHD collision, and resistive diffusion. The solution of the zeroth-order equations results in the force-free equilibria and ideal Ohm's law. The solution of the first-order equations leads under the assumption of a weak collisional plasma to the ideal MHD equations. On the MHD-collision time scale, not only the full set of the MHD transport equations is obtained, but also turbulent terms, where the related transport quantities are one order in the expansion parameter larger than those of classical transport. Finally, at the resistive diffusion time scale the known transport equations are arrived at including, however, also turbulent contributions. copyright 1995 American Institute of Physics

Modeling of amorphous pocket formation in silicon by numerical solution of the heat transport equation

International Nuclear Information System (INIS)

Kovac, D.; Otto, G.; Hobler, G.

2005-01-01

In this paper we present a model of amorphous pocket formation that is based on binary collision simulations to generate the distribution of deposited energy, and on numerical solution of the heat transport equation to describe the quenching process. The heat transport equation is modified to consider the heat of melting when the melting temperature is crossed at any point in space. It is discretized with finite differences on grid points that coincide with the crystallographic lattice sites, which allows easy determination of molten atoms. Atoms are considered molten if the average of their energy and the energy of their neighbors meets the melting criterion. The results obtained with this model are in good overall agreement with published experimental data on P, As, Te and Tl implantations in Si and with data on the polyatomic effect at cryogenic temperature

Finite-difference solution of the space-angle-lethargy-dependent slowing-down transport equation

Energy Technology Data Exchange (ETDEWEB)

Matausek, M V [Boris Kidric Vinca Institute of Nuclear Sciences, Vinca, Belgrade (Yugoslavia)

1972-07-01

A procedure has been developed for solving the slowing-down transport equation for a cylindrically symmetric reactor system. The anisotropy of the resonance neutron flux is treated by the spherical harmonics formalism, which reduces the space-angle-Iethargy-dependent transport equation to a matrix integro-differential equation in space and lethargy. Replacing further the lethargy transfer integral by a finite-difference form, a set of matrix ordinary differential equations is obtained, with lethargy-and space dependent coefficients. If the lethargy pivotal points are chosen dense enough so that the difference correction term can be ignored, this set assumes a lower block triangular form and can be solved directly by forward block substitution. As in each step of the finite-difference procedure a boundary value problem has to be solved for a non-homogeneous system of ordinary differential equations with space-dependent coefficients, application of any standard numerical procedure, for example, the finite-difference method or the method of adjoint equations, is too cumbersome and would make the whole procedure practically inapplicable. A simple and efficient approximation is proposed here, allowing analytical solution for the space dependence of the spherical-harmonics flux moments, and hence the derivation of the recurrence relations between the flux moments at successive lethargy pivotal points. According to the procedure indicated above a computer code has been developed for the CDC -3600 computer, which uses the KEDAK nuclear data file. The space and lethargy distribution of the resonance neutrons can be computed in such a detailed fashion as the neutron cross-sections are known for the reactor materials considered. The computing time is relatively short so that the code can be efficiently used, either autonomously, or as part of some complex modular scheme. Typical results will be presented and discussed in order to prove and illustrate the applicability of the

Solution of charged particle transport equation by Monte-Carlo method in the BRANDZ code system

International Nuclear Information System (INIS)

Artamonov, S.N.; Androsenko, P.A.; Androsenko, A.A.

1992-01-01

Consideration is given to the issues of Monte-Carlo employment for the solution of charged particle transport equation and its implementation in the BRANDZ code system under the conditions of real 3D geometry and all the data available on radiation-to-matter interaction in multicomponent and multilayer targets. For the solution of implantation problem the results of BRANDZ data comparison with the experiments and calculations by other codes in complexes systems are presented. The results of direct nuclear pumping process simulation for laser-active media by a proton beam are also included. 4 refs.; 7 figs

Two analytic transport equation solutions for particular cases of particle history

International Nuclear Information System (INIS)

Simovic, R.

1997-01-01

For anisotropic scattering and plane geometry, the linear transport equation of particles generated by a monodirectional unit source A(x,μ) = δ(x-0)δ(μ - μ 0 ) > 0, can be stated in the form of an integral equation

Global Solutions to the Coupled Chemotaxis-Fluid Equations

KAUST Repository

Duan, Renjun

2010-08-10

In this paper, we are concerned with a model arising from biology, which is a coupled system of the chemotaxis equations and the viscous incompressible fluid equations through transport and external forcing. The global existence of solutions to the Cauchy problem is investigated under certain conditions. Precisely, for the Chemotaxis-Navier-Stokes system over three space dimensions, we obtain global existence and rates of convergence on classical solutions near constant states. When the fluid motion is described by the simpler Stokes equations, we prove global existence of weak solutions in two space dimensions for cell density with finite mass, first-order spatial moment and entropy provided that the external forcing is weak or the substrate concentration is small. © Taylor & Francis Group, LLC.

Matrix-oriented implementation for the numerical solution of the partial differential equations governing flows and transport in porous media

KAUST Repository

Sun, Shuyu; Salama, Amgad; El-Amin, Mohamed

2012-01-01

In this paper we introduce a new technique for the numerical solution of the various partial differential equations governing flow and transport phenomena in porous media. This method is proposed to be used in high level programming languages like

A non-linear optimal Discontinuous Petrov-Galerkin method for stabilising the solution of the transport equation

International Nuclear Information System (INIS)

Merton, S. R.; Smedley-Stevenson, R. P.; Pain, C. C.; Buchan, A. G.; Eaton, M. D.

2009-01-01

This paper describes a new Non-Linear Discontinuous Petrov-Galerkin (NDPG) method and application to the one-speed Boltzmann Transport Equation (BTE) for space-time problems. The purpose of the method is to remove unwanted oscillations in the transport solution which occur in the vicinity of sharp flux gradients, while improving computational efficiency and numerical accuracy. This is achieved by applying artificial dissipation in the solution gradient direction, internal to an element using a novel finite element (FE) Riemann approach. The amount of dissipation added acts internal to each element. This is done using a gradient-informed scaling of the advection velocities in the stabilisation term. This makes the method in its most general form non-linear. The method is designed to be independent of angular expansion framework. This is demonstrated for the both discrete ordinates (S N ) and spherical harmonics (P N ) descriptions of the angular variable. Results show the scheme performs consistently well in demanding time dependent and multi-dimensional radiation transport problems. (authors)

Semianalytical Solutions of Radioactive or Reactive Transport in Variably-Fractured Layered Media: 1. Solutes

International Nuclear Information System (INIS)

George J. Moridis

2001-01-01

In this paper, semianalytical solutions are developed for the problem of transport of radioactive or reactive solute tracers through a layered system of heterogeneous fractured media with misaligned fractures. The tracer transport equations in the non-flowing matrix account for (a) diffusion, (b) surface diffusion, (c) mass transfer between the mobile and immobile water fractions, (d) linear kinetic or equilibrium physical, chemical, or combined solute sorption or colloid filtration, and (e) radioactive decay or first-order chemical reactions. The tracer-transport equations in the fractures account for the same processes, in addition to advection and hydrodynamic dispersion. Any number of radioactive decay daughter products (or products of a linear, first-order reaction chain) can be tracked. The solutions, which are analytical in the Laplace space, are numerically inverted to provide the solution in time and can accommodate any number of fractured and/or porous layers. The solutions are verified using analytical solutions for limiting cases of solute and colloid transport through fractured and porous media. The effect of important parameters on the transport of 3 H, 237 Np and 239 Pu (and its daughters) is investigated in several test problems involving layered geological systems of varying complexity

Quadrature with arbitrary weight for the numerical solution of the critical slab Neutron Transport Equation

International Nuclear Information System (INIS)

Sanchez G, J.

2007-01-01

A standard procedure for the solution of singular integral equations is applied to the one-dimensional transport equation for monoenergetic neutrons. The results obtained with two versions of the procedure, differing only in the extent of the basic region to which they are applied, are compared with analytically derived results available for benchmarking. The procedures considered yield consistent results for the calculated neutron densities and eigenvalues. Several approximate expressions of the neutron density are used to render closed-form formulas for the densities which can then be analytically operated on to obtain expressions for extrapolation distances or angular densities or serve other purposes that require an analytical expression of the neutron density. (Author)

A method for solving neutron transport equation

International Nuclear Information System (INIS)

Dimitrijevic, Z.

1993-01-01

The procedure for solving the transport equation by directly integrating for case one-dimensional uniform multigroup medium is shown. The solution is expressed in terms of linear combination of function H n (x,μ), and the coefficient is determined from given conditions. The solution is applied for homogeneous slab of critical thickness. (author)

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.

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

International Nuclear Information System (INIS)

Velarde, G.

1976-01-01

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

Solving the transport equation with quadratic finite elements: Theory and applications

International Nuclear Information System (INIS)

Ferguson, J.M.

1997-01-01

At the 4th Joint Conference on Computational Mathematics, the author presented a paper introducing a new quadratic finite element scheme (QFEM) for solving the transport equation. In the ensuing year the author has obtained considerable experience in the application of this method, including solution of eigenvalue problems, transmission problems, and solution of the adjoint form of the equation as well as the usual forward solution. He will present detailed results, and will also discuss other refinements of his transport codes, particularly for 3-dimensional problems on rectilinear and non-rectilinear grids

A Finite-Difference Solution of Solute Transport through a Membrane Bioreactor

Directory of Open Access Journals (Sweden)

B. Godongwana

2015-01-01

Full Text Available The current paper presents a theoretical analysis of the transport of solutes through a fixed-film membrane bioreactor (MBR, immobilised with an active biocatalyst. The dimensionless convection-diffusion equation with variable coefficients was solved analytically and numerically for concentration profiles of the solutes through the MBR. The analytical solution makes use of regular perturbation and accounts for radial convective flow as well as axial diffusion of the substrate species. The Michaelis-Menten (or Monod rate equation was assumed for the sink term, and the perturbation was extended up to second-order. In the analytical solution only the first-order limit of the Michaelis-Menten equation was considered; hence the linearized equation was solved. In the numerical solution, however, this restriction was lifted. The solution of the nonlinear, elliptic, partial differential equation was based on an implicit finite-difference method (FDM. An upwind scheme was employed for numerical stability. The resulting algebraic equations were solved simultaneously using the multivariate Newton-Raphson iteration method. The solution allows for the evaluation of the effect on the concentration profiles of (i the radial and axial convective velocity, (ii the convective mass transfer rates, (iii the reaction rates, (iv the fraction retentate, and (v the aspect ratio.

Solution of the Boltzmann equation for primary light ions and the transport of their fragments

Directory of Open Access Journals (Sweden)

J. Kempe

2010-10-01

Full Text Available The Boltzmann equation for the transport of pencil beams of light ions in semi-infinite uniform media has been calculated. The equation is solved for the practically important generalized 3D case of Gaussian incident primary light ion beams of arbitrary mean square radius, mean square angular spread, and covariance. The transport of the associated fragments in three dimensions is derived based on the known transport of the primary particles, taking the mean square angular spread of their production processes, as well as their energy loss and multiple scattering, into account. The analytical pencil and broad beam depth fluence and absorbed dose distributions are accurately expressed using recently derived analytical energy and range formulas. The contributions from low and high linear energy transfer (LET dose components were separately identified using analytical expressions. The analytical results are compared with SHIELD-HIT Monte Carlo (MC calculations and found to be in very good agreement. The pencil beam fluence and absorbed dose distributions of the primary particles are mainly influenced by an exponential loss of the primary ions combined with an increasing lateral spread due to multiple scattering and energy loss with increasing penetration depth. The associated fluence of heavy fragments is concentrated at small radii and so is the LET and absorbed dose distribution. Their transport is also characterized by the buildup of a slowing down spectrum which is quite similar to that of the primaries but with a wider energy and angular spread at increasing penetration depths. The range of the fragments is shorter or longer depending on their nuclear mass to charge ratio relative to that of the primary ions. The absorbed dose of the heavier fragments is fairly similar to that of the primary ions and also influenced by a rapidly increasing energy loss towards the end of their ranges. The present analytical solution of the Boltzmann equation

A spatially adaptive grid-refinement approach for the finite element solution of the even-parity Boltzmann transport equation

International Nuclear Information System (INIS)

Mirza, Anwar M.; Iqbal, Shaukat; Rahman, Faizur

2007-01-01

A spatially adaptive grid-refinement approach has been investigated to solve the even-parity Boltzmann transport equation. A residual based a posteriori error estimation scheme has been utilized for checking the approximate solutions for various finite element grids. The local particle balance has been considered as an error assessment criterion. To implement the adaptive approach, a computer program ADAFENT (adaptive finite elements for neutron transport) has been developed to solve the second order even-parity Boltzmann transport equation using K + variational principle for slab geometry. The program has a core K + module which employs Lagrange polynomials as spatial basis functions for the finite element formulation and Legendre polynomials for the directional dependence of the solution. The core module is called in by the adaptive grid generator to determine local gradients and residuals to explore the possibility of grid refinements in appropriate regions of the problem. The a posteriori error estimation scheme has been implemented in the outer grid refining iteration module. Numerical experiments indicate that local errors are large in regions where the flux gradients are large. A comparison of the spatially adaptive grid-refinement approach with that of uniform meshing approach for various benchmark cases confirms its superiority in greatly enhancing the accuracy of the solution without increasing the number of unknown coefficients. A reduction in the local errors of the order of 10 2 has been achieved using the new approach in some cases

A spatially adaptive grid-refinement approach for the finite element solution of the even-parity Boltzmann transport equation

Energy Technology Data Exchange (ETDEWEB)

Mirza, Anwar M. [Department of Computer Science, National University of Computer and Emerging Sciences, NUCES-FAST, A.K. Brohi Road, H-11, Islamabad (Pakistan)], E-mail: anwar.m.mirza@gmail.com; Iqbal, Shaukat [Faculty of Computer Science and Engineering, Ghulam Ishaq Khan (GIK) Institute of Engineering Science and Technology, Topi-23460, Swabi (Pakistan)], E-mail: shaukat@giki.edu.pk; Rahman, Faizur [Department of Physics, Allama Iqbal Open University, H-8 Islamabad (Pakistan)

2007-07-15

A spatially adaptive grid-refinement approach has been investigated to solve the even-parity Boltzmann transport equation. A residual based a posteriori error estimation scheme has been utilized for checking the approximate solutions for various finite element grids. The local particle balance has been considered as an error assessment criterion. To implement the adaptive approach, a computer program ADAFENT (adaptive finite elements for neutron transport) has been developed to solve the second order even-parity Boltzmann transport equation using K{sup +} variational principle for slab geometry. The program has a core K{sup +} module which employs Lagrange polynomials as spatial basis functions for the finite element formulation and Legendre polynomials for the directional dependence of the solution. The core module is called in by the adaptive grid generator to determine local gradients and residuals to explore the possibility of grid refinements in appropriate regions of the problem. The a posteriori error estimation scheme has been implemented in the outer grid refining iteration module. Numerical experiments indicate that local errors are large in regions where the flux gradients are large. A comparison of the spatially adaptive grid-refinement approach with that of uniform meshing approach for various benchmark cases confirms its superiority in greatly enhancing the accuracy of the solution without increasing the number of unknown coefficients. A reduction in the local errors of the order of 10{sup 2} has been achieved using the new approach in some cases.

Adaptive integral equation methods in transport theory

International Nuclear Information System (INIS)

Kelley, C.T.

1992-01-01

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

Transport equation solving methods

International Nuclear Information System (INIS)

Granjean, P.M.

1984-06-01

This work is mainly devoted to Csub(N) and Fsub(N) methods. CN method: starting from a lemma stated by Placzek, an equivalence is established between two problems: the first one is defined in a finite medium bounded by a surface S, the second one is defined in the whole space. In the first problem the angular flux on the surface S is shown to be the solution of an integral equation. This equation is solved by Galerkin's method. The Csub(N) method is applied here to one-velocity problems: in plane geometry, slab albedo and transmission with Rayleigh scattering, calculation of the extrapolation length; in cylindrical geometry, albedo and extrapolation length calculation with linear scattering. Fsub(N) method: the basic integral transport equation of the Csub(N) method is integrated on Case's elementary distributions; another integral transport equation is obtained: this equation is solved by a collocation method. The plane problems solved by the Csub(N) method are also solved by the Fsub(N) method. The Fsub(N) method is extended to any polynomial scattering law. Some simple spherical problems are also studied. Chandrasekhar's method, collision probability method, Case's method are presented for comparison with Csub(N) and Fsub(N) methods. This comparison shows the respective advantages of the two methods: a) fast convergence and possible extension to various geometries for Csub(N) method; b) easy calculations and easy extension to polynomial scattering for Fsub(N) method [fr

Calculation of radiation effects in solids by direct numerical solution of the adjoint transport equation

International Nuclear Information System (INIS)

Matthes, W.K.

1998-01-01

The 'adjoint transport equation in its integro-differential form' is derived for the radiation damage produced by atoms injected into solids. We reduce it to the one-dimensional form and prepare it for a numerical solution by: --discretizing the continuous variables energy, space and direction, --replacing the partial differential quotients by finite differences and --evaluating the collision integral by a double sum. By a proper manipulation of this double sum the adjoint transport equation turns into a (very large) set of linear equations with tridiagonal matrix which can be solved by a special (simple and fast) algorithm. The solution of this set of linear equations contains complete information on a specified damage type (e.g. the energy deposited in a volume V) in terms of the function D(i,E,c,x) which gives the damage produced by all particles generated in a cascade initiated by a particle of type i starting at x with energy E in direction c. It is essential to remark that one calculation gives the damage function D for the complete ranges of the variables {i,E,c and x} (for numerical reasons of course on grid-points in the {E,c,x}-space). This is most useful to applications where a general source-distribution S(i,E,c,x) of particles is given by the experimental setup (e.g. beam-window and and target in proton accelerator work. The beam-protons along their path through the window--or target material generate recoil atoms by elastic collisions or nuclear reactions. These recoil atoms form the particle source S). The total damage produced then is eventually given by: D = (Σ)i ∫ ∫ ∫ S(i, E, c, x)*D(i, E, c, x)*dE*dc*dx A Fortran-77 program running on a PC-486 was written for the overall procedure and applied to some problems

Semianalytical solutions of radioactive or reactive tracer transport in layered fractured media

International Nuclear Information System (INIS)

Moridis, G.J.; Bodvarsson, G.S.

2001-01-01

In this paper, semianalytical solutions are developed for the problem of transport of radioactive or reactive tracers (solutes or colloids) through a layered system of heterogeneous fractured media with misaligned fractures. The tracer transport equations in the matrix account for (a) diffusion, (b) surface diffusion (for solutes only), (c) mass transfer between the mobile and immobile water fractions, (d) linear kinetic or equilibrium physical, chemical, or combined solute sorption or colloid filtration, and (e) radioactive decay or first order chemical reactions. Any number of radioactive decay daughter products (or products of a linear, first-order reaction chain) can be tracked. The tracer-transport equations in the fractures account for the same processes, in addition to advection and hydrodynamic dispersion. Additionally, the colloid transport equations account for straining and velocity adjustments related to the colloidal size. The solutions, which are analytical in the Laplace space, are numerically inverted to provide the solution in time and can accommodate any number of fractured and/or porous layers. The solutions are verified using analytical solutions for limiting cases of solute and colloid transport through fractured and porous media. The effect of important parameters on the transport of 3 H, 237 Np and 239 Pu (and its daughters) is investigated in several test problems involving layered geological systems of varying complexity. 239 Pu colloid transport problems in multilayered systems indicate significant colloid accumulations at straining interfaces but much faster transport of the colloid than the corresponding strongly sorbing solute species

A biomechanical triphasic approach to the transport of nondilute solutions in articular cartilage.

Science.gov (United States)

Abazari, Alireza; Elliott, Janet A W; Law, Garson K; McGann, Locksley E; Jomha, Nadr M

2009-12-16

Biomechanical models for biological tissues such as articular cartilage generally contain an ideal, dilute solution assumption. In this article, a biomechanical triphasic model of cartilage is described that includes nondilute treatment of concentrated solutions such as those applied in vitrification of biological tissues. The chemical potential equations of the triphasic model are modified and the transport equations are adjusted for the volume fraction and frictional coefficients of the solutes that are not negligible in such solutions. Four transport parameters, i.e., water permeability, solute permeability, diffusion coefficient of solute in solvent within the cartilage, and the cartilage stiffness modulus, are defined as four degrees of freedom for the model. Water and solute transport in cartilage were simulated using the model and predictions of average concentration increase and cartilage weight were fit to experimental data to obtain the values of the four transport parameters. As far as we know, this is the first study to formulate the solvent and solute transport equations of nondilute solutions in the cartilage matrix. It is shown that the values obtained for the transport parameters are within the ranges reported in the available literature, which confirms the proposed model approach.

Analytical Solution for Fractional Derivative Gas-Flow Equation in Porous Media

KAUST Repository

El-Amin, Mohamed; Radwan, Ahmed G.; Sun, Shuyu

2017-01-01

In this paper, we introduce an analytical solution of the fractional derivative gas transport equation using the power-series technique. We present a new universal transform, namely, generalized Boltzmann change of variable which depends on the fractional order, time and space. This universal transform is employed to transfer the partial differential equation into an ordinary differential equation. Moreover, the convergence of the solution has been investigated and found that solutions are unconditionally converged. Results are introduced and discussed for the universal variable and other physical parameters such as porosity and permeability of the reservoir; time and space.

Analytical Solution for Fractional Derivative Gas-Flow Equation in Porous Media

KAUST Repository

El-Amin, Mohamed

2017-07-06

In this paper, we introduce an analytical solution of the fractional derivative gas transport equation using the power-series technique. We present a new universal transform, namely, generalized Boltzmann change of variable which depends on the fractional order, time and space. This universal transform is employed to transfer the partial differential equation into an ordinary differential equation. Moreover, the convergence of the solution has been investigated and found that solutions are unconditionally converged. Results are introduced and discussed for the universal variable and other physical parameters such as porosity and permeability of the reservoir; time and space.

Analytic solutions of hydrodynamics equations

International Nuclear Information System (INIS)

Coggeshall, S.V.

1991-01-01

Many similarity solutions have been found for the equations of one-dimensional (1-D) hydrodynamics. These special combinations of variables allow the partial differential equations to be reduced to ordinary differential equations, which must then be solved to determine the physical solutions. Usually, these reduced ordinary differential equations are solved numerically. In some cases it is possible to solve these reduced equations analytically to obtain explicit solutions. In this work a collection of analytic solutions of the 1-D hydrodynamics equations is presented. These can be used for a variety of purposes, including (i) numerical benchmark problems, (ii) as a basis for analytic models, and (iii) to provide insight into more complicated solutions

Application of preconditioned GMRES to the numerical solution of the neutron transport equation

International Nuclear Information System (INIS)

Patton, B.W.; Holloway, J.P.

2002-01-01

The generalized minimal residual (GMRES) method with right preconditioning is examined as an alternative to both standard and accelerated transport sweeps for the iterative solution of the diamond differenced discrete ordinates neutron transport equation. Incomplete factorization (ILU) type preconditioners are used to determine their effectiveness in accelerating GMRES for this application. ILU(τ), which requires the specification of a dropping criteria τ, proves to be a good choice for the types of problems examined in this paper. The combination of ILU(τ) and GMRES is compared with both DSA and unaccelerated transport sweeps for several model problems. It is found that the computational workload of the ILU(τ)-GMRES combination scales nonlinearly with the number of energy groups and quadrature order, making this technique most effective for problems with a small number of groups and discrete ordinates. However, the cost of preconditioner construction can be amortized over several calculations with different source and/or boundary values. Preconditioners built upon standard transport sweep algorithms are also evaluated as to their effectiveness in accelerating the convergence of GMRES. These preconditioners show better scaling with such problem parameters as the scattering ratio, the number of discrete ordinates, and the number of spatial meshes. These sweeps based preconditioners can also be cast in a matrix free form that greatly reduces storage requirements

Solution of linear transport equation using Chebyshev polynomials and Laplace transform

International Nuclear Information System (INIS)

Cardona, A.V.; Vilhena, M.T.M.B. de

1994-01-01

The Chebyshev polynomials and the Laplace transform are combined to solve, analytically, the linear transport equation in planar geometry, considering isotropic scattering and the one-group model. Numerical simulation is presented. (author)

Correspondence Between One- and Two-Equation Models for Solute Transport in Two-Region Heterogeneous Porous Media

KAUST Repository

Davit, Y.

2012-07-26

In this work, we study the transient behavior of homogenized models for solute transport in two-region porous media. We focus on the following three models: (1) a time non-local, two-equation model (2eq-nlt). This model does not rely on time constraints and, therefore, is particularly useful in the short-time regime, when the timescale of interest (t) is smaller than the characteristic time (τ 1) for the relaxation of the effective macroscale parameters (i. e., when t ≤ τ 1); (2) a time local, two-equation model (2eq). This model can be adopted when (t) is significantly larger than (τ 1) (i.e., when t≫τ 1); and (3) a one-equation, time-asymptotic formulation (1eq ∞). This model can be adopted when (t) is significantly larger than the timescale (τ 2) associated with exchange processes between the two regions (i. e., when t≫τ 2). In order to obtain insight into this transient behavior, we combine a theoretical approach based on the analysis of spatial moments with numerical and analytical results in several simple cases. The main result of this paper is to show that there is only a weak asymptotic convergence of the solution of (2eq) towards the solution of (1eq ∞) in terms of standardized moments but, interestingly, not in terms of centered moments. The physical interpretation of this result is that deviations from the Fickian situation persist in the limit of long times but that the spreading of the solute is eventually dominating these higher order effects. © 2012 Springer Science+Business Media B.V.

Nodal collocation approximation for the multidimensional PL equations applied to transport source problems

Energy Technology Data Exchange (ETDEWEB)

Verdu, G. [Departamento de Ingenieria Quimica Y Nuclear, Universitat Politecnica de Valencia, Cami de Vera, 14, 46022. Valencia (Spain); Capilla, M.; Talavera, C. F.; Ginestar, D. [Dept. of Nuclear Engineering, Departamento de Matematica Aplicada, Universitat Politecnica de Valencia, Cami de Vera, 14, 46022. Valencia (Spain)

2012-07-01

PL equations are classical high order approximations to the transport equations which are based on the expansion of the angular dependence of the angular neutron flux and the nuclear cross sections in terms of spherical harmonics. A nodal collocation method is used to discretize the PL equations associated with a neutron source transport problem. The performance of the method is tested solving two 1D problems with analytical solution for the transport equation and a classical 2D problem. (authors)

Nodal collocation approximation for the multidimensional PL equations applied to transport source problems

International Nuclear Information System (INIS)

Verdu, G.; Capilla, M.; Talavera, C. F.; Ginestar, D.

2012-01-01

PL equations are classical high order approximations to the transport equations which are based on the expansion of the angular dependence of the angular neutron flux and the nuclear cross sections in terms of spherical harmonics. A nodal collocation method is used to discretize the PL equations associated with a neutron source transport problem. The performance of the method is tested solving two 1D problems with analytical solution for the transport equation and a classical 2D problem. (authors)

Coupling of neutron transport equations. First results; Couplage d`equations en transport neutronique. premiere approche 1D monocinetique

Energy Technology Data Exchange (ETDEWEB)

Bal, G.

1995-07-01

To achieve whole core calculations of the neutron transport equation, we have to follow this 2 step method: space and energy homogenization of the assemblies; resolution of the homogenized equation on the whole core. However, this is no more valid when accidents occur (for instance depressurization causing locally strong heterogeneous media). One solution consists then in coupling two kinds of resolutions: a fine computation on the damaged cell (fine mesh, high number of energy groups) coupled with a coarse one everywhere else. We only deal here with steady state solutions (which already live in 6D spaces). We present here two such methods: The coupling by transmission of homogenized sections and the coupling by transmission of boundary conditions. To understand what this coupling is, we first restrict ourselves to 1D with respect to space in one energy group. The first two chapters deal with a recall of basic properties of the neutron transport equation. We give at chapter 3 some indications of the behaviour of the flux with respect to the cross sections. We present at chapter 4 some couplings and give some properties. Chapter 5 is devoted to a presentation of some numerical applications. (author). 9 refs., 7 figs.

Second order time evolution of the multigroup diffusion and P1 equations for radiation transport

International Nuclear Information System (INIS)

Olson, Gordon L.

2011-01-01

Highlights: → An existing multigroup transport algorithm is extended to be second-order in time. → A new algorithm is presented that does not require a grey acceleration solution. → The two algorithms are tested with 2D, multi-material problems. → The two algorithms have comparable computational requirements. - Abstract: An existing solution method for solving the multigroup radiation equations, linear multifrequency-grey acceleration, is here extended to be second order in time. This method works for simple diffusion and for flux-limited diffusion, with or without material conduction. A new method is developed that does not require the solution of an averaged grey transport equation. It is effective solving both the diffusion and P 1 forms of the transport equation. Two dimensional, multi-material test problems are used to compare the solution methods.

Numerical Solution of Heun Equation Via Linear Stochastic Differential Equation

Directory of Open Access Journals (Sweden)

Hamidreza Rezazadeh

2014-05-01

Full Text Available In this paper, we intend to solve special kind of ordinary differential equations which is called Heun equations, by converting to a corresponding stochastic differential equation(S.D.E.. So, we construct a stochastic linear equation system from this equation which its solution is based on computing fundamental matrix of this system and then, this S.D.E. is solved by numerically methods. Moreover, its asymptotic stability and statistical concepts like expectation and variance of solutions are discussed. Finally, the attained solutions of these S.D.E.s compared with exact solution of corresponding differential equations.

Transport Visualization for Studying Mass Transfer and Solute Transport in Permeable Media

International Nuclear Information System (INIS)

Roy Haggerty

2004-01-01

Understanding and predicting mass transfer coupled with solute transport in permeable media is central to several energy-related programs at the US Department of Energy (e.g., CO 2 sequestration, nuclear waste disposal, hydrocarbon extraction, and groundwater remediation). Mass transfer is the set of processes that control movement of a chemical between mobile (advection-dominated) domains and immobile (diffusion- or sorption-dominated) domains within a permeable medium. Consequences of mass transfer on solute transport are numerous and may include (1) increased sequestration time within geologic formations; (2) reduction in average solute transport velocity by as much as several orders of magnitude; (3) long ''tails'' in concentration histories during removal of a solute from a permeable medium; (4) poor predictions of solute behavior over long time scales; and (5) changes in reaction rates due to mass transfer influences on pore-scale mixing of solutes. Our work produced four principle contributions: (1) the first comprehensive visualization of solute transport and mass transfer in heterogeneous porous media; (2) the beginnings of a theoretical framework that encompasses both macrodispersion and mass transfer within a single set of equations; (3) experimental and analytical tools necessary for understanding mixing and aqueous reaction in heterogeneous, granular porous media; (4) a clear experimental demonstration that reactive transport is often not accurately described by a simple coupling of the convection-dispersion equation with chemical reaction equations. The work shows that solute transport in heterogeneous media can be divided into 3 regimes--macrodispersion, advective mass transfer, and diffusive mass transfer--and that these regimes can be predicted quantitatively in binary media. We successfully predicted mass transfer in each of these regimes and verified the prediction by completing quantitative visualization experiments in each of the regimes, the

End-Member Formulation of Solid Solutions and Reactive Transport

Energy Technology Data Exchange (ETDEWEB)

Lichtner, Peter C. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2015-09-01

A model for incorporating solid solutions into reactive transport equations is presented based on an end-member representation. Reactive transport equations are solved directly for the composition and bulk concentration of the solid solution. Reactions of a solid solution with an aqueous solution are formulated in terms of an overall stoichiometric reaction corresponding to a time-varying composition and exchange reactions, equivalent to reaction end-members. Reaction rates are treated kinetically using a transition state rate law for the overall reaction and a pseudo-kinetic rate law for exchange reactions. The composition of the solid solution at the onset of precipitation is assumed to correspond to the least soluble composition, equivalent to the composition at equilibrium. The stoichiometric saturation determines if the solid solution is super-saturated with respect to the aqueous solution. The method is implemented for a simple prototype batch reactor using Mathematica for a binary solid solution. Finally, the sensitivity of the results on the kinetic rate constant for a binary solid solution is investigated for reaction of an initially stoichiometric solid phase with an undersaturated aqueous solution.

A variational solution of transport equation based on spherical geometry

International Nuclear Information System (INIS)

Liu Hui; Zhang Ben'ai

2002-01-01

A variational method with differential forms gives better precision for numerical solution of transport critical problem based on spherical geometry, and its computation seems simple than other approximate methods

Evaluation of finite difference and FFT-based solutions of the transport of intensity equation.

Science.gov (United States)

Zhang, Hongbo; Zhou, Wen-Jing; Liu, Ying; Leber, Donald; Banerjee, Partha; Basunia, Mahmudunnabi; Poon, Ting-Chung

2018-01-01

A finite difference method is proposed for solving the transport of intensity equation. Simulation results show that although slower than fast Fourier transform (FFT)-based methods, finite difference methods are able to reconstruct the phase with better accuracy due to relaxed assumptions for solving the transport of intensity equation relative to FFT methods. Finite difference methods are also more flexible than FFT methods in dealing with different boundary conditions.

A family of analytical solutions of a nonlinear diffusion-convection equation

Science.gov (United States)

Hayek, Mohamed

2018-01-01

Despite its popularity in many engineering fields, the nonlinear diffusion-convection equation has no general analytical solutions. This work presents a family of closed-form analytical traveling wave solutions for the nonlinear diffusion-convection equation with power law nonlinearities. This kind of equations typically appears in nonlinear problems of flow and transport in porous media. The solutions that are addressed are simple and fully analytical. Three classes of analytical solutions are presented depending on the type of the nonlinear diffusion coefficient (increasing, decreasing or constant). It has shown that the structure of the traveling wave solution is strongly related to the diffusion term. The main advantage of the proposed solutions is that they are presented in a unified form contrary to existing solutions in the literature where the derivation of each solution depends on the specific values of the diffusion and convection parameters. The proposed closed-form solutions are simple to use, do not require any numerical implementation, and may be implemented in a simple spreadsheet. The analytical expressions are also useful to mathematically analyze the structure and properties of the solutions.

A closed-form solution for the two-dimensional Fokker-Planck equation for electron transport in the range of Compton Effect

Energy Technology Data Exchange (ETDEWEB)

Rodriguez, B.D.A. [Universidade Federal Rio Grande do Sul, Programa de Pos-Graduacao em Engenharia Mecanica, Rua Portuguesa 218/304, 90650-12 Porto Alegre, RS (Brazil)], E-mail: barbara.arodriguez@gmail.com; Vilhena, M.T. [Universidade Federal Rio Grande do Sul, Departamento de Matematica Pura e Aplicada, Porto Alegre, RS (Brazil)], E-mail: vilhena@mat.ufrgs.br; Borges, V. [Universidade Federal Rio Grande do Sul, Programa de Pos-Graduacao em Engenharia Mecanica, Rua Portuguesa 218/304, 90650-12 Porto Alegre, RS (Brazil)], E-mail: borges@ufrgs.br; Hoff, G. [Pontificia Universidade Catolica do Rio Grande do Sul, Faculdade de Fisica, Porto Alegre, RS (Brazil)], E-mail: hoff@pucrs.br

2008-05-15

In this paper we solve the Fokker-Planck (FP) equation, an alternative approach for the Boltzmann transport equation for charged particles in a rectangular domain. To construct the solution we begin applying the P{sub N} approximation in the angular variable and the Laplace Transform in the x-variable, thus obtaining a first order linear differential equation in y-variable, which the solution is straightforward. The angular flux of electrons and the parameters of the medium are used for the calculation of the energy deposited by the secondary electrons generated by Compton Effect. The remaining effects will not be taken into account. The results will be presented under absorbed energy form in several points of interested. We present numerical simulations and comparisons with results obtained by using Geant4 (version 8) program which applies the Monte Carlo's technique to low energy libraries for a two-dimensional problem assuming the screened Rutherford differential scattering cross-section.

A closed-form solution for the two-dimensional Fokker-Planck equation for electron transport in the range of Compton Effect

International Nuclear Information System (INIS)

Rodriguez, B.D.A.; Vilhena, M.T.; Borges, V.; Hoff, G.

2008-01-01

In this paper we solve the Fokker-Planck (FP) equation, an alternative approach for the Boltzmann transport equation for charged particles in a rectangular domain. To construct the solution we begin applying the P N approximation in the angular variable and the Laplace Transform in the x-variable, thus obtaining a first order linear differential equation in y-variable, which the solution is straightforward. The angular flux of electrons and the parameters of the medium are used for the calculation of the energy deposited by the secondary electrons generated by Compton Effect. The remaining effects will not be taken into account. The results will be presented under absorbed energy form in several points of interested. We present numerical simulations and comparisons with results obtained by using Geant4 (version 8) program which applies the Monte Carlo's technique to low energy libraries for a two-dimensional problem assuming the screened Rutherford differential scattering cross-section

Application of the finite element method to the neutron transport equation

International Nuclear Information System (INIS)

Martin, W.R.

1976-01-01

This paper examines the theoretical and practical application of the finite element method to the neutron transport equation. It is shown that in principle the system of equations obtained by application of the finite element method can be solved with certain physical restrictions concerning the criticality of the medium. The convergence of this approximate solution to the exact solution with mesh refinement is examined, and a non-optical estimate of the convergence rate is obtained analytically. It is noted that the numerical results indicate a faster convergence rate and several approaches to obtain this result analytically are outlined. The practical application of the finite element method involved the development of a computer code capable of solving the neutron transport equation in 1-D plane geometry. Vacuum, reflecting, or specified incoming boundary conditions may be analyzed, and all are treated as natural boundary conditions. The time-dependent transport equation is also examined and it is shown that the application of the finite element method in conjunction with the Crank-Nicholson time discretization method results in a system of algebraic equations which is readily solved. Numerical results are given for several critical slab eigenvalue problems, including anisotropic scattering, and the results compare extremely well with benchmark results. It is seen that the finite element code is more efficient than a standard discrete ordinates code for certain problems. A problem with severe heterogeneities is considered and it is shown that the use of discontinuous spatial and angular elements results in a marked improvement in the results. Finally, time-dependent problems are examined and it is seen that the phenomenon of angular mode separation makes the numerical treatment of the transport equation in slab geometry a considerable challenge, with the result that the angular mesh has a dominant effect on obtaining acceptable solutions

An integral equation arising in two group neutron transport theory

International Nuclear Information System (INIS)

Cassell, J S; Williams, M M R

2003-01-01

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

Solution of the Boltzmann-Fokker-Planck transport equation using exponential nodal schemes

International Nuclear Information System (INIS)

Ortega J, R.; Valle G, E. del

2003-01-01

There are carried out charge and energy calculations deposited due to the interaction of electrons with a plate of a certain material, solving numerically the electron transport equation for the Boltzmann-Fokker-Planck approach of first order in plate geometry with a computer program denominated TEOD-NodExp (Transport of Electrons in Discreet Ordinates, Nodal Exponentials), using the proposed method by the Dr. J. E. Morel to carry out the discretization of the variable energy and several spatial discretization schemes, denominated exponentials nodal. It is used the Fokker-Planck equation since it represents an approach of the Boltzmann transport equation that is been worth whenever it is predominant the dispersion of small angles, that is to say, resulting dispersion in small dispersion angles and small losses of energy in the transport of charged particles. Such electrons could be those that they face with a braking plate in a device of thermonuclear fusion. In the present work its are considered electrons of 1 MeV that impact isotropically on an aluminum plate. They were considered three different thickness of plate that its were designated as problems 1, 2 and 3. In the calculations it was used the discrete ordinate method S 4 with expansions of the dispersion cross sections until P 3 order. They were considered 25 energy groups of uniform size between the minimum energy of 0.1 MeV and the maximum of 1.0 MeV; the one spatial intervals number it was considered variable and it was assigned the values of 10, 20 and 30. (Author)

Linear superposition solutions to nonlinear wave equations

International Nuclear Information System (INIS)

Liu Yu

2012-01-01

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

Measurement of Membrane Characteristics Using the Phenomenological Equation and the Overall Mass Transport Equation in Ion-Exchange Membrane Electrodialysis of Saline Water

Directory of Open Access Journals (Sweden)

Yoshinobu Tanaka

2012-01-01

Full Text Available The overall membrane pair characteristics included in the overall mass transport equation are understandable using the phenomenological equations expressed in the irreversible thermodynamics. In this investigation, the overall membrane pair characteristics (overall transport number , overall solute permeability , overall electro-osmotic permeability and overall hydraulic permeability were measured by seawater electrodialysis changing current density, temperature and salt concentration, and it was found that occasionally takes minus value. For understanding the above phenomenon, new concept of the overall concentration reflection coefficient ∗ is introduced from the phenomenological equation. This is the aim of this investigation. ∗ is defined for describing the permselectivity between solutes and water molecules in the electrodialysis system just after an electric current interruption. ∗ is expressed by the function of and . ∗ is generally larger than 1 and is positive, but occasionally ∗ becomes less than 1 and becomes negative. Negative means that ions are transferred with water molecules (solvent from desalting cells toward concentrating cells just after an electric current interruption, indicating up-hill transport or coupled transport between water molecules and solutes.

Estimating the solute transport parameters of the spatial fractional advection-dispersion equation using Bees Algorithm.

Science.gov (United States)

Mehdinejadiani, Behrouz

2017-08-01

This study represents the first attempt to estimate the solute transport parameters of the spatial fractional advection-dispersion equation using Bees Algorithm. The numerical studies as well as the experimental studies were performed to certify the integrity of Bees Algorithm. The experimental ones were conducted in a sandbox for homogeneous and heterogeneous soils. A detailed comparative study was carried out between the results obtained from Bees Algorithm and those from Genetic Algorithm and LSQNONLIN routines in FracFit toolbox. The results indicated that, in general, the Bees Algorithm much more accurately appraised the sFADE parameters in comparison with Genetic Algorithm and LSQNONLIN, especially in the heterogeneous soil and for α values near to 1 in the numerical study. Also, the results obtained from Bees Algorithm were more reliable than those from Genetic Algorithm. The Bees Algorithm showed the relative similar performances for all cases, while the Genetic Algorithm and the LSQNONLIN yielded different performances for various cases. The performance of LSQNONLIN strongly depends on the initial guess values so that, compared to the Genetic Algorithm, it can more accurately estimate the sFADE parameters by taking into consideration the suitable initial guess values. To sum up, the Bees Algorithm was found to be very simple, robust and accurate approach to estimate the transport parameters of the spatial fractional advection-dispersion equation. Copyright © 2017 Elsevier B.V. All rights reserved.

Numerical solution of the Neutron Transport Equation using discontinuous nodal methods at X-Y geometry

International Nuclear Information System (INIS)

Delfin L, A.

1996-01-01

The purpose of this work is to solve the neutron transport equation in discrete-ordinates and X-Y geometry by developing and using the strong discontinuous and strong modified discontinuous nodal finite element schemes. The strong discontinuous and modified strong discontinuous nodal finite element schemes go from two to ten interpolation parameters per cell. They are describing giving a set D c and polynomial space S c corresponding for each scheme BDMO, RTO, BL, BDM1, HdV, BDFM1, RT1, BQ and BDM2. The solution is obtained solving the neutron transport equation moments for each nodal scheme by developing the basis functions defined by Pascal triangle and the Legendre moments giving in the polynomial space S c and, finally, looking for the non singularity of the resulting linear system. The linear system is numerically solved using a computer program for each scheme mentioned . It uses the LU method and forward and backward substitution and makes a partition of the domain in cells. The source terms and angular flux are calculated, using the directions and weights associated to the S N approximation and solving the angular flux moments to find the effective multiplication constant. The programs are written in Fortran language, using the dynamic allocation of memory to increase efficiently the available memory of the computing equipment. (Author)

Estimating the solute transport parameters of the spatial fractional advection-dispersion equation using Bees Algorithm

Science.gov (United States)

Mehdinejadiani, Behrouz

2017-08-01

This study represents the first attempt to estimate the solute transport parameters of the spatial fractional advection-dispersion equation using Bees Algorithm. The numerical studies as well as the experimental studies were performed to certify the integrity of Bees Algorithm. The experimental ones were conducted in a sandbox for homogeneous and heterogeneous soils. A detailed comparative study was carried out between the results obtained from Bees Algorithm and those from Genetic Algorithm and LSQNONLIN routines in FracFit toolbox. The results indicated that, in general, the Bees Algorithm much more accurately appraised the sFADE parameters in comparison with Genetic Algorithm and LSQNONLIN, especially in the heterogeneous soil and for α values near to 1 in the numerical study. Also, the results obtained from Bees Algorithm were more reliable than those from Genetic Algorithm. The Bees Algorithm showed the relative similar performances for all cases, while the Genetic Algorithm and the LSQNONLIN yielded different performances for various cases. The performance of LSQNONLIN strongly depends on the initial guess values so that, compared to the Genetic Algorithm, it can more accurately estimate the sFADE parameters by taking into consideration the suitable initial guess values. To sum up, the Bees Algorithm was found to be very simple, robust and accurate approach to estimate the transport parameters of the spatial fractional advection-dispersion equation.

Solute transport in aquifers: The comeback of the advection dispersion equation and the First Order Approximation

Science.gov (United States)

Fiori, A.; Zarlenga, A.; Jankovic, I.; Dagan, G.

2017-12-01

Natural gradient steady flow of mean velocity U takes place in heterogeneous aquifers of random logconductivity Y = lnK , characterized by the normal univariate PDF f(Y) and autocorrelation ρY, of variance σY2 and horizontal integral scale I. Solute transport is quantified by the Breakthrough Curve (BTC) M at planes at distance x from the injection plane. The study builds on the extensive 3D numerical simulations of flow and transport of Jankovic et al. (2017) for different conductivity structures. The present study further explores the predictive capabilities of the Advection Dispersion Equation (ADE), with macrodispersivity αL given by the First Order Approximation (FOA), by checking in a quantitative manner its applicability. After a discussion on the suitable boundary conditions for ADE, we find that the ADE-FOA solution is a sufficiently accurate predictor for applications, the many other sources of uncertainty prevailing in practice notwithstanding. We checked by least squares and by comparison of travel time of quantiles of M that indeed the analytical Inverse Gaussian M with αL =σY2 I , is able to fit well the bulk of the simulated BTCs. It tends to underestimate the late arrival time of the thin and persistent tail. The tail is better reproduced by the semi-analytical MIMSCA model, which also allows for a physical explanation of the success of the Inverse Gaussian solution. Examination of the pertinent longitudinal mass distribution shows that it is different from the commonly used Gaussian one in the analysis of field experiments, and it captures the main features of the plume measurements of the MADE experiment. The results strengthen the confidence in the applicability of the ADE and the FOA to predicting longitudinal spreading in solute transport through heterogeneous aquifers of stationary random structure.

Adomian decomposition method for solving the telegraph equation in charged particle transport

International Nuclear Information System (INIS)

Abdou, M.A.

2005-01-01

In this paper, the analysis for the telegraph equation in case of isotropic small angle scattering from the Boltzmann transport equation for charged particle is presented. The Adomian decomposition is used to solve the telegraph equation. By means of MAPLE the Adomian polynomials of obtained series (ADM) solution have been calculated. The behaviour of the distribution function are shown graphically. The results reported in this article provide further evidence of the usefulness of Adomain decomposition for obtaining solution of linear and nonlinear problems

Multigroup discrete ordinates solution of Boltzmann-Fokker-Planck equations and cross section library development of ion transport

International Nuclear Information System (INIS)

Prinja, A.K.

1995-08-01

We have developed and successfully implemented a two-dimensional bilinear discontinuous in space and time, used in conjunction with the S N angular approximation, to numerically solve the time dependent, one-dimensional, one-speed, slab geometry, (ion) transport equation. Numerical results and comparison with analytical solutions have shown that the bilinear-discontinuous (BLD) scheme is third-order accurate in the space ad time dimensions independently. Comparison of the BLD results with diamond-difference methods indicate that the BLD method is both quantitavely and qualitatively superior to the DD scheme. We note that the form of the transport operator is such that these conclusions carry over to energy dependent problems that include the constant-slowing-down-approximation term, and to multiple space dimensions or combinations thereof. An optimized marching or inversion scheme or a parallel algorithm should be investigated to determine if the increased accuracy can compensate for the extra overhead required for a BLD solution, and then could be compared to other discretization methods such as nodal or characteristic schemes

CACTUS, a characteristics solution to the neutron transport equations in complicated geometries

International Nuclear Information System (INIS)

Halsall, M.J.

1980-04-01

CACTUS has been written to solve the multigroup neutron transport equation in a general two-dimensional geometry. The method is based upon a characteristics formulation for the problem in which the transport equation is integrated explicitly along straight line tracks that are suitably distributed throughout the problem. Source distributions and scattering are assumed to be isotropic, but the only restriction on geometry is that the outer boundary should be rectangular. Within this rectangular boundary the user is free to build his problem geometry using any combination of intersecting straight lines and circular arcs. The theory of the method is described, followed by some details of a coding, a sensitivity study on the number of tracks required to integrate fluxes in a particular problem, a user's guide, and a few test cases. (author)

Unified solution of the Boltzmann equation for electron and ion velocity distribution functions and transport coefficients in weakly ionized plasmas

Science.gov (United States)

Konovalov, Dmitry A.; Cocks, Daniel G.; White, Ronald D.

2017-10-01

The velocity distribution function and transport coefficients for charged particles in weakly ionized plasmas are calculated via a multi-term solution of Boltzmann's equation and benchmarked using a Monte-Carlo simulation. A unified framework for the solution of the original full Boltzmann's equation is presented which is valid for ions and electrons, avoiding any recourse to approximate forms of the collision operator in various limiting mass ratio cases. This direct method using Lebedev quadratures over the velocity and scattering angles avoids the need to represent the ion mass dependence in the collision operator through an expansion in terms of the charged particle to neutral mass ratio. For the two-temperature Burnett function method considered in this study, this amounts to avoiding the need for the complex Talmi-transformation methods and associated mass-ratio expansions. More generally, we highlight the deficiencies in the two-temperature Burnett function method for heavy ions at high electric fields to calculate the ion velocity distribution function, even though the transport coefficients have converged. Contribution to the Topical Issue "Physics of Ionized Gases (SPIG 2016)", edited by Goran Poparic, Bratislav Obradovic, Dragana Maric and Aleksandar Milosavljevic.

Modeling of water and solute transport under variably saturated conditions: state of the art

International Nuclear Information System (INIS)

Lappala, E.G.

1980-01-01

This paper reviews the equations used in deterministic models of mass and energy transport in variably saturated porous media. Analytic, quasi-analytic, and numerical solution methods to the nonlinear forms of transport equations are discussed with respect to their advantages and limitations. The factors that influence the selection of a modeling method are discussed in this paper; they include the following: (1) the degree of coupling required among the equations describing the transport of liquids, gases, solutes, and energy; (2) the inclusion of an advection term in the equations; (3) the existence of sharp fronts; (4) the degree of nonlinearity and hysteresis in the transport coefficients and boundary conditions; (5) the existence of complex boundaries; and (6) the availability and reliability of data required by the models

Quasi-exact solutions of nonlinear differential equations

OpenAIRE

Kudryashov, Nikolay A.; Kochanov, Mark B.

2014-01-01

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

Spherical harmonics and energy polynomial solution of the Boltzmann equation for neutrons, 1

International Nuclear Information System (INIS)

Toledo, P.S. de

1974-01-01

The approximate solution of the source-free energy-dependent Boltzmann transport equation for neutrons in plane geometry and isotropic scattering case was given by Leonard and Ferziger using a truncated development in a series of energy-polynomials for the energy dependent neutron flux and solving exactly for the angular dependence. The presence in the general solution of eigenfunctions belonging to a continuous spectrum gives rise to difficult analytical problems in the application of their method even to simple problems. To avoid such difficulties, the angular dependence is treated by a spherical harmonics method and a general solution of the energy-dependent transport equation in plane geometry and isotropic scattering is obtained, in spite of the appearance of matrices as argument of the angular polynomials [pt

Approximate method for solving the velocity dependent transport equation in a slab lattice

International Nuclear Information System (INIS)

Ferrari, A.

1966-01-01

A method is described that is intended to provide an approximate solution of the transport equation in a medium simulating a water-moderated plate filled reactor core. This medium is constituted by a periodic array of water channels and absorbing plates. The velocity dependent transport equation in slab geometry is included. The computation is performed in a water channel: the absorbing plates are accounted for by the boundary conditions. The scattering of neutrons in water is assumed isotropic, which allows the use of a double Pn approximation to deal with the angular dependence. This method is able to represent the discontinuity of the angular distribution at the channel boundary. The set of equations thus obtained is dependent only on x and v and the coefficients are independent on x. This solution suggests to try solutions involving Legendre polynomials. This scheme leads to a set of equations v dependent only. To obtain an explicit solution, a thermalization model must now be chosen. Using the secondary model of Cadilhac a solution of this set is easy to get. The numerical computations were performed with a particular secondary model, the well-known model of Wigner and Wilkins. (author) [fr

Stochastic uncertainty analysis for solute transport in randomly heterogeneous media using a Karhunen‐Loève‐based moment equation approach

Science.gov (United States)

Liu, Gaisheng; Lu, Zhiming; Zhang, Dongxiao

2007-01-01

A new approach has been developed for solving solute transport problems in randomly heterogeneous media using the Karhunen‐Loève‐based moment equation (KLME) technique proposed by Zhang and Lu (2004). The KLME approach combines the Karhunen‐Loève decomposition of the underlying random conductivity field and the perturbative and polynomial expansions of dependent variables including the hydraulic head, flow velocity, dispersion coefficient, and solute concentration. The equations obtained in this approach are sequential, and their structure is formulated in the same form as the original governing equations such that any existing simulator, such as Modular Three‐Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems (MT3DMS), can be directly applied as the solver. Through a series of two‐dimensional examples, the validity of the KLME approach is evaluated against the classical Monte Carlo simulations. Results indicate that under the flow and transport conditions examined in this work, the KLME approach provides an accurate representation of the mean concentration. For the concentration variance, the accuracy of the KLME approach is good when the conductivity variance is 0.5. As the conductivity variance increases up to 1.0, the mismatch on the concentration variance becomes large, although the mean concentration can still be accurately reproduced by the KLME approach. Our results also indicate that when the conductivity variance is relatively large, neglecting the effects of the cross terms between velocity fluctuations and local dispersivities, as done in some previous studies, can produce noticeable errors, and a rigorous treatment of the dispersion terms becomes more appropriate.

Coupling between solute transport and chemical reactions models. Acoplamiento de modelos de transporte de solutos y de modelos de reacciones quimicas

Energy Technology Data Exchange (ETDEWEB)

Samper, J.; Ajora, C. (Instituto de Ciencias de la Tierra, CSIC, Barcerlona (Spain))

1993-01-01

During subsurface transport, reactive solutes are subject to a variety of hydrodynamic and chemical processes. The major hydrodynamic processes include advection and convection, dispersion and diffusion. The key chemical processes are complexation including hydrolysis and acid-base reactions, dissolution-precipitation, reduction-oxidation, adsorption and ion exchange. The combined effects of all these processes on solute transport must satisfy the principle of conservation of mass. The statement of conservation of mass for N mobile species leads to N partial differential equations. Traditional solute transport models often incorporate the effects of hydrodynamic processes rigorously but oversimplify chemical interactions among aqueous species. Sophisticated chemical equilibrium models, on the other hand, incorporate a variety of chemical processes but generally assume no-flow systems. In the past decade, coupled models accounting for complex hydrological and chemical processes, with varying degrees of sophistication, have been developed. The existing models of reactive transport employ two basic sets of equations. The transport of solutes is described by a set of partial differential equations, and the chemical processes, under the assumption of equilibrium, are described by a set of nonlinear algebraic equations. An important consideration in any approach is the choice of primary dependent variables. Most existing models cannot account for the complete set of chemical processes, cannot be easily extended to include mixed chemical equilibria and kinetics, and cannot handle practical two and three dimensional problems. The difficulties arise mainly from improper selection of the primary variables in the transport equations. (Author) 38 refs.

Comparison of analytical transport and stochastic solutions for neutron slowing down in an infinite medium

International Nuclear Information System (INIS)

Jahshan, S.N.; Wemple, C.A.; Ganapol, B.D.

1993-01-01

A comparison of the numerical solutions of the transport equation describing the steady neutron slowing down in an infinite medium with constant cross sections is made with stochastic solutions obtained from tracking successive neutron histories in the same medium. The transport equation solution is obtained using a numerical Laplace transform inversion algorithm. The basis for the algorithm is an evaluation of the Bromwich integral without analytical continuation. Neither the transport nor the stochastic solution is limited in the number of scattering species allowed. The medium may contain an absorption component as well. (orig.)

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

International Nuclear Information System (INIS)

Sirendaoreji

2006-01-01

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

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

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)

Cable Connected Spinning Spacecraft, 1. the Canonical Equations, 2. Urban Mass Transportation, 3

Science.gov (United States)

Sitchin, A.

1972-01-01

Work on the dynamics of cable-connected spinning spacecraft was completed by formulating the equations of motion by both the canonical equations and Lagrange's equations and programming them for numerical solution on a digital computer. These energy-based formulations will permit future addition of the effect of cable mass. Comparative runs indicate that the canonical formulation requires less computer time. Available literature on urban mass transportation was surveyed. Areas of the private rapid transit concept of urban transportation are also studied.

Exponentially-convergent Monte Carlo for the 1-D transport equation

International Nuclear Information System (INIS)

Peterson, J. R.; Morel, J. E.; Ragusa, J. C.

2013-01-01

We define a new exponentially-convergent Monte Carlo method for solving the one-speed 1-D slab-geometry transport equation. This method is based upon the use of a linear discontinuous finite-element trial space in space and direction to represent the transport solution. A space-direction h-adaptive algorithm is employed to restore exponential convergence after stagnation occurs due to inadequate trial-space resolution. This methods uses jumps in the solution at cell interfaces as an error indicator. Computational results are presented demonstrating the efficacy of the new approach. (authors)

The gBL transport equations

International Nuclear Information System (INIS)

Mynick, H.E.

1989-05-01

The transport equations arising from the ''generalized Balescu- Lenard'' (gBL) collision operator are obtained, and some of their properties examined. The equations contain neoclassical and turbulent transport as two special cases, having the same structure. The resultant theory offers potential explanation for a number of results not well understood, including the anomalous pinch, observed ratios of Q/ΓT on TFTR, and numerical reproduction of ASDEX profiles by a model for turbulent transport invoked without derivation, but by analogy to neoclassical theory. The general equations are specialized to consideration of a number of particular transport mechanisms of interest. 10 refs

New solutions of Heun's general equation

International Nuclear Information System (INIS)

Ishkhanyan, Artur; Suominen, Kalle-Antti

2003-01-01

We show that in four particular cases the derivative of the solution of Heun's general equation can be expressed in terms of a solution to another Heun's equation. Starting from this property, we use the Gauss hypergeometric functions to construct series solutions to Heun's equation for the mentioned cases. Each of the hypergeometric functions involved has correct singular behaviour at only one of the singular points of the equation; the sum, however, has correct behaviour. (letter to the editor)

Numerical simulation for fractional order stationary neutron transport equation using Haar wavelet collocation method

Energy Technology Data Exchange (ETDEWEB)

Saha Ray, S., E-mail: santanusaharay@yahoo.com; Patra, A.

2014-10-15

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: In this paper the numerical solution for the fractional order stationary neutron transport equation is presented using Haar wavelet Collocation Method (HWCM). Haar wavelet collocation method is efficient and powerful in solving wide class of linear and nonlinear differential equations. This paper intends to provide an application of Haar wavelets to nuclear science problems. This paper describes the application of Haar wavelets for the numerical solution of fractional order stationary neutron transport equation in homogeneous medium with isotropic scattering. The proposed method is mathematically very simple, easy and fast. To demonstrate about the efficiency and applicability of the method, two test problems are discussed.

Radiative transport equation for the Mittag-Leffler path length distribution

Science.gov (United States)

Liemert, André; Kienle, Alwin

2017-05-01

In this paper, we consider the radiative transport equation for infinitely extended scattering media that are characterized by the Mittag-Leffler path length distribution p (ℓ ) =-∂ℓEα(-σtℓα ) , which is a generalization of the usually assumed Lambert-Beer law p (ℓ ) =σtexp(-σtℓ ) . In this context, we derive the infinite-space Green's function of the underlying fractional transport equation for the spherically symmetric medium as well as for the one-dimensional string. Moreover, simple analytical solutions are presented for the prediction of the radiation field in the single-scattering approximation. The resulting equations are compared with Monte Carlo simulations in the steady-state and time domain showing, within the stochastic nature of the simulations, an excellent agreement.

Differential-algebraic solutions of the heat equation

OpenAIRE

Buchstaber, Victor M.; Netay, Elena Yu.

2014-01-01

In this work we introduce the notion of differential-algebraic ansatz for the heat equation and explicitly construct heat equation and Burgers equation solutions given a solution of a homogeneous non-linear ordinary differential equation of a special form. The ansatz for such solutions is called the $n$-ansatz, where $n+1$ is the order of the differential equation.

Stochastic dynamics modeling solute transport in porous media modeling solute transport in porous media

CERN Document Server

Kulasiri, Don

2002-01-01

Most of the natural and biological phenomena such as solute transport in porous media exhibit variability which can not be modeled by using deterministic approaches. There is evidence in natural phenomena to suggest that some of the observations can not be explained by using the models which give deterministic solutions. Stochastic processes have a rich repository of objects which can be used to express the randomness inherent in the system and the evolution of the system over time. The attractiveness of the stochastic differential equations (SDE) and stochastic partial differential equations (SPDE) come from the fact that we can integrate the variability of the system along with the scientific knowledge pertaining to the system. One of the aims of this book is to explaim some useufl concepts in stochastic dynamics so that the scientists and engineers with a background in undergraduate differential calculus could appreciate the applicability and appropriateness of these developments in mathematics. The ideas ...

Transport-constrained extensions of collision and track length estimators for solutions of radiative transport problems

International Nuclear Information System (INIS)

Kong, Rong; Spanier, Jerome

2013-01-01

In this paper we develop novel extensions of collision and track length estimators for the complete space-angle solutions of radiative transport problems. We derive the relevant equations, prove that our new estimators are unbiased, and compare their performance with that of more conventional estimators. Such comparisons based on numerical solutions of simple one dimensional slab problems indicate the the potential superiority of the new estimators for a wide variety of more general transport problems

Transport methods: general. 7. Formulation of a Fourier-Boltzmann Transformation to Solve the Three-Dimensional Transport Equation

International Nuclear Information System (INIS)

Stancic, V.

2001-01-01

This paper presents some elements of a new approach to solve analytically the linearized three-dimensional (3-D) transport equation of neutral particles. Since this task is of such special importance, we present some results of a paper that is still in progress. The most important is that using this transformation, an integro-differential equation with an analytical solution is obtained. For this purpose, a simplest 3-D equation is being considered which describes the transport process in an infinite medium. Until now, this equation has been analytically considered either using the Laplace transform with respect to time parameter t or applying the Fourier transform over the space coordinate. Both of them reduce the number of differential terms in the equation; however, evaluation of the inverse transformation is complicated. In this paper, we introduce for the first time a Fourier transform induced by the Boltzmann operator. For this, we use a complete set of 3-D eigenfunctions of the Boltzmann transport operator defined in a similar way as those that have been already used in 3-D transport theory as a basic set to transform the transport equation. This set consists of a continuous part and a discrete one with spectral measure. The density distribution equation shows the known form asymptotic behavior. Several applications are to be performed using this equation and compared to the benchmark one. Such an analysis certainly would be out of the available space

Polygons of differential equations for finding exact solutions

International Nuclear Information System (INIS)

Kudryashov, Nikolai A.; Demina, Maria V.

2007-01-01

A method for finding exact solutions of nonlinear differential equations is presented. Our method is based on the application of polygons corresponding to nonlinear differential equations. It allows one to express exact solutions of the equation studied through solutions of another equation using properties of the basic equation itself. The ideas of power geometry are used and developed. Our approach has a pictorial interpretation, which is illustrative and effective. The method can be also applied for finding transformations between solutions of differential equations. To demonstrate the method application exact solutions of several equations are found. These equations are: the Korteveg-de Vries-Burgers equation, the generalized Kuramoto-Sivashinsky equation, the fourth-order nonlinear evolution equation, the fifth-order Korteveg-de Vries equation, the fifth-order modified Korteveg-de Vries equation and the sixth-order nonlinear evolution equation describing turbulent processes. Some new exact solutions of nonlinear evolution equations are given

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

On polynomial solutions of the Heun equation

International Nuclear Information System (INIS)

Gurappa, N; Panigrahi, Prasanta K

2004-01-01

By making use of a recently developed method to solve linear differential equations of arbitrary order, we find a wide class of polynomial solutions to the Heun equation. We construct the series solution to the Heun equation before identifying the polynomial solutions. The Heun equation extended by the addition of a term, -σ/x, is also amenable for polynomial solutions. (letter to the editor)

Cusping, transport and variance of solutions to generalized Fokker-Planck equations

Science.gov (United States)

Carnaffan, Sean; Kawai, Reiichiro

2017-06-01

We study properties of solutions to generalized Fokker-Planck equations through the lens of the probability density functions of anomalous diffusion processes. In particular, we examine solutions in terms of their cusping, travelling wave behaviours, and variance, within the framework of stochastic representations of generalized Fokker-Planck equations. We give our analysis in the cases of anomalous diffusion driven by the inverses of the stable, tempered stable and gamma subordinators, demonstrating the impact of changing the distribution of waiting times in the underlying anomalous diffusion model. We also analyse the cases where the underlying anomalous diffusion contains a Lévy jump component in the parent process, and when a diffusion process is time changed by an uninverted Lévy subordinator. On the whole, we present a combination of four criteria which serve as a theoretical basis for model selection, statistical inference and predictions for physical experiments on anomalously diffusing systems. We discuss possible applications in physical experiments, including, with reference to specific examples, the potential for model misclassification and how combinations of our four criteria may be used to overcome this issue.

Explicit finite-difference solution of two-dimensional solute transport with periodic flow in homogenous porous media

Directory of Open Access Journals (Sweden)

Djordjevich Alexandar

2017-12-01

Full Text Available The two-dimensional advection-diffusion equation with variable coefficients is solved by the explicit finitedifference method for the transport of solutes through a homogenous two-dimensional domain that is finite and porous. Retardation by adsorption, periodic seepage velocity, and a dispersion coefficient proportional to this velocity are permitted. The transport is from a pulse-type point source (that ceases after a period of activity. Included are the firstorder decay and zero-order production parameters proportional to the seepage velocity, and periodic boundary conditions at the origin and at the end of the domain. Results agree well with analytical solutions that were reported in the literature for special cases. It is shown that the solute concentration profile is influenced strongly by periodic velocity fluctuations. Solutions for a variety of combinations of unsteadiness of the coefficients in the advection-diffusion equation are obtainable as particular cases of the one demonstrated here. This further attests to the effectiveness of the explicit finite difference method for solving two-dimensional advection-diffusion equation with variable coefficients in finite media, which is especially important when arbitrary initial and boundary conditions are required.

Mixed-hybrid finite element method for the transport equation and diffusion approximation of transport problems

International Nuclear Information System (INIS)

Cartier, J.

2006-04-01

This thesis focuses on mathematical analysis, numerical resolution and modelling of the transport equations. First of all, we deal with numerical approximation of the solution of the transport equations by using a mixed-hybrid scheme. We derive and study a mixed formulation of the transport equation, then we analyse the related variational problem and present the discretization and the main properties of the scheme. We particularly pay attention to the behavior of the scheme and we show its efficiency in the diffusion limit (when the mean free path is small in comparison with the characteristic length of the physical domain). We present academical benchmarks in order to compare our scheme with other methods in many physical configurations and validate our method on analytical test cases. Unstructured and very distorted meshes are used to validate our scheme. The second part of this thesis deals with two transport problems. The first one is devoted to the study of diffusion due to boundary conditions in a transport problem between two plane plates. The second one consists in modelling and simulating radiative transfer phenomenon in case of the industrial context of inertial confinement fusion. (author)

LTSN solution of the adjoint neutron transport equation with arbitrary source for high order of quadrature in a homogeneous slab

International Nuclear Information System (INIS)

Goncalves, Glenio A.; Orengo, Gilberto; Vilhena, Marco Tullio M.B. de; Graca, Claudio O.

2002-01-01

In this work we present the LTS N solution of the adjoint transport equation for an arbitrary source, testing the aptness of this analytical solution for high order of quadrature in transport problems and comparing some preliminary results with the ANISN computations in a homogeneneous slab geometry. In order to do that we apply the new formulation for the LTS N method based on the invariance projection property, becoming possible to handle problems with arbitrary sources and demanding high order of quadrature or deep penetration. This new approach for the LTS N method is important both for direct and adjoint transport calculations and its development was inspired by the necessity of using generalized adjoint sources for important calculations. Although the mathematical convergence has been proved for an arbitrary source, when the quadrature order or deep penetration is required the LTS N method presents computational overflow even for simple sources (sin, cos, exp, polynomial). With the new formulation we eliminate this drawback and in this work we report the numerical simulations testing the new approach

Singularly perturbed Burger-Huxley equation: Analytical solution ...

African Journals Online (AJOL)

user

solutions of singularly perturbed nonlinear differential equations. ... for solving generalized Burgers-Huxley equation but this equation is not singularly ...... Solitary waves solutions of the generalized Burger Huxley equations, Journal of.

Coupled geochemical and solute transport code development

International Nuclear Information System (INIS)

Morrey, J.R.; Hostetler, C.J.

1985-01-01

A number of coupled geochemical hydrologic codes have been reported in the literature. Some of these codes have directly coupled the source-sink term to the solute transport equation. The current consensus seems to be that directly coupling hydrologic transport and chemical models through a series of interdependent differential equations is not feasible for multicomponent problems with complex geochemical processes (e.g., precipitation/dissolution reactions). A two-step process appears to be the required method of coupling codes for problems where a large suite of chemical reactions must be monitored. Two-step structure requires that the source-sink term in the transport equation is supplied by a geochemical code rather than by an analytical expression. We have developed a one-dimensional two-step coupled model designed to calculate relatively complex geochemical equilibria (CTM1D). Our geochemical module implements a Newton-Raphson algorithm to solve heterogeneous geochemical equilibria, involving up to 40 chemical components and 400 aqueous species. The geochemical module was designed to be efficient and compact. A revised version of the MINTEQ Code is used as a parent geochemical code

A three-dimensional neutron transport benchmark solution

International Nuclear Information System (INIS)

Ganapol, B.D.; Kornreich, D.E.

1993-01-01

For one-group neutron transport theory in one dimension, several powerful analytical techniques have been developed to solve the neutron transport equation, including Caseology, Wiener-Hopf factorization, and Fourier and Laplace transform methods. In addition, after a Fourier transform in the transverse plane and formulation of a pseudo problem, two-dimensional (2-D) and three-dimensional (3-D) problems can be solved using the techniques specifically developed for the one-dimensional (1-D) case. Numerical evaluation of the resulting expressions requiring an inversion in the transverse plane have been successful for 2-D problems but becomes exceedingly difficult in the 3-D case. In this paper, we show that by using the symmetry along the beam direction, a 2-D problem can be transformed into a 3-D problem in an infinite medium. The numerical solution to the 3-D problem is then demonstrated. Thus, a true 3-D transport benchmark solution can be obtained from a well-established numerical solution to a 2-D problem

Calculations of reactivity based in the solution of the Neutron transport equation in X Y geometry and Lineal perturbation theory

International Nuclear Information System (INIS)

Valle G, E. del; Mugica R, C.A.

2005-01-01

In our country, in last congresses, Gomez et al carried out reactivity calculations based on the solution of the diffusion equation for an energy group using nodal methods in one dimension and the TPL approach (Lineal Perturbation Theory). Later on, Mugica extended the application to the case of multigroup so much so much in one as in two dimensions (X Y geometry) with excellent results. Presently work is carried out similar calculations but this time based on the solution of the neutron transport equation in X Y geometry using nodal methods and again the TPL approximation. The idea is to provide a calculation method that allows to obtain in quick form the reactivity solving the direct problem as well as the enclosed problem of the not perturbed problem. A test problem for the one that results are provided for the effective multiplication factor is described and its are offered some conclusions. (Author)

Cellular neural network to the spherical harmonics approximation of neutron transport equation in x-y geometry. Part I: Modeling and verification for time-independent solution

International Nuclear Information System (INIS)

Pirouzmand, Ahmad; Hadad, Kamal

2011-01-01

Highlights: → This paper describes the solution of time-independent neutron transport equation. → Using a novel method based on cellular neural networks (CNNs) coupled with P N method. → Utilize the CNN model to simulate spatial scalar flux distribution in steady state. → The accuracy, stability, and capabilities of CNN model are examined in x-y geometry. - Abstract: This paper describes a novel method based on using cellular neural networks (CNN) coupled with spherical harmonics method (P N ) to solve the time-independent neutron transport equation in x-y geometry. To achieve this, an equivalent electrical circuit based on second-order form of neutron transport equation and relevant boundary conditions is obtained using CNN method. We use the CNN model to simulate spatial response of scalar flux distribution in the steady state condition for different order of spherical harmonics approximations. The accuracy, stability, and capabilities of CNN model are examined in 2D Cartesian geometry for fixed source and criticality problems.

Doubly periodic solutions of the modified Kawahara equation

International Nuclear Information System (INIS)

Zhang Dan

2005-01-01

Some doubly periodic (Jacobi elliptic function) solutions of the modified Kawahara equation are presented in closed form. Our approach is to introduce a new auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct doubly periodic solutions of the modified Kawahara equation. When the module m → 1, these solutions degenerate to the exact solitary wave solutions of the equation. Then we reveal the relation of some exact solutions for the modified Kawahara equation obtained by other authors

Spurious solutions in few-body equations

International Nuclear Information System (INIS)

Adhikari, S.K.; Gloeckle, W.

1979-01-01

After Faddeev and Yakubovskii showed how to write connected few-body equations which are free from discrete spurious solutions various authors have proposed different connected few-body scattering equations. Federbush first pointed out that Weinberg's formulation admits the existence of discrete spurious solutions. In this paper we investigate the possibility and consequence of the existence of spurious solutions in some of the few-body formulations. Contrary to a proof by Hahn, Kouri, and Levin and by Bencze and Tandy the channel coupling array scheme of Kouri, Levin, and Tobocman which is also the starting point of a formulation by Hahn is shown to admit spurious solutions. We can show that the set of six coupled four-body equations proposed independently by Mitra, Gillespie, Sugar, and Panchapakesan, by Rosenberg, by Alessandrini, and by Takahashi and Mishima and the seven coupled four-body equations proposed by Sloan and related by matrix multipliers to basic sets which correspond uniquely to the Schroedinger equation. These multipliers are likely to give spurious solutions to these equations. In all these cases spuriosities are shown to have no hazardous consequence if one is interested in studying the scattering problem

Solute transport model for radioisotopes in layered soil

International Nuclear Information System (INIS)

Essel, P.

2010-01-01

The study considered the transport of a radioactive solute in solution from the surface of the earth down through the soil to the ground water when there is an accidental or intentional spillage of a radioactive material on the surface. The finite difference method was used to model the spatial and temporal profile of moisture content in a soil column using the θ-based Richard's equation leading to solution of the convective-dispersive equation for non-adsorbing solutes numerically. A matlab code has been generated to predict the transport of the radioactive contaminant, spilled on the surface of a vertically heterogeneous soil made up of two layers to determine the residence time of the solute in the unsaturated zone, the time it takes the contaminant to reach the groundwater and the amount of the solute entering the groundwater in various times and the levels of pollution in those times. The model predicted that, then there is a spillage of 7.2g of tritium, on the surface of the ground at the study area, it will take two years for the radionuclide to enter the groundwater and fifteen years to totally leave the unsaturated zone. There is therefore the need to try as much as possible to avoid intentional or accidental spillage of the radionuclide since it has long term effect. (au)

Some results on the neutron transport and the coupling of equations; Quelques resultats sur le transport neutronique et le couplage d`equations

Energy Technology Data Exchange (ETDEWEB)

Bal, G. [Electricite de France (EDF), Direction des Etudes et Recherches, 92 - Clamart (France)

1997-12-31

Neutron transport in nuclear reactors is well modeled by the linear Boltzmann transport equation. Its resolution is relatively easy but very expensive. To achieve whole core calculations, one has to consider simpler models, such as diffusion or homogeneous transport equations. However, the solutions may become inaccurate in particular situations (as accidents for instance). That is the reason why we wish to solve the equations on small area accurately and more coarsely on the remaining part of the core. It is than necessary to introduce some links between different discretizations or modelizations. In this note, we give some results on the coupling of different discretizations of all degrees of freedom of the integral-differential neutron transport equation (two degrees for the angular variable, on for the energy component, and two or three degrees for spatial position respectively in 2D (cylindrical symmetry) and 3D). Two chapters are devoted to the coupling of discrete ordinates methods (for angular discretization). The first one is theoretical and shows the well posing of the coupled problem, whereas the second one deals with numerical applications of practical interest (the results have been obtained from the neutron transport code developed at the R and D, which has been modified for introducing the coupling). Next, we present the nodal scheme RTN0, used for the spatial discretization. We show well posing results for the non-coupled and the coupled problems. At the end, we deal with the coupling of energy discretizations for the multigroup equations obtained by homogenization. Some theoretical results of the discretization of the velocity variable (well-posing of problems), which do not deal directly with the purposes of coupling, are presented in the annexes. (author). 34 refs.

Development of interfacial area transport equation

International Nuclear Information System (INIS)

Kim, Seung Jin; Ishii, Mamoru; Kelly, Joseph

2005-01-01

The interfacial area transport equation dynamically models the changes in interfacial structures along the flow field by mechanistically modeling the creation and destruction of dispersed phase. Hence, when employed in the numerical thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Accounting for the substantial differences in the transport mechanism for various sizes of bubbles, the transport equation is formulated for two characteristic groups of bubbles. The group 1 equation describes the transport of small-dispersed bubbles, whereas the group 2 equation describes the transport of large cap, slug or churn-turbulent bubbles. To evaluate the feasibility and reliability of interfacial area transport equation available at present, it is benchmarked by an extensive database established in various two-phase flow configurations spanning from bubbly to churn-turbulent flow regimes. The geometrical effect in interfacial area transport is examined by the data acquired in vertical air-water two-phase flow through round pipes of various sizes and a confined flow duct, and by those acquired in vertical co-current downward air-water two-phase flow through round pipes of two different sizes

Elliptic random-walk equation for suspension and tracer transport in porous media

DEFF Research Database (Denmark)

Shapiro, Alexander; Bedrikovetsky, P. G.

2008-01-01

. The new theory predicts delay of the maximum of the tracer, compared to the velocity of the flow, while its forward "tail" contains much more particles than in the solution of the classical parabolic (advection-dispersion) equation. This is in agreement with the experimental observations and predictions......We propose a new approach to transport of the suspensions and tracers in porous media. The approach is based on a modified version of the continuous time random walk (CTRW) theory. In the framework of this theory we derive an elliptic transport equation. The new equation contains the time...... of the CTRW theory. (C) 2008 Elsevier B.V. All rights reserved....

On an analytical representation of the solution of the one-dimensional transport equation for a multi-group model in planar geometry

Energy Technology Data Exchange (ETDEWEB)

Fernandes, Julio C.L.; Vilhena, Marco T.; Bodmann, Bardo E.J., E-mail: julio.lombaldo@ufrgs.br, E-mail: mtmbvilhena@gmail.com, E-mail: bardo.bodmann@ufrgs.br [Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, RS (Brazil). Dept. de Matematica Pura e Aplicada; Dulla, Sandra; Ravetto, Piero, E-mail: sandra.dulla@polito.it, E-mail: piero.ravetto@polito.it [Dipartimento di Energia, Politecnico di Torino, Piemonte (Italy)

2015-07-01

In this work we generalize the solution of the one-dimensional neutron transport equation to a multi- group approach in planar geometry. The basic idea of this work consists in consider the hierarchical construction of a solution for a generic number G of energy groups, starting from a mono-energetic solution. The hierarchical method follows the reasoning of the decomposition method. More specifically, the additional terms from adding energy groups is incorporated into the recursive scheme as source terms. This procedure leads to an analytical representation for the solution with G energy groups. The recursion depth is related to the accuracy of the solution, that may be evaluated after each recursion step. The authors present a heuristic analysis of stability for the results. Numerical simulations for a specific example with four energy groups and a localized pulsed source. (author)

Solutions for a non-Markovian diffusion equation

International Nuclear Information System (INIS)

Lenzi, E.K.; Evangelista, L.R.; Lenzi, M.K.; Ribeiro, H.V.; Oliveira, E.C. de

2010-01-01

Solutions for a non-Markovian diffusion equation are investigated. For this equation, we consider a spatial and time dependent diffusion coefficient and the presence of an absorbent term. The solutions exhibit an anomalous behavior which may be related to the solutions of fractional diffusion equations and anomalous diffusion.

The Dirac equation and its solutions

CERN Document Server

Bagrov, Vladislav G

2014-01-01

Dirac equations are of fundamental importance for relativistic quantum mechanics and quantum electrodynamics. In relativistic quantum mechanics, the Dirac equation is referred to as one-particle wave equation of motion for electron in an external electromagnetic field. In quantum electrodynamics, exact solutions of this equation are needed to treat the interaction between the electron and the external field exactly.In particular, all propagators of a particle, i.e., the various Green's functions, are constructed in a certain way by using exact solutions of the Dirac equation.

The Dirac equation and its solutions

International Nuclear Information System (INIS)

Bagrov, Vladislav G.; Gitman, Dmitry; P.N. Lebedev Physical Institute, Moscow; Tomsk State Univ., Tomsk

2013-01-01

The Dirac equation is of fundamental importance for relativistic quantum mechanics and quantum electrodynamics. In relativistic quantum mechanics, the Dirac equation is referred to as one-particle wave equation of motion for electron in an external electromagnetic field. In quantum electrodynamics, exact solutions of this equation are needed to treat the interaction between the electron and the external field exactly. In particular, all propagators of a particle, i.e., the various Green's functions, are constructed in a certain way by using exact solutions of the Dirac equation.

Solute transport modelling with the variable temporally dependent ...

Indian Academy of Sciences (India)

Pintu Das

2018-02-07

Feb 7, 2018 ... in a finite domain with time-dependent sources and dis- tance-dependent dispersivities. Also, existing ... solute transport in multi-layered porous media using gen- eralized integral transform technique with .... methods for solving the fractional reaction-–sub-diffusion equation. To solve numerically the Eqs.

The Dirac equation and its solutions

Energy Technology Data Exchange (ETDEWEB)

Bagrov, Vladislav G. [Tomsk State Univ., Tomsk (Russian Federation). Dept. of Quantum Field Theroy; Gitman, Dmitry [Sao Paulo Univ. (Brazil). Inst. de Fisica; P.N. Lebedev Physical Institute, Moscow (Russian Federation); Tomsk State Univ., Tomsk (Russian Federation). Faculty of Physics

2013-07-01

The Dirac equation is of fundamental importance for relativistic quantum mechanics and quantum electrodynamics. In relativistic quantum mechanics, the Dirac equation is referred to as one-particle wave equation of motion for electron in an external electromagnetic field. In quantum electrodynamics, exact solutions of this equation are needed to treat the interaction between the electron and the external field exactly. In particular, all propagators of a particle, i.e., the various Green's functions, are constructed in a certain way by using exact solutions of the Dirac equation.

Solution of the Baxter equation

International Nuclear Information System (INIS)

Janik, R.A.

1996-01-01

We present a method of construction of a family of solutions of the Baxter equation arising in the Generalized Leading Logarithmic Approximation (GLLA) of the QCD pomeron. The details are given for the exchange of N = 2 reggeons but everything can be generalized in a straightforward way to arbitrary N. A specific choice of solutions is shown to reproduce the correct energy levels for half integral conformal weights. It is shown that the Baxter's equation must be supplemented by an additional condition on the solution. (author)

Generalized multivariate Fokker-Planck equations derived from kinetic transport theory and linear nonequilibrium thermodynamics

International Nuclear Information System (INIS)

Frank, T.D.

2002-01-01

We study many particle systems in the context of mean field forces, concentration-dependent diffusion coefficients, generalized equilibrium distributions, and quantum statistics. Using kinetic transport theory and linear nonequilibrium thermodynamics we derive for these systems a generalized multivariate Fokker-Planck equation. It is shown that this Fokker-Planck equation describes relaxation processes, has stationary maximum entropy distributions, can have multiple stationary solutions and stationary solutions that differ from Boltzmann distributions

A compartmentalized solute transport model for redox zones in contaminated aquifers: 1. Theory and development

Science.gov (United States)

Abrams , Robert H.; Loague, Keith

2000-01-01

This paper, the first of two parts [see Abrams and Loague, this issue], takes the compartmentalized approach for the geochemical evolution of redox zones presented by Abrams et al. [1998] and embeds it within a solute transport framework. In this paper the compartmentalized approach is generalized to facilitate the description of its incorporation into a solute transport simulator. An equivalent formulation is developed which removes any discontinuities that may occur when switching compartments. Rate‐limited redox reactions are modeled with a modified Monod relationship that allows either the organic substrate or the electron acceptor to be the rate‐limiting reactant. Thermodynamic constraints are used to inhibit lower‐energy redox reactions from occurring under infeasible geochemical conditions without imposing equilibrium on the lower‐energy reactions. The procedure used allows any redox reaction to be simulated as being kinetically limited or thermodynamically limited, depending on local geochemical conditions. Empirical reaction inhibition methods are not needed. The sequential iteration approach (SIA), a technique which allows the number of solute transport equations to be reduced, is adopted to solve the coupled geochemical/solute transport problem. When the compartmentalized approach is embedded within the SIA, with the total analytical concentration of each component as the dependent variable in the transport equation, it is possible to reduce the number of transport equations even further than with the unmodified SIA. A one‐dimensional, coupled geochemical/solute transport simulation is presented in which redox zones evolve dynamically in time and space. The compartmentalized solute transport (COMPTRAN) model described in this paper enables the development of redox zones to be simulated under both kinetic and thermodynamic constraints. The modular design of COMPTRAN facilitates the use of many different, preexisting solute transport and

Numerical solution of the equation of neutrons transport on plane geometry by analytical schemes using acceleration by synthetic diffusion

International Nuclear Information System (INIS)

Alonso-Vargas, G.

1991-01-01

A computer program has been developed which uses a technique of synthetic acceleration by diffusion by analytical schemes. Both in the diffusion equation as in that of transport, analytical schemes were used which allowed a substantial time saving in the number of iterations required by source iteration method to obtain the K e ff. The program developed ASD (Synthetic Diffusion Acceleration) by diffusion was written in FORTRAN and can be executed on a personal computer with a hard disc and mathematical O-processor. The program is unlimited as to the number of regions and energy groups. The results obtained by the ASD program for K e ff is nearly completely concordant with those of obtained utilizing the ANISN-PC code for different analytical type problems in this work. The ASD program allowed obtention of an approximate solution of the neutron transport equation with a relatively low number of internal reiterations with good precision. One of its applications would be in the direct determinations of axial distribution neutronic flow in a fuel assembly as well as in the obtention of the effective multiplication factor. (Author)

Solution of spatially homogeneous model Boltzmann equations by means of Lie groups of transformations

International Nuclear Information System (INIS)

Foroutan, A.

1992-05-01

The essential mathematical challenge in transport theory is based on the nonlinearity of the integro-differential equations governing classical thermodynamic systems on molecular kinetic level. It is the aim of this thesis to gain exact analytical solutions to the model Boltzmann equation suggested by Tjon and Wu. Such solutions afford a deeper insight into the dynamics of rarefied gases. Tjon and Wu have provided a stochastic model of a Boltzmann equation. Its transition probability depends only on the relative speed of the colliding particles. This assumption leads in the case of two translational degrees of freedom to an integro-differential equation of convolution type. According to this convolution structure the integro-differential equation is Laplace transformed. The result is a nonlinear partial differential equation. The investigation of the symmetries of this differential equation by means of Lie groups of transformation enables us to transform the originally nonlinear partial differential equation into ordinary differential equation into ordinary differential equations of Bernoulli type. (author)

The H-N method for solving linear transport equation: theory and application

International Nuclear Information System (INIS)

Kaskas, A.; Gulecyuz, M.C.; Tezcan, C.

2002-01-01

The system of singular integral equation which is obtained from the integro-differential form of the linear transport equation as a result of Placzec lemma is solved. Application are given using the exit distributions and the infinite medium Green's function. The same theoretical results are also obtained with the use of the singular eigenfunction of the method of elementary solutions

A lattice Boltzmann model for solute transport in open channel flow

Science.gov (United States)

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

2018-01-01

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

Mixed-hybrid finite element method for the transport equation and diffusion approximation of transport problems; Resolution de l'equation du transport par une methode d'elements finis mixtes-hybrides et approximation par la diffusion de problemes de transport

Energy Technology Data Exchange (ETDEWEB)

Cartier, J

2006-04-15

This thesis focuses on mathematical analysis, numerical resolution and modelling of the transport equations. First of all, we deal with numerical approximation of the solution of the transport equations by using a mixed-hybrid scheme. We derive and study a mixed formulation of the transport equation, then we analyse the related variational problem and present the discretization and the main properties of the scheme. We particularly pay attention to the behavior of the scheme and we show its efficiency in the diffusion limit (when the mean free path is small in comparison with the characteristic length of the physical domain). We present academical benchmarks in order to compare our scheme with other methods in many physical configurations and validate our method on analytical test cases. Unstructured and very distorted meshes are used to validate our scheme. The second part of this thesis deals with two transport problems. The first one is devoted to the study of diffusion due to boundary conditions in a transport problem between two plane plates. The second one consists in modelling and simulating radiative transfer phenomenon in case of the industrial context of inertial confinement fusion. (author)

Solution of the within-group multidimensional discrete ordinates transport equations on massively parallel architectures

Science.gov (United States)

Zerr, Robert Joseph

2011-12-01

The integral transport matrix method (ITMM) has been used as the kernel of new parallel solution methods for the discrete ordinates approximation of the within-group neutron transport equation. The ITMM abandons the repetitive mesh sweeps of the traditional source iterations (SI) scheme in favor of constructing stored operators that account for the direct coupling factors among all the cells and between the cells and boundary surfaces. The main goals of this work were to develop the algorithms that construct these operators and employ them in the solution process, determine the most suitable way to parallelize the entire procedure, and evaluate the behavior and performance of the developed methods for increasing number of processes. This project compares the effectiveness of the ITMM with the SI scheme parallelized with the Koch-Baker-Alcouffe (KBA) method. The primary parallel solution method involves a decomposition of the domain into smaller spatial sub-domains, each with their own transport matrices, and coupled together via interface boundary angular fluxes. Each sub-domain has its own set of ITMM operators and represents an independent transport problem. Multiple iterative parallel solution methods have investigated, including parallel block Jacobi (PBJ), parallel red/black Gauss-Seidel (PGS), and parallel GMRES (PGMRES). The fastest observed parallel solution method, PGS, was used in a weak scaling comparison with the PARTISN code. Compared to the state-of-the-art SI-KBA with diffusion synthetic acceleration (DSA), this new method without acceleration/preconditioning is not competitive for any problem parameters considered. The best comparisons occur for problems that are difficult for SI DSA, namely highly scattering and optically thick. SI DSA execution time curves are generally steeper than the PGS ones. However, until further testing is performed it cannot be concluded that SI DSA does not outperform the ITMM with PGS even on several thousand or tens of

Limit Theorems and Their Relation to Solute Transport in Simulated Fractured Media

Science.gov (United States)

Reeves, D. M.; Benson, D. A.; Meerschaert, M. M.

2003-12-01

Solute particles that travel through fracture networks are subject to wide velocity variations along a restricted set of directions. This may result in super-Fickian dispersion along a few primary scaling directions. The fractional advection-dispersion equation (FADE), a modification of the original advection-dispersion equation in which a fractional derivative replaces the integer-order dispersion term, has the ability to model rapid, non-Gaussian solute transport. The FADE assumes that solute particle motions converge to either α -stable or operator stable densities, which are modeled by spatial fractional derivatives. In multiple dimensions, the multi-fractional dispersion derivative dictates the order and weight of differentiation in all directions, which correspond to the statistics of large particle motions in all directions. This study numerically investigates the presence of super- Fickian solute transport through simulated two-dimensional fracture networks. An ensemble of networks is gen

Maximal stochastic transport in the Lorenz equations

Energy Technology Data Exchange (ETDEWEB)

Agarwal, Sahil, E-mail: sahil.agarwal@yale.edu [Program in Applied Mathematics, Yale University, New Haven (United States); Wettlaufer, J.S., E-mail: john.wettlaufer@yale.edu [Program in Applied Mathematics, Yale University, New Haven (United States); Departments of Geology & Geophysics, Mathematics and Physics, Yale University, New Haven (United States); Mathematical Institute, University of Oxford, Oxford (United Kingdom); Nordita, Royal Institute of Technology and Stockholm University, Stockholm (Sweden)

2016-01-08

We calculate the stochastic upper bounds for the Lorenz equations using an extension of the background method. In analogy with Rayleigh–Bénard convection the upper bounds are for heat transport versus Rayleigh number. As might be expected, the stochastic upper bounds are larger than the deterministic counterpart of Souza and Doering [1], but their variation with noise amplitude exhibits interesting behavior. Below the transition to chaotic dynamics the upper bounds increase monotonically with noise amplitude. However, in the chaotic regime this monotonicity depends on the number of realizations in the ensemble; at a particular Rayleigh number the bound may increase or decrease with noise amplitude. The origin of this behavior is the coupling between the noise and unstable periodic orbits, the degree of which depends on the degree to which the ensemble represents the ergodic set. This is confirmed by examining the close returns plots of the full solutions to the stochastic equations and the numerical convergence of the noise correlations. The numerical convergence of both the ensemble and time averages of the noise correlations is sufficiently slow that it is the limiting aspect of the realization of these bounds. Finally, we note that the full solutions of the stochastic equations demonstrate that the effect of noise is equivalent to the effect of chaos.

Methodology for obtaining a solution for the three-dimensional Boltzmann transport equation and an expression for the calculation of the total doses considering Compton scattering simulated by Klein-Nishina

International Nuclear Information System (INIS)

Rodriguez, Barbara A.; Borges, Volnei; Vilhena, Marco Tullio

2005-01-01

In this work we would like to obtain a formulation of an analytic method for the solution of the three dimensional transport equation considering Compton scattering and an expression for total doses due to gamma radiation, where the deposited energy by the free electron will be considered. For that, we will work with two equations: the first one for the photon transport, considering the Klein-Nishina kernel and energy multigroup model, and the second one considering the free electron with the screened Rutherford scattering. (author)

Two-dimensional finite element solution for the simultaneous transport of water and solutes through a nonhomogeneous aquifer under transient saturated unsaturated flow conditions

International Nuclear Information System (INIS)

Gureghian, A.B.

1979-01-01

A mathematical model of ground water transport through an aquifer is presented. The solute of interest is a metal tracer or radioactive material which may undergo decay through a sorbing unconfined aquifer. The subject is developed under the following headings: flow equation, solute equation, boundary conditions, finite element formulation, element formulation, solution scheme (flow equation, solute equation), results and discussions, water movement in a ditch drained aquifer under transient state, water and solute movement in a homogeneous and unsaturated soil, transport of 226 Ra in nonhomogeneous aquifer, tailings pond lined, and tailings pond unlined. It is concluded that this mathematical model may have a wide variety of applications. The uranium milling industry may find it useful to evaluate the hydrogeological suitability of their disposal sites. It may prove suited for the design of clay disposal ponds destined to hold hazardous liquids. It may also provide a means of estimating the long-term impact of radionuclides or other pollutants on the quality of ground water. 31 references, 9 figures, 3 tables

The 'generalized Balescu-Lenard' transport equations

International Nuclear Information System (INIS)

Mynick, H.E.

1990-01-01

The transport equations arising from the 'generalized Balescu-Lenard' collision operator are obtained and some of their properties examined. The equations contain neoclassical and turbulent transport as two special cases having the same structure. The resultant theory offers a possible explanation for a number of results not well understood, including the anomalous pinch, observed ratios of Q/ΓT on TFTR, and numerical reproduction of ASDEX profiles by a model for turbulent transport invoked without derivation, but by analogy with neoclassical theory. The general equations are specialized to consideration of a number of particular transport mechanisms of interest. (author). Letter-to-the-editor. 10 refs

Seed and soliton solutions for Adler's lattice equation

International Nuclear Information System (INIS)

Atkinson, James; Hietarinta, Jarmo; Nijhoff, Frank

2007-01-01

Adler's lattice equation has acquired the status of a master equation among 2D discrete integrable systems. In this paper we derive what we believe are the first explicit solutions of this equation. In particular it turns out to be necessary to establish a non-trivial seed solution from which soliton solutions can subsequently be constructed using the Baecklund transformation. As a corollary we find the corresponding solutions of the Krichever-Novikov equation which is obtained from Adler's equation in a continuum limit. (fast track communication)

Finite moments approach to the time-dependent neutron transport equation

International Nuclear Information System (INIS)

Kim, Sang Hyun

1994-02-01

Currently, nodal techniques are widely used in solving the multidimensional diffusion equation because of savings in computing time and storage. Thanks to the development of computer technology, one can now solve the transport equation instead of the diffusion equation to obtain more accurate solution. The finite moments method, one of the nodal methods, attempts to represent the fluxes in the cell and on cell surfaces more rigorously by retaining additional spatial moments. Generally, there are two finite moments schemes to solve the time-dependent transport equation. In one, the time variable is treated implicitly with finite moments method in space variable (implicit finite moments method), the other method uses finite moments method in both space and time (space-time finite moments method). In this study, these two schemes are applied to two types of time-dependent neutron transport problems. One is a fixed source problem, the other a heterogeneous fast reactor problem with delayed neutrons. From the results, it is observed that the two finite moments methods give almost the same solutions in both benchmark problems. However, the space-time finite moments method requires a little longer computing time than that of the implicit finite moments method. In order to reduce the longer computing time in the space-time finite moments method, a new iteration strategy is exploited, where a few time-stepwise calculation, in which original time steps are grouped into several coarse time divisions, is performed sequentially instead of performing iterations over the entire time steps. This strategy results in significant reduction of the computing time and we observe that 2-or 3-stepwise calculation is preferable. In addition, we propose a new finite moments method which is called mixed finite moments method in this thesis. Asymptotic analysis for the finite moments method shows that accuracy of the solution in a heterogeneous problem mainly depends on the accuracy of the

Discrete Riccati equation solutions: Distributed algorithms

Directory of Open Access Journals (Sweden)

D. G. Lainiotis

1996-01-01

Full Text Available In this paper new distributed algorithms for the solution of the discrete Riccati equation are introduced. The algorithms are used to provide robust and computational efficient solutions to the discrete Riccati equation. The proposed distributed algorithms are theoretically interesting and computationally attractive.

Normal scheme for solving the transport equation independently of spatial discretization

International Nuclear Information System (INIS)

Zamonsky, O.M.

1993-01-01

To solve the discrete ordinates neutron transport equation, a general order nodal scheme is used, where nodes are allowed to have different orders of approximation and the whole system reaches a final order distribution. Independence in the election of system discretization and order of approximation is obtained without loss of accuracy. The final equations and the iterative method to reach a converged order solution were implemented in a two-dimensional computer code to solve monoenergetic, isotropic scattering, external source problems. Two benchmark problems were solved using different automatic selection order methods. Results show accurate solutions without spatial discretization, regardless of the initial selection of distribution order. (author)

EOS9nT: A TOUGH2 module for the simulation of flow and solute/colloid transport

International Nuclear Information System (INIS)

Moridis, G.J.; Wu, Y.S.; Pruess, K.

1998-04-01

EOS9nT is a new TOUGH2 module for the simulation of flow and transport of an arbitrary number n of tracers (solutes and/or colloids) in the subsurface. The module first solves the flow-related equations, which are comprised of (a) the Richards equation and, depending on conditions, may also include (b) the flow equation of a dense brine or aqueous suspension and/or (c) the heat equation. A second set of transport equations, corresponding to the n tracers, are then solved sequentially. The low concentrations of the n tracers are considered to have no effect on the liquid phase, thus making possible the decoupling of their equations. The first set of equations in EOS9nT provides the flow regime and account for fluid density variations due to thermal and/or solute concentration effects. The n tracer transport equations account for sorption, radioactive decay, advection, hydrodynamic dispersion, molecular diffusion, as well as filtration (for colloids only). EOS9nT can handle gridblocks or irregular geometry in three-dimensional domains. Preliminary results from four 1-D verification problems show an excellent agreement between the numerical predictions and the known analytical solutions

Exactly averaged equations for flow and transport in random media

International Nuclear Information System (INIS)

Shvidler, Mark; Karasaki, Kenzi

2001-01-01

It is well known that exact averaging of the equations of flow and transport in random porous media can be realized only for a small number of special, occasionally exotic, fields. On the other hand, the properties of approximate averaging methods are not yet fully understood. For example, the convergence behavior and the accuracy of truncated perturbation series. Furthermore, the calculation of the high-order perturbations is very complicated. These problems for a long time have stimulated attempts to find the answer for the question: Are there in existence some exact general and sufficiently universal forms of averaged equations? If the answer is positive, there arises the problem of the construction of these equations and analyzing them. There exist many publications related to these problems and oriented on different applications: hydrodynamics, flow and transport in porous media, theory of elasticity, acoustic and electromagnetic waves in random fields, etc. We present a method of finding the general form of exactly averaged equations for flow and transport in random fields by using (1) an assumption of the existence of Green's functions for appropriate stochastic problems, (2) some general properties of the Green's functions, and (3) the some basic information about the random fields of the conductivity, porosity and flow velocity. We present a general form of the exactly averaged non-local equations for the following cases. 1. Steady-state flow with sources in porous media with random conductivity. 2. Transient flow with sources in compressible media with random conductivity and porosity. 3. Non-reactive solute transport in random porous media. We discuss the problem of uniqueness and the properties of the non-local averaged equations, for the cases with some types of symmetry (isotropic, transversal isotropic, orthotropic) and we analyze the hypothesis of the structure non-local equations in general case of stochastically homogeneous fields. (author)

Solution of the two dimensional diffusion and transport equations in a rectangular lattice with an elliptical fuel element using Fourier transform methods: One and two group cases

International Nuclear Information System (INIS)

Williams, M.M.R.; Hall, S.K.; Eaton, M.D.

2014-01-01

Highlights: • A rectangular reactor cell with an elliptical fuel element. • Solution of transport and diffusion equations by Fourier expansion. • Numerical examples showing convergence. • Two group cell problems. - Abstract: A method for solving the diffusion and transport equations in a rectangular lattice cell with an elliptical fuel element has been developed using a Fourier expansion of the neutron flux. The method is applied to a one group model with a source in the moderator. The cell flux is obtained and also the associated disadvantage factor. In addition to the one speed case, we also consider the two group equations in the cell which now become an eigenvalue problem for the lattice multiplication factor. The method of solution relies upon an efficient procedure to solve a large set of simultaneous linear equations and for this we use the IMSL library routines. Our method is compared with the results from a finite element code. The main drawback of the problem arises from the very large number of terms required in the Fourier series which taxes the storage and speed of the computer. Nevertheless, useful solutions are obtained in geometries that would normally require the use of finite element or analogous methods, for this reason the Fourier method is useful for comparison with that type of numerical approach. Extension of the method to more intricate fuel shapes, such as stars and cruciforms as well as superpositions of these, is possible

A numerical spectral approach to solve the dislocation density transport equation

International Nuclear Information System (INIS)

Djaka, K S; Taupin, V; Berbenni, S; Fressengeas, C

2015-01-01

A numerical spectral approach is developed to solve in a fast, stable and accurate fashion, the quasi-linear hyperbolic transport equation governing the spatio-temporal evolution of the dislocation density tensor in the mechanics of dislocation fields. The approach relies on using the Fast Fourier Transform algorithm. Low-pass spectral filters are employed to control both the high frequency Gibbs oscillations inherent to the Fourier method and the fast-growing numerical instabilities resulting from the hyperbolic nature of the transport equation. The numerical scheme is validated by comparison with an exact solution in the 1D case corresponding to dislocation dipole annihilation. The expansion and annihilation of dislocation loops in 2D and 3D settings are also produced and compared with finite element approximations. The spectral solutions are shown to be stable, more accurate for low Courant numbers and much less computation time-consuming than the finite element technique based on an explicit Galerkin-least squares scheme. (paper)

External intermittency prediction using AMR solutions of RANS turbulence and transported PDF models

Science.gov (United States)

Olivieri, D. A.; Fairweather, M.; Falle, S. A. E. G.

2011-12-01

External intermittency in turbulent round jets is predicted using a Reynolds-averaged Navier-Stokes modelling approach coupled to solutions of the transported probability density function (pdf) equation for scalar variables. Solutions to the descriptive equations are obtained using a finite-volume method, combined with an adaptive mesh refinement algorithm, applied in both physical and compositional space. This method contrasts with conventional approaches to solving the transported pdf equation which generally employ Monte Carlo techniques. Intermittency-modified eddy viscosity and second-moment turbulence closures are used to accommodate the effects of intermittency on the flow field, with the influence of intermittency also included, through modifications to the mixing model, in the transported pdf equation. Predictions of the overall model are compared with experimental data on the velocity and scalar fields in a round jet, as well as against measurements of intermittency profiles and scalar pdfs in a number of flows, with good agreement obtained. For the cases considered, predictions based on the second-moment turbulence closure are clearly superior, although both turbulence models give realistic predictions of the bimodal scalar pdfs observed experimentally.

Modeling of water flow and solute transport in unsaturated heterogeneous fields

International Nuclear Information System (INIS)

Bresler, E.; Dagan, G.

1982-01-01

A comprehensive model which considers dispersive solute transport, nonsteady moisture flow regimes and complex boundary conditions is described. The main assumptions are: vertical flow; spatial variability which is associated with the saturated hydraulic conductivity K/sub s/ occurs in the horizontal plane, but is constant in the profile, and has a lognormal probability distribution function (PDF); deterministic recharge and solute concentration are applied during infiltration; the soil is at uniform water content and salt concentration prior to infiltration. The problem is to solve, for arbitrary K/sub s/, the Richards' equation of flow simultaneously with the diffusion-convection equation for salt transport, with the boundary and initial conditions appropriate to infiltration-redistribution. Once this is achieved, the expectation and variance of various quantities of interest (solute concentration, moisture content) are obtained by using the statistical averaging procedure and the given PDF of K/sub s/. Since the solution of Richards' equation for the infiltration-redistribution cycle is extremely difficult (for a given K/sub s/), an approxiate solution is derived by using the concept of piston flow type wetting fronts. Similarly, accurate numerical solutions are used as input for the same statistical averaging procedure. The stochastic model is applied to two spatially variable soils by using both accurate numerical solutions and the simplified water and salt transport models. A comparison between the results shows that the approximate simplified models lead to quite accurate values of the expectations and variances of the flow variables for the entire field. It is suggested that in spatially variable fields, stochastic modeling represents the actual flow phenomena realistically, and provides the main statistical moments by using simplified flow models which can be used with confidence in applications

The Laplace transformation of adjoint transport equations

International Nuclear Information System (INIS)

Hoogenboom, J.E.

1985-01-01

A clarification is given of the difference between the equation adjoint to the Laplace-transformed time-dependent transport equation and the Laplace-transformed time-dependent adjoint transport equation. Proper procedures are derived to obtain the Laplace transform of the instantaneous detector response. (author)

New exact solutions of the Dirac equation

International Nuclear Information System (INIS)

Bagrov, V.G.; Gitman, D.M.; Zadorozhnyj, V.N.; Lavrov, P.M.; Shapovalov, V.N.

1980-01-01

Search for new exact solutions of the Dirac and Klein-Gordon equations are in progress. Considered are general properties of the Dirac equation solutions for an electron in a purely magnetic field, in combination with a longitudinal magnetic and transverse electric fields. New solutions for the equations of charge motion in an electromagnetic field of axial symmetry and in a nonstationary field of a special form have been found for potentials selected concretely

Mixed-hybrid finite element method for the transport equation and diffusion approximation of transport problems; Resolution de l'equation du transport par une methode d'elements finis mixtes-hybrides et approximation par la diffusion de problemes de transport

Energy Technology Data Exchange (ETDEWEB)

Cartier, J

2006-04-15

This thesis focuses on mathematical analysis, numerical resolution and modelling of the transport equations. First of all, we deal with numerical approximation of the solution of the transport equations by using a mixed-hybrid scheme. We derive and study a mixed formulation of the transport equation, then we analyse the related variational problem and present the discretization and the main properties of the scheme. We particularly pay attention to the behavior of the scheme and we show its efficiency in the diffusion limit (when the mean free path is small in comparison with the characteristic length of the physical domain). We present academical benchmarks in order to compare our scheme with other methods in many physical configurations and validate our method on analytical test cases. Unstructured and very distorted meshes are used to validate our scheme. The second part of this thesis deals with two transport problems. The first one is devoted to the study of diffusion due to boundary conditions in a transport problem between two plane plates. The second one consists in modelling and simulating radiative transfer phenomenon in case of the industrial context of inertial confinement fusion. (author)

Exact solutions for modified Korteweg-de Vries equation

International Nuclear Information System (INIS)

Sarma, Jnanjyoti

2009-01-01

Using the simple wave or traveling wave solution technique, many different types of solutions are derived for modified Korteweg-de Vries (KdV) equation. The solutions are obtained from the set of nonlinear algebraic equations, which can be derived from the modified Korteweg-de Vries (KdV) equation by using the hyperbolic transformation method. The method can be applicable for similar nonlinear wave equations.

Solution to the transport equation with anisotropic dispersion in a BWR type assembly using the AZTRAN code

International Nuclear Information System (INIS)

Chepe P, M.; Xolocostli M, J. V.; Gomez T, A. M.; Del Valle G, E.

2016-09-01

Due to the current computing power, the deterministic codes for analyzing nuclear reactors that have been used for several years are becoming more relevant, since much more precise solution techniques can be used; the last century would have been very difficult, since memory and processor capacities were very limited or had high prices on the components. In this work we analyze the effect of the anisotropic dispersion of the effective dispersion section, compared to the isotropic dispersion. The anisotropy implementation was carried out in the AZTRAN transport code, which is part of the AZTLAN platform for nuclear reactors analysis (in development). The AZTRAN code solves the Boltzmann transport equation in one, two and three dimensions at steady state, using the multi-group technique for energy discretization, the RTN-0 nodal method in spatial discretization and for angular discretization the discrete ordinates without considering anisotropy originally. The effect of the anisotropy dispersion on the effective multiplication factor and the axial and radial power on a fuel assembly BWR type are analyzed. (Author)

Integrated compartmental model for describing the transport of solute in a fractured porous medium. [FRACPORT

Energy Technology Data Exchange (ETDEWEB)

DeAngelis, D.L.; Yeh, G.T.; Huff, D.D.

1984-10-01

This report documents a model, FRACPORT, that simulates the transport of a solute through a fractured porous matrix. The model should be useful in analyzing the possible transport of radionuclides from shallow-land burial sites in humid environments. The use of the model is restricted to transport through saturated zones. The report first discusses the general modeling approach used, which is based on the Integrated Compartmental Method. The basic equations of solute transport are then presented. The model, which assumes a known water velocity field, solves these equations on two different time scales; one related to rapid transport of solute along fractures and the other related to slower transport through the porous matrix. FRACPORT is validated by application to a simple example of fractured porous medium transport that has previously been analyzed by other methods. Then its utility is demonstrated in analyzing more complex cases of pulses of solute into a fractured matrix. The report serves as a user's guide to FRACPORT. A detailed description of data input, along with a listing of input for a sample problem, is provided. 16 references, 18 figures, 3 tables.

Travelling wave solutions to the Kuramoto-Sivashinsky equation

International Nuclear Information System (INIS)

Nickel, J.

2007-01-01

Combining the approaches given by Baldwin [Baldwin D et al. Symbolic computation of exact solutions expressible in hyperbolic and elliptic functions for nonlinear PDEs. J Symbol Comput 2004;37:669-705], Peng [Peng YZ. A polynomial expansion method and new general solitary wave solutions to KS equation. Comm Theor Phys 2003;39:641-2] and by Schuermann [Schuermann HW, Serov VS. Weierstrass' solutions to certain nonlinear wave and evolution equations. Proc progress electromagnetics research symposium, 28-31 March 2004, Pisa. p. 651-4; Schuermann HW. Traveling-wave solutions to the cubic-quintic nonlinear Schroedinger equation. Phys Rev E 1996;54:4312-20] leads to a method for finding exact travelling wave solutions of nonlinear wave and evolution equations (NLWEE). The first idea is to generalize ansaetze given by Baldwin and Peng to find elliptic solutions of NLWEEs. Secondly, conditions used by Schuermann to find physical (real and bounded) solutions and to discriminate between periodic and solitary wave solutions are used. The method is shown in detail by evaluating new solutions of the Kuramoto-Sivashinsky equation

Special solutions of neutral functional differential equations

Directory of Open Access Journals (Sweden)

Győri István

2001-01-01

Full Text Available For a system of nonlinear neutral functional differential equations we prove the existence of an -parameter family of "special solutions" which characterize the asymptotic behavior of all solutions at infinity. For retarded functional differential equations the special solutions used in this paper were introduced by Ryabov.

Matrix-oriented implementation for the numerical solution of the partial differential equations governing flows and transport in porous media

KAUST Repository

Sun, Shuyu

2012-09-01

In this paper we introduce a new technique for the numerical solution of the various partial differential equations governing flow and transport phenomena in porous media. This method is proposed to be used in high level programming languages like MATLAB, Python, etc., which show to be more efficient for certain mathematical operations than for others. The proposed technique utilizes those operations in which these programming languages are efficient the most and keeps away as much as possible from those inefficient, time-consuming operations. In particular, this technique is based on the minimization of using multiple indices looping operations by reshaping the unknown variables into one-dimensional column vectors and performing the numerical operations using shifting matrices. The cell-centered information as well as the face-centered information are shifted to the adjacent face-center and cell-center, respectively. This enables the difference equations to be done for all the cells at once using matrix operations rather than within loops. Furthermore, for results post-processing, the face-center information can further be mapped to the physical grid nodes for contour plotting and stream lines constructions. In this work we apply this technique to flow and transport phenomena in porous media. © 2012 Elsevier Ltd.

Complex solutions for generalised fitzhughnagumo equation

International Nuclear Information System (INIS)

Neirameh, A.

2014-01-01

During present investigation, a direct algebraic method on complex solutions of nonlinear partial differential equation is developed and tested in the case of generalized Burgers-Huxley equation. The proposed scheme can be used in a wide class of nonlinear reaction-diffusion equations. These calculations demonstrate that the accuracy of the direct algebraic solutions is quite high even in the case of a small number of grid points. This method is a very reliable, simple, small computation costs, flexible, and convenient alternative method. (author)

New exact solutions of the mBBM equation

International Nuclear Information System (INIS)

Zhang Zhe; Li Desheng

2013-01-01

The enhanced modified simple equation method presented in this article is applied to construct the exact solutions of modified Benjamin-Bona-Mahoney equation. Some new exact solutions are derived by using this method. When some parameters are taken as special values, the solitary wave solutions can be got from the exact solutions. It is shown that the method introduced in this paper has general significance in searching for exact solutions to the nonlinear evolution equations. (authors)

Real solutions to equations from geometry

CERN Document Server

Sottile, Frank

2011-01-01

Understanding, finding, or even deciding on the existence of real solutions to a system of equations is a difficult problem with many applications outside of mathematics. While it is hopeless to expect much in general, we know a surprising amount about these questions for systems which possess additional structure often coming from geometry. This book focuses on equations from toric varieties and Grassmannians. Not only is much known about these, but such equations are common in applications. There are three main themes: upper bounds on the number of real solutions, lower bounds on the number of real solutions, and geometric problems that can have all solutions be real. The book begins with an overview, giving background on real solutions to univariate polynomials and the geometry of sparse polynomial systems. The first half of the book concludes with fewnomial upper bounds and with lower bounds to sparse polynomial systems. The second half of the book begins by sampling some geometric problems for which all ...

Abundant Interaction Solutions of Sine-Gordon Equation

Directory of Open Access Journals (Sweden)

DaZhao Lü

2012-01-01

Full Text Available With the help of computer symbolic computation software (e.g., Maple, abundant interaction solutions of sine-Gordon equation are obtained by means of a constructed Wronskian form expansion method. The method is based upon the forms and structures of Wronskian solutions of sine-Gordon equation, and the functions used in the Wronskian determinants do not satisfy linear partial differential equations. Such interaction solutions are difficultly obtained via other methods. And the method can be automatically carried out in computer.

Discontinuous finite element and characteristics methods for neutrons transport equation solution in heterogeneous grids

International Nuclear Information System (INIS)

Masiello, E.

2006-01-01

The principal goal of this manuscript is devoted to the investigation of a new type of heterogeneous mesh adapted to the shape of the fuel pins (fuel-clad-moderator). The new heterogeneous mesh guarantees the spatial modelling of the pin-cell with a minimum of regions. Two methods are investigated for the spatial discretization of the transport equation: the discontinuous finite element method and the method of characteristics for structured cells. These methods together with the new representation of the pin-cell result in an appreciable reduction of calculation points. They allow an exact modelling of the fuel pin-cell without spatial homogenization. A new synthetic acceleration technique based on an angular multigrid is also presented for the speed up of the inner iterations. These methods are good candidates for transport calculations for a nuclear reactor core. A second objective of this work is the application of method of characteristics for non-structured geometries to the study of double heterogeneity problem. The letters is characterized by fuel material with a stochastic dispersion of heterogeneous grains, and until now was solved with a model based on collision probabilities. We propose a new statistical model based on renewal-Markovian theory, which makes possible to take into account the stochastic nature of the problem and to avoid the approximations of the collision probability model. The numerical solution of this model is guaranteed by the method of characteristics. (author)

Exact solutions of a nonpolynomially nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Parwani, R.; Tan, H.S.

2007-01-01

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

A stochastic model of multiple scattering of charged particles: process, transport equation and solutions

International Nuclear Information System (INIS)

Papiez, L.; Moskvin, V.; Tulovsky, V.

2001-01-01

The process of angular-spatial evolution of multiple scattering of charged particles can be described by a special case of Boltzmann integro-differential equation called Lewis equation. The underlying stochastic process for this evolution is the compound Poisson process on the surface of the unit sphere. The significant portion of events that constitute compound Poisson process that describes multiple scattering have diffusional character. This property allows to analyze the process of angular-spatial evolution of multiple scattering of charged particles as combination of soft and hard collision processes and compute appropriately its transition densities. These computations provide a method of the approximate solution to the Lewis equation. (orig.)

Exact solution for the generalized Telegraph Fisher's equation

International Nuclear Information System (INIS)

Abdusalam, H.A.; Fahmy, E.S.

2009-01-01

In this paper, we applied the factorization scheme for the generalized Telegraph Fisher's equation and an exact particular solution has been found. The exact particular solution for the generalized Fisher's equation was obtained as a particular case of the generalized Telegraph Fisher's equation and the two-parameter solution can be obtained when n=2.

Non-autonomous equations with unpredictable solutions

Science.gov (United States)

Akhmet, Marat; Fen, Mehmet Onur

2018-06-01

To make research of chaos more amenable to investigating differential and discrete equations, we introduce the concepts of an unpredictable function and sequence. The topology of uniform convergence on compact sets is applied to define unpredictable functions [1,2]. The unpredictable sequence is defined as a specific unpredictable function on the set of integers. The definitions are convenient to be verified as solutions of differential and discrete equations. The topology is metrizable and easy for applications with integral operators. To demonstrate the effectiveness of the approach, the existence and uniqueness of the unpredictable solution for a delay differential equation are proved as well as for quasilinear discrete systems. As a corollary of the theorem, a similar assertion for a quasilinear ordinary differential equation is formulated. The results are demonstrated numerically, and an application to Hopfield neural networks is provided. In particular, Poincaré chaos near periodic orbits is observed. The completed research contributes to the theory of chaos as well as to the theory of differential and discrete equations, considering unpredictable solutions.

Generalized fluid equations for parallel transport in collisional to weakly collisional plasmas

International Nuclear Information System (INIS)

Zawaideh, E.; Najmabadi, F.; Conn, R.W.

1986-01-01

A new set of two-fluid equations that are valid from collisional to weakly collisional limits is derived. Starting from gyrokinetic equations in flux coordinates with no zero-order drifts, a set of moment equations describing plasma transport along the field lines of a space- and time-dependent magnetic field is derived. No restriction on the anisotropy of the ion distribution function is imposed. In the highly collisional limit, these equations reduce to those of Braginskii, while in the weakly collisional limit they are similar to the double adiabatic or Chew, Goldberger, and Low (CGL) equations [Proc. R. Soc. London, Ser. A 236, 112 (1956)]. The new set of equations also exhibits a physical singularity at the sound speed. This singularity is used to derive and compute the sound speed. Numerical examples comparing these equations with conventional transport equations show that in the limit where the ratio of the mean free path lambda to the scale length of the magnetic field gradient L/sub B/ approaches zero, there is no significant difference between the solution of the new and conventional transport equations. However, conventional fluid equations, ordinarily expected to be correct to the order (lambda/L/sub B/) 2 , are found to have errors of order (lambda/L/sub u/) 2 = (lambda/L/sub B/) 2 /(1-M 2 ) 2 , where L/sub u/ is the scale length of the flow velocity gradient and M is the Mach number. As such, the conventional equations may contain large errors near the sound speed (Mroughly-equal1)

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

International Nuclear Information System (INIS)

Ren Ji; Ruan Hangyu

2008-01-01

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

Explicit solutions of the Camassa-Holm equation

International Nuclear Information System (INIS)

Parkes, E.J.; Vakhnenko, V.O.

2005-01-01

Explicit travelling-wave solutions of the Camassa-Holm equation are sought. The solutions are characterized by two parameters. For propagation in the positive x-direction, both periodic and solitary smooth-hump, peakon, cuspon and inverted-cuspon waves are found. For propagation in the negative x-direction, there are solutions which are just the mirror image in the x-axis of the aforementioned solutions. Some composite wave solutions of the Degasperis-Procesi equation are given in an appendix

Weak solutions of magma equations

International Nuclear Information System (INIS)

Krishnan, E.V.

1999-01-01

Periodic solutions in terms of Jacobian cosine elliptic functions have been obtained for a set of values of two physical parameters for the magma equation which do not reduce to solitary-wave solutions. It was also obtained solitary-wave solutions for another set of these parameters as an infinite period limit of periodic solutions in terms of Weierstrass and Jacobian elliptic functions

Spatially adaptive hp refinement approach for PN neutron transport equation using spectral element method

International Nuclear Information System (INIS)

Nahavandi, N.; Minuchehr, A.; Zolfaghari, A.; Abbasi, M.

2015-01-01

Highlights: • Powerful hp-SEM refinement approach for P N neutron transport equation has been presented. • The method provides great geometrical flexibility and lower computational cost. • There is a capability of using arbitrary high order and non uniform meshes. • Both posteriori and priori local error estimation approaches have been employed. • High accurate results are compared against other common adaptive and uniform grids. - Abstract: In this work we presented the adaptive hp-SEM approach which is obtained from the incorporation of Spectral Element Method (SEM) and adaptive hp refinement. The SEM nodal discretization and hp adaptive grid-refinement for even-parity Boltzmann neutron transport equation creates powerful grid refinement approach with high accuracy solutions. In this regard a computer code has been developed to solve multi-group neutron transport equation in one-dimensional geometry using even-parity transport theory. The spatial dependence of flux has been developed via SEM method with Lobatto orthogonal polynomial. Two commonly error estimation approaches, the posteriori and the priori has been implemented. The incorporation of SEM nodal discretization method and adaptive hp grid refinement leads to high accurate solutions. Coarser meshes efficiency and significant reduction of computer program runtime in comparison with other common refining methods and uniform meshing approaches is tested along several well-known transport benchmarks

Finite element approximation to the even-parity transport equation

International Nuclear Information System (INIS)

Lewis, E.E.

1981-01-01

This paper studies the finite element method, a procedure for reducing partial differential equations to sets of algebraic equations suitable for solution on a digital computer. The differential equation is cast into the form of a variational principle, the resulting domain then subdivided into finite elements. The dependent variable is then approximated by a simple polynomial, and these are linked across inter-element boundaries by continuity conditions. The finite element method is tailored to a variety of transport problems. Angular approximations are formulated, and the extent of ray effect mitigation is examined. Complex trial functions are introduced to enable the inclusion of buckling approximations. The ubiquitous curved interfaces of cell calculations, and coarse mesh methods are also treated. A concluding section discusses limitations of the work to date and suggests possible future directions

Exact traveling wave solutions of the Boussinesq equation

International Nuclear Information System (INIS)

Ding Shuangshuang; Zhao Xiqiang

2006-01-01

The repeated homogeneous balance method is used to construct exact traveling wave solutions of the Boussinesq equation, in which the homogeneous balance method is applied to solve the Riccati equation and the reduced nonlinear ordinary differential equation, respectively. Many new exact traveling wave solutions of the Boussinesq equation are successfully obtained

Analytical Solutions of the KDV-KZK Equation

Science.gov (United States)

Gan, W. S.

The KdV-KZK equation for fluids developed by me was presented at the ICSV 11 in St. Petersburg in July 2004. In this paper, I made an attempt on the analytical solutions of this equation using the perturbation method. Some physical interpretation of the solutions is given. A brief introduction to KdV-KZK equation for solids is given

Particlelike solutions of the Einstein-Dirac equations

Science.gov (United States)

Finster, Felix; Smoller, Joel; Yau, Shing-Tung

1999-05-01

The coupled Einstein-Dirac equations for a static, spherically symmetric system of two fermions in a singlet spinor state are derived. Using numerical methods, we construct an infinite number of solitonlike solutions of these equations. The stability of the solutions is analyzed. For weak coupling (i.e., small rest mass of the fermions), all the solutions are linearly stable (with respect to spherically symmetric perturbations), whereas for stronger coupling, both stable and unstable solutions exist. For the physical interpretation, we discuss how the energy of the fermions and the (ADM) mass behave as functions of the rest mass of the fermions. Although gravitation is not renormalizable, our solutions of the Einstein-Dirac equations are regular and well behaved even for strong coupling.

A single continuum approximation of the solute transport in fractured porous media

International Nuclear Information System (INIS)

Jeong, J.T.; Lee, K.J.

1992-01-01

Solute transport in fractured porous media is described by the single continuum model, i.e., equivalent porous medium model. In this model, one-dimensional solute transport in the fracture and two-dimensional solute transport in the porous rock matrix is considered. The network of fractures embedded in the porous rock matrix is idealized as two orthogonally intersecting families of equally spaced, parallel fractures directed at 45 o to the regional groundwater flow direction. Governing equations are solved by the finite element method, and an upstream weighting technique is used in order to prevent the oscillation of the solution in the case of highly advection dominated transport. Breakthrough curves, similar to those of the one-dimensional solute transport problem in ordinary porous media, are obtained as a function of time according to volume or flux averaging of the concentration profile across the width of the flow region. The equivalent parameters, i.e., porosity and overall coefficient of longitudinal dispersivity, are obtained by a trial-and-error method. Analyses for the non-sorbing solute transport case show that within the range of considered parameters, and except for the region very close to the source, application of the single continuum model in the idealized fracture system is sufficient for modeling solute transport in fractured porous media. This numerical scheme is shown to be applicable to a sorbing solute and radionuclide transport. (author)

Solute transport in aggregated and layered porous media

International Nuclear Information System (INIS)

Koch, S.

1993-01-01

This work is a contribution to research in soil physics dealing with solute transport in porous media. The influence of structural inhomogeneities on solute transport is investigated. Detailed experiments at the laboratory scale are used to enlighten distinct processes which cannot be studied separately at field scale. Two main aspects are followed up: (i) to show the influence of aggregation of a porous medium on breakthrough time and spreading of an inert tracer and consequences on the estimation of parameter values of models describing solute transport in aggregated systems, (ii) to investigate the influences on the dispersion process when stratification is perpendicular to the direction of flow. Several concepts of modelling solute transport in soil are discussed. Models based on the convection-dispersion equation (CDE) are emphasized because they are used here to model solute transport experiments conducted with aggregated porous media. Stochastic concepts are introduced to show the limitations of the deterministic CDE approaches. Experiments are done in columns containing two kinds of solid phases and were saturated with water. The solid phases are porous and solid glass beads exhibiting a distinctly unimodal or bimodal pore size distribution. Experimental breakthrough curves (BTCs) are modelled with the CDE, a bicontinuum model with a phenomenological mass transfer rate and a bicontinuum spherical diffusion model. Experiments are also done in columns that are unsaturated containing porous materials that are layered. Flow is made at a steady rate. It is shown that layer boundaries have a severe influence on lateral mixing. They may force streamlines to converge or cause a lateral redistribution of solutes. (author) figs., tabs., 122 refs

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

OpenAIRE

Efimova, Olga Yu.

2010-01-01

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

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

International Nuclear Information System (INIS)

Schwenk, G.A. Jr.

1980-01-01

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

Green's function method for the monoenergetic transport equation in heterogeneous plane geometry

International Nuclear Information System (INIS)

Ganapol, B.D.

1995-01-01

For the past several years, a series of papers by the transport group at the University of Arizona dealing with benchmark solutions of the monoenergetic transport equation has appeared. The approach has been to take advantage of highly successful numerical Laplace Fourier transform inversions to provide benchmark quality solutions in infinite media, half-space in one and two dimensions and in homogeneous slabs. This paper extends the set of solutions to include heterogeneous slab geometry by using the recently established Green's Function Method (GFM). Analytical benchmark solutions are an essential part of the quality control of computational algorithms developed for particle transport. In addition, benchmarking methods have applications in the classroom by providing examples of how computational mathematics is used to solve physical problems to obtain meaningful answers. In a structural context, monoenergetic solutions are directly applicable to the investigation of the microlight environment within a leaf. The leaf is considered to be a composition of alternating layers of highly absorbing pigments and water superimposed on a refractively scattering background

Symmetry and exact solutions of nonlinear spinor equations

International Nuclear Information System (INIS)

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

1989-01-01

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

New compacton solutions and solitary wave solutions of fully nonlinear generalized Camassa-Holm equations

International Nuclear Information System (INIS)

Tian Lixin; Yin Jiuli

2004-01-01

In this paper, we introduce the fully nonlinear generalized Camassa-Holm equation C(m,n,p) and by using four direct ansatzs, we obtain abundant solutions: compactons (solutions with the absence of infinite wings), solitary patterns solutions having infinite slopes or cups, solitary waves and singular periodic wave solutions and obtain kink compacton solutions and nonsymmetry compacton solutions. We also study other forms of fully nonlinear generalized Camassa-Holm equation, and their compacton solutions are governed by linear equations

New traveling wave solutions to AKNS and SKdV equations

International Nuclear Information System (INIS)

Ozer, Teoman

2009-01-01

We analyze the traveling wave solutions of Ablowitz-Kaup-Newell-Segur (AKNS) and Schwarz-Korteweg-de Vries (SKdV) equations. As the solution method for differential equations we consider the improved tanh approach. This approach provides to transform the partial differential equation into the ordinary differential equation and then obtain the new families of exact solutions based on the solutions of the Riccati equation. The different values of the coefficients of the Riccati equation allow us to obtain new type of traveling wave solutions to AKNS and SKdV equations.

Mathematical description of adsorption and transport of reactive solutes in soil: a review of selected literature

International Nuclear Information System (INIS)

Travis, C.C.

1978-10-01

This report reviews selected literature related to the mathematical description of the transport of reactive solutes through soil. The primary areas of the literature reviewed are (1) mathematical models in current use for description of the adsorption-desorption interaction between the soil solution and the soil matrix and (2) analytic solutions of the differential equations describing the convective-dispersive transport of reactive solutes through soil

New exact solutions of the Dirac equation. 11

International Nuclear Information System (INIS)

Bagrov, V.G.; Noskov, M.D.

1984-01-01

Investigations into determining new exact solutions of relativistic wave equations started in another paper were continued. Exact solutions of the Dirac, Klein-Gordon equations and classical relativistic equations of motion in four new types of external electromagnetic fields were found

Inelastic Quantum Transport in Superlattices: Success and Failure of the Boltzmann Equation

DEFF Research Database (Denmark)

Wacker, Andreas; Jauho, Antti-Pekka; Rott, Stephan

1999-01-01

the whole held range from linear response to negative differential conductivity. The quantum results are compared with the respective results obtained from a Monte Carlo solution of the Boltzmann equation. Our analysis thus sets the limits of validity for the semiclassical theory in a nonlinear transport...

Development of a polynomial nodal model to the multigroup transport equation in one dimension

International Nuclear Information System (INIS)

Feiz, M.

1986-01-01

A polynomial nodal model that uses Legendre polynomial expansions was developed for the multigroup transport equation in one dimension. The development depends upon the least-squares minimization of the residuals using the approximate functions over the node. Analytical expressions were developed for the polynomial coefficients. The odd moments of the angular neutron flux over the half ranges were used at the internal interfaces, and the Marshak boundary condition was used at the external boundaries. Sample problems with fine-mesh finite-difference solutions of the diffusion and transport equations were used for comparison with the model

Adjoint P1 equations solution for neutron slowing down; Solucao das equacoes P1 adjuntas para moderacao de neutrons

Energy Technology Data Exchange (ETDEWEB)

Cardoso, Carlos Eduardo Santos; Martinez, Aquilino Senra; Silva, Fernando Carvalho da [Universidade Federal, Rio de Janeiro, RJ (Brazil). Coordenacao dos Programas de Pos-graduacao de Engenharia. Programa de Engenharia Nuclear

2002-07-01

In some applications of perturbation theory, it is necessary know the adjoint neutron flux, which is obtained by the solution of adjoint neutron diffusion equation. However, the multigroup constants used for this are weighted in only the direct neutron flux, from the solution of direct P1 equations. In this work, the adjoint P1 equations are derived by the neutron transport equation, the reversion operators rules and analogies between direct and adjoint parameters. The direct and adjoint neutron fluxes resulting from the solution of P{sub 1} equations were used to three different weighting processes, to obtain the macrogroup macroscopic cross sections. It was found out noticeable differences among them. (author)

Quadratic inner element subgrid scale discretisation of the Boltzmann transport equation

International Nuclear Information System (INIS)

Baker, C.M.J.; Buchan, A.G.; Pain, C.C.; Tollit, B.; Eaton, M.D.; Warner, P.

2012-01-01

This paper explores the application of the inner element subgrid scale method to the Boltzmann transport equation using quadratic basis functions. Previously, only linear basis functions for both the coarse scale and the fine scale were considered. This paper, therefore, analyses the advantages of using different coarse and subgrid basis functions for increasing the accuracy of the subgrid scale method. The transport of neutral particle radiation may be described by the Boltzmann transport equation (BTE) which, due to its 7 dimensional phase space, is computationally expensive to resolve. Multi-scale methods offer an approach to efficiently resolve the spatial dimensions of the BTE by separating the solution into its coarse and fine scales and formulating a solution whereby only the computationally efficient coarse scales need to be solved. In previous work an inner element subgrid scale method was developed that applied a linear continuous and discontinuous finite element method to represent the solution’s coarse and fine scale components. This approach was shown to generate efficient and stable solutions, and so this article continues its development by formulating higher order quadratic finite element expansions over the continuous and discontinuous scales. Here it is shown that a solution’s convergence can be improved significantly using higher order basis functions. Furthermore, by using linear finite elements to represent coarse scales in combination with quadratic fine scales, convergence can also be improved with only a modest increase in computational expense.

Semianalytic Solution of Space-Time Fractional Diffusion Equation

Directory of Open Access Journals (Sweden)

A. Elsaid

2016-01-01

Full Text Available We study the space-time fractional diffusion equation with spatial Riesz-Feller fractional derivative and Caputo fractional time derivative. The continuation of the solution of this fractional equation to the solution of the corresponding integer order equation is proved. The series solution of this problem is obtained via the optimal homotopy analysis method (OHAM. Numerical simulations are presented to validate the method and to show the effect of changing the fractional derivative parameters on the solution behavior.

New solutions of Heun's general equation

Energy Technology Data Exchange (ETDEWEB)

Ishkhanyan, Artur [Engineering Center of Armenian National Academy of Sciences, Ashtarak (Armenia); Suominen, Kalle-Antti [Helsinki Institute of Physics, PL 64, Helsinki (Finland)

2003-02-07

We show that in four particular cases the derivative of the solution of Heun's general equation can be expressed in terms of a solution to another Heun's equation. Starting from this property, we use the Gauss hypergeometric functions to construct series solutions to Heun's equation for the mentioned cases. Each of the hypergeometric functions involved has correct singular behaviour at only one of the singular points of the equation; the sum, however, has correct behaviour. (letter to the editor)

Exact solutions to some modified sine-Gordon equations

International Nuclear Information System (INIS)

Saermark, K.

1983-01-01

Exact, translational solutions to a number of modified sine-Gordon equations are presented. In deriving the equations and the solutions use is made of results from the theory of ordinary differential equations without moving critical points as given by Ince. It is found that kink-like solutions exist also in cases where the coefficients of the trigonometric terms are space- and time-dependent. (Auth.)

Assessing lateral flows and solute transport during floods in a conduit-flow-dominated karst system using the inverse problem for the advection-diffusion equation

Science.gov (United States)

Cholet, Cybèle; Charlier, Jean-Baptiste; Moussa, Roger; Steinmann, Marc; Denimal, Sophie

2017-07-01

The aim of this study is to present a framework that provides new ways to characterize the spatio-temporal variability of lateral exchanges for water flow and solute transport in a karst conduit network during flood events, treating both the diffusive wave equation and the advection-diffusion equation with the same mathematical approach, assuming uniform lateral flow and solute transport. A solution to the inverse problem for the advection-diffusion equations is then applied to data from two successive gauging stations to simulate flows and solute exchange dynamics after recharge. The study site is the karst conduit network of the Fourbanne aquifer in the French Jura Mountains, which includes two reaches characterizing the network from sinkhole to cave stream to the spring. The model is applied, after separation of the base from the flood components, on discharge and total dissolved solids (TDSs) in order to assess lateral flows and solute concentrations and compare them to help identify water origin. The results showed various lateral contributions in space - between the two reaches located in the unsaturated zone (R1), and in the zone that is both unsaturated and saturated (R2) - as well as in time, according to hydrological conditions. Globally, the two reaches show a distinct response to flood routing, with important lateral inflows on R1 and large outflows on R2. By combining these results with solute exchanges and the analysis of flood routing parameters distribution, we showed that lateral inflows on R1 are the addition of diffuse infiltration (observed whatever the hydrological conditions) and localized infiltration in the secondary conduit network (tributaries) in the unsaturated zone, except in extreme dry periods. On R2, despite inflows on the base component, lateral outflows are observed during floods. This pattern was attributed to the concept of reversal flows of conduit-matrix exchanges, inducing a complex water mixing effect in the saturated zone

Assessing lateral flows and solute transport during floods in a conduit-flow-dominated karst system using the inverse problem for the advection–diffusion equation

Directory of Open Access Journals (Sweden)

C. Cholet

2017-07-01

Full Text Available The aim of this study is to present a framework that provides new ways to characterize the spatio-temporal variability of lateral exchanges for water flow and solute transport in a karst conduit network during flood events, treating both the diffusive wave equation and the advection–diffusion equation with the same mathematical approach, assuming uniform lateral flow and solute transport. A solution to the inverse problem for the advection–diffusion equations is then applied to data from two successive gauging stations to simulate flows and solute exchange dynamics after recharge. The study site is the karst conduit network of the Fourbanne aquifer in the French Jura Mountains, which includes two reaches characterizing the network from sinkhole to cave stream to the spring. The model is applied, after separation of the base from the flood components, on discharge and total dissolved solids (TDSs in order to assess lateral flows and solute concentrations and compare them to help identify water origin. The results showed various lateral contributions in space – between the two reaches located in the unsaturated zone (R1, and in the zone that is both unsaturated and saturated (R2 – as well as in time, according to hydrological conditions. Globally, the two reaches show a distinct response to flood routing, with important lateral inflows on R1 and large outflows on R2. By combining these results with solute exchanges and the analysis of flood routing parameters distribution, we showed that lateral inflows on R1 are the addition of diffuse infiltration (observed whatever the hydrological conditions and localized infiltration in the secondary conduit network (tributaries in the unsaturated zone, except in extreme dry periods. On R2, despite inflows on the base component, lateral outflows are observed during floods. This pattern was attributed to the concept of reversal flows of conduit–matrix exchanges, inducing a complex water mixing effect

Lump solutions to the Kadomtsev–Petviashvili equation

Energy Technology Data Exchange (ETDEWEB)

Ma, Wen-Xiu, E-mail: mawx@cas.usf.edu

2015-09-25

Through symbolic computation with Maple, a class of lump solutions, rationally localized in all directions in the space, to the (2 + 1)-dimensional Kadomtsev–Petviashvili (KP) equation is presented, making use of its Hirota bilinear form. The resulting lump solutions contain six free parameters, two of which are due to the translation invariance of the KP equation and the other four of which satisfy a non-zero determinant condition guaranteeing analyticity and rational localization of the solutions. Three contour plots with different determinant values are sequentially made to show that the corresponding lump solution tends to zero when the determinant approaches zero. Two particular lump solutions with specific values of the involved parameters are plotted, as illustrative examples. - Highlights: • Positive quadratic function solutions. • Solitons rationally-localized in all directions in the space. • Solving systems of nonlinear algebraic equations by symbolic computation with Maple.

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)

Deterministic methods to solve the integral transport equation in neutronic

International Nuclear Information System (INIS)

Warin, X.

1993-11-01

We present a synthesis of the methods used to solve the integral transport equation in neutronic. This formulation is above all used to compute solutions in 2D in heterogeneous assemblies. Three kinds of methods are described: - the collision probability method; - the interface current method; - the current coupling collision probability method. These methods don't seem to be the most effective in 3D. (author). 9 figs

New multiple soliton solutions to the general Burgers-Fisher equation and the Kuramoto-Sivashinsky equation

International Nuclear Information System (INIS)

Chen Huaitang; Zhang Hongqing

2004-01-01

A generalized tanh function method is used for constructing exact travelling wave solutions of nonlinear partial differential equations in a unified way. The main idea of this method is to take full advantage of the Riccati equation which has more new solutions. More new multiple soliton solutions are obtained for the general Burgers-Fisher equation and the Kuramoto-Sivashinsky equation

New exact solutions of the Dirac equation. 8

International Nuclear Information System (INIS)

Bagrov, V.G.; Gitman, D.M.; Zadorozhnyj, V.N.; Sukhomlin, N.B.; Shapovalov, V.N.

1978-01-01

The paper continues the investigation into the exact solutions of the Dirac, Klein-Gordon, and Lorentz equations for a charge in an external electromagnetic field. The fields studied do not allow for separation of variables in the Dirac equation, but solutions to the Dirac equation are obtained

Rational Solutions and Lump Solutions of the Potential YTSF Equation

Science.gov (United States)

Sun, Hong-Qian; Chen, Ai-Hua

2017-07-01

By using of the bilinear form, rational solutions and lump solutions of the potential Yu-Toda-Sasa-Fukuyama (YTSF) equation are derived. Dynamics of the fundamental lump solution, n1-order lump solutions, and N-lump solutions are studied for some special cases. We also find some interaction behaviours of solitary waves and one lump of rational solutions.

Asymptotic formulae for solutions of the two-group integral neutron-transport equation

International Nuclear Information System (INIS)

Duracz, T.

1976-01-01

The steady-state, two-group integral neutron-transport equation is considered for two cases. First, for plane geometry, formulae for the asymptotic flux are obtained, under assumptions of homogeneous medium with isotropic scattering, extended to infinity (whole space and half-space), with sources vanishing at infinity as 0(esup(-IXI)). Next, for spherical geometry, the Milne problem is considered and formulae for the asymptotic flux are obtained. These formulae have the form of asymptotic expansions for small and large radii of the black sphere. (orig.) [de

Solution of neutron transport equation using Daubechies' wavelet expansion in the angular discretization

International Nuclear Information System (INIS)

Cao Liangzhi; Wu Hongchun; Zheng Youqi

2008-01-01

Daubechies' wavelet expansion is introduced to discretize the angular variables of the neutron transport equation when the neutron angular flux varies very acutely with the angular directions. An improvement is made by coupling one-dimensional wavelet expansion and discrete ordinate method to make two-dimensional angular discretization efficient and stable. The angular domain is divided into several subdomains for treating the vacuum boundary condition exactly in the unstructured geometry. A set of wavelet equations coupled with each other is obtained in each subdomain. An iterative method is utilized to decouple the wavelet moments. The numerical results of several benchmark problems demonstrate that the wavelet expansion method can provide more accurate results by lower-order expansion than other angular discretization methods

New analytic solutions of stochastic coupled KdV equations

International Nuclear Information System (INIS)

Dai Chaoqing; Chen Junlang

2009-01-01

In this paper, firstly, we use the exp-function method to seek new exact solutions of the Riccati equation. Then, with the help of Hermit transformation, we employ the Riccati equation and its new exact solutions to find new analytic solutions of the stochastic coupled KdV equation in the white noise environment. As some special examples, some analytic solutions can degenerate into these solutions reported in open literatures.

An Experimental Study on Solute Transport in One-Dimensional Clay Soil Columns

Directory of Open Access Journals (Sweden)

Muhammad Zaheer

2017-01-01

Full Text Available Solute transport in low-permeability media such as clay has not been studied carefully up to present, and we are often unclear what the proper governing law is for describing the transport process in such media. In this study, we composed and analyzed the breakthrough curve (BTC data and the development of leaching in one-dimensional solute transport experiments in low-permeability homogeneous and saturated media at small scale, to identify key parameters controlling the transport process. Sodium chloride (NaCl was chosen to be the tracer. A number of tracer tests were conducted to inspect the transport process under different conditions. The observed velocity-time behavior for different columns indicated the decline of soil permeability when switching from tracer introducing to tracer flushing. The modeling approaches considered were the Advection-Dispersion Equation (ADE, Two-Region Model (TRM, Continuous Time Random Walk (CTRW, and Fractional Advection-Dispersion Equation (FADE. It was found that all the models can fit the transport process very well; however, ADE and TRM were somewhat unable to characterize the transport behavior in leaching. The CTRW and FADE models were better in capturing the full evaluation of tracer-breakthrough curve and late-time tailing in leaching.

Some results on the neutron transport and the coupling of equations

International Nuclear Information System (INIS)

Bal, G.

1997-01-01

Neutron transport in nuclear reactors is well modeled by the linear Boltzmann transport equation. Its resolution is relatively easy but very expensive. To achieve whole core calculations, one has to consider simpler models, such as diffusion or homogeneous transport equations. However, the solutions may become inaccurate in particular situations (as accidents for instance). That is the reason why we wish to solve the equations on small area accurately and more coarsely on the remaining part of the core. It is than necessary to introduce some links between different discretizations or modelizations. In this note, we give some results on the coupling of different discretizations of all degrees of freedom of the integral-differential neutron transport equation (two degrees for the angular variable, on for the energy component, and two or three degrees for spatial position respectively in 2D (cylindrical symmetry) and 3D). Two chapters are devoted to the coupling of discrete ordinates methods (for angular discretization). The first one is theoretical and shows the well posing of the coupled problem, whereas the second one deals with numerical applications of practical interest (the results have been obtained from the neutron transport code developed at the R and D, which has been modified for introducing the coupling). Next, we present the nodal scheme RTN0, used for the spatial discretization. We show well posing results for the non-coupled and the coupled problems. At the end, we deal with the coupling of energy discretizations for the multigroup equations obtained by homogenization. Some theoretical results of the discretization of the velocity variable (well-posing of problems), which do not deal directly with the purposes of coupling, are presented in the annexes. (author)

Neutron transport equation - indications on homogenization and neutron diffusion

International Nuclear Information System (INIS)

Argaud, J.P.

1992-06-01

In PWR nuclear reactor, the practical study of the neutrons in the core uses diffusion equation to describe the problem. On the other hand, the most correct method to describe these neutrons is to use the Boltzmann equation, or neutron transport equation. In this paper, we give some theoretical indications to obtain a diffusion equation from the general transport equation, with some simplifying hypothesis. The work is organised as follows: (a) the most general formulations of the transport equation are presented: integro-differential equation and integral equation; (b) the theoretical approximation of this Boltzmann equation by a diffusion equation is introduced, by the way of asymptotic developments; (c) practical homogenization methods of transport equation is then presented. In particular, the relationships with some general and useful methods in neutronic are shown, and some homogenization methods in energy and space are indicated. A lot of other points of view or complements are detailed in the text or the remarks

Bayesian estimation of the hydraulic and solute transport properties of a small-scale unsaturated soil column

Directory of Open Access Journals (Sweden)

Moreira Paulo H. S.

2016-03-01

Full Text Available In this study the hydraulic and solute transport properties of an unsaturated soil were estimated simultaneously from a relatively simple small-scale laboratory column infiltration/outflow experiment. As governing equations we used the Richards equation for variably saturated flow and a physical non-equilibrium dual-porosity type formulation for solute transport. A Bayesian parameter estimation approach was used in which the unknown parameters were estimated with the Markov Chain Monte Carlo (MCMC method through implementation of the Metropolis-Hastings algorithm. Sensitivity coefficients were examined in order to determine the most meaningful measurements for identifying the unknown hydraulic and transport parameters. Results obtained using the measured pressure head and solute concentration data collected during the unsaturated soil column experiment revealed the robustness of the proposed approach.

Spurious solutions in few-body equations. II. Numerical investigations

International Nuclear Information System (INIS)

Adhikari, S.K.

1979-01-01

A recent analytic study of spurious solutions in few-body equations by Adhikari and Gloeckle is here complemented by numerical investigations. As proposed by Adhikari and Gloeckle we study numerically the spurious solutions in the three-body Weinberg type equations and draw some general conclusions about the existence of spurious solutions in three-body equations with the Weinberg kernel and in other few-body formulations. In particular we conclude that for most of the potentials we encounter in problems of nuclear physics the three-body Weinberg type equation will not have a spurious solution which may interfere with the bound state or scattering calculation. Hence, if proven convenient, the three-body Weinberg type equation can be used in practical calculations. The same conclusion is true for the three-body channel coupling array scheme of Kouri, Levin, and Tobocman. In the case of the set of six coupled four-body equations proposed by Rosenberg et al. and the set of the Bencze-Redish-Sloan equations a careful study of the possible spurious solutions is needed before using these equations in practical calculations

Exact, multiple soliton solutions of the double sine Gordon equation

International Nuclear Information System (INIS)

Burt, P.B.

1978-01-01

Exact, particular solutions of the double sine Gordon equation in n dimensional space are constructed. Under certain restrictions these solutions are N solitons, where N <= 2q - 1 and q is the dimensionality of space-time. The method of solution, known as the base equation technique, relates solutions of nonlinear partial differential equations to solutions of linear partial differential equations. This method is reviewed and its applicability to the double sine Gordon equation shown explicitly. The N soliton solutions have the remarkable property that they collapse to a single soliton when the wave vectors are parallel. (author)

Exact Solutions to a Combined sinh-cosh-Gordon Equation

International Nuclear Information System (INIS)

Wei Long

2010-01-01

Based on a transformed Painleve property and the variable separated ODE method, a function transformation method is proposed to search for exact solutions of some partial differential equations (PDEs) with hyperbolic or exponential functions. This approach provides a more systematical and convenient handling of the solution process of this kind of nonlinear equations. Its key point is to eradicate the hyperbolic or exponential terms by a transformed Painleve property and reduce the given PDEs to a variable-coefficient ordinary differential equations, then we seek for solutions to the resulting equations by some methods. As an application, exact solutions for the combined sinh-cosh-Gordon equation are formally derived. (general)

Finite element method solution of simplified P3 equation for flexible geometry handling

International Nuclear Information System (INIS)

Ryu, Eun Hyun; Joo, Han Gyu

2011-01-01

In order to obtain efficiently core flux solutions which would be much closer to the transport solution than the diffusion solution is, not being limited by the geometry of the core, the simplified P 3 (SP 3 ) equation is solved with the finite element method (FEM). A generic mesh generator, GMSH, is used to generate linear and quadratic mesh data. The linear system resulting from the SP 3 FEM discretization is solved by Krylov subspace methods (KSM). A symmetric form of the SP 3 equation is derived to apply the conjugate gradient method rather than the KSMs for nonsymmetric linear systems. An optional iso-parametric quadratic mapping scheme, which is to selectively model nonlinear shapes with a quadratic mapping to prevent significant mismatch in local domain volume, is also implemented for efficient application of arbitrary geometry handling. The gain in the accuracy attainable by the SP 3 solution over the diffusion solution is assessed by solving numerous benchmark problems having various core geometries including the IAEA PWR problems involving rectangular fuels and the Takeda fast reactor problems involving hexagonal fuels. The reference transport solution is produced by the McCARD Monte Carlo code and the multiplication factor and power distribution errors are assessed. In addition, the effect of quadratic mapping is examined for circular cell problems. It is shown that significant accuracy gain is possible with the SP 3 solution for the fast reactor problems whereas only marginal improvement is noted for thermal reactor problems. The quadratic mapping is also quite effective handling geometries with curvature. (author)

Analytical Solution of Pantograph Equation with Incommensurate Delay

Science.gov (United States)

Patade, Jayvant; Bhalekar, Sachin

2017-08-01

Pantograph equation is a delay differential equation (DDE) arising in electrodynamics. This paper studies the pantograph equation with two delays. The existence, uniqueness, stability and convergence results for DDEs are presented. The series solution of the proposed equation is obtained by using Daftardar-Gejji and Jafari method and given in terms of a special function. This new special function has several properties and relations with other functions. Further, we generalize the proposed equation to fractional-order case and obtain its solution.

Solution of the Neutron transport equation in hexagonal geometry using strongly discontinuous nodal schemes

International Nuclear Information System (INIS)

Mugica R, C.A.; Valle G, E. del

2005-01-01

In 2002, E. del Valle and Ernest H. Mund developed a technique to solve numerically the Neutron transport equations in discrete ordinates and hexagonal geometry using two nodal schemes type finite element weakly discontinuous denominated WD 5,3 and WD 12,8 (of their initials in english Weakly Discontinuous). The technique consists on representing each hexagon in the union of three rhombuses each one of which it is transformed in a square in the one that the methods WD 5,3 and WD 12,8 were applied. In this work they are solved the mentioned equations of transport using the same discretization technique by hexagon but using two nodal schemes type finite element strongly discontinuous denominated SD 3 and SD 8 (of their initials in english Strongly Discontinuous). The application in each case as well as a reference problem for those that results are provided for the effective multiplication factor is described. It is carried out a comparison with the obtained results by del Valle and Mund for different discretization meshes so much angular as spatial. (Author)

Intuitive Understanding of Solutions of Partially Differential Equations

Science.gov (United States)

Kobayashi, Y.

2008-01-01

This article uses diagrams that help the observer see how solutions of the wave equation and heat conduction equation are obtained. The analytical approach cannot necessarily show the mechanisms of the key to the solution without transforming the differential equation into a more convenient form by separation of variables. The visual clues based…

Positive solution of a time and energy dependent neutron transport problem

International Nuclear Information System (INIS)

Pao, C.V.

1975-01-01

A constructive method is given for the determination of a solution and an existence--uniqueness theorem for some nonlinear time and energy dependent neutron transport problems, including the linear transport system. The geometry of the medium under consideration is allowed to be either bounded or unbounded which includes the geometry of a finite or infinite cylinder, a half-space and the whole space R/subm/ (m=1,2,center-dotcenter-dotcenter-dot). Our approach to the problem is by successive approximation which leads to various recursion formulas for the approximations in terms of explicit integrations. It is shown under some Lipschitz conditions on the nonlinear functions, which describe the process of neutrons absorption, fission, and scattering, that the sequence of approximations converges to a unique positive solution. Since these conditions are satisfied by the linear transport equation, all the results for the nonlinear system are valid for the linear transport problem. In the general nonlinear problem, the existence of both local and global solutions are discussed, and an iterative process for the construction of the solution is given

SU-G-TeP1-15: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation

Energy Technology Data Exchange (ETDEWEB)

Yang, R [University of Alberta, Edmonton, AB (Canada); Fallone, B [University of Alberta, Edmonton, AB (Canada); Cross Cancer Institute, Edmonton, AB (Canada); MagnetTx Oncology Solutions, Edmonton, AB (Canada); St Aubin, J [University of Alberta, Edmonton, AB (Canada); Cross Cancer Institute, Edmonton, AB (Canada)

2016-06-15

Purpose: To develop a Graphic Processor Unit (GPU) accelerated deterministic solution to the Linear Boltzmann Transport Equation (LBTE) for accurate dose calculations in radiotherapy (RT). A deterministic solution yields the potential for major speed improvements due to the sparse matrix-vector and vector-vector multiplications and would thus be of benefit to RT. Methods: In order to leverage the massively parallel architecture of GPUs, the first order LBTE was reformulated as a second order self-adjoint equation using the Least Squares Finite Element Method (LSFEM). This produces a symmetric positive-definite matrix which is efficiently solved using a parallelized conjugate gradient (CG) solver. The LSFEM formalism is applied in space, discrete ordinates is applied in angle, and the Multigroup method is applied in energy. The final linear system of equations produced is tightly coupled in space and angle. Our code written in CUDA-C was benchmarked on an Nvidia GeForce TITAN-X GPU against an Intel i7-6700K CPU. A spatial mesh of 30,950 tetrahedral elements was used with an S4 angular approximation. Results: To avoid repeating a full computationally intensive finite element matrix assembly at each Multigroup energy, a novel mapping algorithm was developed which minimized the operations required at each energy. Additionally, a parallelized memory mapping for the kronecker product between the sparse spatial and angular matrices, including Dirichlet boundary conditions, was created. Atomicity is preserved by graph-coloring overlapping nodes into separate kernel launches. The one-time mapping calculations for matrix assembly, kronecker product, and boundary condition application took 452±1ms on GPU. Matrix assembly for 16 energy groups took 556±3s on CPU, and 358±2ms on GPU using the mappings developed. The CG solver took 93±1s on CPU, and 468±2ms on GPU. Conclusion: Three computationally intensive subroutines in deterministically solving the LBTE have been

A network thermodynamic method for numerical solution of the Nernst-Planck and Poisson equation system with application to ionic transport through membranes.

Science.gov (United States)

Horno, J; González-Caballero, F; González-Fernández, C F

1990-01-01

Simple techniques of network thermodynamics are used to obtain the numerical solution of the Nernst-Planck and Poisson equation system. A network model for a particular physical situation, namely ionic transport through a thin membrane with simultaneous diffusion, convection and electric current, is proposed. Concentration and electric field profiles across the membrane, as well as diffusion potential, have been simulated using the electric circuit simulation program, SPICE. The method is quite general and extremely efficient, permitting treatments of multi-ion systems whatever the boundary and experimental conditions may be.

Accuracy analysis of automodel solutions for Lévy flight-based transport: from resonance radiative transfer to a simple general model

Science.gov (United States)

Kukushkin, A. B.; Sdvizhenskii, P. A.

2017-12-01

The results of accuracy analysis of automodel solutions for Lévy flight-based transport on a uniform background are presented. These approximate solutions have been obtained for Green’s function of the following equations: the non-stationary Biberman-Holstein equation for three-dimensional (3D) radiative transfer in plasma and gases, for various (Doppler, Lorentz, Voigt and Holtsmark) spectral line shapes, and the 1D transport equation with a simple longtailed step-length probability distribution function with various power-law exponents. The results suggest the possibility of substantial extension of the developed method of automodel solution to other fields far beyond physics.

New Poisson–Boltzmann type equations: one-dimensional solutions

International Nuclear Information System (INIS)

Lee, Chiun-Chang; Lee, Hijin; Hyon, YunKyong; Lin, Tai-Chia; Liu, Chun

2011-01-01

The Poisson–Boltzmann (PB) equation is conventionally used to model the equilibrium of bulk ionic species in different media and solvents. In this paper we study a new Poisson–Boltzmann type (PB n ) equation with a small dielectric parameter ε 2 and non-local nonlinearity which takes into consideration the preservation of the total amount of each individual ion. This equation can be derived from the original Poisson–Nernst–Planck system. Under Robin-type boundary conditions with various coefficient scales, we demonstrate the asymptotic behaviours of one-dimensional solutions of PB n equations as the parameter ε approaches zero. In particular, we show that in case of electroneutrality, i.e. α = β, solutions of 1D PB n equations have a similar asymptotic behaviour as those of 1D PB equations. However, as α ≠ β (non-electroneutrality), solutions of 1D PB n equations may have blow-up behaviour which cannot be found in 1D PB equations. Such a difference between 1D PB and PB n equations can also be verified by numerical simulations

Exact Solutions of the Harry-Dym Equation

International Nuclear Information System (INIS)

Mokhtari, Reza

2011-01-01

The aim of this paper is to generate exact travelling wave solutions of the Harry-Dym equation through the methods of Adomian decomposition, He's variational iteration, direct integration, and power series. We show that the two later methods are more successful than the two former to obtain more solutions of the equation. (general)

Cellular neural network to the spherical harmonics approximation of neutron transport equation in x–y geometry

International Nuclear Information System (INIS)

Pirouzmand, Ahmad; Hadad, Kamal

2012-01-01

Highlights: ► This paper describes the solution of time-dependent neutron transport equation. ► We use a novel method based on cellular neural networks (CNNs) coupled with the spherical harmonics method. ► We apply the CNN model to simulate step and ramp perturbation transients in a core. ► The accuracy and capabilities of the CNN model are examined for x–y geometry. - Abstract: In an earlier paper we utilized a novel method using cellular neural networks (CNNs) coupled with spherical harmonics method to solve the steady state neutron transport equation in x–y geometry. Here, the previous work is extended to the study of time-dependent neutron transport equation. To achieve this goal, an equivalent electrical circuit based on a second-order form of time-dependent neutron transport equation and one equivalent group of neutron precursor density is obtained by the CNN method. The CNN model is used to simulate step and ramp perturbation transients in a typical 2D core.

Analytical solutions of time-fractional models for homogeneous Gardner equation and non-homogeneous differential equations

Directory of Open Access Journals (Sweden)

Olaniyi Samuel Iyiola

2014-09-01

Full Text Available In this paper, we obtain analytical solutions of homogeneous time-fractional Gardner equation and non-homogeneous time-fractional models (including Buck-master equation using q-Homotopy Analysis Method (q-HAM. Our work displays the elegant nature of the application of q-HAM not only to solve homogeneous non-linear fractional differential equations but also to solve the non-homogeneous fractional differential equations. The presence of the auxiliary parameter h helps in an effective way to obtain better approximation comparable to exact solutions. The fraction-factor in this method gives it an edge over other existing analytical methods for non-linear differential equations. Comparisons are made upon the existence of exact solutions to these models. The analysis shows that our analytical solutions converge very rapidly to the exact solutions.

Rational approximations to solutions of linear differential equations.

Science.gov (United States)

Chudnovsky, D V; Chudnovsky, G V

1983-08-01

Rational approximations of Padé and Padé type to solutions of differential equations are considered. One of the main results is a theorem stating that a simultaneous approximation to arbitrary solutions of linear differential equations over C(x) cannot be "better" than trivial ones implied by the Dirichlet box principle. This constitutes, in particular, the solution in the linear case of Kolchin's problem that the "Roth's theorem" holds for arbitrary solutions of algebraic differential equations. Complete effective proofs for several valuations are presented based on the Wronskian methods and graded subrings of Picard-Vessiot extensions.

Numerical solution of the time dependent neutron transport equation by the method of the characteristics

International Nuclear Information System (INIS)

Talamo, Alberto

2013-01-01

This study presents three numerical algorithms to solve the time dependent neutron transport equation by the method of the characteristics. The algorithms have been developed taking into account delayed neutrons and they have been implemented into the novel MCART code, which solves the neutron transport equation for two-dimensional geometry and an arbitrary number of energy groups. The MCART code uses regular mesh for the representation of the spatial domain, it models up-scattering, and takes advantage of OPENMP and OPENGL algorithms for parallel computing and plotting, respectively. The code has been benchmarked with the multiplication factor results of a Boiling Water Reactor, with the analytical results for a prompt jump transient in an infinite medium, and with PARTISN and TDTORT results for cross section and source transients. The numerical simulations have shown that only two numerical algorithms are stable for small time steps

Numerical solution of the time dependent neutron transport equation by the method of the characteristics

Energy Technology Data Exchange (ETDEWEB)

Talamo, Alberto, E-mail: alby@anl.gov [Nuclear Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439 (United States)

2013-05-01

This study presents three numerical algorithms to solve the time dependent neutron transport equation by the method of the characteristics. The algorithms have been developed taking into account delayed neutrons and they have been implemented into the novel MCART code, which solves the neutron transport equation for two-dimensional geometry and an arbitrary number of energy groups. The MCART code uses regular mesh for the representation of the spatial domain, it models up-scattering, and takes advantage of OPENMP and OPENGL algorithms for parallel computing and plotting, respectively. The code has been benchmarked with the multiplication factor results of a Boiling Water Reactor, with the analytical results for a prompt jump transient in an infinite medium, and with PARTISN and TDTORT results for cross section and source transients. The numerical simulations have shown that only two numerical algorithms are stable for small time steps.

Classification of exact solutions to the generalized Kadomtsev-Petviashvili equation

International Nuclear Information System (INIS)

Pandir, Yusuf; Gurefe, Yusuf; Misirli, Emine

2013-01-01

In this paper, we study the Kadomtsev-Petviashvili equation with generalized evolution and derive some new results using the approach called the trial equation method. The obtained results can be expressed by the soliton solutions, rational function solutions, elliptic function solutions and Jacobi elliptic function solutions. In the discussion, we give a new version of the trial equation method for nonlinear differential equations.

On the history of a stochastic ansatz for solving the transport equation

International Nuclear Information System (INIS)

Williams, M.M.R.

2010-01-01

A very useful approximate tool for understanding the role of random material properties on solutions of the transport equation is described and its historical derivation given. The development of this stochastic tool, from its introduction by Randall, to its use in describing current problems involving dichotomic or pseudo-dichotomic Markov processes is discussed.

Swarm analysis by using transport equations

International Nuclear Information System (INIS)

Dote, Toshihiko.

1985-01-01

As the basis of weak ionization plasma phenomena, the motion, i.e. swarm, of charged particles in the gas is analyzed by use of the transport equations, from which basic nature of the swarm is discussed. The present report is an overview of the studies made in the past several years. Described are principally the most basic aspects concerning behaviors of the electrons and positive ions, that is, the basic equations and their significance, characteristics of the behaviors of the electron and positive ion swarms as revealed by solving the equations, and various characteristics of the swarm parameters. Contents are: Maxwell-Boltzmann's transport equations, behavior of the electron swarm, energy loss of the electrons, and behavior of the positive ion swarm. (Mori, K.)

Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS

Science.gov (United States)

Morway, Eric D.; Niswonger, Richard G.; Langevin, Christian D.; Bailey, Ryan T.; Healy, Richard W.

2013-01-01

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship.

Crossing-symmetric solutions to low equations

International Nuclear Information System (INIS)

McLeod, R.J.; Ernst, D.J.

1985-01-01

Crossing symmetric models of the pion-nucleon interaction in which crossing symmetry is kept to lowest order in msub(π)/msub(N) are investigated. Two iterative techniques are developed to solve the crossing-symmetric Low equation. The techniques are used to solve the original Chew-Low equations and their generalizations to include the coupling to the pion-production channels. Small changes are found in comparison with earlier results which used an iterative technique proposed by Chew and Low and which did not produce crossing-symmetric results. The iterative technique of Chew and Low is shown to fail because of its inability to produce zeroes in the amplitude at complex energies while physical solutions to the model require such zeroes. We also prove that, within the class of solutions such that phase shifts approach zero for infinite energy, the solution to the Low equation is unique. (orig.)

A time dependent mixing model to close PDF equations for transport in heterogeneous aquifers

Science.gov (United States)

Schüler, L.; Suciu, N.; Knabner, P.; Attinger, S.

2016-10-01

Probability density function (PDF) methods are a promising alternative to predicting the transport of solutes in groundwater under uncertainty. They make it possible to derive the evolution equations of the mean concentration and the concentration variance, used in moment methods. The mixing model, describing the transport of the PDF in concentration space, is essential for both methods. Finding a satisfactory mixing model is still an open question and due to the rather elaborate PDF methods, a difficult undertaking. Both the PDF equation and the concentration variance equation depend on the same mixing model. This connection is used to find and test an improved mixing model for the much easier to handle concentration variance. Subsequently, this mixing model is transferred to the PDF equation and tested. The newly proposed mixing model yields significantly improved results for both variance modelling and PDF modelling.

Singularly perturbed Burger-Huxley equation: Analytical solution ...

African Journals Online (AJOL)

user

numbers, Navier-Stokes flows with large Reynolds numbers, chemical reactor ... It is to observe the layer behavior of the solution for smaller values of ε leading to singular ...... Burger equation, momentum gas equation and heat equation.

Exact solutions to sine-Gordon-type equations

International Nuclear Information System (INIS)

Liu Shikuo; Fu Zuntao; Liu Shida

2006-01-01

In this Letter, sine-Gordon-type equations, including single sine-Gordon equation, double sine-Gordon equation and triple sine-Gordon equation, are systematically solved by Jacobi elliptic function expansion method. It is shown that different transformations for these three sine-Gordon-type equations play different roles in obtaining exact solutions, some transformations may not work for a specific sine-Gordon equation, while work for other sine-Gordon equations

Soliton-like solutions to the ordinary Schroedinger equation

International Nuclear Information System (INIS)

Zamboni-Rached, Michel; Recami, Erasmo

2011-01-01

In recent times it has been paid attention to the fact that (linear) wave equations admit of soliton-like solutions, known as Localized Waves or Non-diffracting Waves, which propagate without distortion in one direction. Such Localized Solutions (existing also for K-G or Dirac equations) are a priori suitable, more than Gaussian's, for describing elementary particle motion. In this paper we show that, mutatis mutandis, Localized Solutions exist even for the ordinary Schroedinger equation within standard Quantum Mechanics; and we obtain both approximate and exact solutions, also setting forth for them particular examples. In the ideal case such solutions bear infinite energy, as well as plane or spherical waves: we show therefore how to obtain nite-energy solutions. At last, we briefly consider solutions for a particle moving in the presence of a potential. (author)

Soliton-like solutions to the ordinary Schroedinger equation

Energy Technology Data Exchange (ETDEWEB)

Zamboni-Rached, Michel [Universidade Estadual de Campinas (DMO/FEEC/UNICAMP), Campinas, SP (Brazil). Fac. de Engenharia Eletrica e de Computacao. Dept. de Microondas e Optica; Recami, Erasmo, E-mail: recami@mi.infn.i [Universita Statale di Bergamo, Bergamo (Italy). Facolta di Ingegneria

2011-07-01

In recent times it has been paid attention to the fact that (linear) wave equations admit of soliton-like solutions, known as Localized Waves or Non-diffracting Waves, which propagate without distortion in one direction. Such Localized Solutions (existing also for K-G or Dirac equations) are a priori suitable, more than Gaussian's, for describing elementary particle motion. In this paper we show that, mutatis mutandis, Localized Solutions exist even for the ordinary Schroedinger equation within standard Quantum Mechanics; and we obtain both approximate and exact solutions, also setting forth for them particular examples. In the ideal case such solutions bear infinite energy, as well as plane or spherical waves: we show therefore how to obtain nite-energy solutions. At last, we briefly consider solutions for a particle moving in the presence of a potential. (author)

Exact analytical solutions for nonlinear reaction-diffusion equations

International Nuclear Information System (INIS)

Liu Chunping

2003-01-01

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

Range of validity of transport equations

International Nuclear Information System (INIS)

Berges, Juergen; Borsanyi, Szabolcs

2006-01-01

Transport equations can be derived from quantum field theory assuming a loss of information about the details of the initial state and a gradient expansion. While the latter can be systematically improved, the assumption about a memory loss is not known to be controlled by a small expansion parameter. We determine the range of validity of transport equations for the example of a scalar g 2 Φ 4 theory. We solve the nonequilibrium time evolution using the three-loop 2PI effective action. The approximation includes off-shell and memory effects and assumes no gradient expansion. This is compared to transport equations to lowest order (LO) and beyond (NLO). We find that the earliest time for the validity of transport equations is set by the characteristic relaxation time scale t damp =-2ω/Σ ρ (eq) , where -Σ ρ (eq) /2 denotes the on-shell imaginary-part of the self-energy. This time scale agrees with the characteristic time for partial memory loss, but is much shorter than thermal equilibration times. For times larger than about t damp the gradient expansion to NLO is found to describe the full results rather well for g 2 (less-or-similar sign)1

Molecular representation of molar domain (volume), evolution equations, and linear constitutive relations for volume transport.

Science.gov (United States)

Eu, Byung Chan

2008-09-07

In the traditional theories of irreversible thermodynamics and fluid mechanics, the specific volume and molar volume have been interchangeably used for pure fluids, but in this work we show that they should be distinguished from each other and given distinctive statistical mechanical representations. In this paper, we present a general formula for the statistical mechanical representation of molecular domain (volume or space) by using the Voronoi volume and its mean value that may be regarded as molar domain (volume) and also the statistical mechanical representation of volume flux. By using their statistical mechanical formulas, the evolution equations of volume transport are derived from the generalized Boltzmann equation of fluids. Approximate solutions of the evolution equations of volume transport provides kinetic theory formulas for the molecular domain, the constitutive equations for molar domain (volume) and volume flux, and the dissipation of energy associated with volume transport. Together with the constitutive equation for the mean velocity of the fluid obtained in a previous paper, the evolution equations for volume transport not only shed a fresh light on, and insight into, irreversible phenomena in fluids but also can be applied to study fluid flow problems in a manner hitherto unavailable in fluid dynamics and irreversible thermodynamics. Their roles in the generalized hydrodynamics will be considered in the sequel.

Solution of some types of differential equations: operational calculus and inverse differential operators.

Science.gov (United States)

Zhukovsky, K

2014-01-01

We present a general method of operational nature to analyze and obtain solutions for a variety of equations of mathematical physics and related mathematical problems. We construct inverse differential operators and produce operational identities, involving inverse derivatives and families of generalised orthogonal polynomials, such as Hermite and Laguerre polynomial families. We develop the methodology of inverse and exponential operators, employing them for the study of partial differential equations. Advantages of the operational technique, combined with the use of integral transforms, generating functions with exponentials and their integrals, for solving a wide class of partial derivative equations, related to heat, wave, and transport problems, are demonstrated.

Special function solutions of the free particle Dirac equation

International Nuclear Information System (INIS)

Strange, P

2012-01-01

The Dirac equation is one of the fundamental equations in physics. Here we present and discuss two novel solutions of the free particle Dirac equation. These solutions have an exact analytical form in terms of Airy or Mathieu functions and exhibit unexpected properties including an enhanced Doppler effect, accelerating wavefronts and solutions with a degree of localization. (paper)

New types of exact solutions for a breaking soliton equation

International Nuclear Information System (INIS)

Mei Jianqin; Zhang Hongqing

2004-01-01

In this paper based on a system of Riccati equations, we present a newly generally projective Riccati equation expansion method and its algorithm, which can be used to construct more new exact solutions of nonlinear differential equations in mathematical physics. A typical breaking soliton equation is chosen to illustrate our algorithm such that more families of new exact solutions are obtained, which contain soliton-like solutions and periodic solutions. This algorithm can also be applied to other nonlinear differential equations

Manufactured solutions for the three-dimensional Euler equations with relevance to Inertial Confinement Fusion

International Nuclear Information System (INIS)

Waltz, J.; Canfield, T.R.; Morgan, N.R.; Risinger, L.D.; Wohlbier, J.G.

2014-01-01

We present a set of manufactured solutions for the three-dimensional (3D) Euler equations. The purpose of these solutions is to allow for code verification against true 3D flows with physical relevance, as opposed to 3D simulations of lower-dimensional problems or manufactured solutions that lack physical relevance. Of particular interest are solutions with relevance to Inertial Confinement Fusion (ICF) capsules. While ICF capsules are designed for spherical symmetry, they are hypothesized to become highly 3D at late time due to phenomena such as Rayleigh–Taylor instability, drive asymmetry, and vortex decay. ICF capsules also involve highly nonlinear coupling between the fluid dynamics and other physics, such as radiation transport and thermonuclear fusion. The manufactured solutions we present are specifically designed to test the terms and couplings in the Euler equations that are relevant to these phenomena. Example numerical results generated with a 3D Finite Element hydrodynamics code are presented, including mesh convergence studies

Analytic solution of integral equations for molecular fluids

International Nuclear Information System (INIS)

Cummings, P.T.

1984-01-01

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

An extended step characteristic method for solving the transport equation in general geometries

International Nuclear Information System (INIS)

DeHart, M.D.; Pevey, R.E.; Parish, T.A.

1994-01-01

A method for applying the discrete ordinates method to solve the Boltzmann transport equation on arbitrary two-dimensional meshes has been developed. The finite difference approach normally used to approximate spatial derivatives in extrapolating angular fluxes across a cell is replaced by direct solution of the characteristic form of the transport equation for each discrete direction. Thus, computational cells are not restricted to the geometrical shape of a mesh element characteristic of a given coordinate system. However, in terms of the treatment of energy and angular dependencies, this method resembles traditional discrete ordinates techniques. By using the method developed here, a general two-dimensional space can be approximated by an irregular mesh comprised of arbitrary polygons. Results for a number of test problems have been compared with solutions obtained from traditional methods, with good agreement. Comparisons include benchmarks against analytical results for problems with simple geometry, as well as numerical results obtained from traditional discrete ordinates methods by applying the ANISN and TWOTRAN-II computer programs

Pure soliton solutions of some nonlinear partial differential equations

International Nuclear Information System (INIS)

Fuchssteiner, B.

1977-01-01

A general approach is given to obtain the system of ordinary differential equations which determines the pure soliton solutions for the class of generalized Korteweg-de Vries equations. This approach also leads to a system of ordinary differential equations for the pure soliton solutions of the sine-Gordon equation. (orig.) [de

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)

A linear multiple balance method for discrete ordinates neutron transport equations

International Nuclear Information System (INIS)

Park, Chang Je; Cho, Nam Zin

2000-01-01

A linear multiple balance method (LMB) is developed to provide more accurate and positive solutions for the discrete ordinates neutron transport equations. In this multiple balance approach, one mesh cell is divided into two subcells with quadratic approximation of angular flux distribution. Four multiple balance equations are used to relate center angular flux with average angular flux by Simpson's rule. From the analysis of spatial truncation error, the accuracy of the linear multiple balance scheme is ο(Δ 4 ) whereas that of diamond differencing is ο(Δ 2 ). To accelerate the linear multiple balance method, we also describe a simplified additive angular dependent rebalance factor scheme which combines a modified boundary projection acceleration scheme and the angular dependent rebalance factor acceleration schme. It is demonstrated, via fourier analysis of a simple model problem as well as numerical calculations, that the additive angular dependent rebalance factor acceleration scheme is unconditionally stable with spectral radius < 0.2069c (c being the scattering ration). The numerical results tested so far on slab-geometry discrete ordinates transport problems show that the solution method of linear multiple balance is effective and sufficiently efficient

Thin-Layer Solutions of the Helmholtz and Related Equations

KAUST Repository

Ockendon, J. R.

2012-01-01

This paper concerns a certain class of two-dimensional solutions to four generic partial differential equations-the Helmholtz, modified Helmholtz, and convection-diffusion equations, and the heat conduction equation in the frequency domain-and the connections between these equations for this particular class of solutions.S pecifically, we consider thin-layer solutions, valid in narrow regions across which there is rapid variation, in the singularly perturbed limit as the coefficient of the Laplacian tends to zero.F or the wellstudied Helmholtz equation, this is the high-frequency limit and the solutions in question underpin the conventional ray theory/WKB approach in that they provide descriptions valid in some of the regions where these classical techniques fail.E xamples are caustics, shadow boundaries, whispering gallery, and creeping waves and focusing and bouncing ball modes.It transpires that virtually all such thin-layer models reduce to a class of generalized parabolic wave equations, of which the heat conduction equation is a special case. Moreover, in most situations, we will find that the appropriate parabolic wave equation solutions can be derived as limits of exact solutions of the Helmholtz equation.W e also show how reasonably well-understood thin-layer phenomena associated with any one of the four generic equations may translate into less well-known effects associated with the others.In addition, our considerations also shed some light on the relationship between the methods of matched asymptotic, WKB, and multiple-scales expansions. © 2012 Society for Industrial and Applied Mathematics.

Effects of coal gangue content on water movement and solute transport in a China loess plateau soil

Energy Technology Data Exchange (ETDEWEB)

Beibei, Zhou; Quanjiu, Wang [Institute of Water Resources and Hydro-electric Engineering, Xi' an University of Technology, Xi' an (China); State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Northwest A and F University, Yangling, Shaanxi (China); Ming' an, Shao; Mingxia, Wen [State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Northwest A and F University, Yangling, Shaanxi (China); College of Resources and Environment, Northwest A and F University, Yangling, Shaanxi (China); Horton, Robert [Department of Agronomy, Iowa State University, Ames, Iowa (United States)

2010-11-15

The mining industry has grown strongly in China in recent decades, resulting in large amounts of coal gangues, which cause water and soil pollution, soil erosion, and various other environmental problems. They are often used in reclamation projects in attempts to restore land damaged by mining, hence they are frequently present (in widely varying proportions) in the topsoil in areas around mines. Their presence can strongly affect key soil variables, including its bulk density, structure, water retention, water movement, and solute transport rates. In the study presented here, the effects of gangue contents on infiltration, saturated hydraulic conductivity, and solute transport parameters of a Chinese Loess plateau soil were examined. The results show that infiltration rates and saturated hydraulic conductivity decreased with increasing gangue content. The Peck-Watson equation modeled these relationships well, but Bouwer-Rice equations provided poorer matches with the acquired data. Cumulative infiltration over time was described well by both the Philip equation and Kostiakov equation. Both the simplified convection-dispersion equation and a two-region model described the solute transport processes well. In addition, the dispersion increased, while both the Peclet number and mobile water fraction decreased, with increases in gangue contents. (Copyright copyright 2010 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim)

Quantifying the relative contributions of different solute carriers to aggregate substrate transport

Science.gov (United States)

Taslimifar, Mehdi; Oparija, Lalita; Verrey, Francois; Kurtcuoglu, Vartan; Olgac, Ufuk; Makrides, Victoria

2017-01-01

Determining the contributions of different transporter species to overall cellular transport is fundamental for understanding the physiological regulation of solutes. We calculated the relative activities of Solute Carrier (SLC) transporters using the Michaelis-Menten equation and global fitting to estimate the normalized maximum transport rate for each transporter (Vmax). Data input were the normalized measured uptake of the essential neutral amino acid (AA) L-leucine (Leu) from concentration-dependence assays performed using Xenopus laevis oocytes. Our methodology was verified by calculating Leu and L-phenylalanine (Phe) data in the presence of competitive substrates and/or inhibitors. Among 9 potentially expressed endogenous X. laevis oocyte Leu transporter species, activities of only the uniporters SLC43A2/LAT4 (and/or SLC43A1/LAT3) and the sodium symporter SLC6A19/B0AT1 were required to account for total uptake. Furthermore, Leu and Phe uptake by heterologously expressed human SLC6A14/ATB0,+ and SLC43A2/LAT4 was accurately calculated. This versatile systems biology approach is useful for analyses where the kinetics of each active protein species can be represented by the Hill equation. Furthermore, its applicable even in the absence of protein expression data. It could potentially be applied, for example, to quantify drug transporter activities in target cells to improve specificity. PMID:28091567

Nature of complex time eigenvalues of the one speed transport equation in a homogeneous sphere

International Nuclear Information System (INIS)

Dahl, E.B.; Sahni, D.C.

1990-01-01

The complex time eigenvalues of the transport equation have been studied for one speed neutrons, scattered isotropically in a homogeneous sphere with vacuum boundary conditions. It is shown that the complex decay constants vary continuously with the radius of the sphere. Our earlier conjecture (Dahl and Sahni (1983-84)) regarding disjoint arcs is thus shown to be true. We also indicate that complex decay constants exist even for large assemblies, though with rapid oscillations in the corresponding eigenvectors. These modes cannot be predicted by the diffusion equation as this behaviour of the eigenvectors contradicts the assumption of 'slowly varying flux' needed to derive the diffusion approximation from the transport equation. For an infinite system, the existence of complex modes is related to the solution of a homogeneous equation. (author)

Exact solutions for the cubic-quintic nonlinear Schroedinger equation

International Nuclear Information System (INIS)

Zhu Jiamin; Ma Zhengyi

2007-01-01

In this paper, the cubic-quintic nonlinear Schroedinger equation is solved through the extended elliptic sub-equation method. As a consequence, many types of exact travelling wave solutions are obtained which including bell and kink profile solitary wave solutions, triangular periodic wave solutions and singular solutions

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

International Nuclear Information System (INIS)

Xiao Yafeng; Xue Haili; Zhang Hongqing

2011-01-01

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

Approximate solutions to Mathieu's equation

Science.gov (United States)

Wilkinson, Samuel A.; Vogt, Nicolas; Golubev, Dmitry S.; Cole, Jared H.

2018-06-01

Mathieu's equation has many applications throughout theoretical physics. It is especially important to the theory of Josephson junctions, where it is equivalent to Schrödinger's equation. Mathieu's equation can be easily solved numerically, however there exists no closed-form analytic solution. Here we collect various approximations which appear throughout the physics and mathematics literature and examine their accuracy and regimes of applicability. Particular attention is paid to quantities relevant to the physics of Josephson junctions, but the arguments and notation are kept general so as to be of use to the broader physics community.

Approximate variational solutions of the Grad-Shafranov equation

International Nuclear Information System (INIS)

Ludwig, G.O.

2001-01-01

Approximate solutions of the Grad-Schlueter-Shafranov equation based on variational methods are developed. The power series solutions of the Euler-Lagrange equations for equilibrium are compared with direct variational results for a low aspect ratio tokamak equilibrium. (author)

The quasi-diffusive approximation in transport theory: Local solutions

International Nuclear Information System (INIS)

Celaschi, M.; Montagnini, B.

1995-01-01

The one velocity, plane geometry integral neutron transport equation is transformed into a system of two equations, one of them being the equation of continuity and the other a generalized Fick's law, in which the usual diffusion coefficient is replaced by a self-adjoint integral operator. As the kernel of this operator is very close to the Green function of a diffusion equation, an approximate inversion by means of a second order differential operator allows to transform these equations into a purely differential system which is shown to be equivalent, in the simplest case, to a diffusion-like equation. The method, the principles of which have been exposed in a previous paper, is here extended and applied to a variety of problems. If the inversion is properly performed, the quasi-diffusive solutions turn out to be quite accurate, even in the vicinity of the interface between different material regions, where elementary diffusion theory usually fails. 16 refs., 3 tabs

Exact traveling wave solutions of modified KdV-Zakharov-Kuznetsov equation and viscous Burgers equation.

Science.gov (United States)

Islam, Md Hamidul; Khan, Kamruzzaman; Akbar, M Ali; Salam, Md Abdus

2014-01-01

Mathematical modeling of many physical systems leads to nonlinear evolution equations because most physical systems are inherently nonlinear in nature. The investigation of traveling wave solutions of nonlinear partial differential equations (NPDEs) plays a significant role in the study of nonlinear physical phenomena. In this article, we construct the traveling wave solutions of modified KDV-ZK equation and viscous Burgers equation by using an enhanced (G '/G) -expansion method. A number of traveling wave solutions in terms of unknown parameters are obtained. Derived traveling wave solutions exhibit solitary waves when special values are given to its unknown parameters. 35C07; 35C08; 35P99.

A fractional Dirac equation and its solution

International Nuclear Information System (INIS)

Muslih, Sami I; Agrawal, Om P; Baleanu, Dumitru

2010-01-01

This paper presents a fractional Dirac equation and its solution. The fractional Dirac equation may be obtained using a fractional variational principle and a fractional Klein-Gordon equation; both methods are considered here. We extend the variational formulations for fractional discrete systems to fractional field systems defined in terms of Caputo derivatives. By applying the variational principle to a fractional action S, we obtain the fractional Euler-Lagrange equations of motion. We present a Lagrangian and a Hamiltonian for the fractional Dirac equation of order α. We also use a fractional Klein-Gordon equation to obtain the fractional Dirac equation which is the same as that obtained using the fractional variational principle. Eigensolutions of this equation are presented which follow the same approach as that for the solution of the standard Dirac equation. We also provide expressions for the path integral quantization for the fractional Dirac field which, in the limit α → 1, approaches to the path integral for the regular Dirac field. It is hoped that the fractional Dirac equation and the path integral quantization of the fractional field will allow further development of fractional relativistic quantum mechanics.

Bounds and maximum principles for the solution of the linear transport equation

International Nuclear Information System (INIS)

Larsen, E.W.

1981-01-01

Pointwise bounds are derived for the solution of time-independent linear transport problems with surface sources in convex spatial domains. Under specified conditions, upper bounds are derived which, as a function of position, decrease with distance from the boundary. Also, sufficient conditions are obtained for the existence of maximum and minimum principles, and a counterexample is given which shows that such principles do not always exist

Numerical solutions of the monoenergetic neutron transport equation with anisotropic scattering

International Nuclear Information System (INIS)

Dahl, B.

1985-01-01

The Boltzmann equation for monoenergetic neutrons has been solved numerically with high accuracy for homogeneous slabs and spheres with various degree of linear anisotropy. Vacuum boundary conditions are used. The numerical method is based on previous work by Carlvik. Benchmark values of the criticality factor and higher order eigenvalues are given for multiplying systems of thickness or diameter from 10 -5 to 20 mean free paths and with anisotropy coefficients from 0.0 to 0.3. For slab geometry, both even and odd mode eigenvalues are treated. With increasing anisotropy, an increasing number of complex eigenvalues is observer. The total flux is calculated from the eigenvector and tables of the fundamental mode flux are given. Accurate extrapolation distances are derived for various dimensions and anisotropy coefficients from our eigenvalue results on slabs and spheres and from the work by Sanchez on infinite cylinders.The time eigenvalue spectrum in subcritical systems has also been studied. First, the connection between the eigenvalues arising from the time dependent and stationary transport equation is established. Based on this, the spectrum of real time eigenvalues in slabs and spheres is calculated. For spheres, the existence of complex time eigenvalues in the region beyond the value corresponding to the Corngold limit is numerically established. The presence of such eigenvalues has earlier not been proved. It is further shown that the Boltzmann equation for a sphere is significantly simplified when the decay constant is at the Corngold limit. The spectrum of sphere diameters corresponding to this decay constant is calculated for various linear anisotropies, and detailed numerical results are given. (Author)

New exact travelling wave solutions for the Ostrovsky equation

International Nuclear Information System (INIS)

Kangalgil, Figen; Ayaz, Fatma

2008-01-01

In this Letter, auxiliary equation method is proposed for constructing more general exact solutions of nonlinear partial differential equation with the aid of symbolic computation. In order to illustrate the validity and the advantages of the method we choose the Ostrovsky equation. As a result, many new and more general exact solutions have been obtained for the equation

Decay Mode Solutions for Kadomtsev-Petviashvili Equation

International Nuclear Information System (INIS)

Fan Guohao; Deng Shufang; Zhang Meng

2012-01-01

The decay mode solutions for the Kadomtsev-Petviashvili (KP) equation are derived by Hirota method (direct method). The decay mode solution is a new set of analytical solutions with Airy function. (general)

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

International Nuclear Information System (INIS)

Zhang Huiqun

2009-01-01

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

Elliptic equation rational expansion method and new exact travelling solutions for Whitham-Broer-Kaup equations

International Nuclear Information System (INIS)

Chen Yong; Wang Qi; Li Biao

2005-01-01

Based on a new general ansatz and a general subepuation, a new general algebraic method named elliptic equation rational expansion method is devised for constructing multiple travelling wave solutions in terms of rational special function for nonlinear evolution equations (NEEs). We apply the proposed method to solve Whitham-Broer-Kaup equation and explicitly construct a series of exact solutions which include rational form solitary wave solution, rational form triangular periodic wave solutions and rational wave solutions as special cases. In addition, the links among our proposed method with the method by Fan [Chaos, Solitons and Fractals 2004;20:609], are also clarified generally

Shallow water equations: viscous solutions and inviscid limit

Science.gov (United States)

Chen, Gui-Qiang; Perepelitsa, Mikhail

2012-12-01

We establish the inviscid limit of the viscous shallow water equations to the Saint-Venant system. For the viscous equations, the viscosity terms are more degenerate when the shallow water is close to the bottom, in comparison with the classical Navier-Stokes equations for barotropic gases; thus, the analysis in our earlier work for the classical Navier-Stokes equations does not apply directly, which require new estimates to deal with the additional degeneracy. We first introduce a notion of entropy solutions to the viscous shallow water equations and develop an approach to establish the global existence of such solutions and their uniform energy-type estimates with respect to the viscosity coefficient. These uniform estimates yield the existence of measure-valued solutions to the Saint-Venant system generated by the viscous solutions. Based on the uniform energy-type estimates and the features of the Saint-Venant system, we further establish that the entropy dissipation measures of the viscous solutions for weak entropy-entropy flux pairs, generated by compactly supported C 2 test-functions, are confined in a compact set in H -1, which yields that the measure-valued solutions are confined by the Tartar-Murat commutator relation. Then, the reduction theorem established in Chen and Perepelitsa [5] for the measure-valued solutions with unbounded support leads to the convergence of the viscous solutions to a finite-energy entropy solution of the Saint-Venant system with finite-energy initial data, which is relative with respect to the different end-states of the bottom topography of the shallow water at infinity. The analysis also applies to the inviscid limit problem for the Saint-Venant system in the presence of friction.

An Algebraic Method for Constructing Exact Solutions to Difference-Differential Equations

International Nuclear Information System (INIS)

Wang Zhen; Zhang Hongqing

2006-01-01

In this paper, we present a method to solve difference differential equation(s). As an example, we apply this method to discrete KdV equation and Ablowitz-Ladik lattice equation. As a result, many exact solutions are obtained with the help of Maple including soliton solutions presented by hyperbolic functions sinh and cosh, periodic solutions presented by sin and cos and rational solutions. This method can also be used to other nonlinear difference-differential equation(s).

Analytical exact solution of the non-linear Schroedinger equation

International Nuclear Information System (INIS)

Martins, Alisson Xavier; Rocha Filho, Tarcisio Marciano da

2011-01-01

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

Analytical solution for multi-species contaminant transport in finite media with time-varying boundary conditions

Science.gov (United States)

Most analytical solutions available for the equations governing the advective-dispersive transport of multiple solutes undergoing sequential first-order decay reactions have been developed for infinite or semi-infinite spatial domains and steady-state boundary conditions. In this work we present an ...

Numerical solutions of diffusive logistic equation

International Nuclear Information System (INIS)

Afrouzi, G.A.; Khademloo, S.

2007-01-01

In this paper we investigate numerically positive solutions of a superlinear Elliptic equation on bounded domains. The study of Diffusive logistic equation continues to be an active field of research. The subject has important applications to population migration as well as many other branches of science and engineering. In this paper the 'finite difference scheme' will be developed and compared for solving the one- and three-dimensional Diffusive logistic equation. The basis of the analysis of the finite difference equations considered here is the modified equivalent partial differential equation approach, developed from many authors these years

Preservation of support and positivity for solutions of degenerate evolution equations

International Nuclear Information System (INIS)

Ambrose, David M; Wright, J Douglas

2010-01-01

We prove that sufficiently smooth solutions of equations of a certain class have two interesting properties. These evolution equations are in a sense degenerate, in that every term on the right-hand side of the evolution equation has either the unknown or its first spatial derivative as a factor. We first find a conserved quantity for the equation: the measure of the set on which the solution is non-zero. Second, we show that solutions which are initially non-negative remain non-negative for all times. These properties rely heavily upon the degeneracy of the leading order term. When the equation is more degenerate, we are able to prove that there are additional conserved quantities: the measure of the set on which the solution is positive and the measure of the set on which the solution is negative. To illustrate these results, we give examples of equations with nonlinear dispersion which have solutions in spaces with sufficient regularity to satisfy the hypotheses of the support and positivity theorems. An important family of equations with nonlinear dispersion are the Rosenau–Hyman compacton equations; there is no existence theory yet for these equations, but the known solutions of the compacton equations are of lower regularity than is needed for the preceding theorems. We prove an additional positivity theorem which applies to solutions of the same family of equations in a function space which includes some solutions of compacton equations

Differential and difference equations a comparison of methods of solution

CERN Document Server

Maximon, Leonard C

2016-01-01

This book, intended for researchers and graduate students in physics, applied mathematics and engineering, presents a detailed comparison of the important methods of solution for linear differential and difference equations - variation of constants, reduction of order, Laplace transforms and generating functions - bringing out the similarities as well as the significant differences in the respective analyses. Equations of arbitrary order are studied, followed by a detailed analysis for equations of first and second order. Equations with polynomial coefficients are considered and explicit solutions for equations with linear coefficients are given, showing significant differences in the functional form of solutions of differential equations from those of difference equations. An alternative method of solution involving transformation of both the dependent and independent variables is given for both differential and difference equations. A comprehensive, detailed treatment of Green’s functions and the associat...

A class of exact solutions to the Einstein field equations

International Nuclear Information System (INIS)

Goyal, Nisha; Gupta, R K

2012-01-01

The Einstein-Rosen metric is considered and a class of new exact solutions of the field equations for stationary axisymmetric Einstein-Maxwell fields is obtained. The Lie classical approach is applied to obtain exact solutions. By using the Lie classical method, we are able to derive symmetries that are used for reducing the coupled system of partial differential equations into ordinary differential equations. From reduced differential equations we have derived some new exact solutions of Einstein-Maxwell equations. (paper)

On the Solution of the Rational Matrix Equation

Directory of Open Access Journals (Sweden)

Faßbender Heike

2007-01-01

Full Text Available We study numerical methods for finding the maximal symmetric positive definite solution of the nonlinear matrix equation , where is symmetric positive definite and is nonsingular. Such equations arise for instance in the analysis of stationary Gaussian reciprocal processes over a finite interval. Its unique largest positive definite solution coincides with the unique positive definite solution of a related discrete-time algebraic Riccati equation (DARE. We discuss how to use the butterfly algorithm to solve the DARE. This approach is compared to several fixed-point and doubling-type iterative methods suggested in the literature.

On the solutions of the second heavenly and Pavlov equations

Science.gov (United States)

Manakov, S. V.; Santini, P. M.

2009-10-01

We have recently solved the inverse scattering problem for one-parameter families of vector fields, and used this result to construct the formal solution of the Cauchy problem for a class of integrable nonlinear partial differential equations connected with the commutation of multidimensional vector fields, such as the heavenly equation of Plebanski, the dispersionless Kadomtsev-Petviashvili (dKP) equation and the two-dimensional dispersionless Toda (2ddT) equation, as well as with the commutation of one-dimensional vector fields, such as the Pavlov equation. We also showed that the associated Riemann-Hilbert inverse problems are powerful tools to establish if the solutions of the Cauchy problem break at finite time, to construct their long-time behaviour and characterize classes of implicit solutions. In this paper, using the above theory, we concentrate on the heavenly and Pavlov equations, (i) establishing that their localized solutions evolve without breaking, unlike the cases of dKP and 2ddT; (ii) constructing the long-time behaviour of the solutions of their Cauchy problems; (iii) characterizing a distinguished class of implicit solutions of the heavenly equation.

On the solutions of the second heavenly and Pavlov equations

International Nuclear Information System (INIS)

Manakov, S V; Santini, P M

2009-01-01

We have recently solved the inverse scattering problem for one-parameter families of vector fields, and used this result to construct the formal solution of the Cauchy problem for a class of integrable nonlinear partial differential equations connected with the commutation of multidimensional vector fields, such as the heavenly equation of Plebanski, the dispersionless Kadomtsev-Petviashvili (dKP) equation and the two-dimensional dispersionless Toda (2ddT) equation, as well as with the commutation of one-dimensional vector fields, such as the Pavlov equation. We also showed that the associated Riemann-Hilbert inverse problems are powerful tools to establish if the solutions of the Cauchy problem break at finite time, to construct their long-time behaviour and characterize classes of implicit solutions. In this paper, using the above theory, we concentrate on the heavenly and Pavlov equations, (i) establishing that their localized solutions evolve without breaking, unlike the cases of dKP and 2ddT; (ii) constructing the long-time behaviour of the solutions of their Cauchy problems; (iii) characterizing a distinguished class of implicit solutions of the heavenly equation.

Positive Solutions for Coupled Nonlinear Fractional Differential Equations

Directory of Open Access Journals (Sweden)

Wenning Liu

2014-01-01

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

Bayesian estimation of the hydraulic and solute transport properties of a small-scale unsaturated soil column

NARCIS (Netherlands)

Moreira, Paulo H S; Van Genuchten, Martinus Th; Orlande, Helcio R B; Cotta, Renato M.

2016-01-01

In this study the hydraulic and solute transport properties of an unsaturated soil were estimated simultaneously from a relatively simple small-scale laboratory column infiltration/outflow experiment. As governing equations we used the Richards equation for variably saturated flow and a physical

Growth of meromorphic solutions of delay differential equations

OpenAIRE

Halburd, Rod; Korhonen, Risto

2016-01-01

Necessary conditions are obtained for certain types of rational delay differential equations to admit a non-rational meromorphic solution of hyper-order less than one. The equations obtained include delay Painlev\\'e equations and equations solved by elliptic functions.

Unconditionally stable diffusion-acceleration of the transport equation

International Nuclear Information System (INIS)

Larsen, E.W.

1982-01-01

The standard iterative procedure for solving fixed-source discrete-ordinates problems converges very slowly for problems in optically large regions with scattering ratios c near unity. The diffusion-synthetic acceleration method has been proposed to make use of the fact that for this class of problems the diffusion equation is often an accurate approximation to the transport equation. However, stability difficulties have historically hampered the implementation of this method for general transport differencing schemes. In this article we discuss a recently developed procedure for obtaining unconditionally stable diffusion-synthetic acceleration methods for various transport differencing schemes. We motivate the analysis by first discussing the exact transport equation; then we illustrate the procedure by deriving a new stable acceleration method for the linear discontinuous transport differencing scheme. We also provide some numerical results

Unconditionally stable diffusion-acceleration of the transport equation

International Nuclear Information System (INIS)

Larson, E.W.

1982-01-01

The standard iterative procedure for solving fixed-source discrete-ordinates problems converges very slowly for problems in optically thick regions with scattering ratios c near unity. The diffusion-synthetic acceleration method has been proposed to make use of the fact that for this class of problems, the diffusion equation is often an accurate approximation to the transport equation. However, stability difficulties have historically hampered the implementation of this method for general transport differencing schemes. In this article we discuss a recently developed procedure for obtaining unconditionally stable diffusion-synthetic acceleration methods for various transport differencing schemes. We motivate the analysis by first discussing the exact transport equation; then we illustrate the procedure by deriving a new stable acceleration method for the linear discontinuous transport differencing scheme. We also provide some numerical results

Solution of the neutron transport equation by means of Hermite-Ssub(infinity)-theory

International Nuclear Information System (INIS)

Brandt, D.; Haelg, W.; Mennig, J.

1979-01-01

A stable numerical approximation Hsub(α)-Ssub(infinity) is obtained through the use of Hermite's method of order α(Hsub(α)) in the spatial integration of the ID neutron transport equation. The theory for α = 1 is applied to a one-group shielding problem. Numerical calculations show the new method to converge much faster than earlier versions of Ssub(infinity)-theory. Comparison of H 1 - Ssub(infinity) with the well-known Ssub(N)-code ANISN indicates a large gain in computing time for the former. (Auth.)

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

International Nuclear Information System (INIS)

Burt, P.B.

1978-01-01

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

Exact solution of nonsteady thermal boundary layer equation

International Nuclear Information System (INIS)

Dorfman, A.S.

1995-01-01

There are only a few exact solutions of the thermal boundary layer equation. Most of them are derived for a specific surface temperature distribution. The first exact solution of the steady-state boundary layer equation was given for a plate with constant surface temperature and free-stream velocity. The same problem for a plate with polynomial surface temperature distribution was solved by Chapmen and Rubesin. Levy gave the exact solution for the case of a power law distribution of both surface temperature and free-stream velocity. The exact solution of the steady-state boundary layer equation for an arbitrary surface temperature and a power law free-stream velocity distribution was given by the author in two forms: of series and of the integral with an influence function of unheated zone. A similar solution of the nonsteady thermal boundary layer equation for an arbitrary surface temperature and a power law free-stream velocity distribution is presented here. In this case, the coefficients of series depend on time, and in the limit t → ∞ they become the constant coefficients of a similar solution published before. This solution, unlike the one presented here, does not satisfy the initial conditions at t = 0, and, hence, can be used only in time after the beginning of the process. The solution in the form of a series becomes a closed-form exact solution for polynomial surface temperature and a power law free-stream velocity distribution. 7 refs., 2 figs

EXACT TRAVELLING WAVE SOLUTIONS TO BBM EQUATION

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

Abundant new travelling wave solutions to the BBM (Benjamin-Bona-Mahoni) equation are obtained by the generalized Jacobian elliptic function method. This method can be applied to other nonlinear evolution equations.

Methodology for solving the equation of transport ordered discrete TORT code in the reactor IPEN/MB-01; Metodologia para resolver la ecuacion del transporte con el codigo de Ordenadas Discretas TORT en el reactor IPEN/MB-01

Energy Technology Data Exchange (ETDEWEB)

Bernal, A.; Abarca, A.; Barrachina, T.; Miro, R.; Verdu, G.

2013-07-01

The resolution of the neutron transport equation in steady state in pool-type nuclear reactors, is normally achieved through 2 different numerical methods: Monte Carlo (stochastic) and discrete ordinates (deterministic). The discrete ordinates method solves the neutron transport equation for a set of specific addresses, obtaining a set of equations and solutions for each direction, where the solution for each direction is the angular flux. With the aim of treating energy dependence, used energy multigroup approximation, thus obtaining a set of equations that depends on the number of energy groups considered.

Similarity Solutions for Multiterm Time-Fractional Diffusion Equation

OpenAIRE

Elsaid, A.; Abdel Latif, M. S.; Maneea, M.

2016-01-01

Similarity method is employed to solve multiterm time-fractional diffusion equation. The orders of the fractional derivatives belong to the interval (0,1] and are defined in the Caputo sense. We illustrate how the problem is reduced from a multiterm two-variable fractional partial differential equation to a multiterm ordinary fractional differential equation. Power series solution is obtained for the resulting ordinary problem and the convergence of the series solution is discussed. Based on ...

Properties of meromorphic solutions to certain differential-difference equations

Directory of Open Access Journals (Sweden)

Xiaoguang Qi

2013-06-01

Full Text Available We consider the properties of meromorphic solutions to certain type of non-linear difference equations. Also we show the existence of meromorphic solutions with finite order for differential-difference equations related to the Fermat type functional equation. This article extends earlier results by Liu et al [12].

Some new exact solutions of Jacobian elliptic function about the generalized Boussinesq equation and Boussinesq-Burgers equation

International Nuclear Information System (INIS)

Zhang Liang; Zhang Lifeng; Li Chongyin

2008-01-01

By using the modified mapping method, we find some new exact solutions of the generalized Boussinesq equation and the Boussinesq-Burgers equation. The solutions obtained in this paper include Jacobian elliptic function solutions, combined Jacobian elliptic function solutions, soliton solutions, triangular function solutions

Exact solutions of some nonlinear partial differential equations using ...

Indian Academy of Sciences (India)

The functional variable method is a powerful solution method for obtaining exact solutions of some nonlinear partial differential equations. In this paper, the functional variable method is used to establish exact solutions of the generalized forms of Klein–Gordon equation, the (2 + 1)-dimensional Camassa–Holm ...

Darboux Transformation and Explicit Solutions for Drinfel'd-Sokolov-Wilson Equation

International Nuclear Information System (INIS)

Geng Xianguo; Wu Lihua

2010-01-01

A generalized Drinfel'd-Sokolov-Wilson (DSW) equation and its Lax pair are proposed. A Darboux transformation for the generalized DSW equation is constructed with the help of the gauge transformation between spectral problems, from which a Darboux transformation for the DSW equation is obtained through a reduction technique. As an application of the Darboux transformations, we give some explicit solutions of the generalized DSW equation and DSW equation such as rational solutions, soliton solutions, periodic solutions. (general)

Homoclinic solutions for Davey-Stewartson equation

International Nuclear Information System (INIS)

Huang Jian; Dai Zhengde

2008-01-01

In this paper, we firstly prove the existence of homoclinic solutions for Davey-Stewartson I equation (DSI) with the periodic boundary condition. Then we obtain a set of exact homoclinic solutions by the novel method-Hirota's method. Moreover, the structure of homoclinic solutions has been investigated. At the same time, we give some numerical simulations which validate these theoretical results

Diffusion equation and spin drag in spin-polarized transport

DEFF Research Database (Denmark)

Flensberg, Karsten; Jensen, Thomas Stibius; Mortensen, Asger

2001-01-01

We study the role of electron-electron interactions for spin-polarized transport using the Boltzmann equation, and derive a set of coupled transport equations. For spin-polarized transport the electron-electron interactions are important, because they tend to equilibrate the momentum of the two-s...

Parallel computing solution of Boltzmann neutron transport equation

International Nuclear Information System (INIS)

Ansah-Narh, T.

2010-01-01

The focus of the research was on developing parallel computing algorithm for solving Eigen-values of the Boltzmam Neutron Transport Equation (BNTE) in a slab geometry using multi-grid approach. In response to the problem of slow execution of serial computing when solving large problems, such as BNTE, the study was focused on the design of parallel computing systems which was an evolution of serial computing that used multiple processing elements simultaneously to solve complex physical and mathematical problems. Finite element method (FEM) was used for the spatial discretization scheme, while angular discretization was accomplished by expanding the angular dependence in terms of Legendre polynomials. The eigenvalues representing the multiplication factors in the BNTE were determined by the power method. MATLAB Compiler Version 4.1 (R2009a) was used to compile the MATLAB codes of BNTE. The implemented parallel algorithms were enabled with matlabpool, a Parallel Computing Toolbox function. The option UseParallel was set to 'always' and the default value of the option was 'never'. When those conditions held, the solvers computed estimated gradients in parallel. The parallel computing system was used to handle all the bottlenecks in the matrix generated from the finite element scheme and each domain of the power method generated. The parallel algorithm was implemented on a Symmetric Multi Processor (SMP) cluster machine, which had Intel 32 bit quad-core x 86 processors. Convergence rates and timings for the algorithm on the SMP cluster machine were obtained. Numerical experiments indicated the designed parallel algorithm could reach perfect speedup and had good stability and scalability. (au)

Classes of general axisymmetric solutions of Einstein-Maxwell equations

International Nuclear Information System (INIS)

Krori, K.D.; Choudhury, T.

1981-01-01

An exact solution of the Einstein equations for a stationary axially symmetric distribution of mass composed of all types of multipoles is obtained. Following Ernst (1968), from this vacuum solution the corresponding solution of the coupled Einstein-Maxwell equations is derived. A solution of Einstein-Maxwell fields for a static axially symmetric system composed of all types of multipoles is also obtained. (author)

On genus-two solutions for the ILW equation

Science.gov (United States)

Tutiya, Y.

2018-02-01

The existence of theta function solutions of genus two for the intermediate long-wave equation is established. A numerical example is also presented. The method basically goes along with Krichever's construction of theta function solutions for soliton equations, such as the Kronecker product equation. This idea leads us to a question whether a Riemann surface exists which allows a peculiar abelian integral of the third kind. The answer is affirmative at least for genus-two curves.

Calculation of similarity solutions of partial differential equations

International Nuclear Information System (INIS)

Dresner, L.

1980-08-01

When a partial differential equation in two independent variables is invariant to a group G of stretching transformations, it has similarity solutions that can be found by solving an ordinary differential equation. Under broad conditions, this ordinary differential equation is also invariant to another stretching group G', related to G. The invariance of the ordinary differential equation to G' can be used to simplify its solution, particularly if it is of second order. Then a method of Lie's can be used to reduce it to a first-order equation, the study of which is greatly facilitated by analysis of its direction field. The method developed here is applied to three examples: Blasius's equation for boundary layer flow over a flat plate and two nonlinear diffusion equations, cc/sub t/ = c/sub zz/ and c/sub t/ = (cc/sub z/)/sub z/

On solutions of variable-order fractional differential equations

Directory of Open Access Journals (Sweden)

Ali Akgül

2017-01-01

solutions to fractional differential equations are compelling to get in real applications, due to the nonlocality and complexity of the fractional differential operators, especially for variable-order fractional differential equations. Therefore, it is significant to enhanced numerical methods for fractional differential equations. In this work, we consider variable-order fractional differential equations by reproducing kernel method. There has been much attention in the use of reproducing kernels for the solutions to many problems in the recent years. We give two examples to demonstrate how efficiently our theory can be implemented in practice.

General solutions of second-order linear difference equations of Euler type

Directory of Open Access Journals (Sweden)

Akane Hongyo

2017-01-01

Full Text Available The purpose of this paper is to give general solutions of linear difference equations which are related to the Euler-Cauchy differential equation \\(y^{\\prime\\prime}+(\\lambda/t^2y=0\\ or more general linear differential equations. We also show that the asymptotic behavior of solutions of the linear difference equations are similar to solutions of the linear differential equations.

Exact Solutions for Nonlinear Differential Difference Equations in Mathematical Physics

Directory of Open Access Journals (Sweden)

Khaled A. Gepreel

2013-01-01

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

White noise solutions to the stochastic mKdV equation

International Nuclear Information System (INIS)

Zhang Zhongjun; Wei Caimin

2009-01-01

In this paper, we present the white noise solutions of the stochastic mKdV equation via the Hermite transformation and variable-coefficient generalized projected Ricatti equation expansion method. These solutions include white noise solitary wave solutions, white noise soliton-like solutions and white noise trigonometric function solutions.

Analytical solution and simplified analysis of coupled parent-daughter steady-state transport with multirate mass transfer

Science.gov (United States)

R. Haggerty

2013-01-01

In this technical note, a steady-state analytical solution of concentrations of a parent solute reacting to a daughter solute, both of which are undergoing transport and multirate mass transfer, is presented. Although the governing equations are complicated, the resulting solution can be expressed in simple terms. A function of the ratio of concentrations, In (daughter...

Extremal solutions of measure differential equations

Czech Academy of Sciences Publication Activity Database

Monteiro, Giselle Antunes; Slavík, A.

2016-01-01

Roč. 444, č. 1 (2016), s. 568-597 ISSN 0022-247X Institutional support: RVO:67985840 Keywords : measure differential equations * extremal solution * lower solution Subject RIV: BA - General Mathematics Impact factor: 1.064, year: 2016 http://www.sciencedirect.com/science/article/pii/S0022247X16302724

Theory of contributon transport

International Nuclear Information System (INIS)

Painter, J.W.; Gerstl, S.A.W.; Pomraning, G.C.

1980-10-01

A general discussion of the physics of contributon transport is presented. To facilitate this discussion, a Boltzmann-like transport equation for contributons is obtained, and special contributon cross sections are defined. However, the main goal of this study is to identify contributon transport equations and investigate possible deterministic solution techniques. Four approaches to the deterministic solution of the contributon transport problem are investigated. These approaches are an attempt to exploit certain attractive properties of the contributon flux, psi = phi phi + , where phi and phi + are the solutions to the forward and adjoint Boltzmann transport equations

Conformally flat spaces and solutions to Yang-Mills equations

International Nuclear Information System (INIS)

Chaohao, G.

1980-01-01

Using the conformal invariance of Yang-Mills equations in four-dimensional manifolds, it is proved that in a simply connected space of negative constant curvature Yang-Mills equations admit solutions with any real number as their Pontryagin number. It is also shown that the space S 3 x S 1 which is the regular counterpart of the meron solution is one example of a class of solutions to Yang-Mills equations on compact manifolds that are neither self-dual nor anti-self-dual

Multiloop soliton and multibreather solutions of the short pulse model equation

International Nuclear Information System (INIS)

Matsuno, Yoshimasa

2007-01-01

We develop a systematic procedure for constructing the multisoliton solutions of the short pulse (SP) model equation which describes the propagation of ultra-short pulses in nonlinear medica. We first introduce a novel hodograph transformation to convert the SP equation into the sine-Gordon (sG) equation. With the soliton solutions of the sG equation, the system of linear partial differential equations governing the inverse mapping can be integrated analytically to obtain the soliton solutions of the SP equation in the form of the parametric representation. By specifying the soliton parameters, we obtain the multiloop and multibreather solutions. We investigate the asymptotic behavior of both solutions and confirm their solitonic feature. The nonsingular breather solutions may play an important role in studying the propagation of ultra-short pulses in an optical fibre. (author)

Solute transport in a well under slow-purge and no-purge conditions

Science.gov (United States)

Plummer, M. A.; Britt, S. L.; Martin-Hayden, J. M.

2010-12-01

Non-purge sampling techniques, such as diffusion bags and in-situ sealed samplers, offer reliable and cost-effective groundwater monitoring methods that are a step closer to the goal of real-time monitoring without pumping or sample collection. Non-purge methods are, however, not yet completely accepted because questions remain about how solute concentrations in an unpurged well relate to concentrations in the adjacent formation. To answer questions about how undisturbed well water samples compare to formation concentrations, and to provide the information necessary to interpret results from non-purge monitoring systems, we have conducted a variety of physical experiments and numerical simulations of flow and transport in and through monitoring wells under low-flow and ambient flow conditions. Previous studies of flow and transport in wells used a Darcy’s law - based continuity equation for flow, which is often justified under the strong, forced-convection flow caused by pumping or large vertical hydraulic potential gradients. In our study, we focus on systems with weakly forced convection, where density-driven free convection may be of similar strength. We therefore solved Darcy’s law for porous media domains and the Navier Stokes equations for flow in the well, and coupled solution of the flow equations to that of solute transport. To illustrate expected in-well transport behavior under low-flow conditions, we present results of three particular studies: (1) time-dependent effluent concentrations from a well purged at low-flow pumping rates, (2) solute-driven density effects in a well under ambient horizontal flow and (3) temperature-driven mixing in a shallow well subject to seasonal temperature variations. Results of the first study illustrate that assumptions about the nature of in-well flow have a significant impact on effluent concentration curves even during pumping, with Poiseuille-type flow producing more rapid removal of concentration differences

Exact solutions to the Lienard equation and its applications

International Nuclear Information System (INIS)

Feng Zhaosheng

2004-01-01

In this paper, a kind of explicit exact solutions to the Lienard equation is obtained, and the applications of the result in seeking traveling solitary wave solution of the nonlinear Schroedinger equation are presented

New solutions of the confluent Heun equation

Directory of Open Access Journals (Sweden)

Harold Exton

1998-05-01

Full Text Available New compact triple series solutions of the confluent Heun equation (CHE are obtained by the appropriate applications of the Laplace transform and its inverse to a suitably constructed system of soluble differential equations. The computer-algebra package MAPLE V is used to tackle an auxiliary system of non-linear algebraic equations. This study is partly motivated by the relationship between the CHE and certain Schrödininger equations.

Nontrivial Periodic Solutions for Nonlinear Second-Order Difference Equations

Directory of Open Access Journals (Sweden)

Tieshan He

2011-01-01

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

Multi-component Wronskian solution to the Kadomtsev-Petviashvili equation

Science.gov (United States)

Xu, Tao; Sun, Fu-Wei; Zhang, Yi; Li, Juan

2014-01-01

It is known that the Kadomtsev-Petviashvili (KP) equation can be decomposed into the first two members of the coupled Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy by the binary non-linearization of Lax pairs. In this paper, we construct the N-th iterated Darboux transformation (DT) for the second- and third-order m-coupled AKNS systems. By using together the N-th iterated DT and Cramer's rule, we find that the KPII equation has the unreduced multi-component Wronskian solution and the KPI equation admits a reduced multi-component Wronskian solution. In particular, based on the unreduced and reduced two-component Wronskians, we obtain two families of fully-resonant line-soliton solutions which contain arbitrary numbers of asymptotic solitons as y → ∓∞ to the KPII equation, and the ordinary N-soliton solution to the KPI equation. In addition, we find that the KPI line solitons propagating in parallel can exhibit the bound state at the moment of collision.

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

Bifurcations of traveling wave solutions for an integrable equation

International Nuclear Information System (INIS)

Li Jibin; Qiao Zhijun

2010-01-01

This paper deals with the following equation m t =(1/2)(1/m k ) xxx -(1/2)(1/m k ) x , which is proposed by Z. J. Qiao [J. Math. Phys. 48, 082701 (2007)] and Qiao and Liu [Chaos, Solitons Fractals 41, 587 (2009)]. By adopting the phase analysis method of planar dynamical systems and the theory of the singular traveling wave systems to the traveling wave solutions of the equation, it is shown that for different k, the equation may have infinitely many solitary wave solutions, periodic wave solutions, kink/antikink wave solutions, cusped solitary wave solutions, and breaking loop solutions. We discuss in a detail the cases of k=-2,-(1/2),(1/2),2, and parametric representations of all possible bounded traveling wave solutions are given in the different (c,g)-parameter regions.

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

International Nuclear Information System (INIS)

Sun Yuhuai; Ma Zhimin; Li Yan

2010-01-01

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

Periodic and solitary-wave solutions of the Degasperis-Procesi equation

International Nuclear Information System (INIS)

Vakhnenko, V.O.; Parkes, E.J.

2004-01-01

Travelling-wave solutions of the Degasperis-Procesi equation are investigated. The solutions are characterized by two parameters. For propagation in the positive x-direction, hump-like, inverted loop-like and coshoidal periodic-wave solutions are found; hump-like, inverted loop-like and peakon solitary-wave solutions are obtained as well. For propagation in the negative x-direction, there are solutions which are just the mirror image in the x-axis of the aforementioned solutions. A transformed version of the Degasperis-Procesi equation, which is a generalization of the Vakhnenko equation, is also considered. For propagation in the positive x-direction, hump-like, loop-like, inverted loop-like, bell-like and coshoidal periodic-wave solutions are found; loop-like, inverted loop-like and kink-like solitary-wave solutions are obtained as well. For propagation in the negative x-direction, well-like and inverted coshoidal periodic-wave solutions are found; well-like and inverted peakon solitary-wave solutions are obtained as well. In an appropriate limit, the previously known solutions of the Vakhnenko equation are recovered

FORSIM-6, Automatic Solution of Coupled Differential Equation System

International Nuclear Information System (INIS)

Carver, M.B.; Stewart, D.G.; Blair, J.M.; Selander, W.N.

1983-01-01

1 - Description of problem or function: The FORSIM program is a versatile package which automates the solution of coupled differential equation systems. The independent variables are time, and up to three space coordinates, and the equations may be any mixture of partial and/or ordinary differential equations. The philosophy of the program is to provide a tool which will solve a system of differential equations for a user who has basic but unspecialized knowledge of numerical analysis and FORTRAN. The equations to be solved, together with the initial conditions and any special instructions, may be specified by the user in a single FORTRAN subroutine, although he may write a number of routines if this is more suitable. These are then loaded with the control routines, which perform the solution and any requested input and output. 2 - Method of solution: Partial differential equations are automatically converted into sets of coupled ordinary differential equations by variable order discretization in the spatial dimensions. These and other ordinary differential equations are integrated continuously in time using efficient variable order, variable step, error-controlled algorithms

Null solution of the Yang-Mills equations

International Nuclear Information System (INIS)

Tafel, J.

1986-05-01

We investigate the correspondence between null solutions of the Yang-Mills equations and shearfree geodesic null congruences. We give an example of a non-Abelian null solution with twisting rays. (orig.)

Solution of radial spin-1 field equation in Robertson-Walker space-time via Heun's equation

International Nuclear Information System (INIS)

Zecca, A.

2010-01-01

The spin-1 field equation is considered in Robertson-Walker spacetime. The problem of the solution of the separated radial equations, previously discussed in the flat space-time case, is solved also for both the closed and open curvature case. The radial equation is reduced to Heun's differential equation that recently has been widely reconsidered. It is shown that the solution of the present Heun equation does not fall into the class of polynomial-like or hypergeometric functions. Heun's operator results also non-factorisable. The properties follow from application of general theorems and power series expansion. In the positive curvature case of the universe a discrete energy spectrum of the system is found. The result follows by requiring a polynomial-like behaviour of at least one component of the spinor field. Developments and applications of the theory suggest further study of the solution of Heun's equation.

TRAVELING WAVE SOLUTIONS OF SOME FRACTIONAL DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

SERIFE MUGE EGE

2016-07-01

Full Text Available The modified Kudryashov method is powerful, efficient and can be used as an alternative to establish new solutions of different type of fractional differential equations applied in mathematical physics. In this article, we’ve constructed new traveling wave solutions including symmetrical Fibonacci function solutions, hyperbolic function solutions and rational solutions of the space-time fractional Cahn Hillihard equation D_t^α u − γD_x^α u − 6u(D_x^α u^2 − (3u^2 − 1D_x^α (D_x^α u + D_x^α(D_x^α(D_x^α(D_x^α u = 0 and the space-time fractional symmetric regularized long wave (SRLW equation D_t^α(D_t^α u + D_x^α(D_x^α u + uD_t^α(D_x^α u + D_x^α u D_t^α u + D_t^α(D_t^α(D_x^α(D_x^α u = 0 via modified Kudryashov method. In addition, some of the solutions are described in the figures with the help of Mathematica.

New correct solutions of the Dirac equation. 5

International Nuclear Information System (INIS)

Bagrov, V.G.; Byzov, N.N.; Gitman, D.M.; Klimenko, Yu.I.; Meshkov, A.G.; Shapovalov, V.N.; Shakhmatov, V.M.

1975-01-01

Some exact solutions for the Dirac equation, Klein-Gordon equation and classical relativistic equations of motion of an electron in external electromagnetic fields of a special type are considered. When fields E vector and H vector are related by the expression H vector=[n vector E vector]+n vector H 3 , where n vector is a constant unit vector, it turns out that among fields permitting the separation of variables in the Klein-Gordon equation more than half satisfy this relationship. For such fields the solution of the Dirac equation may be simplified considerably. Four specific kinds of fields are examined. The character of electron motion in such fields is peculiar but in the mathematical aspect, part of the problem is reduced to those considered previously

Solution of Moving Boundary Space-Time Fractional Burger’s Equation

Directory of Open Access Journals (Sweden)

E. A.-B. Abdel-Salam

2014-01-01

Full Text Available The fractional Riccati expansion method is used to solve fractional differential equations with variable coefficients. To illustrate the effectiveness of the method, the moving boundary space-time fractional Burger’s equation is studied. The obtained solutions include generalized trigonometric and hyperbolic function solutions. Among these solutions, some are found for the first time. The linear and periodic moving boundaries for the kink solution of the Burger’s equation are presented graphically and discussed.

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

Remarks about singular solutions to the Dirac equation

International Nuclear Information System (INIS)

Uhlir, M.

1975-01-01

In the paper singular solutions of the Dirac equation are investigated. They are derived in the Lorentz-covariant way of functions proportional to static multipole fields of scalar and (or) electromagnetic fields and of regular solutions of the Dirac equations. The regularization procedure excluding divergences of total energy, momentum and angular momentum of the spinor field considered is proposed

Exact solution of some linear matrix equations using algebraic methods

Science.gov (United States)

Djaferis, T. E.; Mitter, S. K.

1977-01-01

A study is done of solution methods for Linear Matrix Equations including Lyapunov's equation, using methods of modern algebra. The emphasis is on the use of finite algebraic procedures which are easily implemented on a digital computer and which lead to an explicit solution to the problem. The action f sub BA is introduced a Basic Lemma is proven. The equation PA + BP = -C as well as the Lyapunov equation are analyzed. Algorithms are given for the solution of the Lyapunov and comment is given on its arithmetic complexity. The equation P - A'PA = Q is studied and numerical examples are given.

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

Numerical solution of Boltzmann's equation

International Nuclear Information System (INIS)

Sod, G.A.

1976-04-01

The numerical solution of Boltzmann's equation is considered for a gas model consisting of rigid spheres by means of Hilbert's expansion. If only the first two terms of the expansion are retained, Boltzmann's equation reduces to the Boltzmann-Hilbert integral equation. Successive terms in the Hilbert expansion are obtained by solving the same integral equation with a different source term. The Boltzmann-Hilbert integral equation is solved by a new very fast numerical method. The success of the method rests upon the simultaneous use of four judiciously chosen expansions; Hilbert's expansion for the distribution function, another expansion of the distribution function in terms of Hermite polynomials, the expansion of the kernel in terms of the eigenvalues and eigenfunctions of the Hilbert operator, and an expansion involved in solving a system of linear equations through a singular value decomposition. The numerical method is applied to the study of the shock structure in one space dimension. Numerical results are presented for Mach numbers of 1.1 and 1.6. 94 refs, 7 tables, 1 fig

On nonlinear differential equation with exact solutions having various pole orders

International Nuclear Information System (INIS)

Kudryashov, N.A.

2015-01-01

We consider a nonlinear ordinary differential equation having solutions with various movable pole order on the complex plane. We show that the pole order of exact solution is determined by values of parameters of the equation. Exact solutions in the form of the solitary waves for the second order nonlinear differential equation are found taking into account the method of the logistic function. Exact solutions of differential equations are discussed and analyzed

Transport equation and shock waves

International Nuclear Information System (INIS)

Besnard, D.

1981-04-01

A multi-group method is derived from a one dimensional transport equation for the slowing down and spatial transport of energetic positive ions in a plasma. This method is used to calculate the behaviour of energetic charged particles in non homogeneous and non stationary plasma, and the effect of energy deposition of the particles on the heating of the plasma. In that purpose, an equation for the density of fast ions is obtained from the Fokker-Planck equation, and a closure condition for the second moment of this equation is deduced from phenomenological considerations. This method leads to a numerical method, simple and very efficient, which doesn't require much computer storage. Two types of numerical results are obtained. First, results on the slowing down of 3.5 MeV alpha particles in a 50 keV plasma plublished by Corman and al and Moses are compared with the results obtained with both our method and a Monte Carlo type method. Good agreement was obtained, even for energy deposition on the ions of the plasma. Secondly, we have calculated propagation of alpha particles heating a cold plasma. These results are in very good agreement with those given by an accurate Monte Carlo method, for both the thermal velocity, and the energy deposition in the plasma

Efficient traveltime solutions of the acoustic TI eikonal equation

KAUST Repository

Waheed, Umair bin

2015-02-01

Numerical solutions of the eikonal (Hamilton-Jacobi) equation for transversely isotropic (TI) media are essential for imaging and traveltime tomography applications. Such solutions, however, suffer from the inherent higher-order nonlinearity of the TI eikonal equation, which requires solving a quartic polynomial for every grid point. Analytical solutions of the quartic polynomial yield numerically unstable formulations. Thus, it requires a numerical root finding algorithm, adding significantly to the computational load. Using perturbation theory we approximate, in a first order discretized form, the TI eikonal equation with a series of simpler equations for the coefficients of a polynomial expansion of the eikonal solution, in terms of the anellipticity anisotropy parameter. Such perturbation, applied to the discretized form of the eikonal equation, does not impose any restrictions on the complexity of the perturbed parameter field. Therefore, it provides accurate traveltime solutions even for models with complex distribution of velocity and anisotropic anellipticity parameter, such as that for the complicated Marmousi model. The formulation allows for large cost reduction compared to using the direct TI eikonal solver. Furthermore, comparative tests with previously developed approximations illustrate remarkable gain in accuracy in the proposed algorithm, without any addition to the computational cost.

Elliptic and solitary wave solutions for Bogoyavlenskii equations system, couple Boiti-Leon-Pempinelli equations system and Time-fractional Cahn-Allen equation

Directory of Open Access Journals (Sweden)

Mostafa M.A. Khater

Full Text Available In this article and for the first time, we introduce and describe Khater method which is a new technique for solving nonlinear partial differential equations (PDEs.. We apply this method for each of the following models Bogoyavlenskii equation, couple Boiti-Leon-Pempinelli system and Time-fractional Cahn-Allen equation. Khater method is very powerful, Effective, felicitous and fabulous method to get exact and solitary wave solution of (PDEs.. Not only just like that but it considers too one of the general methods for solving that kind of equations since it involves some methods as we will see in our discuss of the results. We make a comparison between the results of this new method and another method. Keywords: Bogoyavlenskii equations system, Couple Boiti-Leon-Pempinelli equations system, Time-fractional Cahn-Allen equation, Khater method, Traveling wave solutions, Solitary wave solutions

ALMOST PERIODIC SOLUTIONS TO SOME NONLINEAR DELAY DIFFERENTIAL EQUATION

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

The existence of an almost periodic solutions to a nonlinear delay diffierential equation is considered in this paper. A set of sufficient conditions for the existence and uniqueness of almost periodic solutions to some delay diffierential equations is obtained.

Wronskians, generalized Wronskians and solutions to the Korteweg-de Vries equation

International Nuclear Information System (INIS)

Ma Wenxiu

2004-01-01

A bridge going from Wronskian solutions to generalized Wronskian solutions of the Korteweg-de Vries (KdV) equation is built. It is then shown that generalized Wronskian solutions can be viewed as Wronskian solutions. The idea is used to generate positons, negatons and their interaction solutions to the KdV equation. Moreover, general positons and negatons are constructed through the Wronskian formulation. A few new exact solutions to the KdV equation are explicitly presented as examples of Wronskian solutions

The relationship among the solutions of two auxiliary ordinary differential equations

International Nuclear Information System (INIS)

Liu Xiaoping; Liu Chunping

2009-01-01

In a recent article [Phys. Lett. A 356 (2006) 124], Sirendaoreji extended their auxiliary equation method by introducing a new auxiliary ordinary differential equation (NAODE) and its 14 solutions. Then the author studied some nonlinear evolution equations (NLEEs) and got more exact travelling wave solutions. In this paper, we will show that the 14 solutions of the NAODE are actually the same as the solutions obtained by original auxiliary equation method, and they are only different in the form.

Minimal solution for inconsistent singular fuzzy matrix equations

Directory of Open Access Journals (Sweden)

M. Nikuie

2013-10-01

Full Text Available The fuzzy matrix equations $Ailde{X}=ilde{Y}$ is called a singular fuzzy matrix equations while the coefficients matrix of its equivalent crisp matrix equations be a singular matrix. The singular fuzzy matrix equations are divided into two parts: consistent singular matrix equations and inconsistent fuzzy matrix equations. In this paper, the inconsistent singular fuzzy matrix equations is studied and the effect of generalized inverses in finding minimal solution of an inconsistent singular fuzzy matrix equations are investigated.

Numerical solution of a reaction-diffusion equation

International Nuclear Information System (INIS)

Moyano, Edgardo A.; Scarpettini, Alberto F.

2000-01-01

The purpose of the present work to continue the observations and the numerical experiences on a reaction-diffusion model, that is a simplified form of the neutronic flux equation. The model is parabolic, nonlinear, with Dirichlet boundary conditions. The purpose is to approximate non trivial solutions, asymptotically stables for t → ∞, that is solutions that tend to the elliptic problem, in the Lyapunov sense. It belongs to the so-called reaction-diffusion equations of semi linear kind, that is, linear equations in the heat operator and they have a nonlinear reaction function, in this case f (u, a, b) = u (a - b u), being u concentration, a and b parameters. The study of the incidence of these parameters take an interest to the neutronic flux physics. So that we search non trivial, positive and bounded solutions. The used algorithm is based on the concept of monotone and ordered sequences, and on the existence theorem of Amann and Sattinger. (author)

A numerical model for the solution of the Shallow Water equations in composite channels with movable bed

Science.gov (United States)

minatti, L.

2013-12-01

A finite volume model solving the shallow water equations coupled with the sediments continuity equation in composite channels with irregular geometry is presented. The model is essentially 1D but can handle composite cross-sections in which bedload transport is considered to occur inside the main channel only. This assumption is coherent with the observed behavior of rivers on short time scales where main channel areas exhibit more relevant morphological variations than overbanks. Furthermore, such a model allows a more precise prediction of thalweg elevation and cross section shape variations than fully 1D models where bedload transport is considered to occur uniformly over the entire cross section. The coupling of the equations describing water and sediments dynamics results in a hyperbolic non-conservative system that cannot be solved numerically with the use of a conservative scheme. Therefore, a path-conservative scheme, based on the approach proposed by Pares and Castro (2004) has been devised in order to account for the coupling with the sediments continuity equation and for the concurrent presence of bottom elevation and breadth variations of the cross section. In order to correctly compute numerical fluxes related to bedload transport in main channel areas, a special treatment of the equations is employed in the model. The resulting scheme is well balanced and fully coupled and can accurately model abrupt time variations of flow and bedload transport conditions in wide rivers, characterized by the presence of overbank areas that are less active than the main channel. The accuracy of the model has been first tested in fixed bed conditions by solving problems with a known analytical solution: in these tests the model proved to be able to handle shocks and supercritical flow conditions properly(see Fig. 01). A practical application of the model to the Ombrone river, southern Tuscany (Italy) is shown. The river has shown relevant morphological changes during

Exact solutions to the Mo-Papas and Landau-Lifshitz equations

Science.gov (United States)

Rivera, R.; Villarroel, D.

2002-10-01

Two exact solutions of the Mo-Papas and Landau-Lifshitz equations for a point charge in classical electrodynamics are presented here. Both equations admit as an exact solution the motion of a charge rotating with constant speed in a circular orbit. These equations also admit as an exact solution the motion of two identical charges rotating with constant speed at the opposite ends of a diameter. These exact solutions allow one to obtain, starting from the equation of motion, a definite formula for the rate of radiation. In both cases the rate of radiation can also be obtained, with independence of the equation of motion, from the well known fields of a point charge, that is, from the Maxwell equations. The rate of radiation obtained from the Mo-Papas equation in the one-charge case coincides with the rate of radiation that comes from the Maxwell equations; but in the two-charge case the results do not coincide. On the other hand, the rate of radiation obtained from the Landau-Lifshitz equation differs from the one that follows from the Maxwell equations in both the one-charge and two-charge cases. This last result does not support a recent statement by Rohrlich in favor of considering the Landau-Lifshitz equation as the correct and exact equation of motion for a point charge in classical electrodynamics.

Exact solutions to the Mo-Papas and Landau-Lifshitz equations

International Nuclear Information System (INIS)

Rivera, R.; Villarroel, D.

2002-01-01

Two exact solutions of the Mo-Papas and Landau-Lifshitz equations for a point charge in classical electrodynamics are presented here. Both equations admit as an exact solution the motion of a charge rotating with constant speed in a circular orbit. These equations also admit as an exact solution the motion of two identical charges rotating with constant speed at the opposite ends of a diameter. These exact solutions allow one to obtain, starting from the equation of motion, a definite formula for the rate of radiation. In both cases the rate of radiation can also be obtained, with independence of the equation of motion, from the well known fields of a point charge, that is, from the Maxwell equations. The rate of radiation obtained from the Mo-Papas equation in the one-charge case coincides with the rate of radiation that comes from the Maxwell equations; but in the two-charge case the results do not coincide. On the other hand, the rate of radiation obtained from the Landau-Lifshitz equation differs from the one that follows from the Maxwell equations in both the one-charge and two-charge cases. This last result does not support a recent statement by Rohrlich in favor of considering the Landau-Lifshitz equation as the correct and exact equation of motion for a point charge in classical electrodynamics

Exact solutions of space-time fractional EW and modified EW equations

International Nuclear Information System (INIS)

Korkmaz, Alper

2017-01-01

The bright soliton solutions and singular solutions are constructed for the space-time fractional EW and the space-time fractional modified EW (MEW) equations. Both equations are reduced to ordinary differential equations by the use of fractional complex transform (FCT) and properties of modified Riemann–Liouville derivative. Then, various ansatz method are implemented to construct the solutions for both equations.

Solutions manual to accompany Ordinary differential equations

CERN Document Server

Greenberg, Michael D

2014-01-01

Features a balance between theory, proofs, and examples and provides applications across diverse fields of study Ordinary Differential Equations presents a thorough discussion of first-order differential equations and progresses to equations of higher order. The book transitions smoothly from first-order to higher-order equations, allowing readers to develop a complete understanding of the related theory. Featuring diverse and interesting applications from engineering, bioengineering, ecology, and biology, the book anticipates potential difficulties in understanding the various solution steps

PARALLEL SOLUTION METHODS OF PARTIAL DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

Korhan KARABULUT

1998-03-01

Full Text Available Partial differential equations arise in almost all fields of science and engineering. Computer time spent in solving partial differential equations is much more than that of in any other problem class. For this reason, partial differential equations are suitable to be solved on parallel computers that offer great computation power. In this study, parallel solution to partial differential equations with Jacobi, Gauss-Siedel, SOR (Succesive OverRelaxation and SSOR (Symmetric SOR algorithms is studied.

Self-similar solutions of the modified nonlinear schrodinger equation

International Nuclear Information System (INIS)

Kitaev, A.V.

1986-01-01

This paper considers a 2 x 2 matrix linear ordinary differential equation with large parameter t and irregular singular point of fourth order at infinity. The leading order of the monodromy data of this equation is calculated in terms of its coefficients. Isomonodromic deformations of the equation are self-similar solutions of the modified nonlinear Schrodinger equation, and therefore inversion of the expressions obtained for the monodromy data gives the leading term in the time-asymptotic behavior of the self-similar solution. The application of these results to the type IV Painleve equation is considered in detail

Similarity Solutions for Multiterm Time-Fractional Diffusion Equation

Directory of Open Access Journals (Sweden)

A. Elsaid

2016-01-01

Full Text Available Similarity method is employed to solve multiterm time-fractional diffusion equation. The orders of the fractional derivatives belong to the interval (0,1] and are defined in the Caputo sense. We illustrate how the problem is reduced from a multiterm two-variable fractional partial differential equation to a multiterm ordinary fractional differential equation. Power series solution is obtained for the resulting ordinary problem and the convergence of the series solution is discussed. Based on the obtained results, we propose a definition for a multiterm error function with generalized coefficients.

Periodic Solutions and S-Asymptotically Periodic Solutions to Fractional Evolution Equations

Directory of Open Access Journals (Sweden)

Jia Mu

2017-01-01

Full Text Available This paper deals with the existence and uniqueness of periodic solutions, S-asymptotically periodic solutions, and other types of bounded solutions for some fractional evolution equations with the Weyl-Liouville fractional derivative defined for periodic functions. Applying Fourier transform we give reasonable definitions of mild solutions. Then we accurately estimate the spectral radius of resolvent operator and obtain some existence and uniqueness results.

Dynamic modeling of interfacial structures via interfacial area transport equation

International Nuclear Information System (INIS)