Linear and quasi-linear equations of parabolic type

CERN Document Server

Ladyženskaja, O A; Ural′ceva, N N; Uralceva, N N

1968-01-01

Equations of parabolic type are encountered in many areas of mathematics and mathematical physics, and those encountered most frequently are linear and quasi-linear parabolic equations of the second order. In this volume, boundary value problems for such equations are studied from two points of view: solvability, unique or otherwise, and the effect of smoothness properties of the functions entering the initial and boundary conditions on the smoothness of the solutions.

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

Linear determining equations for differential constraints

International Nuclear Information System (INIS)

Kaptsov, O V

1998-01-01

A construction of differential constraints compatible with partial differential equations is considered. Certain linear determining equations with parameters are used to find such differential constraints. They generalize the classical determining equations used in the search for admissible Lie operators. As applications of this approach equations of an ideal incompressible fluid and non-linear heat equations are discussed

Quantum Non-Markovian Langevin Equations and Transport Coefficients

International Nuclear Information System (INIS)

Sargsyan, V.V.; Antonenko, N.V.; Kanokov, Z.; Adamian, G.G.

2005-01-01

Quantum diffusion equations featuring explicitly time-dependent transport coefficients are derived from generalized non-Markovian Langevin equations. Generalized fluctuation-dissipation relations and analytic expressions for calculating the friction and diffusion coefficients in nuclear processes are obtained. The asymptotic behavior of the transport coefficients and correlation functions for a damped harmonic oscillator that is linearly coupled in momentum to a heat bath is studied. The coupling to a heat bath in momentum is responsible for the appearance of the diffusion coefficient in coordinate. The problem of regression of correlations in quantum dissipative systems is analyzed

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.

Analysis of discrete reaction-diffusion equations for autocatalysis and continuum diffusion equations for transport

Energy Technology Data Exchange (ETDEWEB)

Wang, Chi-Jen [Iowa State Univ., Ames, IA (United States)

2013-01-01

In this thesis, we analyze both the spatiotemporal behavior of: (A) non-linear “reaction” models utilizing (discrete) reaction-diffusion equations; and (B) spatial transport problems on surfaces and in nanopores utilizing the relevant (continuum) diffusion or Fokker-Planck equations. Thus, there are some common themes in these studies, as they all involve partial differential equations or their discrete analogues which incorporate a description of diffusion-type processes. However, there are also some qualitative differences, as shall be discussed below.

Comparison of neutronic transport equation resolution nodal methods

International Nuclear Information System (INIS)

Zamonsky, O.M.; Gho, C.J.

1990-01-01

In this work, some transport equation resolution nodal methods are comparatively studied: the constant-constant (CC), linear-nodal (LN) and the constant-quadratic (CQ). A nodal scheme equivalent to finite differences has been used for its programming, permitting its inclusion in existing codes. Some bidimensional problems have been solved, showing that linear-nodal (LN) are, in general, obtained with accuracy in CPU shorter times. (Author) [es

Isomorphism of Intransitive Linear Lie Equations

Directory of Open Access Journals (Sweden)

Jose Miguel Martins Veloso

2009-11-01

Full Text Available We show that formal isomorphism of intransitive linear Lie equations along transversal to the orbits can be extended to neighborhoods of these transversal. In analytic cases, the word formal is dropped from theorems. Also, we associate an intransitive Lie algebra with each intransitive linear Lie equation, and from the intransitive Lie algebra we recover the linear Lie equation, unless of formal isomorphism. The intransitive Lie algebra gives the structure functions introduced by É. Cartan.

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

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)

Solution of the linear transport equation, monoenergetic in multiregions with anisotopic scattering by the method F sub(N)

International Nuclear Information System (INIS)

Pontedeiro, E.M.B.D.; Maiorino, J.R.

1982-01-01

The linear equation transport, monoenergetic, with anysotropic scattering, in multiregions, by F sub(N) method, is resolved. The mathematical analysis used for this method consists in to use parcially the expansion method in singular autofunctions, or Case's method, aiming to derive a set of integral equations coupled to the angular distribution in the boundaries and interfaces, and then to approximate these distributions by polynomics of N order, aiming to derive, with the use of these boundary and continuity conditions in the interfaces, a set of algebric equations for the coef. of polynomical approximation. With the goal to obtain numerical results, a computer code (FNAM-1) with options for the number of regions, boundary conditions, F sub(N) approx order, were developed. Numerical results were then obtained for various sample problems and compared with the results published in the literature with the objective to demonstrate the precision and applicability of the F sub(N) method. (E.G.) [pt

Students’ difficulties in solving linear equation problems

Science.gov (United States)

Wati, S.; Fitriana, L.; Mardiyana

2018-03-01

A linear equation is an algebra material that exists in junior high school to university. It is a very important material for students in order to learn more advanced mathematics topics. Therefore, linear equation material is essential to be mastered. However, the result of 2016 national examination in Indonesia showed that students’ achievement in solving linear equation problem was low. This fact became a background to investigate students’ difficulties in solving linear equation problems. This study used qualitative descriptive method. An individual written test on linear equation tasks was administered, followed by interviews. Twenty-one sample students of grade VIII of SMPIT Insan Kamil Karanganyar did the written test, and 6 of them were interviewed afterward. The result showed that students with high mathematics achievement donot have difficulties, students with medium mathematics achievement have factual difficulties, and students with low mathematics achievement have factual, conceptual, operational, and principle difficulties. Based on the result there is a need of meaningfulness teaching strategy to help students to overcome difficulties in solving linear equation problems.

TLC scheme for numerical solution of the transport equation on equilateral triangular meshes

International Nuclear Information System (INIS)

Walters, W.F.

1983-01-01

A new triangular linear characteristic TLC scheme for numerically solving the transport equation on equilateral triangular meshes has been developed. This scheme uses the analytic solution of the transport equation in the triangle as its basis. The data on edges of the triangle are assumed linear as is the source representation. A characteristic approach or nodal approach is used to obtain the analytic solution. Test problems indicate that the new TLC is superior to the widely used DITRI scheme for accuracy

Computing with linear equations and matrices

International Nuclear Information System (INIS)

Churchhouse, R.F.

1983-01-01

Systems of linear equations and matrices arise in many disciplines. The equations may accurately represent conditions satisfied by a system or, more likely, provide an approximation to a more complex system of non-linear or differential equations. The system may involve a few or many thousand unknowns and each individual equation may involve few or many of them. Over the past 50 years a vast literature on methods for solving systems of linear equations and the associated problems of finding the inverse or eigenvalues of a matrix has been produced. These lectures cover those methods which have been found to be most useful for dealing with such types of problem. References are given where appropriate and attention is drawn to the possibility of improved methods for use on vector and parallel processors. (orig.)

Correct Linearization of Einstein's Equations

Directory of Open Access Journals (Sweden)

Rabounski D.

2006-06-01

Full Text Available Regularly Einstein's equations can be reduced to a wave form (linearly dependent from the second derivatives of the space metric in the absence of gravitation, the space rotation and Christoffel's symbols. As shown here, the origin of the problem is that one uses the general covariant theory of measurement. Here the wave form of Einstein's equations is obtained in the terms of Zelmanov's chronometric invariants (physically observable projections on the observer's time line and spatial section. The obtained equations depend on solely the second derivatives even if gravitation, the space rotation and Christoffel's symbols. The correct linearization proves: the Einstein equations are completely compatible with weak waves of the metric.

Inverse problems in linear transport theory

International Nuclear Information System (INIS)

Dressler, K.

1988-01-01

Inverse problems for a class of linear kinetic equations are investigated. The aim is to identify the scattering kernel of a transport equation (corresponding to the structure of a background medium) by observing the 'albedo' part of the solution operator for the corresponding direct initial boundary value problem. This means to get information on some integral operator in an integrodifferential equation through on overdetermined boundary value problem. We first derive a constructive method for solving direct halfspace problems and prove a new factorization theorem for the solutions. Using this result we investigate stationary inverse problems with respect to well posedness (e.g. reduce them to classical ill-posed problems, such as integral equations of first kind). In the time-dependent case we show that a quite general inverse problem is well posed and solve it constructively. (orig.)

Spectral theories for linear differential equations

International Nuclear Information System (INIS)

Sell, G.R.

1976-01-01

The use of spectral analysis in the study of linear differential equations with constant coefficients is not only a fundamental technique but also leads to far-reaching consequences in describing the qualitative behaviour of the solutions. The spectral analysis, via the Jordan canonical form, will not only lead to a representation theorem for a basis of solutions, but will also give a rather precise statement of the (exponential) growth rates of various solutions. Various attempts have been made to extend this analysis to linear differential equations with time-varying coefficients. The most complete such extensions is the Floquet theory for equations with periodic coefficients. For time-varying linear differential equations with aperiodic coefficients several authors have attempted to ''extend'' the Foquet theory. The precise meaning of such an extension is itself a problem, and we present here several attempts in this direction that are related to the general problem of extending the spectral analysis of equations with constant coefficients. The main purpose of this paper is to introduce some problems of current research. The primary problem we shall examine occurs in the context of linear differential equations with almost periodic coefficients. We call it ''the Floquet problem''. (author)

Lie algebras and linear differential equations.

Science.gov (United States)

Brockett, R. W.; Rahimi, A.

1972-01-01

Certain symmetry properties possessed by the solutions of linear differential equations are examined. For this purpose, some basic ideas from the theory of finite dimensional linear systems are used together with the work of Wei and Norman on the use of Lie algebraic methods in differential equation theory.

Alternative formulation of the monokinetic transport equation

International Nuclear Information System (INIS)

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

1985-01-01

After recalling a technique already exploited in stationary neutron transport, the dynamic linear monokinetic equation for general geometry is cast into an integro-differential form where a second order space Laplace operator and both a second and first time derivatives appear. The introduced unknowns are given a physical interpretation for plane geometry and their relations with the total flux and current are derived

New nonlinear methods for linear transport calculations

International Nuclear Information System (INIS)

Adams, M.L.

1993-01-01

We present a new family of methods for the numerical solution of the linear transport equation. With these methods an iteration consists of an 'S N sweep' followed by an 'S 2 -like' calculation. We show, by analysis as well as numerical results, that iterative convergence is always rapid. We show that this rapid convergence does not depend on a consistent discretization of the S 2 -like equations - they can be discretized independently from the S N equations. We show further that independent discretizations can offer significant advantages over consistent ones. In particular, we find that in a wide range of problems, an accurate discretization of the S 2 -like equation can be combined with a crude discretization of the S N equations to produce an accurate S N answer. We demonstrate this by analysis as well as numerical results. (orig.)

Linear q-nonuniform difference equations

International Nuclear Information System (INIS)

Bangerezako, Gaspard

2010-01-01

We introduce basic concepts of q-nonuniform differentiation and integration and study linear q-nonuniform difference equations and systems, as well as their application in q-nonuniform difference linear control systems. (author)

Schwarz maps of algebraic linear ordinary differential equations

Science.gov (United States)

Sanabria Malagón, Camilo

2017-12-01

A linear ordinary differential equation is called algebraic if all its solution are algebraic over its field of definition. In this paper we solve the problem of finding closed form solution to algebraic linear ordinary differential equations in terms of standard equations. Furthermore, we obtain a method to compute all algebraic linear ordinary differential equations with rational coefficients by studying their associated Schwarz map through the Picard-Vessiot Theory.

Basic linear partial differential equations

CERN Document Server

Treves, Francois

1975-01-01

Focusing on the archetypes of linear partial differential equations, this text for upper-level undergraduates and graduate students features most of the basic classical results. The methods, however, are decidedly nontraditional: in practically every instance, they tend toward a high level of abstraction. This approach recalls classical material to contemporary analysts in a language they can understand, as well as exploiting the field's wealth of examples as an introduction to modern theories.The four-part treatment covers the basic examples of linear partial differential equations and their

Non-local quasi-linear parabolic equations

International Nuclear Information System (INIS)

Amann, H

2005-01-01

This is a survey of the most common approaches to quasi-linear parabolic evolution equations, a discussion of their advantages and drawbacks, and a presentation of an entirely new approach based on maximal L p regularity. The general results here apply, above all, to parabolic initial-boundary value problems that are non-local in time. This is illustrated by indicating their relevance for quasi-linear parabolic equations with memory and, in particular, for time-regularized versions of the Perona-Malik equation of image processing

Relativistic transport equation for a discontinuity wave of multiplicity one

Energy Technology Data Exchange (ETDEWEB)

Giambo, S; Palumbo, A [Istituto di Matematica, Universita degli Studi, Messina (Italy)

1980-04-14

In the framework of the theory of the singular hypersurfaces, the transport equation for the amplitude of a discontinuity wave, corresponding to a simple characteristic of a quasi-linear hyperbolic system, is established in the context of special relativity.

The recovery of a time-dependent point source in a linear transport equation: application to surface water pollution

International Nuclear Information System (INIS)

Hamdi, Adel

2009-01-01

The aim of this paper is to localize the position of a point source and recover the history of its time-dependent intensity function that is both unknown and constitutes the right-hand side of a 1D linear transport equation. Assuming that the source intensity function vanishes before reaching the final control time, we prove that recording the state with respect to the time at two observation points framing the source region leads to the identification of the source position and the recovery of its intensity function in a unique manner. Note that at least one of the two observation points should be strategic. We establish an identification method that determines quasi-explicitly the source position and transforms the task of recovering its intensity function into solving directly a well-conditioned linear system. Some numerical experiments done on a variant of the water pollution BOD model are presented

The Cauchy problem for non-linear Klein-Gordon equations

International Nuclear Information System (INIS)

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

1993-01-01

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

Hamiltonian structures of some non-linear evolution equations

International Nuclear Information System (INIS)

Tu, G.Z.

1983-06-01

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

Solving polynomial differential equations by transforming them to linear functional-differential equations

OpenAIRE

Nahay, John Michael

2008-01-01

We present a new approach to solving polynomial ordinary differential equations by transforming them to linear functional equations and then solving the linear functional equations. We will focus most of our attention upon the first-order Abel differential equation with two nonlinear terms in order to demonstrate in as much detail as possible the computations necessary for a complete solution. We mention in our section on further developments that the basic transformation idea can be generali...

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

OpenAIRE

Nair, C. Radhakrishnan

2004-01-01

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

Solving the equation of neutron transport

International Nuclear Information System (INIS)

Nasfi, Rim

2009-01-01

This work is devoted to the study of some numerical methods of resolution of the problem of transport of the neutrons. We started by introducing the equation integro-differential transport of the neutrons. Then we applied the finite element method traditional for stationary and nonstationary linear problems in 2D. A great part is reserved for the presentation of the mixed numerical diagram and mixed hybrid with two types of uniform grids: triangular and rectangular. Thereafter we treated some numerical examples by implementations in Matlab in order to test the convergence of each method. To finish, we had results of simulation by the Monte Carlo method on a problem of two-dimensional transport with an aim of comparing them with the results resulting from the finite element method mixed hybrids. Some remarks and prospects conclude this work.

Linear measure functional differential equations with infinite delay

OpenAIRE

Monteiro, G. (Giselle Antunes); Slavík, A.

2014-01-01

We use the theory of generalized linear ordinary differential equations in Banach spaces to study linear measure functional differential equations with infinite delay. We obtain new results concerning the existence, uniqueness, and continuous dependence of solutions. Even for equations with a finite delay, our results are stronger than the existing ones. Finally, we present an application to functional differential equations with impulses.

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

Science.gov (United States)

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

1972-01-01

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

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

Linear Discontinuous Expansion Method using the Subcell Balances for Unstructured Geometry SN Transport

International Nuclear Information System (INIS)

Hong, Ser Gi; Kim, Jong Woon; Lee, Young Ouk; Kim, Kyo Youn

2010-01-01

The subcell balance methods have been developed for one- and two-dimensional SN transport calculations. In this paper, a linear discontinuous expansion method using sub-cell balances (LDEM-SCB) is developed for neutral particle S N transport calculations in 3D unstructured geometrical problems. At present, this method is applied to the tetrahedral meshes. As the name means, this method assumes the linear distribution of the particle flux in each tetrahedral mesh and uses the balance equations for four sub-cells of each tetrahedral mesh to obtain the equations for the four sub-cell average fluxes which are unknowns. This method was implemented in the computer code MUST (Multi-group Unstructured geometry S N Transport). The numerical tests show that this method gives more robust solution than DFEM (Discontinuous Finite Element Method)

Linearized pseudo-Einstein equations on the Heisenberg group

Science.gov (United States)

Barletta, Elisabetta; Dragomir, Sorin; Jacobowitz, Howard

2017-02-01

We study the pseudo-Einstein equation R11bar = 0 on the Heisenberg group H1 = C × R. We consider first order perturbations θɛ =θ0 + ɛ θ and linearize the pseudo-Einstein equation about θ0 (the canonical Tanaka-Webster flat contact form on H1 thought of as a strictly pseudoconvex CR manifold). If θ =e2uθ0 the linearized pseudo-Einstein equation is Δb u - 4 | Lu|2 = 0 where Δb is the sublaplacian of (H1 ,θ0) and L bar is the Lewy operator. We solve the linearized pseudo-Einstein equation on a bounded domain Ω ⊂H1 by applying subelliptic theory i.e. existence and regularity results for weak subelliptic harmonic maps. We determine a solution u to the linearized pseudo-Einstein equation, possessing Heisenberg spherical symmetry, and such that u(x) → - ∞ as | x | → + ∞.

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

Two-scale approach to oscillatory singularly perturbed transport equations

CERN Document Server

Frénod, Emmanuel

2017-01-01

This book presents the classical results of the two-scale convergence theory and explains – using several figures – why it works. It then shows how to use this theory to homogenize ordinary differential equations with oscillating coefficients as well as oscillatory singularly perturbed ordinary differential equations. In addition, it explores the homogenization of hyperbolic partial differential equations with oscillating coefficients and linear oscillatory singularly perturbed hyperbolic partial differential equations. Further, it introduces readers to the two-scale numerical methods that can be built from the previous approaches to solve oscillatory singularly perturbed transport equations (ODE and hyperbolic PDE) and demonstrates how they can be used efficiently. This book appeals to master’s and PhD students interested in homogenization and numerics, as well as to the Iter community.

Linear causal modeling with structural equations

CERN Document Server

Mulaik, Stanley A

2009-01-01

Emphasizing causation as a functional relationship between variables that describe objects, Linear Causal Modeling with Structural Equations integrates a general philosophical theory of causation with structural equation modeling (SEM) that concerns the special case of linear causal relations. In addition to describing how the functional relation concept may be generalized to treat probabilistic causation, the book reviews historical treatments of causation and explores recent developments in experimental psychology on studies of the perception of causation. It looks at how to perceive causal

Transport equations, Level Set and Eulerian mechanics. Application to fluid-structure coupling

International Nuclear Information System (INIS)

Maitre, E.

2008-11-01

My works were devoted to numerical analysis of non-linear elliptic-parabolic equations, to neutron transport equation and to the simulation of fabrics draping. More recently I developed an Eulerian method based on a level set formulation of the immersed boundary method to deal with fluid-structure coupling problems arising in bio-mechanics. Some of the more efficient algorithms to solve the neutron transport equation make use of the splitting of the transport operator taking into account its characteristics. In the present work we introduced a new algorithm based on this splitting and an adaptation of minimal residual methods to infinite dimensional case. We present the case where the velocity space is of dimension 1 (slab geometry) and 2 (plane geometry) because the splitting is simpler in the former

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.

Systems of Inhomogeneous Linear Equations

Science.gov (United States)

Scherer, Philipp O. J.

Many problems in physics and especially computational physics involve systems of linear equations which arise e.g. from linearization of a general nonlinear problem or from discretization of differential equations. If the dimension of the system is not too large standard methods like Gaussian elimination or QR decomposition are sufficient. Systems with a tridiagonal matrix are important for cubic spline interpolation and numerical second derivatives. They can be solved very efficiently with a specialized Gaussian elimination method. Practical applications often involve very large dimensions and require iterative methods. Convergence of Jacobi and Gauss-Seidel methods is slow and can be improved by relaxation or over-relaxation. An alternative for large systems is the method of conjugate gradients.

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.

Numerical solution of the radionuclide transport equation

International Nuclear Information System (INIS)

Hadermann, J.; Roesel, F.

1983-11-01

A numerical solution of the one-dimensional geospheric radionuclide chain transport equation based on the pseudospectral method is developed. The advantages of this approach are flexibility in incorporating space and time dependent migration parameters, arbitrary boundary conditions and solute rock interactions as well as efficiency and reliability. As an application the authors investigate the impact of non-linear sorption isotherms on migration in crystalline rock. It is shown that non-linear sorption, in the present case a Freundlich isotherm, may reduce concentration at the geosphere outlet by orders of magnitude provided the migration time is comparable or larger than the half-life of the nuclide in question. The importance of fixing dispersivity within the continuum approach is stressed. (Auth.)

Diffusion phenomenon for linear dissipative wave equations

KAUST Repository

Said-Houari, Belkacem

2012-01-01

In this paper we prove the diffusion phenomenon for the linear wave equation. To derive the diffusion phenomenon, a new method is used. In fact, for initial data in some weighted spaces, we prove that for {equation presented} decays with the rate {equation presented} [0,1] faster than that of either u or v, where u is the solution of the linear wave equation with initial data {equation presented} [0,1], and v is the solution of the related heat equation with initial data v 0 = u 0 + u 1. This result improves the result in H. Yang and A. Milani [Bull. Sci. Math. 124 (2000), 415-433] in the sense that, under the above restriction on the initial data, the decay rate given in that paper can be improved by t -γ/2. © European Mathematical Society.

A new linearized equation for servo valve in hydraulic control systems

International Nuclear Information System (INIS)

Kim, Tae Hyung; Lee, Ill Yeong

2002-01-01

In the procedure of the hydraulic control system analysis, a linearized approximate equation described by the first order term of Taylor's series has been widely used. Such a linearized equation is effective just near the operating point. And, as of now, there are no general standards on how to determine the operating point of a servo valve in the process of applying the linearized equation. So, in this study, a new linearized equation for valve characteristics is proposed as a modified form of the existing linearized equation. And, a method for selecting an optimal operating point is proposed for the new linearized equation. The effectiveness of the new linearized equation is confirmed through numerical simulations and experiments for a model hydraulic control system

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.

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

OpenAIRE

Nakamura, Gen; Vashisth, Manmohan

2017-01-01

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

The numerical solution of linear multi-term fractional differential equations: systems of equations

Science.gov (United States)

Edwards, John T.; Ford, Neville J.; Simpson, A. Charles

2002-11-01

In this paper, we show how the numerical approximation of the solution of a linear multi-term fractional differential equation can be calculated by reduction of the problem to a system of ordinary and fractional differential equations each of order at most unity. We begin by showing how our method applies to a simple class of problems and we give a convergence result. We solve the Bagley Torvik equation as an example. We show how the method can be applied to a general linear multi-term equation and give two further examples.

Infinite sets of conservation laws for linear and non-linear field equations

International Nuclear Information System (INIS)

Niederle, J.

1984-01-01

The work was motivated by a desire to understand group theoretically the existence of an infinite set of conservation laws for non-interacting fields and to carry over these conservation laws to the case of interacting fields. The relation between an infinite set of conservation laws of a linear field equation and the enveloping algebra of its space-time symmetry group was established. It is shown that in the case of the Korteweg-de Vries (KdV) equation to each symmetry of the corresponding linear equation delta sub(o)uxxx=u sub() determined by an element of the enveloping algebra of the space translation algebra, there corresponds a symmetry of the full KdV equation

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

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

On the relativistic transport equation for a discontinuity wave of multiplicity one

International Nuclear Information System (INIS)

Giambo, Sebastiano; Palumbo, Annunziata

1980-01-01

In the framework of the theory of the singular hypersurfaces, the transport equation for the amplitude of a discontinuity wave, corresponding to a simple characteristic of a quasi-linear hyperbolic system, is established in the context of special relativity [fr

Invariant imbedding equations for linear scattering problems

International Nuclear Information System (INIS)

Apresyan, L.

1988-01-01

A general form of the invariant imbedding equations is investigated for the linear problem of scattering by a bounded scattering volume. The conditions for the derivability of such equations are described. It is noted that the possibility of the explicit representation of these equations for a sphere and for a layer involves the separation of variables in the unperturbed wave equation

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

Linear Einstein equations and Kerr-Schild maps

International Nuclear Information System (INIS)

Gergely, Laszlo A

2002-01-01

We prove that given a solution of the Einstein equations g ab for the matter field T ab , an autoparallel null vector field l a and a solution (l a l c , T ac ) of the linearized Einstein equation on the given background, the Kerr-Schild metric g ac + λl a l c (λ arbitrary constant) is an exact solution of the Einstein equation for the energy-momentum tensor T ac + λT ac + λ 2 l (a T c)b l b . The mixed form of the Einstein equation for Kerr-Schild metrics with autoparallel null congruence is also linear. Some more technical conditions hold when the null congruence is not autoparallel. These results generalize previous theorems for vacuum due to Xanthopoulos and for flat seed spacetime due to Guerses and Guersey

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)

Optimal traffic control in highway transportation networks using linear programming

KAUST Repository

Li, Yanning

2014-06-01

This article presents a framework for the optimal control of boundary flows on transportation networks. The state of the system is modeled by a first order scalar conservation law (Lighthill-Whitham-Richards PDE). Based on an equivalent formulation of the Hamilton-Jacobi PDE, the problem of controlling the state of the system on a network link in a finite horizon can be posed as a Linear Program. Assuming all intersections in the network are controllable, we show that the optimization approach can be extended to an arbitrary transportation network, preserving linear constraints. Unlike previously investigated transportation network control schemes, this framework leverages the intrinsic properties of the Halmilton-Jacobi equation, and does not require any discretization or boolean variables on the link. Hence this framework is very computational efficient and provides the globally optimal solution. The feasibility of this framework is illustrated by an on-ramp metering control example.

Non-linear wave equations:Mathematical techniques

International Nuclear Information System (INIS)

1978-01-01

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

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

International Nuclear Information System (INIS)

Chen, G.S.

1997-01-01

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

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)

Solution of the neutron transport equation by the Method of Characteristics using a linear representation of the source within a mesh

International Nuclear Information System (INIS)

Mazumdar, Tanay; Degweker, S.B.

2017-01-01

Highlights: • In Method of Characteristics, the neutron source within a mesh is expanded up to linear term. • This expansion reduces the number of meshes as compared to flat source assumption. • Poor representation of circular geometry with coarser meshes is corrected. • Few benchmark problems are solved to show the advantages of linear expansion of source. • The advantage of the present formalism is quite visible in problems with large flux gradient. - Abstract: A common assumption in the solution of the neutron transport equation by the Method of Characteristics (MOC) is that the source (or flux) is constant within a mesh. This assumption is adequate provided the meshes are small enough so that the spatial variation of flux within a mesh may be ignored. Whether a mesh is small enough or not depends upon the flux gradient across a mesh, which in turn depends on factors like the presence of strong absorbers, localized sources or vacuum boundaries. The flat flux assumption often requires a very large number of meshes for solving the neutron transport equation with acceptable accuracy as was observed in our earlier work on the subject. A significant reduction in the required number of meshes is attainable by using a higher order representation of the flux within a mesh. In this paper, we expand the source within a mesh up to first order (linear) terms, which permits the use of larger sized (and therefore fewer) meshes and thereby reduces the computation time without compromising the accuracy of calculation. Since the division of the geometry into meshes is through an automatic triangulation procedure using the Bowyer-Watson algorithm, representation of circular objects (cylindrical fuel rods) with coarse meshes is poorer and causes geometry related errors. A numerical recipe is presented to make a correction to the automatic triangulation process and thereby eliminate this source of error. A number of benchmark problems are analyzed to emphasize the

Dynamical symmetries of semi-linear Schrodinger and diffusion equations

International Nuclear Information System (INIS)

Stoimenov, Stoimen; Henkel, Malte

2005-01-01

Conditional and Lie symmetries of semi-linear 1D Schrodinger and diffusion equations are studied if the mass (or the diffusion constant) is considered as an additional variable. In this way, dynamical symmetries of semi-linear Schrodinger equations become related to the parabolic and almost-parabolic subalgebras of a three-dimensional conformal Lie algebra (conf 3 ) C . We consider non-hermitian representations and also include a dimensionful coupling constant of the non-linearity. The corresponding representations of the parabolic and almost-parabolic subalgebras of (conf 3 ) C are classified and the complete list of conditionally invariant semi-linear Schrodinger equations is obtained. Possible applications to the dynamical scaling behaviour of phase-ordering kinetics are discussed

Half-trek criterion for generic identifiability of linear structural equation models

NARCIS (Netherlands)

Foygel, R.; Draisma, J.; Drton, M.

2012-01-01

A linear structural equation model relates random variables of interest and corresponding Gaussian noise terms via a linear equation system. Each such model can be represented by a mixed graph in which directed edges encode the linear equations, and bidirected edges indicate possible correlations

Half-trek criterion for generic identifiability of linear structural equation models

NARCIS (Netherlands)

Foygel, R.; Draisma, J.; Drton, M.

2011-01-01

A linear structural equation model relates random variables of interest and corresponding Gaussian noise terms via a linear equation system. Each such model can be represented by a mixed graph in which directed edges encode the linear equations, and bidirected edges indicate possible correlations

Introduction to linear systems of differential equations

CERN Document Server

Adrianova, L Ya

1995-01-01

The theory of linear systems of differential equations is one of the cornerstones of the whole theory of differential equations. At its root is the concept of the Lyapunov characteristic exponent. In this book, Adrianova presents introductory material and further detailed discussions of Lyapunov exponents. She also discusses the structure of the space of solutions of linear systems. Classes of linear systems examined are from the narrowest to widest: 1)�autonomous, 2)�periodic, 3)�reducible to autonomous, 4)�nearly reducible to autonomous, 5)�regular. In addition, Adrianova considers the following: stability of linear systems and the influence of perturbations of the coefficients on the stability the criteria of uniform stability and of uniform asymptotic stability in terms of properties of the solutions several estimates of the growth rate of solutions of a linear system in terms of its coefficients How perturbations of the coefficients change all the elements of the spectrum of the system is defin...

GLOBAL LINEARIZATION OF DIFFERENTIAL EQUATIONS WITH SPECIAL STRUCTURES

Institute of Scientific and Technical Information of China (English)

无

2011-01-01

This paper introduces the global linearization of the differential equations with special structures.The function in the differential equation is unbounded.We prove that the differential equation with unbounded function can be topologically linearlized if it has a special structure.

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

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

Geometric Insight into Scalar Combination of Linear Equations

Indian Academy of Sciences (India)

... Journals; Resonance – Journal of Science Education; Volume 14; Issue 11. Geometric Insight into Scalar Combination of Linear Equations. Ranjit Konkar. Classroom Volume 14 Issue 11 November 2009 pp 1092-1097 ... Keywords. Linear algebra; linear dependence; linear combination; family of lines; family of planes.

A parallel version of a multigrid algorithm for isotropic transport equations

International Nuclear Information System (INIS)

Manteuffel, T.; McCormick, S.; Yang, G.; Morel, J.; Oliveira, S.

1994-01-01

The focus of this paper is on a parallel algorithm for solving the transport equations in a slab geometry using multigrid. The spatial discretization scheme used is a finite element method called the modified linear discontinuous (MLD) scheme. The MLD scheme represents a lumped version of the standard linear discontinuous (LD) scheme. The parallel algorithm was implemented on the Connection Machine 2 (CM2). Convergence rates and timings for this algorithm on the CM2 and Cray-YMP are shown

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.

What happens to linear properties as we move from the Klein-Gordon equation to the sine-Gordon equation

International Nuclear Information System (INIS)

Kovalyov, Mikhail

2010-01-01

In this article the sets of solutions of the sine-Gordon equation and its linearization the Klein-Gordon equation are discussed and compared. It is shown that the set of solutions of the sine-Gordon equation possesses a richer structure which partly disappears during linearization. Just like the solutions of the Klein-Gordon equation satisfy the linear superposition principle, the solutions of the sine-Gordon equation satisfy a nonlinear superposition principle.

Emmy Noether and Linear Evolution Equations

Directory of Open Access Journals (Sweden)

P. G. L. Leach

2013-01-01

Full Text Available Noether’s Theorem relates the Action Integral of a Lagrangian with symmetries which leave it invariant and the first integrals consequent upon the variational principle and the existence of the symmetries. These each have an equivalent in the Schrödinger Equation corresponding to the Lagrangian and by extension to linear evolution equations in general. The implications of these connections are investigated.

Simplified Linear Equation Solvers users manual

Energy Technology Data Exchange (ETDEWEB)

Gropp, W. [Argonne National Lab., IL (United States); Smith, B. [California Univ., Los Angeles, CA (United States)

1993-02-01

The solution of large sparse systems of linear equations is at the heart of many algorithms in scientific computing. The SLES package is a set of easy-to-use yet powerful and extensible routines for solving large sparse linear systems. The design of the package allows new techniques to be used in existing applications without any source code changes in the applications.

A linearizing transformation for the Korteweg-de Vries equation; generalizations to higher-dimensional nonlinear partial differential equations

NARCIS (Netherlands)

Dorren, H.J.S.

1998-01-01

It is shown that the Korteweg–de Vries (KdV) equation can be transformed into an ordinary linear partial differential equation in the wave number domain. Explicit solutions of the KdV equation can be obtained by subsequently solving this linear differential equation and by applying a cascade of

Construction of a Roe linearization for the ideal MHD equations

International Nuclear Information System (INIS)

Cargo, P.; Gallice, G.; Raviart, P.A.

1996-01-01

In [3], Munz has constructed a Roe linearization for the equations of gas dynamics in Lagrangian coordinates. We extend this construction to the case of the ideal magnetohydrodynamics equations again in Lagrangian coordinates. As a consequence we obtain a Roe linearization for the MHD equations in Eulerian coordinates. (author)

A general analytical approach to the one-group, one-dimensional transport equation

International Nuclear Information System (INIS)

Barichello, L.B.; Vilhena, M.T.

1993-01-01

The main feature of the presented approach to solve the neutron transport equation consists in the application of the Laplace transform to the discrete ordinates equations, which yields a linear system of order N to be solved (LTS N method). In this paper this system is solved analytically and the inversion is performed using the Heaviside expansion technique. The general formulation achieved by this procedure is then applied to homogeneous and heterogeneous one-group slab-geometry problems. (orig.) [de

Variational linear algebraic equations method

International Nuclear Information System (INIS)

Moiseiwitsch, B.L.

1982-01-01

A modification of the linear algebraic equations method is described which ensures a variational bound on the phaseshifts for potentials having a definite sign at all points. The method is illustrated by the elastic scattering of s-wave electrons by the static field of atomic hydrogen. (author)

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.

Non-linear effects in the Boltzmann equation

International Nuclear Information System (INIS)

Barrachina, R.O.

1985-01-01

The Boltzmann equation is studied by defining an integral transformation of the energy distribution function for an isotropic and homogeneous gas. This transformation may be interpreted as a linear superposition of equilibrium states with variable temperatures. It is shown that the temporal evolution features of the distribution function are determined by the singularities of said transformation. This method is applied to Maxwell and Very Hard Particle interaction models. For the latter, the solution of the Boltzmann equation with the solution of its linearized version is compared, finding out many basic discrepancies and non-linear effects. This gives a hint to propose a new rational approximation method with a clear physical meaning. Applying this technique, the relaxation features of the BKW (Bobylev, Krook anf Wu) mode is analyzed, finding a conclusive counter-example for the Krook and Wu conjecture. The anisotropic Boltzmann equation for Maxwell models is solved as an expansion in terms of the eigenfunctions of the corresponding linearized collision operator, finding interesting transient overpopulation and underpopulation effects at thermal energies as well as a new preferential spreading effect. By analyzing the initial collision, a criterion is established to deduce the general features of the final approach to equilibrium. Finally, it is shown how to improve the convergence of the eigenfunction expansion for high energy underpopulated distribution functions. As an application of this theory, the linear cascade model for sputtering is analyzed, thus finding out that many differences experimentally observed are due to non-linear effects. (M.E.L.) [es

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)

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.

Asymptotic properties for half-linear difference equations

Czech Academy of Sciences Publication Activity Database

Cecchi, M.; Došlá, Z.; Marini, M.; Vrkoč, Ivo

2006-01-01

Roč. 131, č. 4 (2006), s. 347-363 ISSN 0862-7959 R&D Projects: GA ČR(CZ) GA201/04/0580 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear second order difference equation * nonoscillatory solutions * Riccati difference equation Subject RIV: BA - General Mathematics

Students' errors in solving linear equation word problems: Case ...

African Journals Online (AJOL)

The study examined errors students make in solving linear equation word problems with a view to expose the nature of these errors and to make suggestions for classroom teaching. A diagnostic test comprising 10 linear equation word problems, was administered to a sample (n=130) of senior high school first year Home ...

New Equating Methods and Their Relationships with Levine Observed Score Linear Equating under the Kernel Equating Framework

Science.gov (United States)

Chen, Haiwen; Holland, Paul

2010-01-01

In this paper, we develop a new curvilinear equating for the nonequivalent groups with anchor test (NEAT) design under the assumption of the classical test theory model, that we name curvilinear Levine observed score equating. In fact, by applying both the kernel equating framework and the mean preserving linear transformation of…

Complex eigenvalues for neutron transport equation with quadratically anisotropic scattering

International Nuclear Information System (INIS)

Sjoestrand, N.G.

1981-01-01

Complex eigenvalues for the monoenergetic neutron transport equation in the buckling approximation have been calculated for various combinations of linearly and quadratically anisotropic scattering. The results are discussed in terms of the time-dependent case. Tables are given of complex bucklings for real decay constants and of complex decay constants for real bucklings. The results fit nicely into the pattern of real and purely imaginary eigenvalues obtained earlier. (author)

High-order quantum algorithm for solving linear differential equations

International Nuclear Information System (INIS)

Berry, Dominic W

2014-01-01

Linear differential equations are ubiquitous in science and engineering. Quantum computers can simulate quantum systems, which are described by a restricted type of linear differential equations. Here we extend quantum simulation algorithms to general inhomogeneous sparse linear differential equations, which describe many classical physical systems. We examine the use of high-order methods (where the error over a time step is a high power of the size of the time step) to improve the efficiency. These provide scaling close to Δt 2 in the evolution time Δt. As with other algorithms of this type, the solution is encoded in amplitudes of the quantum state, and it is possible to extract global features of the solution. (paper)

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.

Computer codes for three dimensional mass transport with non-linear sorption

International Nuclear Information System (INIS)

Noy, D.J.

1985-03-01

The report describes the mathematical background and data input to finite element programs for three dimensional mass transport in a porous medium. The transport equations are developed and sorption processes are included in a general way so that non-linear equilibrium relations can be introduced. The programs are described and a guide given to the construction of the required input data sets. Concluding remarks indicate that the calculations require substantial computer resources and suggest that comprehensive preliminary analysis with lower dimensional codes would be important in the assessment of field data. (author)

Periodic feedback stabilization for linear periodic evolution equations

CERN Document Server

Wang, Gengsheng

2016-01-01

This book introduces a number of recent advances regarding periodic feedback stabilization for linear and time periodic evolution equations. First, it presents selected connections between linear quadratic optimal control theory and feedback stabilization theory for linear periodic evolution equations. Secondly, it identifies several criteria for the periodic feedback stabilization from the perspective of geometry, algebra and analyses respectively. Next, it describes several ways to design periodic feedback laws. Lastly, the book introduces readers to key methods for designing the control machines. Given its coverage and scope, it offers a helpful guide for graduate students and researchers in the areas of control theory and applied mathematics.

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)

Evaluation of an analytic linear Boltzmann transport equation solver for high-density inhomogeneities

Energy Technology Data Exchange (ETDEWEB)

Lloyd, S. A. M.; Ansbacher, W. [Department of Physics and Astronomy, University of Victoria, Victoria, British Columbia V8W 3P6 (Canada); Department of Physics and Astronomy, University of Victoria, Victoria, British Columbia V8W 3P6 (Canada) and Department of Medical Physics, British Columbia Cancer Agency-Vancouver Island Centre, Victoria, British Columbia V8R 6V5 (Canada)

2013-01-15

Purpose: Acuros external beam (Acuros XB) is a novel dose calculation algorithm implemented through the ECLIPSE treatment planning system. The algorithm finds a deterministic solution to the linear Boltzmann transport equation, the same equation commonly solved stochastically by Monte Carlo methods. This work is an evaluation of Acuros XB, by comparison with Monte Carlo, for dose calculation applications involving high-density materials. Existing non-Monte Carlo clinical dose calculation algorithms, such as the analytic anisotropic algorithm (AAA), do not accurately model dose perturbations due to increased electron scatter within high-density volumes. Methods: Acuros XB, AAA, and EGSnrc based Monte Carlo are used to calculate dose distributions from 18 MV and 6 MV photon beams delivered to a cubic water phantom containing a rectangular high density (4.0-8.0 g/cm{sup 3}) volume at its center. The algorithms are also used to recalculate a clinical prostate treatment plan involving a unilateral hip prosthesis, originally evaluated using AAA. These results are compared graphically and numerically using gamma-index analysis. Radio-chromic film measurements are presented to augment Monte Carlo and Acuros XB dose perturbation data. Results: Using a 2% and 1 mm gamma-analysis, between 91.3% and 96.8% of Acuros XB dose voxels containing greater than 50% the normalized dose were in agreement with Monte Carlo data for virtual phantoms involving 18 MV and 6 MV photons, stainless steel and titanium alloy implants and for on-axis and oblique field delivery. A similar gamma-analysis of AAA against Monte Carlo data showed between 80.8% and 87.3% agreement. Comparing Acuros XB and AAA evaluations of a clinical prostate patient plan involving a unilateral hip prosthesis, Acuros XB showed good overall agreement with Monte Carlo while AAA underestimated dose on the upstream medial surface of the prosthesis due to electron scatter from the high-density material. Film measurements

Infinite sets of conservation laws for linear and nonlinear field equations

International Nuclear Information System (INIS)

Mickelsson, J.

1984-01-01

The relation between an infinite set of conservation laws of a linear field equation and the enveloping algebra of the space-time symmetry group is established. It is shown that each symmetric element of the enveloping algebra of the space-time symmetry group of a linear field equation generates a one-parameter group of symmetries of the field equation. The cases of the Maxwell and Dirac equations are studied in detail. Then it is shown that (at least in the sense of a power series in the 'coupling constant') the conservation laws of the linear case can be deformed to conservation laws of a nonlinear field equation which is obtained from the linear one by adding a nonlinear term invariant under the group of space-time symmetries. As an example, our method is applied to the Korteweg-de Vries equation and to the massless Thirring model. (orig.)

A Comparison between Linear IRT Observed-Score Equating and Levine Observed-Score Equating under the Generalized Kernel Equating Framework

Science.gov (United States)

Chen, Haiwen

2012-01-01

In this article, linear item response theory (IRT) observed-score equating is compared under a generalized kernel equating framework with Levine observed-score equating for nonequivalent groups with anchor test design. Interestingly, these two equating methods are closely related despite being based on different methodologies. Specifically, when…

Dual exponential polynomials and linear differential equations

Science.gov (United States)

Wen, Zhi-Tao; Gundersen, Gary G.; Heittokangas, Janne

2018-01-01

We study linear differential equations with exponential polynomial coefficients, where exactly one coefficient is of order greater than all the others. The main result shows that a nontrivial exponential polynomial solution of such an equation has a certain dual relationship with the maximum order coefficient. Several examples illustrate our results and exhibit possibilities that can occur.

Finite element analysis of the neutron transport equation in spherical geometry

International Nuclear Information System (INIS)

Kim, Yong Ill; Kim, Jong Kyung; Suk, Soo Dong

1992-01-01

The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation. (Author)

Exact non-linear equations for cosmological perturbations

Energy Technology Data Exchange (ETDEWEB)

Gong, Jinn-Ouk [Asia Pacific Center for Theoretical Physics, Pohang 37673 (Korea, Republic of); Hwang, Jai-chan [Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Daegu 41566 (Korea, Republic of); Noh, Hyerim [Korea Astronomy and Space Science Institute, Daejeon 34055 (Korea, Republic of); Wu, David Chan Lon; Yoo, Jaiyul, E-mail: jinn-ouk.gong@apctp.org, E-mail: jchan@knu.ac.kr, E-mail: hr@kasi.re.kr, E-mail: clwu@physik.uzh.ch, E-mail: jyoo@physik.uzh.ch [Center for Theoretical Astrophysics and Cosmology, Institute for Computational Science, Universität Zürich, CH-8057 Zürich (Switzerland)

2017-10-01

We present a complete set of exact and fully non-linear equations describing all three types of cosmological perturbations—scalar, vector and tensor perturbations. We derive the equations in a thoroughly gauge-ready manner, so that any spatial and temporal gauge conditions can be employed. The equations are completely general without any physical restriction except that we assume a flat homogeneous and isotropic universe as a background. We also comment briefly on the application of our formulation to the non-expanding Minkowski background.

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

International Nuclear Information System (INIS)

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

2002-01-01

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

Inverse scattering solution of non-linear evolution equations in one space dimension: an introduction

International Nuclear Information System (INIS)

Alvarez-Estrada, R.F.

1979-01-01

A comprehensive review of the inverse scattering solution of certain non-linear evolution equations of physical interest in one space dimension is presented. We explain in some detail the interrelated techniques which allow to linearize exactly the following equations: (1) the Korteweg and de Vries equation; (2) the non-linear Schrodinger equation; (3) the modified Korteweg and de Vries equation; (4) the Sine-Gordon equation. We concentrate in discussing the pairs of linear operators which accomplish such an exact linearization and the solution of the associated initial value problem. The application of the method to other non-linear evolution equations is reviewed very briefly

HESS Opinions: Linking Darcy's equation to the linear reservoir

Science.gov (United States)

Savenije, Hubert H. G.

2018-03-01

In groundwater hydrology, two simple linear equations exist describing the relation between groundwater flow and the gradient driving it: Darcy's equation and the linear reservoir. Both equations are empirical and straightforward, but work at different scales: Darcy's equation at the laboratory scale and the linear reservoir at the watershed scale. Although at first sight they appear similar, it is not trivial to upscale Darcy's equation to the watershed scale without detailed knowledge of the structure or shape of the underlying aquifers. This paper shows that these two equations, combined by the water balance, are indeed identical provided there is equal resistance in space for water entering the subsurface network. This implies that groundwater systems make use of an efficient drainage network, a mostly invisible pattern that has evolved over geological timescales. This drainage network provides equally distributed resistance for water to access the system, connecting the active groundwater body to the stream, much like a leaf is organized to provide all stomata access to moisture at equal resistance. As a result, the timescale of the linear reservoir appears to be inversely proportional to Darcy's conductance, the proportionality being the product of the porosity and the resistance to entering the drainage network. The main question remaining is which physical law lies behind pattern formation in groundwater systems, evolving in a way that resistance to drainage is constant in space. But that is a fundamental question that is equally relevant for understanding the hydraulic properties of leaf veins in plants or of blood veins in animals.

ADM For Solving Linear Second-Order Fredholm Integro-Differential Equations

Science.gov (United States)

Karim, Mohd F.; Mohamad, Mahathir; Saifullah Rusiman, Mohd; Che-Him, Norziha; Roslan, Rozaini; Khalid, Kamil

2018-04-01

In this paper, we apply Adomian Decomposition Method (ADM) as numerically analyse linear second-order Fredholm Integro-differential Equations. The approximate solutions of the problems are calculated by Maple package. Some numerical examples have been considered to illustrate the ADM for solving this equation. The results are compared with the existing exact solution. Thus, the Adomian decomposition method can be the best alternative method for solving linear second-order Fredholm Integro-Differential equation. It converges to the exact solution quickly and in the same time reduces computational work for solving the equation. The result obtained by ADM shows the ability and efficiency for solving these equations.

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

International Nuclear Information System (INIS)

Bonalumi, Riccardo A.

1998-01-01

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

Linear System of Equations, Matrix Inversion, and Linear Programming Using MS Excel

Science.gov (United States)

El-Gebeily, M.; Yushau, B.

2008-01-01

In this note, we demonstrate with illustrations two different ways that MS Excel can be used to solve Linear Systems of Equation, Linear Programming Problems, and Matrix Inversion Problems. The advantage of using MS Excel is its availability and transparency (the user is responsible for most of the details of how a problem is solved). Further, we…

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

Turbulence Spreading into Linearly Stable Zone and Transport Scaling

International Nuclear Information System (INIS)

Hahm, T.S.; Diamond, P.H.; Lin, Z.; Itoh, K.; Itoh, S.-I.

2003-01-01

We study the simplest problem of turbulence spreading corresponding to the spatio-temporal propagation of a patch of turbulence from a region where it is locally excited to a region of weaker excitation, or even local damping. A single model equation for the local turbulence intensity I(x, t) includes the effects of local linear growth and damping, spatially local nonlinear coupling to dissipation and spatial scattering of turbulence energy induced by nonlinear coupling. In the absence of dissipation, the front propagation into the linearly stable zone occurs with the property of rapid progression at small t, followed by slower subdiffusive progression at late times. The turbulence radial spreading into the linearly stable zone reduces the turbulent intensity in the linearly unstable zone, and introduces an additional dependence on the rho* is always equal to rho i/a to the turbulent intensity and the transport scaling. These are in broad, semi-quantitative agreements with a number of global gyrokinetic simulation results with zonal flows and without zonal flows. The front propagation stops when the radial flux of fluctuation energy from the linearly unstable region is balanced by local dissipation in the linearly stable region

Application of the Galerkin's method to the solution of the one-dimensional integral transport equation: generalized collision probabilities taken in account the flux gradient and the linearly anisotropic scattering

International Nuclear Information System (INIS)

Sanchez, Richard.

1975-04-01

For the one-dimensional geometries, the transport equation with linearly anisotropic scattering can be reduced to a single integral equation; this is a singular-kernel FREDHOLM equation of the second kind. When applying a conventional projective method that of GALERKIN, to the solution of this equation the well-known collision probability algorithm is obtained. Piecewise polynomial expansions are used to represent the flux. In the ANILINE code, the flux is supposed to be linear in plane geometry and parabolic in both cylindrical and spherical geometries. An integral relationship was found between the one-dimensional isotropic and anisotropic kernels; this allows to reduce the new matrix elements (issuing from the anisotropic kernel) to classic collision probabilities of the isotropic scattering equation. For cylindrical and spherical geometries used an approximate representation of the current was used to avoid an additional numerical integration. Reflective boundary conditions were considered; in plane geometry the reflection is supposed specular, for the other geometries the isotropic reflection hypothesis has been adopted. Further, the ANILINE code enables to deal with an incoming isotropic current. Numerous checks were performed in monokinetic theory. Critical radii and albedos were calculated for homogeneous slabs, cylinders and spheres. For heterogeneous media, the thermal utilization factor obtained by this method was compared with the theoretical result based upon a formula by BENOIST. Finally, ANILINE was incorporated into the multigroup APOLLO code, which enabled to analyse the MINERVA experimental reactor in transport theory with 99 groups. The ANILINE method is particularly suited to the treatment of strongly anisotropic media with considerable flux gradients. It is also well adapted to the calculation of reflectors, and in general, to the exact analysis of anisotropic effects in large-sized media [fr

Linear matrix differential equations of higher-order and applications

Directory of Open Access Journals (Sweden)

Mustapha Rachidi

2008-07-01

Full Text Available In this article, we study linear differential equations of higher-order whose coefficients are square matrices. The combinatorial method for computing the matrix powers and exponential is adopted. New formulas representing auxiliary results are obtained. This allows us to prove properties of a large class of linear matrix differential equations of higher-order, in particular results of Apostol and Kolodner are recovered. Also illustrative examples and applications are presented.

Equations of motion for a (non-linear) scalar field model as derived from the field equations

International Nuclear Information System (INIS)

Kaniel, S.; Itin, Y.

2006-01-01

The problem of derivation of the equations of motion from the field equations is considered. Einstein's field equations have a specific analytical form: They are linear in the second order derivatives and quadratic in the first order derivatives of the field variables. We utilize this particular form and propose a novel algorithm for the derivation of the equations of motion from the field equations. It is based on the condition of the balance between the singular terms of the field equation. We apply the algorithm to a non-linear Lorentz invariant scalar field model. We show that it results in the Newton law of attraction between the singularities of the field moved on approximately geodesic curves. The algorithm is applicable to the N-body problem of the Lorentz invariant field equations. (Abstract Copyright [2006], Wiley Periodicals, Inc.)

On linear transport problems

International Nuclear Information System (INIS)

Ignatovich, V.K.

1989-01-01

The equations. governing the transport of radiation in plane media of finite thickness are formulated and solved in terms reflection and extintion of radiation inthe case of semi infinite media. 13 refs

Large-time asymptotic behaviour of solutions of non-linear Sobolev-type equations

International Nuclear Information System (INIS)

Kaikina, Elena I; Naumkin, Pavel I; Shishmarev, Il'ya A

2009-01-01

The large-time asymptotic behaviour of solutions of the Cauchy problem is investigated for a non-linear Sobolev-type equation with dissipation. For small initial data the approach taken is based on a detailed analysis of the Green's function of the linear problem and the use of the contraction mapping method. The case of large initial data is also closely considered. In the supercritical case the asymptotic formulae are quasi-linear. The asymptotic behaviour of solutions of a non-linear Sobolev-type equation with a critical non-linearity of the non-convective kind differs by a logarithmic correction term from the behaviour of solutions of the corresponding linear equation. For a critical convective non-linearity, as well as for a subcritical non-convective non-linearity it is proved that the leading term of the asymptotic expression for large times is a self-similar solution. For Sobolev equations with convective non-linearity the asymptotic behaviour of solutions in the subcritical case is the product of a rarefaction wave and a shock wave. Bibliography: 84 titles.

A General Linear Method for Equating with Small Samples

Science.gov (United States)

Albano, Anthony D.

2015-01-01

Research on equating with small samples has shown that methods with stronger assumptions and fewer statistical estimates can lead to decreased error in the estimated equating function. This article introduces a new approach to linear observed-score equating, one which provides flexible control over how form difficulty is assumed versus estimated…

A local-global problem for linear differential equations

NARCIS (Netherlands)

Put, Marius van der; Reversat, Marc

An inhomogeneous linear differential equation Ly = f over a global differential field can have a formal solution for each place without having a global solution. The vector space lgl(L) measures this phenomenon. This space is interpreted in terms of cohomology of linear algebraic groups and is

A local-global problem for linear differential equations

NARCIS (Netherlands)

Put, Marius van der; Reversat, Marc

2008-01-01

An inhomogeneous linear differential equation Ly = f over a global differential field can have a formal solution for each place without having a global solution. The vector space lgl(L) measures this phenomenon. This space is interpreted in terms of cohomology of linear algebraic groups and is

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

Iterative solution of linear equations in ODE codes. [Krylov subspaces

Energy Technology Data Exchange (ETDEWEB)

Gear, C. W.; Saad, Y.

1981-01-01

Each integration step of a stiff equation involves the solution of a nonlinear equation, usually by a quasi-Newton method that leads to a set of linear problems. Iterative methods for these linear equations are studied. Of particular interest are methods that do not require an explicit Jacobian, but can work directly with differences of function values using J congruent to f(x + delta) - f(x). Some numerical experiments using a modification of LSODE are reported. 1 figure, 2 tables.

Runge-Kutta Methods for Linear Ordinary Differential Equations

Science.gov (United States)

Zingg, David W.; Chisholm, Todd T.

1997-01-01

Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.

Projective-Dual Method for Solving Systems of Linear Equations with Nonnegative Variables

Science.gov (United States)

Ganin, B. V.; Golikov, A. I.; Evtushenko, Yu. G.

2018-02-01

In order to solve an underdetermined system of linear equations with nonnegative variables, the projection of a given point onto its solutions set is sought. The dual of this problem—the problem of unconstrained maximization of a piecewise-quadratic function—is solved by Newton's method. The problem of unconstrained optimization dual of the regularized problem of finding the projection onto the solution set of the system is considered. A connection of duality theory and Newton's method with some known algorithms of projecting onto a standard simplex is shown. On the example of taking into account the specifics of the constraints of the transport linear programming problem, the possibility to increase the efficiency of calculating the generalized Hessian matrix is demonstrated. Some examples of numerical calculations using MATLAB are presented.

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

Dissipative behavior of some fully non-linear KdV-type equations

Science.gov (United States)

Brenier, Yann; Levy, Doron

2000-03-01

The KdV equation can be considered as a special case of the general equation u t+f(u) x-δg(u xx) x=0, δ>0, where f is non-linear and g is linear, namely f( u)= u2/2 and g( v)= v. As the parameter δ tends to 0, the dispersive behavior of the KdV equation has been throughly investigated (see, e.g., [P.G. Drazin, Solitons, London Math. Soc. Lect. Note Ser. 85, Cambridge University Press, Cambridge, 1983; P.D. Lax, C.D. Levermore, The small dispersion limit of the Korteweg-de Vries equation, III, Commun. Pure Appl. Math. 36 (1983) 809-829; G.B. Whitham, Linear and Nonlinear Waves, Wiley/Interscience, New York, 1974] and the references therein). We show through numerical evidence that a completely different, dissipative behavior occurs when g is non-linear, namely when g is an even concave function such as g( v)=-∣ v∣ or g( v)=- v2. In particular, our numerical results hint that as δ→0 the solutions strongly converge to the unique entropy solution of the formal limit equation, in total contrast with the solutions of the KdV equation.

Darboux transformations and linear parabolic partial differential equations

International Nuclear Information System (INIS)

Arrigo, Daniel J.; Hickling, Fred

2002-01-01

Solutions for a class of linear parabolic partial differential equation are provided. These solutions are obtained by first solving a system of (n+1) nonlinear partial differential equations. This system arises as the coefficients of a Darboux transformation and is equivalent to a matrix Burgers' equation. This matrix equation is solved using a generalized Hopf-Cole transformation. The solutions for the original equation are given in terms of solutions of the heat equation. These results are applied to the (1+1)-dimensional Schroedinger equation where all bound state solutions are obtained for a 2n-parameter family of potentials. As a special case, the solutions for integral members of the regular and modified Poeschl-Teller potentials are recovered. (author). Letter-to-the-editor

Piecewise-linear and bilinear approaches to nonlinear differential equations approximation problem of computational structural mechanics

OpenAIRE

Leibov Roman

2017-01-01

This paper presents a bilinear approach to nonlinear differential equations system approximation problem. Sometimes the nonlinear differential equations right-hand sides linearization is extremely difficult or even impossible. Then piecewise-linear approximation of nonlinear differential equations can be used. The bilinear differential equations allow to improve piecewise-linear differential equations behavior and reduce errors on the border of different linear differential equations systems ...

Solvable linear potentials in the Dirac equation

International Nuclear Information System (INIS)

Dominguez-Adame, F.; Gonzalez, M.A.

1990-01-01

The Dirac equation for some linear potentials leading to Schroedinger-like oscillator equations for the upper and lower components of the Dirac spinor have been solved. Energy levels for the bound states appear in pairs, so that both particles and antiparticles may be bound with the same energy. For weak coupling, the spacing between levels is proportional to the coupling constant while in the strong limit those levels are depressed compared to the nonrelativistic ones

Linear stochastic neutron transport theory

International Nuclear Information System (INIS)

Lewins, J.

1978-01-01

A new and direct derivation of the Bell-Pal fundamental equation for (low power) neutron stochastic behaviour in the Boltzmann continuum model is given. The development includes correlation of particle emission direction in induced and spontaneous fission. This leads to generalizations of the backward and forward equations for the mean and variance of neutron behaviour. The stochastic importance for neutron transport theory is introduced and related to the conventional deterministic importance. Defining equations and moment equations are derived and shown to be related to the backward fundamental equation with the detector distribution of the operational definition of stochastic importance playing the role of an adjoint source. (author)

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

Technological pedagogical content knowledge of junior high school mathematics teachers in teaching linear equation

Science.gov (United States)

Wati, S.; Fitriana, L.; Mardiyana

2018-04-01

Linear equation is one of the topics in mathematics that are considered difficult. Student difficulties of understanding linear equation can be caused by lack of understanding this concept and the way of teachers teach. TPACK is a way to understand the complex relationships between teaching and content taught through the use of specific teaching approaches and supported by the right technology tools. This study aims to identify TPACK of junior high school mathematics teachers in teaching linear equation. The method used in the study was descriptive. In the first phase, a survey using a questionnaire was carried out on 45 junior high school mathematics teachers in teaching linear equation. While in the second phase, the interview involved three teachers. The analysis of data used were quantitative and qualitative technique. The result PCK revealed teachers emphasized developing procedural and conceptual knowledge through reliance on traditional in teaching linear equation. The result of TPK revealed teachers’ lower capacity to deal with the general information and communications technologies goals across the curriculum in teaching linear equation. The result indicated that PowerPoint constitutes TCK modal technological capability in teaching linear equation. The result of TPACK seems to suggest a low standard in teachers’ technological skills across a variety of mathematics education goals in teaching linear equation. This means that the ability of teachers’ TPACK in teaching linear equation still needs to be improved.

On some perturbation techniques for quasi-linear parabolic equations

Directory of Open Access Journals (Sweden)

Igor Malyshev

1990-01-01

Full Text Available We study a nonhomogeneous quasi-linear parabolic equation and introduce a method that allows us to find the solution of a nonlinear boundary value problem in “explicit” form. This task is accomplished by perturbing the original equation with a source function, which is then found as a solution of some nonlinear operator equation.

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)

Linearized gyro-kinetic equation

International Nuclear Information System (INIS)

Catto, P.J.; Tsang, K.T.

1976-01-01

An ordering of the linearized Fokker-Planck equation is performed in which gyroradius corrections are retained to lowest order and the radial dependence appropriate for sheared magnetic fields is treated without resorting to a WKB technique. This description is shown to be necessary to obtain the proper radial dependence when the product of the poloidal wavenumber and the gyroradius is large (k rho much greater than 1). A like particle collision operator valid for arbitrary k rho also has been derived. In addition, neoclassical, drift, finite β (plasma pressure/magnetic pressure), and unperturbed toroidal electric field modifications are treated

Experimental quantum computing to solve systems of linear equations.

Science.gov (United States)

Cai, X-D; Weedbrook, C; Su, Z-E; Chen, M-C; Gu, Mile; Zhu, M-J; Li, Li; Liu, Nai-Le; Lu, Chao-Yang; Pan, Jian-Wei

2013-06-07

Solving linear systems of equations is ubiquitous in all areas of science and engineering. With rapidly growing data sets, such a task can be intractable for classical computers, as the best known classical algorithms require a time proportional to the number of variables N. A recently proposed quantum algorithm shows that quantum computers could solve linear systems in a time scale of order log(N), giving an exponential speedup over classical computers. Here we realize the simplest instance of this algorithm, solving 2×2 linear equations for various input vectors on a quantum computer. We use four quantum bits and four controlled logic gates to implement every subroutine required, demonstrating the working principle of this algorithm.

Oscillation and non-oscillation criterion for Riemann–Weber type half-linear differential equations

Directory of Open Access Journals (Sweden)

Petr Hasil

2016-08-01

Full Text Available By the combination of the modified half-linear Prüfer method and the Riccati technique, we study oscillatory properties of half-linear differential equations. Taking into account the transformation theory of half-linear equations and using some known results, we show that the analysed equations in the Riemann–Weber form with perturbations in both terms are conditionally oscillatory. Within the process, we identify the critical oscillation values of their coefficients and, consequently, we decide when the considered equations are oscillatory and when they are non-oscillatory. As a direct corollary of our main result, we solve the so-called critical case for a certain type of half-linear non-perturbed equations.

Solution methods for large systems of linear equations in BACCHUS

International Nuclear Information System (INIS)

Homann, C.; Dorr, B.

1993-05-01

The computer programme BACCHUS is used to describe steady state and transient thermal-hydraulic behaviour of a coolant in a fuel element with intact geometry in a fast breeder reactor. In such computer programmes generally large systems of linear equations with sparse matrices of coefficients, resulting from discretization of coolant conservation equations, must be solved thousands of times giving rise to large demands of main storage and CPU time. Direct and iterative solution methods of the systems of linear equations, available in BACCHUS, are described, giving theoretical details and experience with their use in the programme. Besides use of a method of lines, a Runge-Kutta-method, for solution of the partial differential equation is outlined. (orig.) [de

Perturbations of linear delay differential equations at the verge of instability.

Science.gov (United States)

Lingala, N; Namachchivaya, N Sri

2016-06-01

The characteristic equation for a linear delay differential equation (DDE) has countably infinite roots on the complex plane. This paper considers linear DDEs that are on the verge of instability, i.e., a pair of roots of the characteristic equation lies on the imaginary axis of the complex plane and all other roots have negative real parts. It is shown that when small noise perturbations are present, the probability distribution of the dynamics can be approximated by the probability distribution of a certain one-dimensional stochastic differential equation (SDE) without delay. This is advantageous because equations without delay are easier to simulate and one-dimensional SDEs are analytically tractable. When the perturbations are also linear, it is shown that the stability depends on a specific complex number. The theory is applied to study oscillators with delayed feedback. Some errors in other articles that use multiscale approach are pointed out.

Local energy decay for linear wave equations with variable coefficients

Science.gov (United States)

Ikehata, Ryo

2005-06-01

A uniform local energy decay result is derived to the linear wave equation with spatial variable coefficients. We deal with this equation in an exterior domain with a star-shaped complement. Our advantage is that we do not assume any compactness of the support on the initial data, and its proof is quite simple. This generalizes a previous famous result due to Morawetz [The decay of solutions of the exterior initial-boundary value problem for the wave equation, Comm. Pure Appl. Math. 14 (1961) 561-568]. In order to prove local energy decay, we mainly apply two types of ideas due to Ikehata-Matsuyama [L2-behaviour of solutions to the linear heat and wave equations in exterior domains, Sci. Math. Japon. 55 (2002) 33-42] and Todorova-Yordanov [Critical exponent for a nonlinear wave equation with damping, J. Differential Equations 174 (2001) 464-489].

Constructive Development of the Solutions of Linear Equations in Introductory Ordinary Differential Equations

Science.gov (United States)

Mallet, D. G.; McCue, S. W.

2009-01-01

The solution of linear ordinary differential equations (ODEs) is commonly taught in first-year undergraduate mathematics classrooms, but the understanding of the concept of a solution is not always grasped by students until much later. Recognizing what it is to be a solution of a linear ODE and how to postulate such solutions, without resorting to…

Linear orbit parameters for the exact equations of motion

International Nuclear Information System (INIS)

Parzen, G.

1995-01-01

This paper defines the beta function and other linear orbit parameters using the exact equations of motion. The β, α and ψ functions are redefined using the exact equations. Expressions are found for the transfer matrix and the emittance. The differential equations for η = x/β 1/2 is found. New relationships between α, β, ψ and ν are derived

Approximate Method for Solving the Linear Fuzzy Delay Differential Equations

Directory of Open Access Journals (Sweden)

S. Narayanamoorthy

2015-01-01

Full Text Available We propose an algorithm of the approximate method to solve linear fuzzy delay differential equations using Adomian decomposition method. The detailed algorithm of the approach is provided. The approximate solution is compared with the exact solution to confirm the validity and efficiency of the method to handle linear fuzzy delay differential equation. To show this proper features of this proposed method, numerical example is illustrated.

A canonical form of the equation of motion of linear dynamical systems

Science.gov (United States)

Kawano, Daniel T.; Salsa, Rubens Goncalves; Ma, Fai; Morzfeld, Matthias

2018-03-01

The equation of motion of a discrete linear system has the form of a second-order ordinary differential equation with three real and square coefficient matrices. It is shown that, for almost all linear systems, such an equation can always be converted by an invertible transformation into a canonical form specified by two diagonal coefficient matrices associated with the generalized acceleration and displacement. This canonical form of the equation of motion is unique up to an equivalence class for non-defective systems. As an important by-product, a damped linear system that possesses three symmetric and positive definite coefficients can always be recast as an undamped and decoupled system.

LOCFES-B: Solving the one-dimensional transport equation with user-selected spatial approximations

International Nuclear Information System (INIS)

Jarvis, R.D.; Nelson, P.

1993-01-01

Closed linear one-cell functional (CLOF) methods constitute an abstractly defined class of spatial approximations to the one-dimensional discrete ordinates equations of linear particle transport that encompass, as specific instances, the vast majority of the spatial approximations that have been either used or suggested in the computational solution of these equations. A specific instance of the class of CLOF methods is defined by a (typically small) number of functions of the cell width, total cross section, and direction cosine of particle motion. The LOCFES code takes advantage of the latter observation by permitting the use, within a more-or-less standard source iteration solution process, of an arbitrary CLOF method as defined by a user-supplied subroutine. The design objective of LOCFES was to provide automated determination of the order of accuracy (i.e., order of the discretization error) in the fine-mesh limit for an arbitrary user-selected CLOF method. This asymptotic order of accuracy is one widely used measure of the merit of a spatial approximation. This paper discusses LOCFES-B, which is a code that uses methods developed in LOCFES to solve one-dimensional linear particle transport problems with any user-selected CLOF method. LOCFES-B provides automatic solution of a given problem to within an accuracy specified by user input and provides comparison of the computational results against results from externally provided benchmark results

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

An implicit spectral formula for generalized linear Schroedinger equations

International Nuclear Information System (INIS)

Schulze-Halberg, A.; Garcia-Ravelo, J.; Pena Gil, Jose Juan

2009-01-01

We generalize the semiclassical Bohr–Sommerfeld quantization rule to an exact, implicit spectral formula for linear, generalized Schroedinger equations admitting a discrete spectrum. Special cases include the position-dependent mass Schroedinger equation or the Schroedinger equation for weighted energy. Requiring knowledge of the potential and the solution associated with the lowest spectral value, our formula predicts the complete spectrum in its exact form. (author)

Linear Equating for the NEAT Design: Parameter Substitution Models and Chained Linear Relationship Models

Science.gov (United States)

Kane, Michael T.; Mroch, Andrew A.; Suh, Youngsuk; Ripkey, Douglas R.

2009-01-01

This paper analyzes five linear equating models for the "nonequivalent groups with anchor test" (NEAT) design with internal anchors (i.e., the anchor test is part of the full test). The analysis employs a two-dimensional framework. The first dimension contrasts two general approaches to developing the equating relationship. Under a "parameter…

An Evaluation of Five Linear Equating Methods for the NEAT Design

Science.gov (United States)

Mroch, Andrew A.; Suh, Youngsuk; Kane, Michael T.; Ripkey, Douglas R.

2009-01-01

This study uses the results of two previous papers (Kane, Mroch, Suh, & Ripkey, this issue; Suh, Mroch, Kane, & Ripkey, this issue) and the literature on linear equating to evaluate five linear equating methods along several dimensions, including the plausibility of their assumptions and their levels of bias and root mean squared difference…

Exact solution of the neutron transport equation in spherical geometry

Energy Technology Data Exchange (ETDEWEB)

Anli, Fikret; Akkurt, Abdullah; Yildirim, Hueseyin; Ates, Kemal [Kahramanmaras Suetcue Imam Univ. (Turkey). Faculty of Sciences and Letters

2017-03-15

Solution of the neutron transport equation in one dimensional slab geometry construct a basis for the solution of neutron transport equation in a curvilinear geometry. Therefore, in this work, we attempt to derive an exact analytical benchmark solution for both neutron transport equations in slab and spherical medium by using P{sub N} approximation which is widely used in neutron transport theory.

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

SUPPORTING STUDENTS’ UNDERSTANDING OF LINEAR EQUATIONS WITH ONE VARIABLE USING ALGEBRA TILES

Directory of Open Access Journals (Sweden)

Sari Saraswati

2016-01-01

Full Text Available This research aimed to describe how algebra tiles can support students’ understanding of linear equations with one variable. This article is a part of a larger research on learning design of linear equations with one variable using algebra tiles combined with balancing method. Therefore, it will merely discuss one activity focused on how students use the algebra tiles to find a method to solve linear equations with one variable. Design research was used as an approach in this study. It consists of three phases, namely preliminary design, teaching experiment and retrospective analysis. Video registrations, students’ written works, pre-test, post-test, field notes, and interview are technic to collect data. The data were analyzed by comparing the hypothetical learning trajectory (HLT and the actual learning process. The result shows that algebra tiles could supports students’ understanding to find the formal solution of linear equation with one variable.Keywords: linear equation with one variable, algebra tiles, design research, balancing method, HLT DOI: http://dx.doi.org/10.22342/jme.7.1.2814.19-30

New approach to solve fully fuzzy system of linear equations using ...

Indian Academy of Sciences (India)

This paper proposes two new methods to solve fully fuzzy system of linear equations. The fuzzy system has been converted to a crisp system of linear equations by using single and double parametric form of fuzzy numbers to obtain the non-negative solution. Double parametric form of fuzzy numbers is defined and applied ...

Growth of meromorphic solutions of higher-order linear differential equations

Directory of Open Access Journals (Sweden)

Wenjuan Chen

2009-01-01

Full Text Available In this paper, we investigate the higher-order linear differential equations with meromorphic coefficients. We improve and extend a result of M.S. Liu and C.L. Yuan, by using the estimates for the logarithmic derivative of a transcendental meromorphic function due to Gundersen, and the extended Winman-Valiron theory which proved by J. Wang and H.X. Yi. In addition, we also consider the nonhomogeneous linear differential equations.

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

Transport methods: general. 3. An Additive Angular-Dependent Re-balance Acceleration Method for Neutron Transport Equations

International Nuclear Information System (INIS)

Cho, Nam Zin; Park, Chang Je

2001-01-01

An additive angular-dependent re-balance (AADR) factor acceleration method is described to accelerate the source iteration of discrete ordinates transport calculation. The formulation of the AADR method follows that of the angular-dependent re-balance (ADR) method in that the re-balance factor is defined only on the cell interface and in that the low-order equation is derived by integrating the transport equation (high-order equation) over angular subspaces. But, the re-balance factor is applied additively. While the AADR method is similar to the boundary projection acceleration and the alpha-weighted linear acceleration, it is more general and does have distinct features. The method is easily extendible to DP N and low-order S N re-balancing, and it does not require consistent discretizations between the high- and low-order equations as in diffusion synthetic acceleration. We find by Fourier analysis and numerical results that the AADR method with a chosen form of weighting functions is unconditionally stable and very effective. There also exists an optimal weighting parameter that leads to the smallest spectral radius. The AADR acceleration method described in this paper is simple to implement, unconditionally stable, and very effective. It uses a physically based weighting function with an optimal parameter, leading to the best spectral radius of ρ<0.1865, compared to ρ<0.2247 of DSA. The application of the AADR acceleration method with the LMB scheme on a test problem shows encouraging results

SUPPORTING STUDENTS’ UNDERSTANDING OF LINEAR EQUATIONS WITH ONE VARIABLE USING ALGEBRA TILES

Directory of Open Access Journals (Sweden)

Sari Saraswati

2016-01-01

Full Text Available This research aimed to describe how algebra tiles can support students’ understanding of linear equations with one variable. This article is a part of a larger research on learning design of linear equations with one variable using algebra tiles combined with balancing method. Therefore, it will merely discuss one activity focused on how students use the algebra tiles to find a method to solve linear equations with one variable. Design research was used as an approach in this study. It consists of three phases, namely preliminary design, teaching experiment and retrospective analysis. Video registrations, students’ written works, pre-test, post-test, field notes, and interview are technic to collect data. The data were analyzed by comparing the hypothetical learning trajectory (HLT and the actual learning process. The result shows that algebra tiles could supports students’ understanding to find the formal solution of linear equation with one variable.

Optimal Control of Scalar Conservation Laws Using Linear/Quadratic Programming: Application to Transportation Networks

KAUST Repository

Li, Yanning

2014-03-01

This article presents a new optimal control framework for transportation networks in which the state is modeled by a first order scalar conservation law. Using an equivalent formulation based on a Hamilton-Jacobi (H-J) equation and the commonly used triangular fundamental diagram, we pose the problem of controlling the state of the system on a network link, in a finite horizon, as a Linear Program (LP). We then show that this framework can be extended to an arbitrary transportation network, resulting in an LP or a Quadratic Program. Unlike many previously investigated transportation network control schemes, this method yields a globally optimal solution and is capable of handling shocks (i.e., discontinuities in the state of the system). As it leverages the intrinsic properties of the H-J equation used to model the state of the system, it does not require any approximation, unlike classical methods that are based on discretizations of the model. The computational efficiency of the method is illustrated on a transportation network. © 2014 IEEE.

Optimal Control of Scalar Conservation Laws Using Linear/Quadratic Programming: Application to Transportation Networks

KAUST Repository

Li, Yanning; Canepa, Edward S.; Claudel, Christian

2014-01-01

This article presents a new optimal control framework for transportation networks in which the state is modeled by a first order scalar conservation law. Using an equivalent formulation based on a Hamilton-Jacobi (H-J) equation and the commonly used triangular fundamental diagram, we pose the problem of controlling the state of the system on a network link, in a finite horizon, as a Linear Program (LP). We then show that this framework can be extended to an arbitrary transportation network, resulting in an LP or a Quadratic Program. Unlike many previously investigated transportation network control schemes, this method yields a globally optimal solution and is capable of handling shocks (i.e., discontinuities in the state of the system). As it leverages the intrinsic properties of the H-J equation used to model the state of the system, it does not require any approximation, unlike classical methods that are based on discretizations of the model. The computational efficiency of the method is illustrated on a transportation network. © 2014 IEEE.

New non-linear modified massless Klein-Gordon equation

Energy Technology Data Exchange (ETDEWEB)

Asenjo, Felipe A. [Universidad Adolfo Ibanez, UAI Physics Center, Santiago (Chile); Universidad Adolfo Ibanez, Facultad de Ingenieria y Ciencias, Santiago (Chile); Hojman, Sergio A. [Universidad Adolfo Ibanez, UAI Physics Center, Santiago (Chile); Universidad Adolfo Ibanez, Departamento de Ciencias, Facultad de Artes Liberales, Santiago (Chile); Universidad de Chile, Departamento de Fisica, Facultad de Ciencias, Santiago (Chile); Centro de Recursos Educativos Avanzados, CREA, Santiago (Chile)

2017-11-15

The massless Klein-Gordon equation on arbitrary curved backgrounds allows for solutions which develop ''tails'' inside the light cone and, therefore, do not strictly follow null geodesics as discovered by DeWitt and Brehme almost 60 years ago. A modification of the massless Klein-Gordon equation is presented, which always exhibits null geodesic propagation of waves on arbitrary curved spacetimes. This new equation is derived from a Lagrangian which exhibits current-current interaction. Its non-linearity is due to a self-coupling term which is related to the quantum mechanical Bohm potential. (orig.)

Non-linear absorption for concentrated solar energy transport

Energy Technology Data Exchange (ETDEWEB)

Jaramillo, O. A; Del Rio, J.A; Huelsz, G [Centro de Investigacion de Energia, UNAM, Temixco, Morelos (Mexico)

2000-07-01

In order to determine the maximum solar energy that can be transported using SiO{sub 2} optical fibers, analysis of non-linear absorption is required. In this work, we model the interaction between solar radiation and the SiO{sub 2} optical fiber core to determine the dependence of the absorption of the radioactive intensity. Using Maxwell's equations we obtain the relation between the refractive index and the electric susceptibility up to second order in terms of the electric field intensity. This is not enough to obtain an explicit expression for the non-linear absorption. Thus, to obtain the non-linear optical response, we develop a microscopic model of an harmonic driven oscillators with damp ing, based on the Drude-Lorentz theory. We solve this model using experimental information for the SiO{sub 2} optical fiber, and we determine the frequency-dependence of the non-linear absorption and the non-linear extinction of SiO{sub 2} optical fibers. Our results estimate that the average value over the solar spectrum for the non-linear extinction coefficient for SiO{sub 2} is k{sub 2}=10{sup -}29m{sup 2}V{sup -}2. With this result we conclude that the non-linear part of the absorption coefficient of SiO{sub 2} optical fibers during the transport of concentrated solar energy achieved by a circular concentrator is negligible, and therefore the use of optical fibers for solar applications is an actual option. [Spanish] Con el objeto de determinar la maxima energia solar que puede transportarse usando fibras opticas de SiO{sub 2} se requiere el analisis de absorcion no linear. En este trabajo modelamos la interaccion entre la radiacion solar y el nucleo de la fibra optica de SiO{sub 2} para determinar la dependencia de la absorcion de la intensidad radioactiva. Mediante el uso de las ecuaciones de Maxwell obtenemos la relacion entre el indice de refraccion y la susceptibilidad electrica hasta el segundo orden en terminos de intensidad del campo electrico. Esto no es

A discrete homotopy perturbation method for non-linear Schrodinger equation

Directory of Open Access Journals (Sweden)

H. A. Wahab

2015-12-01

Full Text Available A general analysis is made by homotopy perturbation method while taking the advantages of the initial guess, appearance of the embedding parameter, different choices of the linear operator to the approximated solution to the non-linear Schrodinger equation. We are not dependent upon the Adomian polynomials and find the linear forms of the components without these calculations. The discretised forms of the nonlinear Schrodinger equation allow us whether to apply any numerical technique on the discritisation forms or proceed for perturbation solution of the problem. The discretised forms obtained by constructed homotopy provide the linear parts of the components of the solution series and hence a new discretised form is obtained. The general discretised form for the NLSE allows us to choose any initial guess and the solution in the closed form.

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

On index-2 linear implicit difference equations

NARCIS (Netherlands)

Nguyen Huu Du, [No Value; Le Cong Loi, [No Value; Trinh Khanh Duy, [No Value; Vu Tien Viet, [No Value

2011-01-01

This paper deals with an index-2 notion for linear implicit difference equations (LIDEs) and with the solvability of initial value problems (IVPs) for index-2 LIDEs. Besides, the cocycle property as well as the multiplicative ergodic theorem of Oseledets type are also proved. (C) 2010 Elsevier Inc.

Singular Linear Differential Equations in Two Variables

NARCIS (Netherlands)

Braaksma, B.L.J.; Put, M. van der

2008-01-01

The formal and analytic classification of integrable singular linear differential equations has been studied among others by R. Gerard and Y. Sibuya. We provide a simple proof of their main result, namely: For certain irregular systems in two variables there is no Stokes phenomenon, i.e. there is no

Visual construction of characteristic equations of linear electric circuits

Directory of Open Access Journals (Sweden)

V.V. Kostyukov

2013-12-01

Full Text Available A visual identification method with application of partial circuits is developed for characteristic equation coefficients of transients in linear electric circuits. The method is based on interrelationship between the roots of algebraic polynomial and its coefficients. The method is illustrated with an example of a third-order linear electric circuit.

Dynamic modeling of interfacial structures via interfacial area transport equation

International Nuclear Information System (INIS)

Seungjin, Kim; Mamoru, Ishii

2005-01-01

The interfacial area transport equation dynamically models the two-phase flow regime transitions and predicts continuous change of the interfacial area concentration along the flow field. 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 phenomena of various sizes of bubbles, the two-group interfacial area transport equations have been developed. The group 1 equation describes the transport of small-dispersed bubbles that are either distorted or spherical in shapes, and the group 2 equation describes the transport of large cap, slug or churn-turbulent bubbles. The source and sink terms in the right-hand-side of the transport equations have been established by mechanistically modeling the creation and destruction of bubbles due to major bubble interaction mechanisms. In the present paper, the interfacial area transport equations currently available are reviewed to address the feasibility and reliability of the model along with extensive experimental results. These include the data from adiabatic upward air-water two-phase flow in round tubes of various sizes, from a rectangular duct, and from adiabatic co-current downward air-water two-phase flow in round pipes of two sizes. (authors)

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

Energy Technology Data Exchange (ETDEWEB)

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

2007-01-15

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

Nonoscillation of half-linear dynamic equations

Czech Academy of Sciences Publication Activity Database

Matucci, S.; Řehák, Pavel

2010-01-01

Roč. 60, č. 5 (2010), s. 1421-1429 ISSN 0898-1221 R&D Projects: GA AV ČR KJB100190701 Grant - others:GA ČR(CZ) GA201/07/0145 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear dynamic equation * time scale * (non)oscillation * Riccati technique Subject RIV: BA - General Mathematics Impact factor: 1.472, year: 2010 http://www.sciencedirect.com/science/article/pii/S0898122110004384

Subroutine for series solutions of linear differential equations

International Nuclear Information System (INIS)

Tasso, H.; Steuerwald, J.

1976-02-01

A subroutine for Taylor series solutions of systems of ordinary linear differential equations is descriebed. It uses the old idea of Lie series but allows simple implementation and is time-saving for symbolic manipulations. (orig.) [de

A fast iterative scheme for the linearized Boltzmann equation

Science.gov (United States)

Wu, Lei; Zhang, Jun; Liu, Haihu; Zhang, Yonghao; Reese, Jason M.

2017-06-01

Iterative schemes to find steady-state solutions to the Boltzmann equation are efficient for highly rarefied gas flows, but can be very slow to converge in the near-continuum flow regime. In this paper, a synthetic iterative scheme is developed to speed up the solution of the linearized Boltzmann equation by penalizing the collision operator L into the form L = (L + Nδh) - Nδh, where δ is the gas rarefaction parameter, h is the velocity distribution function, and N is a tuning parameter controlling the convergence rate. The velocity distribution function is first solved by the conventional iterative scheme, then it is corrected such that the macroscopic flow velocity is governed by a diffusion-type equation that is asymptotic-preserving into the Navier-Stokes limit. The efficiency of this new scheme is assessed by calculating the eigenvalue of the iteration, as well as solving for Poiseuille and thermal transpiration flows. We find that the fastest convergence of our synthetic scheme for the linearized Boltzmann equation is achieved when Nδ is close to the average collision frequency. The synthetic iterative scheme is significantly faster than the conventional iterative scheme in both the transition and the near-continuum gas flow regimes. Moreover, due to its asymptotic-preserving properties, the synthetic iterative scheme does not need high spatial resolution in the near-continuum flow regime, which makes it even faster than the conventional iterative scheme. Using this synthetic scheme, with the fast spectral approximation of the linearized Boltzmann collision operator, Poiseuille and thermal transpiration flows between two parallel plates, through channels of circular/rectangular cross sections and various porous media are calculated over the whole range of gas rarefaction. Finally, the flow of a Ne-Ar gas mixture is solved based on the linearized Boltzmann equation with the Lennard-Jones intermolecular potential for the first time, and the difference

Asymptotic solutions and spectral theory of linear wave equations

International Nuclear Information System (INIS)

Adam, J.A.

1982-01-01

This review contains two closely related strands. Firstly the asymptotic solution of systems of linear partial differential equations is discussed, with particular reference to Lighthill's method for obtaining the asymptotic functional form of the solution of a scalar wave equation with constant coefficients. Many of the applications of this technique are highlighted. Secondly, the methods and applications of the theory of the reduced (one-dimensional) wave equation - particularly spectral theory - are discussed. While the breadth of application and power of the techniques is emphasised throughout, the opportunity is taken to present to a wider readership, developments of the methods which have occured in some aspects of astrophysical (particularly solar) and geophysical fluid dynamics. It is believed that the topics contained herein may be of relevance to the applied mathematician or theoretical physicist interest in problems of linear wave propagation in these areas. (orig./HSI)

Transport methods: general. 8. Formulation of Transport Equation in a Split Form

International Nuclear Information System (INIS)

Stancic, V.

2001-01-01

The singular eigenfunction expansion method has enabled the application of functional analysis methods in transport theory. However, when applying it, the users were discouraged, since in most problems, including slab problems, an extra problem has occurred. It appears necessary to solve the Fredholm integral equation in order to determine the expansion coefficients. There are several reasons for this difficulty. One reason might be the use of the full-range expansion techniques even in the regions where the function is singular. Such an example is the free boundary condition that requires the distribution to be equal to zero. Moreover, at μ = 0, the transport equation becomes an integral one. Both reasons motivated us to redefine the transport equation in a more natural way. Similar to scattering theory, here we define the flux distribution as a direct sum of forward- and backward-directed neutrons, e.g., μ ≥ 0 and μ < 0, respectively. As a result, the plane geometry transport equation is being split into coupled-pair equations. Further, using an appropriate transformation, this pair of equations reduces to a self-adjoint one having the same form as the known full-range single flux. It is interesting that all the methods of full-range theory are applicable here provided the flux as well as the transformed transport operator are two-dimensional matrices. Applying this to the slab problem, we find explicit expressions for reflected and transmitted particles caused by an arbitrary plane source. That is the news in this paper. Because of space constraints, only fundamentals of this approach will be presented here. We assume that the reader is familiar with this field; therefore, the applications are noted only at the end. (author)

Appearance of eigen modes for the linearized Vlasov-Poisson equation

International Nuclear Information System (INIS)

Degond, P.

1983-01-01

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

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

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

Linear algebra a first course with applications to differential equations

CERN Document Server

Apostol, Tom M

2014-01-01

Developed from the author's successful two-volume Calculus text this book presents Linear Algebra without emphasis on abstraction or formalization. To accommodate a variety of backgrounds, the text begins with a review of prerequisites divided into precalculus and calculus prerequisites. It continues to cover vector algebra, analytic geometry, linear spaces, determinants, linear differential equations and more.

Approximate Controllability for Linear Stochastic Differential Equations in Infinite Dimensions

International Nuclear Information System (INIS)

Goreac, D.

2009-01-01

The objective of the paper is to investigate the approximate controllability property of a linear stochastic control system with values in a separable real Hilbert space. In a first step we prove the existence and uniqueness for the solution of the dual linear backward stochastic differential equation. This equation has the particularity that in addition to an unbounded operator acting on the Y-component of the solution there is still another one acting on the Z-component. With the help of this dual equation we then deduce the duality between approximate controllability and observability. Finally, under the assumption that the unbounded operator acting on the state process of the forward equation is an infinitesimal generator of an exponentially stable semigroup, we show that the generalized Hautus test provides a necessary condition for the approximate controllability. The paper generalizes former results by Buckdahn, Quincampoix and Tessitore (Stochastic Partial Differential Equations and Applications, Series of Lecture Notes in Pure and Appl. Math., vol. 245, pp. 253-260, Chapman and Hall, London, 2006) and Goreac (Applied Analysis and Differential Equations, pp. 153-164, World Scientific, Singapore, 2007) from the finite dimensional to the infinite dimensional case

Spectrum of the linearized operator for the Ginzburg-Landau equation

Directory of Open Access Journals (Sweden)

Tai-Chia Lin

2000-06-01

Full Text Available We study the spectrum of the linearized operator for the Ginzburg-Landau equation about a symmetric vortex solution with degree one. We show that the smallest eigenvalue of the linearized operator has multiplicity two, and then we describe its behavior as a small parameter approaches zero. We also find a positive lower bound for all the other eigenvalues, and find estimates of the first eigenfunction. Then using these results, we give partial results on the dynamics of vortices in the nonlinear heat and Schrodinger equations.

Stability of numerical method for semi-linear stochastic pantograph differential equations

Directory of Open Access Journals (Sweden)

Yu Zhang

2016-01-01

Full Text Available Abstract As a particular expression of stochastic delay differential equations, stochastic pantograph differential equations have been widely used in nonlinear dynamics, quantum mechanics, and electrodynamics. In this paper, we mainly study the stability of analytical solutions and numerical solutions of semi-linear stochastic pantograph differential equations. Some suitable conditions for the mean-square stability of an analytical solution are obtained. Then we proved the general mean-square stability of the exponential Euler method for a numerical solution of semi-linear stochastic pantograph differential equations, that is, if an analytical solution is stable, then the exponential Euler method applied to the system is mean-square stable for arbitrary step-size h > 0 $h>0$ . Numerical examples further illustrate the obtained theoretical results.

Quantum osp-invariant non-linear Schroedinger equation

International Nuclear Information System (INIS)

Kulish, P.P.

1985-04-01

The generalizations of the non-linear Schroedinger equation (NS) associated with the orthosymplectic superalgebras are formulated. The simplest osp(1/2)-NS model is solved by the quantum inverse scattering method on a finite interval under periodic boundary conditions as well as on the wholeline in the case of a finite number of excitations. (author)

Supporting Students' Understanding of Linear Equations with One Variable Using Algebra Tiles

Science.gov (United States)

Saraswati, Sari; Putri, Ratu Ilma Indra; Somakim

2016-01-01

This research aimed to describe how algebra tiles can support students' understanding of linear equations with one variable. This article is a part of a larger research on learning design of linear equations with one variable using algebra tiles combined with balancing method. Therefore, it will merely discuss one activity focused on how students…

Novel algorithm of large-scale simultaneous linear equations

International Nuclear Information System (INIS)

Fujiwara, T; Hoshi, T; Yamamoto, S; Sogabe, T; Zhang, S-L

2010-01-01

We review our recently developed methods of solving large-scale simultaneous linear equations and applications to electronic structure calculations both in one-electron theory and many-electron theory. This is the shifted COCG (conjugate orthogonal conjugate gradient) method based on the Krylov subspace, and the most important issue for applications is the shift equation and the seed switching method, which greatly reduce the computational cost. The applications to nano-scale Si crystals and the double orbital extended Hubbard model are presented.

Stability problems for linear hyperbolic systems

International Nuclear Information System (INIS)

Eckhoff, K.S.

1975-05-01

The stability properties for the trivial solution of a general linear hyperbolic system of partial differential equations of the first order are studied. It is shown that results may be obtained by studying the stability properties of certain systems of ordinary differential equations which can be constructed from the hyperbolic system (the so-called transport equations). In some cases the associated stability problem for the transport equations can in fact be shown to be equivalent to the stability problem for the hyperbolic system, but in general the transport equations will only give the necessary conditions for stability. (Auth.)

Dark energy cosmology with generalized linear equation of state

International Nuclear Information System (INIS)

Babichev, E; Dokuchaev, V; Eroshenko, Yu

2005-01-01

Dark energy with the usually used equation of state p = wρ, where w const 0 ), where the constants α and ρ 0 are free parameters. This non-homogeneous linear equation of state provides the description of both hydrodynamically stable (α > 0) and unstable (α < 0) fluids. In particular, the considered cosmological model describes the hydrodynamically stable dark (and phantom) energy. The possible types of cosmological scenarios in this model are determined and classified in terms of attractors and unstable points by using phase trajectories analysis. For the dark energy case, some distinctive types of cosmological scenarios are possible: (i) the universe with the de Sitter attractor at late times, (ii) the bouncing universe, (iii) the universe with the big rip and with the anti-big rip. In the framework of a linear equation of state the universe filled with a phantom energy, w < -1, may have either the de Sitter attractor or the big rip

POSITIVE SOLUTIONS TO SEMI-LINEAR SECOND-ORDER ORDINARY DIFFERENTIAL EQUATIONS IN BANACH SPACE

Institute of Scientific and Technical Information of China (English)

无

2008-01-01

In this paper,we study the existence of positive periodic solution to some second- order semi-linear differential equation in Banach space.By the fixed point index theory, we prove that the semi-linear differential equation has two positive periodic solutions.

Focal decompositions for linear differential equations of the second order

Directory of Open Access Journals (Sweden)

L. Birbrair

2003-01-01

two-points problems to itself such that the image of the focal decomposition associated to the first equation is a focal decomposition associated to the second one. In this paper, we present a complete classification for linear second-order equations with respect to this equivalence relation.

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

International Nuclear Information System (INIS)

Zawaideh, E.S.

1985-01-01

A new set of two-fluid equations which are valid from collisional to weakly collisional limits are derived. Starting from gyrokinetic equations in flux coordinates with no zeroth order drifts, a set of moment equations describing plasma transport along the field lines of a space and time dependent magnetic field are 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. The new transport equations are used to study the effects of collisionality, magnetic field structure, and plasma anisotropy on plasma parallel transport. Numerical examples comparing these equations with conventional transport equations show that the conventional equations may contain large errors near the sound speed (M approx. = 1). It is also found that plasma anisotropy, which is not included in the conventional equations, is a critical parameter in determining plasma transport in varying magnetic field. The new transport equations are also used to study axial confinement in multiple mirror devices from the strongly to weakly collisional regime. A new ion conduction model was worked out to extend the regime of validity of the transport equations to the low density multiple mirror regime

Approximate solution of the transport equation by methods of Galerkin type

International Nuclear Information System (INIS)

Pitkaranta, J.

1977-01-01

Questions of the existence, uniqueness, and convergence of approximate solutions of transport equations by methods of the Galerkin type (where trial and weighting functions are the same) are discussed. The results presented do not exclude the infinite-dimensional case. Two strategies can be followed in the variational approximation of the transport operator: one proceeds from the original form of the transport equation, while the other is based on the partially symmetrized equation. Both principles are discussed in this paper. The transport equation is assumed in a discretized multigroup form

Analytical solution to the hybrid diffusion-transport equation

International Nuclear Information System (INIS)

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

1986-01-01

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

Prolongation structure and linear eigenvalue equations for Einstein-Maxwell fields

International Nuclear Information System (INIS)

Kramer, D.; Neugebauer, G.

1981-01-01

The Einstein-Maxwell equations for stationary axisymmetric exterior fields are shown to be the integrability conditions of a set of linear eigenvalue equations for pseudopotentials. Using the method of Wahlquist and Estabrook (J. Math Phys.; 16:1 (1975)) it is shown that the prolongation structure of the Einstein-Maxwell equations contains the SU(2,1) Lie algebra. A new mapping of known solutions to other solutions has been found. (author)

Swarm analysis by using transport equations, 1

International Nuclear Information System (INIS)

Dote, Toshihiko; Shimada, Masatoshi

1980-01-01

By evolving Maxwell-Boltzmann transport equations, various quantities on swarm of charged particles have been analyzed. Although this treatment is properly general, and common transport equations for charged particles ought to be given, in particular, equations only for electrons were presented here. The relation between the random energy and the drift energy was first derived and the general expression of the electron velocity was deduced too. For a simple example, one dimensional steady-state electron swarm in a uniform medium was treated. Electron swarm characteristics numerically calculated in He, Ne or Ar exhibited some interesting properties, which were physically clearly elucidated. These results were also compared with several data already published. Agreements between them were qualitatively rather well in detailed structures. (author)

From the hypergeometric differential equation to a non-linear Schrödinger one

International Nuclear Information System (INIS)

Plastino, A.; Rocca, M.C.

2015-01-01

We show that the q-exponential function is a hypergeometric function. Accordingly, it obeys the hypergeometric differential equation. We demonstrate that this differential equation can be transformed into a non-linear Schrödinger equation (NLSE). This NLSE exhibits both similarities and differences vis-a-vis the Nobre–Rego-Monteiro–Tsallis one. - Highlights: • We show that the q-exponential is a hypergeometric function. • It thus obeys the hypergeometric differential equation (HDE). • We show that the HDE can be cast as a non-linear Schrödinger equation. • This is different from the Nobre, Rego-Monteiro, Tsallis one.

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.

Mean Transit Time and Mean Residence Time for Linear Diffusion–Convection–Reaction Transport System

Directory of Open Access Journals (Sweden)

Jacek Waniewski

2007-01-01

Full Text Available Characteristic times for transport processes in biological systems may be evaluated as mean transit times (MTTs (for transit states or mean residence times (MRT (for steady states. It is shown in a general framework of a (linear reaction–diffusion–convection equation that these two times are related. Analytical formulas are also derived to calculate moments of exit time distribution using solutions for a stationary state of the system.

Construction of local and non-local conservation laws for non-linear field equations

International Nuclear Information System (INIS)

Vladimirov, V.S.; Volovich, I.V.

1984-08-01

A method of constructing conserved currents for non-linear field equations is presented. More explicitly for non-linear equations, which can be derived from compatibility conditions of some linear system with a parameter, a procedure of obtaining explicit expressions for local and non-local currents is developed. Some examples such as the classical Heisenberg spin chain and supersymmetric Yang-Mills theory are considered. (author)

Fundaments of transport equation splitting and the eigenvalue problem

International Nuclear Information System (INIS)

Stancic, V.

2000-01-01

In order to remove some singularities concerning the boundary conditions of one dimensional transport equation, a split form of transport equation describing the forward i.e. μ≥0, and a backward μ<0 directed neutrons is being proposed here. The eigenvalue problem has also been considered here (author)

A Hamiltonian structure for the linearized Einstein vacuum field equations

International Nuclear Information System (INIS)

Torres del Castillo, G.F.

1991-01-01

By considering the Einstein vacuum field equations linearized about the Minkowski metric, the evolution equations for the gauge-invariant quantities characterizing the gravitational field are written in a Hamiltonian form. A Poisson bracket between functionals of the field, compatible with the constraints satisfied by the field variables, is obtained (Author)

Collective spin by linearization of the Schrodinger equation for nuclear collective motion

International Nuclear Information System (INIS)

Greiner, M.; Scheid, W.; Herrmann, R.

1988-01-01

The free Schrodinger equation for multipole degrees of freedom is linearized so that energy and momentum operators appear only in first order. As an example, the authors demonstrate the linearization procedure for quadrupole degrees of freedom. The wave function solving this equation carries a spin. The authors derive the operator of the collective spin and its eigen values depending on multipolarity

Dynamic modeling of interfacial structures via interfacial area transport equation

International Nuclear Information System (INIS)

Seungjin, Kim; Mamoru, Ishii

2004-01-01

Full text of publication follows:In the current thermal-hydraulic system analysis codes using the two-fluid model, the empirical correlations that are based on the two-phase flow regimes and regime transition criteria are being employed as closure relations for the interfacial transfer terms. Due to its inherent shortcomings, however, such static correlations are inaccurate and present serious problems in the numerical analysis. In view of this, a new dynamic approach employing the interfacial area transport equation has been studied. The interfacial area transport equation dynamically models the two-phase flow regime transitions and predicts continuous change of the interfacial area concentration along the flow field. Hence, when employed in the thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Therefore, the interfacial area transport equation can make a leapfrog improvement in the current capability of the two-fluid model from both scientific and practical point of view. Accounting for the substantial differences in the transport phenomena of various sizes of bubbles, the two-group interfacial area transport equations have been developed. The group 1 equation describes the transport of small-dispersed bubbles that are either distorted or spherical in shapes, and the group 2 equation describes the transport of large cap, slug or churn-turbulent bubbles. The source and sink terms in the right hand-side of the transport equations have been established by mechanistically modeling the creation and destruction of bubbles due to major bubble interaction mechanisms. The coalescence mechanisms include the random collision driven by turbulence, and the entrainment of trailing bubbles in the wake region of the preceding bubble. The disintegration mechanisms include the break-up by turbulence impact, shearing-off at the rim of large cap bubbles and the break-up of large cap

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

International Nuclear Information System (INIS)

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

2009-01-01

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

Solutions of the linearized Bach-Einstein equation in the static spherically symmetric case

International Nuclear Information System (INIS)

Schmidt, H.J.

1985-01-01

The Bach-Einstein equation linearized around Minkowski space-time is completely solved. The set of solutions depends on three parameters; a two-parameter subset of it becomes asymptotically flat. In that region the gravitational potential is of the type phi = -m/r + epsilon exp (-r/l). Because of the different asymptotic behaviour of both terms, it became necessary to linearize also around the Schwarzschild solution phi = -m/r. The linearized equation resulting in this case is discussed using qualitative methods. The result is that for m = 2l phi = -m/r + epsilon r -2 exp (-r/l) u, where u is some bounded function; m is arbitrary and epsilon again small. Further, the relation between the solution of the linearized and the full equation is discussed. (author)

On a representation of linear differential equations

Czech Academy of Sciences Publication Activity Database

Neuman, František

2010-01-01

Roč. 52, 1-2 (2010), s. 355-360 ISSN 0895-7177 Grant - others:GA ČR(CZ) GA201/08/0469 Institutional research plan: CEZ:AV0Z10190503 Keywords : Brandt and Ehresmann groupoinds * transformations * canonical forms * linear differential equations Subject RIV: BA - General Mathematics Impact factor: 1.066, year: 2010 http://www.sciencedirect.com/science/article/pii/S0895717710001184

CFORM- LINEAR CONTROL SYSTEM DESIGN AND ANALYSIS: CLOSED FORM SOLUTION AND TRANSIENT RESPONSE OF THE LINEAR DIFFERENTIAL EQUATION

Science.gov (United States)

Jamison, J. W.

1994-01-01

CFORM was developed by the Kennedy Space Center Robotics Lab to assist in linear control system design and analysis using closed form and transient response mechanisms. The program computes the closed form solution and transient response of a linear (constant coefficient) differential equation. CFORM allows a choice of three input functions: the Unit Step (a unit change in displacement); the Ramp function (step velocity); and the Parabolic function (step acceleration). It is only accurate in cases where the differential equation has distinct roots, and does not handle the case for roots at the origin (s=0). Initial conditions must be zero. Differential equations may be input to CFORM in two forms - polynomial and product of factors. In some linear control analyses, it may be more appropriate to use a related program, Linear Control System Design and Analysis (KSC-11376), which uses root locus and frequency response methods. CFORM was written in VAX FORTRAN for a VAX 11/780 under VAX VMS 4.7. It has a central memory requirement of 30K. CFORM was developed in 1987.

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

Energy Technology Data Exchange (ETDEWEB)

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

2016-07-05

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

GDTM-Padé technique for the non-linear differential-difference equation

Directory of Open Access Journals (Sweden)

Lu Jun-Feng

2013-01-01

Full Text Available This paper focuses on applying the GDTM-Padé technique to solve the non-linear differential-difference equation. The bell-shaped solitary wave solution of Belov-Chaltikian lattice equation is considered. Comparison between the approximate solutions and the exact ones shows that this technique is an efficient and attractive method for solving the differential-difference equations.

On the Liouvillian solution of second-order linear differential equations and algebraic invariant curves

International Nuclear Information System (INIS)

Man, Yiu-Kwong

2010-01-01

In this communication, we present a method for computing the Liouvillian solution of second-order linear differential equations via algebraic invariant curves. The main idea is to integrate Kovacic's results on second-order linear differential equations with the Prelle-Singer method for computing first integrals of differential equations. Some examples on using this approach are provided. (fast track communication)

On the economical solution method for a system of linear algebraic equations

Directory of Open Access Journals (Sweden)

Jan Awrejcewicz

2004-01-01

Full Text Available The present work proposes a novel optimal and exact method of solving large systems of linear algebraic equations. In the approach under consideration, the solution of a system of algebraic linear equations is found as a point of intersection of hyperplanes, which needs a minimal amount of computer operating storage. Two examples are given. In the first example, the boundary value problem for a three-dimensional stationary heat transfer equation in a parallelepiped in ℝ3 is considered, where boundary value problems of first, second, or third order, or their combinations, are taken into account. The governing differential equations are reduced to algebraic ones with the help of the finite element and boundary element methods for different meshes applied. The obtained results are compared with known analytical solutions. The second example concerns computation of a nonhomogeneous shallow physically and geometrically nonlinear shell subject to transversal uniformly distributed load. The partial differential equations are reduced to a system of nonlinear algebraic equations with the error of O(hx12+hx22. The linearization process is realized through either Newton method or differentiation with respect to a parameter. In consequence, the relations of the boundary condition variations along the shell side and the conditions for the solution matching are reported.

FUNDAMENTAL MATRIX OF LINEAR CONTINUOUS SYSTEM IN THE PROBLEM OF ESTIMATING ITS TRANSPORT DELAY

Directory of Open Access Journals (Sweden)

N. A. Dudarenko

2014-09-01

Full Text Available The paper deals with the problem of quantitative estimation for transport delay of linear continuous systems. The main result is received by means of fundamental matrix of linear differential equations solutions specified in the normal Cauchy form for the cases of SISO and MIMO systems. Fundamental matrix has the dual property. It means that the weight function of the system can be formed as a free motion of systems. Last one is generated by the vector of initial system conditions, which coincides with the matrix input of the system being researched. Thus, using the properties of the system- solving for fundamental matrix has given the possibility to solve the problem of estimating transport linear continuous system delay without the use of derivation procedure in hardware environment and without formation of exogenous Dirac delta function. The paper is illustrated by examples. The obtained results make it possible to solve the problem of modeling the pure delay links using consecutive chain of aperiodic links of the first order with the equal time constants. Modeling results have proved the correctness of obtained computations. Knowledge of transport delay can be used when configuring multi- component technological complexes and in the diagnosis of their possible functional degeneration.

A modular spherical harmonics approach to the neutron transport equation

International Nuclear Information System (INIS)

Inanc, F.; Rohach, A.F.

1989-01-01

A modular nodal method was developed for solving the neutron transport equation in 2-D xy coordinates. The spherical harmonic expansion was used for approximating the second-order even-parity form of the neutron transport equation. The boundary conditions of the spherical harmonics approximation were derived in a form to have forms analogous to the partial currents in the neutron diffusion equation. Relations were developed for generating both the second-order spherical harmonic equations and the boundary conditions in an automated computational algorithm. Nodes using different orders of the spherical harmonics approximation to the transport equation were interfaced through mixed-type boundary conditions. The determination of spherical harmonic orders implemented in the nodes were determined by the scheme in an automated manner. Results of the method compared favorably to benchmark problems. (author)

Generalized heat-transport equations: parabolic and hyperbolic models

Science.gov (United States)

Rogolino, Patrizia; Kovács, Robert; Ván, Peter; Cimmelli, Vito Antonio

2018-03-01

We derive two different generalized heat-transport equations: the most general one, of the first order in time and second order in space, encompasses some well-known heat equations and describes the hyperbolic regime in the absence of nonlocal effects. Another, less general, of the second order in time and fourth order in space, is able to describe hyperbolic heat conduction also in the presence of nonlocal effects. We investigate the thermodynamic compatibility of both models by applying some generalizations of the classical Liu and Coleman-Noll procedures. In both cases, constitutive equations for the entropy and for the entropy flux are obtained. For the second model, we consider a heat-transport equation which includes nonlocal terms and study the resulting set of balance laws, proving that the corresponding thermal perturbations propagate with finite speed.

On the stability, the periodic solutions and the resolution of certain types of non linear equations, and of non linearly coupled systems of these equations, appearing in betatronic oscillations

International Nuclear Information System (INIS)

Valat, J.

1960-12-01

Universal stability diagrams have been calculated and experimentally checked for Hill-Meissner type equations with square-wave coefficients. The study of these equations in the phase-plane has then made it possible to extend the periodic solution calculations to the case of non-linear differential equations with periodic square-wave coefficients. This theory has been checked experimentally. For non-linear coupled systems with constant coefficients, a search was first made for solutions giving an algebraic motion. The elliptical and Fuchs's functions solve such motions. The study of non-algebraic motions is more delicate, apart from the study of nonlinear Lissajous's motions. A functional analysis shows that it is possible however in certain cases to decouple the system and to find general solutions. For non-linear coupled systems with periodic square-wave coefficients it is then possible to calculate the conditions leading to periodic solutions, if the two non-linear associated systems with constant coefficients fall into one of the categories of the above paragraph. (author) [fr

Inhomogeneous linear equation in Rota-Baxter algebra

OpenAIRE

Pietrzkowski, Gabriel

2014-01-01

We consider a complete filtered Rota-Baxter algebra of weight $\\lambda$ over a commutative ring. Finding the unique solution of a non-homogeneous linear algebraic equation in this algebra, we generalize Spitzer's identity in both commutative and non-commutative cases. As an application, considering the Rota-Baxter algebra of power series in one variable with q-integral as the Rota-Baxter operator, we show certain Eulerian identities.

Solution of systems of linear algebraic equations by the method of summation of divergent series

International Nuclear Information System (INIS)

Kirichenko, G.A.; Korovin, Ya.S.; Khisamutdinov, M.V.; Shmojlov, V.I.

2015-01-01

A method for solving systems of linear algebraic equations has been proposed on the basis on the summation of the corresponding continued fractions. The proposed algorithm for solving systems of linear algebraic equations is classified as direct algorithms providing an exact solution in a finite number of operations. Examples of solving systems of linear algebraic equations have been presented and the effectiveness of the algorithm has been estimated [ru

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

International Nuclear Information System (INIS)

Grigoriu, Mircea; Samorodnitsky, Gennady

2004-01-01

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

Variance analysis of the Monte-Carlo perturbation source method in inhomogeneous linear particle transport problems

International Nuclear Information System (INIS)

Noack, K.

1982-01-01

The perturbation source method may be a powerful Monte-Carlo means to calculate small effects in a particle field. In a preceding paper we have formulated this methos in inhomogeneous linear particle transport problems describing the particle fields by solutions of Fredholm integral equations and have derived formulae for the second moment of the difference event point estimator. In the present paper we analyse the general structure of its variance, point out the variance peculiarities, discuss the dependence on certain transport games and on generation procedures of the auxiliary particles and draw conclusions to improve this method

Factorization of a class of almost linear second-order differential equations

International Nuclear Information System (INIS)

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

2007-01-01

A general type of almost linear second-order differential equations, which are directly related to several interesting physical problems, is characterized. The solutions of these equations are obtained using the factorization technique, and their non-autonomous invariants are also found by means of scale transformations

Variations in the Solution of Linear First-Order Differential Equations. Classroom Notes

Science.gov (United States)

Seaman, Brian; Osler, Thomas J.

2004-01-01

A special project which can be given to students of ordinary differential equations is described in detail. Students create new differential equations by changing the dependent variable in the familiar linear first-order equation (dv/dx)+p(x)v=q(x) by means of a substitution v=f(y). The student then creates a table of the new equations and…

Problems of linear electron (polaron) transport theory in semiconductors

CERN Document Server

Klinger, M I

1979-01-01

Problems of Linear Electron (Polaron) Transport Theory in Semiconductors summarizes and discusses the development of areas in electron transport theory in semiconductors, with emphasis on the fundamental aspects of the theory and the essential physical nature of the transport processes. The book is organized into three parts. Part I focuses on some general topics in the theory of transport phenomena: the general dynamical theory of linear transport in dissipative systems (Kubo formulae) and the phenomenological theory. Part II deals with the theory of polaron transport in a crystalline semicon

Homogenization of the critically spectral equation in neutron transport

International Nuclear Information System (INIS)

Allaire, G.; Paris-6 Univ., 75; Bal, G.

1998-01-01

We address the homogenization of an eigenvalue problem for the neutron transport equation in a periodic heterogeneous domain, modeling the criticality study of nuclear reactor cores. We prove that the neutron flux, corresponding to the first and unique positive eigenvector, can be factorized in the product of two terms, up to a remainder which goes strongly to zero with the period. On terms is the first eigenvector of the transport equation in the periodicity cell. The other term is the first eigenvector of a diffusion equation in the homogenized domain. Furthermore, the corresponding eigenvalue gives a second order corrector for the eigenvalue of the heterogeneous transport problem. This result justifies and improves the engineering procedure used in practice for nuclear reactor cores computations. (author)

Homogenization of the critically spectral equation in neutron transport

Energy Technology Data Exchange (ETDEWEB)

Allaire, G. [CEA Saclay, 91 - Gif-sur-Yvette (France). Dept. de Mecanique et de Technologie]|[Paris-6 Univ., 75 (France). Lab. d' Analyse Numerique; Bal, G. [Electricite de France (EDF), 92 - Clamart (France). Direction des Etudes et Recherches

1998-07-01

We address the homogenization of an eigenvalue problem for the neutron transport equation in a periodic heterogeneous domain, modeling the criticality study of nuclear reactor cores. We prove that the neutron flux, corresponding to the first and unique positive eigenvector, can be factorized in the product of two terms, up to a remainder which goes strongly to zero with the period. On terms is the first eigenvector of the transport equation in the periodicity cell. The other term is the first eigenvector of a diffusion equation in the homogenized domain. Furthermore, the corresponding eigenvalue gives a second order corrector for the eigenvalue of the heterogeneous transport problem. This result justifies and improves the engineering procedure used in practice for nuclear reactor cores computations. (author)

From statistic mechanic outside equilibrium to transport equations

International Nuclear Information System (INIS)

Balian, R.

1995-01-01

This lecture notes give a synthetic view on the foundations of non-equilibrium statistical mechanics. The purpose is to establish the transport equations satisfied by the relevant variables, starting from the microscopic dynamics. The Liouville representation is introduced, and a projection associates with any density operator , for given choice of relevant observables, a reduced density operator. An exact integral-differential equation for the relevant variables is thereby derived. A short-memory approximation then yields the transport equations. A relevant entropy which characterizes the coarseness of the description is associated with each level of description. As an illustration, the classical gas, with its three levels of description and with the Chapman-Enskog method, is discussed. (author). 3 figs., 5 refs

A Proposed Method for Solving Fuzzy System of Linear Equations

Directory of Open Access Journals (Sweden)

Reza Kargar

2014-01-01

Full Text Available This paper proposes a new method for solving fuzzy system of linear equations with crisp coefficients matrix and fuzzy or interval right hand side. Some conditions for the existence of a fuzzy or interval solution of m×n linear system are derived and also a practical algorithm is introduced in detail. The method is based on linear programming problem. Finally the applicability of the proposed method is illustrated by some numerical examples.

On a class of fourth order linear recurrence equations

Directory of Open Access Journals (Sweden)

Sui-Sun Cheng

1984-01-01

Full Text Available This paper is concerned with sequences that satisfy a class of fourth order linear recurrence equations. Basic properties of such sequences are derived. In addition, we discuss the oscillatory and nonoscillatory behavior of such sequences.

Could solitons be adiabatic invariants attached to certain non linear equations

International Nuclear Information System (INIS)

Lochak, P.

1984-01-01

Arguments are given to support the claim that solitons should be the adiabatic invariants associated to certain non linear partial differential equations; a precise mathematical form of this conjecture is then stated. As a particular case of the conjecture, the Korteweg-de Vries equation is studied. (Auth.)

On a Linear Equation Arising in Isometric Embedding of Torus-like Surface

Institute of Scientific and Technical Information of China (English)

Chunhe LI

2009-01-01

The solvability of a linear equation and the regularity of the solution are discussed.The equation is arising in a geometric problem which is concerned with the realization of Alexandroff's positive annul in R3.

Basic equations of interfacial area transport in gas-liquid two-phase flow

International Nuclear Information System (INIS)

Kataoka, I.; Yoshida, K.; Naitoh, M.; Okada, H.; Morii, T.

2011-01-01

The rigorous and consistent formulations of basic equations of interfacial area transport were derived using correlation functions of characteristic function of each phase and velocities of each phase. Turbulent transport term of interfacial area concentration was consistently derived and related to the difference between interfacial velocity and averaged velocity of each phase. Constitutive equations of turbulent transport terms of interfacial area concentration were proposed for bubbly flow. New transport model and constitutive equations were developed for churn flow. These models and constitutive equations are validated by experimental data of radial distributions of interfacial area concentration in bubbly and churn flow. (author)

Stochastic modeling of mode interactions via linear parabolized stability equations

Science.gov (United States)

Ran, Wei; Zare, Armin; Hack, M. J. Philipp; Jovanovic, Mihailo

2017-11-01

Low-complexity approximations of the Navier-Stokes equations have been widely used in the analysis of wall-bounded shear flows. In particular, the parabolized stability equations (PSE) and Floquet theory have been employed to capture the evolution of primary and secondary instabilities in spatially-evolving flows. We augment linear PSE with Floquet analysis to formally treat modal interactions and the evolution of secondary instabilities in the transitional boundary layer via a linear progression. To this end, we leverage Floquet theory by incorporating the primary instability into the base flow and accounting for different harmonics in the flow state. A stochastic forcing is introduced into the resulting linear dynamics to model the effect of nonlinear interactions on the evolution of modes. We examine the H-type transition scenario to demonstrate how our approach can be used to model nonlinear effects and capture the growth of the fundamental and subharmonic modes observed in direct numerical simulations and experiments.

KAM for the non-linear Schroedinger equation

CERN Document Server

Eliasson, L H

2006-01-01

We consider the $d$-dimensional nonlinear Schr\\"o\\-dinger equation under periodic boundary conditions:-i\\dot u=\\Delta u+V(x)*u+\\ep|u|^2u;\\quad u=u(t,x),\\;x\\in\\T^dwhere $V(x)=\\sum \\hat V(a)e^{i\\sc{a,x}}$ is an analytic function with $\\hat V$ real. (This equation is a popular model for the `real' NLS equation, where instead of the convolution term $V*u$ we have the potential term $Vu$.) For $\\ep=0$ the equation is linear and has time--quasi-periodic solutions $u$,u(t,x)=\\sum_{s\\in \\AA}\\hat u_0(a)e^{i(|a|^2+\\hat V(a))t}e^{i\\sc{a,x}}, \\quad 0<|\\hat u_0(a)|\\le1,where $\\AA$ is any finite subset of $\\Z^d$. We shall treat $\\omega_a=|a|^2+\\hat V(a)$, $a\\in\\AA$, as free parameters in some domain $U\\subset\\R^{\\AA}$. This is a Hamiltonian system in infinite degrees of freedom, degenerate but with external parameters, and we shall describe a KAM-theory which, in particular, will have the following consequence: \\smallskip {\\it If $|\\ep|$ is sufficiently small, then there is a large subset $U'$ of $U$ such that for all $...

Minimal solution of linear formed fuzzy matrix equations

Directory of Open Access Journals (Sweden)

Maryam Mosleh

2012-10-01

Full Text Available In this paper according to the structured element method, the $mimes n$ inconsistent fuzzy matrix equation $Ailde{X}=ilde{B},$ which are linear formed by fuzzy structured element, is investigated. The necessary and sufficient condition for the existence of a fuzzy solution is also discussed. some examples are presented to illustrate the proposed method.

General particle transport equation. Final report

International Nuclear Information System (INIS)

Lafi, A.Y.; Reyes, J.N. Jr.

1994-12-01

The general objectives of this research are as follows: (1) To develop fundamental models for fluid particle coalescence and breakage rates for incorporation into statistically based (Population Balance Approach or Monte Carlo Approach) two-phase thermal hydraulics codes. (2) To develop fundamental models for flow structure transitions based on stability theory and fluid particle interaction rates. This report details the derivation of the mass, momentum and energy conservation equations for a distribution of spherical, chemically non-reacting fluid particles of variable size and velocity. To study the effects of fluid particle interactions on interfacial transfer and flow structure requires detailed particulate flow conservation equations. The equations are derived using a particle continuity equation analogous to Boltzmann's transport equation. When coupled with the appropriate closure equations, the conservation equations can be used to model nonequilibrium, two-phase, dispersed, fluid flow behavior. Unlike the Eulerian volume and time averaged conservation equations, the statistically averaged conservation equations contain additional terms that take into account the change due to fluid particle interfacial acceleration and fluid particle dynamics. Two types of particle dynamics are considered; coalescence and breakage. Therefore, the rate of change due to particle dynamics will consider the gain and loss involved in these processes and implement phenomenological models for fluid particle breakage and coalescence

A New Theory of Non-Linear Thermo-Elastic Constitutive Equation of Isotropic Hyperelastic Materials

Science.gov (United States)

Li, Chen; Liao, Yufei

2018-03-01

Considering the influence of temperature and strain variables on materials. According to the relationship of conjugate stress-strain, a complete and irreducible non-linear constitutive equation of isotropic hyperelastic materials is derived and the constitutive equations of 16 types of isotropic hyperelastic materials are given we study the transformation methods and routes of 16 kinds of constitutive equations and the study proves that transformation of two forms of constitutive equation. As an example of application, the non-linear thermo-elastic constitutive equation of isotropic hyperelastic materials is combined with the natural vulcanized rubber experimental data in the existing literature base on MATLAB, The results show that the fitting accuracy is satisfactory.

Symmetry groups of integro-differential equations for linear thermoviscoelastic materials with memory

Science.gov (United States)

Zhou, L.-Q.; Meleshko, S. V.

2017-07-01

The group analysis method is applied to a system of integro-differential equations corresponding to a linear thermoviscoelastic model. A recently developed approach for calculating the symmetry groups of such equations is used. The general solution of the determining equations for the system is obtained. Using subalgebras of the admitted Lie algebra, two classes of partially invariant solutions of the considered system of integro-differential equations are studied.

A simple Boltzmann transport equation for ballistic to diffusive transient heat transport

International Nuclear Information System (INIS)

Maassen, Jesse; Lundstrom, Mark

2015-01-01

Developing simplified, but accurate, theoretical approaches to treat heat transport on all length and time scales is needed to further enable scientific insight and technology innovation. Using a simplified form of the Boltzmann transport equation (BTE), originally developed for electron transport, we demonstrate how ballistic phonon effects and finite-velocity propagation are easily and naturally captured. We show how this approach compares well to the phonon BTE, and readily handles a full phonon dispersion and energy-dependent mean-free-path. This study of transient heat transport shows (i) how fundamental temperature jumps at the contacts depend simply on the ballistic thermal resistance, (ii) that phonon transport at early times approach the ballistic limit in samples of any length, and (iii) perceived reductions in heat conduction, when ballistic effects are present, originate from reductions in temperature gradient. Importantly, this framework can be recast exactly as the Cattaneo and hyperbolic heat equations, and we discuss how the key to capturing ballistic heat effects is to use the correct physical boundary conditions

Some Additional Remarks on the Cumulant Expansion for Linear Stochastic Differential Equations

NARCIS (Netherlands)

Roerdink, J.B.T.M.

1984-01-01

We summarize our previous results on cumulant expansions for linear stochastic differential equations with correlated multipliclative and additive noise. The application of the general formulas to equations with statistically independent multiplicative and additive noise is reconsidered in detail,

Some additional remarks on the cumulant expansion for linear stochastic differential equations

NARCIS (Netherlands)

Roerdink, J.B.T.M.

1984-01-01

We summarize our previous results on cumular expasions for linear stochastic differential equations with correlated multipliclative and additive noise. The application of the general formulas to equations with statistically independent multiplicative and additive noise is reconsidered in detail,

Hyers-Ulam stability for second-order linear differential equations with boundary conditions

Directory of Open Access Journals (Sweden)

Pasc Gavruta

2011-06-01

Full Text Available We prove the Hyers-Ulam stability of linear differential equations of second-order with boundary conditions or with initial conditions. That is, if y is an approximate solution of the differential equation $y''+ eta (x y = 0$ with $y(a = y(b =0$, then there exists an exact solution of the differential equation, near y.

Asymptotic integration of a linear fourth order differential equation of Poincaré type

Directory of Open Access Journals (Sweden)

Anibal Coronel

2015-11-01

Full Text Available This article deals with the asymptotic behavior of nonoscillatory solutions of fourth order linear differential equation where the coefficients are perturbations of constants. We define a change of variable and deduce that the new variable satisfies a third order nonlinear differential equation. We assume three hypotheses. The first hypothesis is related to the constant coefficients and set up that the characteristic polynomial associated with the fourth order linear equation has simple and real roots. The other two hypotheses are related to the behavior of the perturbation functions and establish asymptotic integral smallness conditions of the perturbations. Under these general hypotheses, we obtain four main results. The first two results are related to the application of a fixed point argument to prove that the nonlinear third order equation has a unique solution. The next result concerns with the asymptotic behavior of the solutions of the nonlinear third order equation. The fourth main theorem is introduced to establish the existence of a fundamental system of solutions and to precise the formulas for the asymptotic behavior of the linear fourth order differential equation. In addition, we present an example to show that the results introduced in this paper can be applied in situations where the assumptions of some classical theorems are not satisfied.

The Modified Enskog Equation for Mixtures

NARCIS (Netherlands)

Beijeren, H. van; Ernst, M.H.

1973-01-01

In a previous paper it was shown that a modified form of the Enskog equation, applied to mixtures of hard spheres, should be considered as the correct extension of the usual Enskog equation to the case of mixtures. The main argument was that the modified Enskog equation leads to linear transport

A Photon Free Method to Solve Radiation Transport Equations

International Nuclear Information System (INIS)

Chang, B

2006-01-01

The multi-group discrete-ordinate equations of radiation transfer is solved for the first time by Newton's method. It is a photon free method because the photon variables are eliminated from the radiation equations to yield a N group XN direction smaller but equivalent system of equations. The smaller set of equations can be solved more efficiently than the original set of equations. Newton's method is more stable than the Semi-implicit Linear method currently used by conventional radiation codes

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

On oscillation of second-order linear ordinary differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, A.; Šremr, Jiří

2011-01-01

Roč. 54, - (2011), s. 69-81 ISSN 1512-0015 Institutional research plan: CEZ:AV0Z10190503 Keywords : linear second-order ordinary differential equation * Kamenev theorem * oscillation Subject RIV: BA - General Mathematics http://www.rmi.ge/jeomj/memoirs/vol54/abs54-4.htm

An inhomogeneous wave equation and non-linear Diophantine approximation

DEFF Research Database (Denmark)

Beresnevich, V.; Dodson, M. M.; Kristensen, S.

2008-01-01

A non-linear Diophantine condition involving perfect squares and arising from an inhomogeneous wave equation on the torus guarantees the existence of a smooth solution. The exceptional set associated with the failure of the Diophantine condition and hence of the existence of a smooth solution...

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

Improved harmonic balance approach to periodic solutions of non-linear jerk equations

International Nuclear Information System (INIS)

Wu, B.S.; Lim, C.W.; Sun, W.P.

2006-01-01

An analytical approximate approach for determining periodic solutions of non-linear jerk equations involving third-order time-derivative is presented. This approach incorporates salient features of both Newton's method and the method of harmonic balance. By appropriately imposing the method of harmonic balance to the linearized equation, the approach requires only one or two iterations to predict very accurate analytical approximate solutions for a large range of initial velocity amplitude. One typical example is used to verify and illustrate the usefulness and effectiveness of the proposed approach

Supporting second grade lower secondary school students’ understanding of linear equation system in two variables using ethnomathematics

Science.gov (United States)

Nursyahidah, F.; Saputro, B. A.; Rubowo, M. R.

2018-03-01

The aim of this research is to know the students’ understanding of linear equation system in two variables using Ethnomathematics and to acquire learning trajectory of linear equation system in two variables for the second grade of lower secondary school students. This research used methodology of design research that consists of three phases, there are preliminary design, teaching experiment, and retrospective analysis. Subject of this study is 28 second grade students of Sekolah Menengah Pertama (SMP) 37 Semarang. The result of this research shows that the students’ understanding in linear equation system in two variables can be stimulated by using Ethnomathematics in selling buying tradition in Peterongan traditional market in Central Java as a context. All of strategies and model that was applied by students and also their result discussion shows how construction and contribution of students can help them to understand concept of linear equation system in two variables. All the activities that were done by students produce learning trajectory to gain the goal of learning. Each steps of learning trajectory of students have an important role in understanding the concept from informal to the formal level. Learning trajectory using Ethnomathematics that is produced consist of watching video of selling buying activity in Peterongan traditional market to construct linear equation in two variables, determine the solution of linear equation in two variables, construct model of linear equation system in two variables from contextual problem, and solving a contextual problem related to linear equation system in two variables.

Nonoscillation criteria for half-linear second order difference equations

Czech Academy of Sciences Publication Activity Database

Došlý, Ondřej; Řehák, Pavel

2001-01-01

Roč. 42, - (2001), s. 453-464 ISSN 0898-1221 R&D Projects: GA ČR GA201/98/0677; GA ČR GA201/99/0295 Keywords : half-linear difference equation%nonoscillation criteria%variational principle Subject RIV: BA - General Mathematics Impact factor: 0.383, year: 2001

Lie symmetries and differential galois groups of linear equations

NARCIS (Netherlands)

Oudshoorn, W.R.; Put, M. van der

2002-01-01

For a linear ordinary differential equation the Lie algebra of its infinitesimal Lie symmetries is compared with its differential Galois group. For this purpose an algebraic formulation of Lie symmetries is developed. It turns out that there is no direct relation between the two above objects. In

Coupling of discrete ordinates methods by transmission of boundary conditions in solving the neutron transport equation in slab geometry; Couplage de discretisations aux ordonnees discretes d`equations de transport 1D par passage de conditions frontieres

Energy Technology Data Exchange (ETDEWEB)

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

1995-10-01

Neutron transport in nuclear reactors is quite well modelled by the linear Boltzmann transport equation. Its solution is relatively easy, but unfortunately too expensive to achieve whole core computations. Thus, we have to simplify it, for example by homogenizing some physical characteristics. However, the solution may then be inaccurate. Moreover, in strongly homogeneous areas, the error may be too big. Then we would like to deal with such an inconvenient by solving the equation accurately on this area, but more coarsely away from it, so that the computation is not too expensive. This problem is the subject of a thesis. We present here some results obtained for slab geometry. The couplings between the fine and coarse discretization regions could be conceived in a number of approaches. Here, we only deal with the coupling at crossing the interface between two sub-domains. In the first section, we present the coupling of discrete ordinate methods for solving the homogeneous, isotropic and mono-kinetic equation. Coupling operators are defined and shown to be optimal. The second and the third sections are devoted to an extension of the previous results when the equation is non-homogeneous, anisotropic and multigroup (under some restrictive assumptions). Some numerical results are given in the case of isotropic and mono-kinetic equations. (author) 15 refs.

Solving Fully Fuzzy Linear System of Equations in General Form

Directory of Open Access Journals (Sweden)

A. Yousefzadeh

2012-06-01

Full Text Available In this work, we propose an approach for computing the positive solution of a fully fuzzy linear system where the coefficient matrix is a fuzzy $nimes n$ matrix. To do this, we use arithmetic operations on fuzzy numbers that introduced by Kaffman in and convert the fully fuzzy linear system into two $nimes n$ and $2nimes 2n$ crisp linear systems. If the solutions of these linear systems don't satisfy in positive fuzzy solution condition, we introduce the constrained least squares problem to obtain optimal fuzzy vector solution by applying the ranking function in given fully fuzzy linear system. Using our proposed method, the fully fuzzy linear system of equations always has a solution. Finally, we illustrate the efficiency of proposed method by solving some numerical examples.

Moment equation approach to neoclassical transport theory

International Nuclear Information System (INIS)

Hirshman, S.P.

1978-01-01

The neoclassical cross-field fluxes for a toroidally confined, axisymmetric plasma are calculated in terms of the thermodynamic forces from the fluid continuity and momentum balance equations. This macroscopic formulation of neoclassical transport theory unifies the numerous complex expressions for the transport coefficients, previously obtained by solving the Fokker--Planck equation, and elucidates their physical basis. In the large aspect ratio limit, the continuous transition in the scaling of the diffusion coefficient throughout various collisionality regimes is shown to depend on the ratio of parallel viscosity coefficients of the plasma species. Comparison of the present results with the kinetic theory expressions for the neoclassical fluxes determines the parallel viscosity coefficients for a multispecies plasma in the long-mean-free-path regime

Local linearization methods for the numerical integration of ordinary differential equations: An overview

International Nuclear Information System (INIS)

Jimenez, J.C.

2009-06-01

Local Linearization (LL) methods conform a class of one-step explicit integrators for ODEs derived from the following primary and common strategy: the vector field of the differential equation is locally (piecewise) approximated through a first-order Taylor expansion at each time step, thus obtaining successive linear equations that are explicitly integrated. Hereafter, the LL approach may include some additional strategies to improve that basic affine approximation. Theoretical and practical results have shown that the LL integrators have a number of convenient properties. These include arbitrary order of convergence, A-stability, linearization preserving, regularity under quite general conditions, preservation of the dynamics of the exact solution around hyperbolic equilibrium points and periodic orbits, integration of stiff and high-dimensional equations, low computational cost, and others. In this paper, a review of the LL methods and their properties is presented. (author)

On the Linearized Darboux Equation Arising in Isometric Embedding of the Alexandrov Positive Annulus

Institute of Scientific and Technical Information of China (English)

Chunhe LI

2013-01-01

In the present paper,the solvability condition of the linearized Gauss-Codazzi system and the solutions to the homogenous system are given.In the meantime,the Solvability of a relevant linearized Darboux equation is given.The equations are arising in a geometric problem which is concerned with the realization of the Alexandrov's positive annulus in R3.

Variational formulation and projectional methods for the second order transport equation

International Nuclear Information System (INIS)

Borysiewicz, M.; Stankiewicz, R.

1979-01-01

Herein the variational problem for a second-order boundary value problem for the neutron transport equation is formulated. The projectional methods solving the problem are examined. The approach is compared with that based on the original untransformed form of the neutron transport equation

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

CERN Document Server

Rosinger, Elemer E

1990-01-01

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

Anti-symmetrically fused model and non-linear integral equations in the three-state Uimin-Sutherland model

International Nuclear Information System (INIS)

Fujii, Akira; Kluemper, Andreas

1999-01-01

We derive the non-linear integral equations determining the free energy of the three-state pure bosonic Uimin-Sutherland model. In order to find a complete set of auxiliary functions, the anti-symmetric fusion procedure is utilized. We solve the non-linear integral equations numerically and see that the low-temperature behavior coincides with that predicted by conformal field theory. The magnetization and magnetic susceptibility are also calculated by means of the non-linear integral equation

Asymptotic formulae for solutions of half-linear differential equations

Czech Academy of Sciences Publication Activity Database

Řehák, Pavel

2017-01-01

Roč. 292, January (2017), s. 165-177 ISSN 0096-3003 Institutional support: RVO:67985840 Keywords : half-linear differential equation * nonoscillatory solution * regular variation Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 1.738, year: 2016 http://www.sciencedirect.com/science/article/pii/S0096300316304581

Calculations of stationary solutions for the non linear viscous resistive MHD equations in slab geometry

International Nuclear Information System (INIS)

Edery, D.

1983-11-01

The reduced system of the non linear resistive MHD equations is used in the 2-D one helicity approximation in the numerical computations of stationary tearing modes. The critical magnetic Raynolds number S (S=tausub(r)/tausub(H) where tausub(R) and tausub(H) are respectively the characteristic resistive and hydro magnetic times) and the corresponding linear solution are computed as a starting approximation for the full non linear equations. These equations are then treated numerically by an iterative procedure which is shown to be rapidly convergent. A numerical application is given in the last part of this paper

Quasi-linear landau kinetic equations for magnetized plasmas: compact propagator formalism, rotation matrices and interaction

International Nuclear Information System (INIS)

Misguich, J.H.

2004-04-01

As a first step toward a nonlinear renormalized description of turbulence phenomena in magnetized plasmas, the lowest order quasi-linear description is presented here from a unified point of view for collisionless as well as for collisional plasmas in a constant magnetic field. The quasi-linear approximation is applied to a general kinetic equation obtained previously from the Klimontovich exact equation, by means of a generalised Dupree-Weinstock method. The so-obtained quasi-linear description of electromagnetic turbulence in a magnetoplasma is applied to three separate physical cases: -) weak electrostatic turbulence, -) purely magnetic field fluctuations (the classical quasi-linear results are obtained for cosmic ray diffusion in the 'slab model' of magnetostatic turbulence in the solar wind), and -) collisional kinetic equations of magnetized plasmas. This mathematical technique has allowed us to derive basic kinetic equations for turbulent plasmas and collisional plasmas, respectively in the quasi-linear and Landau approximation. In presence of a magnetic field we have shown that the systematic use of rotation matrices describing the helical particle motion allows for a much more compact derivation than usually performed. Moreover, from the formal analogy between turbulent and collisional plasmas, the results derived here in detail for the turbulent plasmas, can be immediately translated to obtain explicit results for the Landau kinetic equation

Quasi-linear landau kinetic equations for magnetized plasmas: compact propagator formalism, rotation matrices and interaction

Energy Technology Data Exchange (ETDEWEB)

Misguich, J.H

2004-04-01

As a first step toward a nonlinear renormalized description of turbulence phenomena in magnetized plasmas, the lowest order quasi-linear description is presented here from a unified point of view for collisionless as well as for collisional plasmas in a constant magnetic field. The quasi-linear approximation is applied to a general kinetic equation obtained previously from the Klimontovich exact equation, by means of a generalised Dupree-Weinstock method. The so-obtained quasi-linear description of electromagnetic turbulence in a magnetoplasma is applied to three separate physical cases: -) weak electrostatic turbulence, -) purely magnetic field fluctuations (the classical quasi-linear results are obtained for cosmic ray diffusion in the 'slab model' of magnetostatic turbulence in the solar wind), and -) collisional kinetic equations of magnetized plasmas. This mathematical technique has allowed us to derive basic kinetic equations for turbulent plasmas and collisional plasmas, respectively in the quasi-linear and Landau approximation. In presence of a magnetic field we have shown that the systematic use of rotation matrices describing the helical particle motion allows for a much more compact derivation than usually performed. Moreover, from the formal analogy between turbulent and collisional plasmas, the results derived here in detail for the turbulent plasmas, can be immediately translated to obtain explicit results for the Landau kinetic equation.

Numerical solution of the transport equation describing the radon transport from subsurface soil to buildings

International Nuclear Information System (INIS)

Savovic, S.; Djordjevich, A.; Ristic, G.

2012-01-01

A theoretical evaluation of the properties and processes affecting the radon transport from subsurface soil into buildings is presented in this work. The solution of the relevant transport equation is obtained using the explicit finite difference method (EFDM). Results are compared with analytical steady-state solution reported in the literature. Good agreement is found. It is shown that EFDM is effective and accurate for solving the equation that describes radon diffusion, advection and decay during its transport from subsurface to buildings, which is especially important when arbitrary initial and boundary conditions are required. (authors)

A symbiotic approach to fluid equations and non-linear flux-driven simulations of plasma dynamics

Science.gov (United States)

Halpern, Federico

2017-10-01

The fluid framework is ubiquitous in studies of plasma transport and stability. Typical forms of the fluid equations are motivated by analytical work dating several decades ago, before computer simulations were indispensable, and can be, therefore, not optimal for numerical computation. We demonstrate a new first-principles approach to obtaining manifestly consistent, skew-symmetric fluid models, ensuring internal consistency and conservation properties even in discrete form. Mass, kinetic, and internal energy become quadratic (and always positive) invariants of the system. The model lends itself to a robust, straightforward discretization scheme with inherent non-linear stability. A simpler, drift-ordered form of the equations is obtained, and first results of their numerical implementation as a binary framework for bulk-fluid global plasma simulations are demonstrated. This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, Theory Program, under Award No. DE-FG02-95ER54309.

Diffusion phenomenon for linear dissipative wave equations in an exterior domain

Science.gov (United States)

Ikehata, Ryo

Under the general condition of the initial data, we will derive the crucial estimates which imply the diffusion phenomenon for the dissipative linear wave equations in an exterior domain. In order to derive the diffusion phenomenon for dissipative wave equations, the time integral method which was developed by Ikehata and Matsuyama (Sci. Math. Japon. 55 (2002) 33) plays an effective role.

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)

Exponential estimates for solutions of half-linear differential equations

Czech Academy of Sciences Publication Activity Database

Řehák, Pavel

2015-01-01

Roč. 147, č. 1 (2015), s. 158-171 ISSN 0236-5294 Institutional support: RVO:67985840 Keywords : half-linear differential equation * decreasing solution * increasing solution * asymptotic behavior Subject RIV: BA - General Mathematics Impact factor: 0.469, year: 2015 http://link.springer.com/article/10.1007%2Fs10474-015-0522-9

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)

A non-Linear transport model for determining shale rock characteristics

Science.gov (United States)

Ali, Iftikhar; Malik, Nadeem

2016-04-01

Unconventional hydrocarbon reservoirs consist of tight porous rocks which are characterised by nano-scale size porous networks with ultra-low permeability [1,2]. Transport of gas through them is not well understood at the present time, and realistic transport models are needed in order to determine rock properties and for estimating future gas pressure distribution in the reservoirs. Here, we consider a recently developed non-linear gas transport equation [3], ∂p-+ U ∂p- = D ∂2p-, t > 0, (1) ∂t ∂x ∂x2 complimented with suitable initial and boundary conditions, in order to determine shale rock properties such as the permeability K, the porosity φ and the tortuosity, τ. In our new model, the apparent convection velocity, U = U(p,px), and the apparent diffusivity D = D(p), are both highly non-linear functions of the pressure. The model incorporate various flow regimes (slip, surface diffusion, transition, continuum) based upon the Knudsen number Kn, and also includes Forchchiemers turbulence correction terms. In application, the model parameters and associated compressibility factors are fully pressure dependent, giving the model more realism than previous models. See [4]. Rock properties are determined by solving an inverse problem, with model parameters adjustment to minimise the error between the model simulation and available data. It is has been found that the proposed model performs better than previous models. Results and details of the model will be presented at the conference. Corresponding author: namalik@kfupm.edu.sa and nadeem_malik@cantab.net References [1] Cui, X., Bustin, A.M. and Bustin, R., "Measurements of gas permeability and diffusivity of tight reservoir rocks: different approaches and their applications", Geofluids 9, 208-223 (2009). [2] Chiba R., Fomin S., Chugunov V., Niibori Y. and Hashida T., "Numerical Simulation of Non Fickian Diffusion and Advection in a Fractured Porous Aquifer", AIP Conference Proceedings 898, 75 (2007

On the equivalence between particular types of Navier-Stokes and non-linear Schroedinger equations

International Nuclear Information System (INIS)

Dietrich, K.; Vautherin, D.

1985-01-01

We derive a Schroedinger equation equivalent to the Navier-Stokes equation in the special case of constant kinematic viscosities. This equation contains a non-linear term similar to that proposed by Kostin for a quantum description of friction [fr

SHADOK-3-6, Transport Equation with Anisotropic Diffusion in P1 Approximation for Spherical and Cylindrical Geometry

International Nuclear Information System (INIS)

Ligou, J.; Thomi, P.A.

1973-01-01

1 - Nature of physical problem solved: Integral transport equation, anisotropy of diffusion in P1 approximation. SHADOK3 - cylindrical geometry; direct solution of the linear system. SHADOK4 - cylindrical geometry; Thermalization iteration; solution of the linear system with inverse matrix calculation. SHADOK5 - like SHADOK3 for spherical geometry. SHADOK6 - like SHADOK4 for spherical geometry. 2 - Method of solution: Analysis in terms of annuli for each of which polynomial approximation is applied. Dynamic allocation (for formulas see report TM(10)). 3 - Restrictions on the complexity of the problem: Relative accuracy of the Bickley functions about 1.0E-13

Discontinuous Galerkin for the Radiative Transport Equation

KAUST Repository

Guermond, Jean-Luc

2013-10-11

This note presents some recent results regarding the approximation of the linear radiative transfer equation using discontinuous Galerkin methods. The locking effect occurring in the diffusion limit with the upwind numerical flux is investigated and a correction technique is proposed.

Discontinuous Galerkin for the Radiative Transport Equation

KAUST Repository

Guermond, Jean-Luc; Kanschat, Guido; Ragusa, Jean C.

2013-01-01

This note presents some recent results regarding the approximation of the linear radiative transfer equation using discontinuous Galerkin methods. The locking effect occurring in the diffusion limit with the upwind numerical flux is investigated and a correction technique is proposed.

Shifted Legendre method with residual error estimation for delay linear Fredholm integro-differential equations

Directory of Open Access Journals (Sweden)

Şuayip Yüzbaşı

2017-03-01

Full Text Available In this paper, we suggest a matrix method for obtaining the approximate solutions of the delay linear Fredholm integro-differential equations with constant coefficients using the shifted Legendre polynomials. The problem is considered with mixed conditions. Using the required matrix operations, the delay linear Fredholm integro-differential equation is transformed into a matrix equation. Additionally, error analysis for the method is presented using the residual function. Illustrative examples are given to demonstrate the efficiency of the method. The results obtained in this study are compared with the known results.

An Implementation of Interfacial Transport Equation into the CUPID code

Energy Technology Data Exchange (ETDEWEB)

Park, Ik Kyu; Cho, Heong Kyu; Yoon, Han Young; Jeong, Jae Jun

2009-11-15

A component scale thermal hydraulic analysis code, CUPID (Component Unstructured Program for Interfacial Dynamics), is being developed for the analysis of components for a nuclear reactor, such as reactor vessel, steam generator, containment, etc. It adopted a three-dimensional, transient, two phase and three-field model. In order to develop the numerical schemes for the three-field model, various numerical schemes have been examined including the SMAS, semi-implicit ICE, SIMPLE. The governing equations for a 2-phase flow are composed of mass, momentum, and energy conservation equations for each phase. These equation sets are closed by the interfacial transfer rate of mass, momentum, and energy. The interfacial transfer of mass, momentum, and energy occurs through the interfacial area, and this area plays an important role in the transfer rate. The flow regime based correlations are used for calculating the interracial area in the traditional style 2-phase flow model. This is dependent upon the flow regime and is limited to the fully developed 2-phase flow region. Its application to the multi-dimensional 2-phase flow has some limitation because it adopts the measured results of 2-phase flow in the 1-dimensional tube. The interfacial area concentration transport equation had been suggested in order to calculate the interfacial area without the interfacial area correlations. The source terms to close the interfacial area transport equation should be further developed for a wide ranger usage of it. In this study, the one group interfacial area concentration transport equation has been implemented into the CUPID code. This interfacial area concentration transport equation can be used instead of the interfacial area concentration correlations for the bubbly flow region.

An Implementation of Interfacial Transport Equation into the CUPID code

International Nuclear Information System (INIS)

Park, Ik Kyu; Cho, Heong Kyu; Yoon, Han Young; Jeong, Jae Jun

2009-11-01

A component scale thermal hydraulic analysis code, CUPID (Component Unstructured Program for Interfacial Dynamics), is being developed for the analysis of components for a nuclear reactor, such as reactor vessel, steam generator, containment, etc. It adopted a three-dimensional, transient, two phase and three-field model. In order to develop the numerical schemes for the three-field model, various numerical schemes have been examined including the SMAS, semi-implicit ICE, SIMPLE. The governing equations for a 2-phase flow are composed of mass, momentum, and energy conservation equations for each phase. These equation sets are closed by the interfacial transfer rate of mass, momentum, and energy. The interfacial transfer of mass, momentum, and energy occurs through the interfacial area, and this area plays an important role in the transfer rate. The flow regime based correlations are used for calculating the interracial area in the traditional style 2-phase flow model. This is dependent upon the flow regime and is limited to the fully developed 2-phase flow region. Its application to the multi-dimensional 2-phase flow has some limitation because it adopts the measured results of 2-phase flow in the 1-dimensional tube. The interfacial area concentration transport equation had been suggested in order to calculate the interfacial area without the interfacial area correlations. The source terms to close the interfacial area transport equation should be further developed for a wide ranger usage of it. In this study, the one group interfacial area concentration transport equation has been implemented into the CUPID code. This interfacial area concentration transport equation can be used instead of the interfacial area concentration correlations for the bubbly flow region

Asymptotic behavior of solutions of linear multi-order fractional differential equation systems

OpenAIRE

Diethelm, Kai; Siegmund, Stefan; Tuan, H. T.

2017-01-01

In this paper, we investigate some aspects of the qualitative theory for multi-order fractional differential equation systems. First, we obtain a fundamental result on the existence and uniqueness for multi-order fractional differential equation systems. Next, a representation of solutions of homogeneous linear multi-order fractional differential equation systems in series form is provided. Finally, we give characteristics regarding the asymptotic behavior of solutions to some classes of line...

Non-linear analysis of wave progagation using transform methods and plates and shells using integral equations

Science.gov (United States)

Pipkins, Daniel Scott

Two diverse topics of relevance in modern computational mechanics are treated. The first involves the modeling of linear and non-linear wave propagation in flexible, lattice structures. The technique used combines the Laplace Transform with the Finite Element Method (FEM). The procedure is to transform the governing differential equations and boundary conditions into the transform domain where the FEM formulation is carried out. For linear problems, the transformed differential equations can be solved exactly, hence the method is exact. As a result, each member of the lattice structure is modeled using only one element. In the non-linear problem, the method is no longer exact. The approximation introduced is a spatial discretization of the transformed non-linear terms. The non-linear terms are represented in the transform domain by making use of the complex convolution theorem. A weak formulation of the resulting transformed non-linear equations yields a set of element level matrix equations. The trial and test functions used in the weak formulation correspond to the exact solution of the linear part of the transformed governing differential equation. Numerical results are presented for both linear and non-linear systems. The linear systems modeled are longitudinal and torsional rods and Bernoulli-Euler and Timoshenko beams. For non-linear systems, a viscoelastic rod and Von Karman type beam are modeled. The second topic is the analysis of plates and shallow shells under-going finite deflections by the Field/Boundary Element Method. Numerical results are presented for two plate problems. The first is the bifurcation problem associated with a square plate having free boundaries which is loaded by four, self equilibrating corner forces. The results are compared to two existing numerical solutions of the problem which differ substantially. linear model are compared to those

Modified two-fluid model for the two-group interfacial area transport equation

International Nuclear Information System (INIS)

Sun Xiaodong; Ishii, Mamoru; Kelly, Joseph M.

2003-01-01

This paper presents a modified two-fluid model that is ready to be applied in the approach of the two-group interfacial area transport equation. The two-group interfacial area transport equation was developed to provide a mechanistic constitutive relation for the interfacial area concentration in the two-fluid model. In the two-group transport equation, bubbles are categorized into two groups: spherical/distorted bubbles as Group 1 while cap/slug/churn-turbulent bubbles as Group 2. Therefore, this transport equation can be employed in the flow regimes spanning from bubbly, cap bubbly, slug to churn-turbulent flows. However, the introduction of the two groups of bubbles requires two gas velocity fields. Yet it is not practical to solve two momentum equations for the gas phase alone. In the current modified two-fluid model, a simplified approach is proposed. The momentum equation for the averaged velocity of both Group-1 and Group-2 bubbles is retained. By doing so, the velocity difference between Group-1 and Group-2 bubbles needs to be determined. This may be made either based on simplified momentum equations for both Group-1 and Group-2 bubbles or by a modified drift-flux model

Students' errors in solving linear equation word problems: Case ...

African Journals Online (AJOL)

kofi.mereku

Development in most areas of life is based on effective knowledge of science and ... Problem solving, as used in mathematics education literature, refers ... word problems, on the other hand, are those linear equation tasks or ... taught LEWPs in the junior high school, many of them reach the senior high school without a.

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)

On nonnegative solutions of second order linear functional differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander; Vodstrčil, Petr

2004-01-01

Roč. 32, č. 1 (2004), s. 59-88 ISSN 1512-0015 Institutional research plan: CEZ:AV0Z1019905 Keywords : second order linear functional differential equations * nonnegative solution * two-point boundary value problem Subject RIV: BA - General Mathematics

A Hamiltonian functional for the linearized Einstein vacuum field equations

International Nuclear Information System (INIS)

Rosas-RodrIguez, R

2005-01-01

By considering the Einstein vacuum field equations linearized about the Minkowski metric, the evolution equations for the gauge-invariant quantities characterizing the gravitational field are written in a Hamiltonian form by using a conserved functional as Hamiltonian; this Hamiltonian is not the analog of the energy of the field. A Poisson bracket between functionals of the field, compatible with the constraints satisfied by the field variables, is obtained. The generator of spatial translations associated with such bracket is also obtained

New exact solutions of the Tzitzéica-type equations in non-linear optics using the expa function method

Science.gov (United States)

Hosseini, K.; Ayati, Z.; Ansari, R.

2018-04-01

One specific class of non-linear evolution equations, known as the Tzitzéica-type equations, has received great attention from a group of researchers involved in non-linear science. In this article, new exact solutions of the Tzitzéica-type equations arising in non-linear optics, including the Tzitzéica, Dodd-Bullough-Mikhailov and Tzitzéica-Dodd-Bullough equations, are obtained using the expa function method. The integration technique actually suggests a useful and reliable method to extract new exact solutions of a wide range of non-linear evolution equations.

An introduction to linear ordinary differential equations using the impulsive response method and factorization

CERN Document Server

Camporesi, Roberto

2016-01-01

This book presents a method for solving linear ordinary differential equations based on the factorization of the differential operator. The approach for the case of constant coefficients is elementary, and only requires a basic knowledge of calculus and linear algebra. In particular, the book avoids the use of distribution theory, as well as the other more advanced approaches: Laplace transform, linear systems, the general theory of linear equations with variable coefficients and variation of parameters. The case of variable coefficients is addressed using Mammana’s result for the factorization of a real linear ordinary differential operator into a product of first-order (complex) factors, as well as a recent generalization of this result to the case of complex-valued coefficients.

Systems of evolution equations and the singular perturbation method

International Nuclear Information System (INIS)

Mika, J.

Several fundamental theorems are presented important for the solution of linear evolution equations in the Banach space. The algorithm is deduced extending the solution of the system of singularly perturbed evolution equations into an asymptotic series with respect to a small positive parameter. The asymptotic convergence is shown of an approximate solution to the accurate solution. Singularly perturbed evolution equations of the resonance type were analysed. The special role is considered of the asymptotic equivalence of P1 equations obtained as the first order approximation if the spherical harmonics method is applied to the linear Boltzmann equation, and the diffusion equations of the linear transport theory where the small parameter approaches zero. (J.B.)

Fuzzy linear programming approach for solving transportation

Indian Academy of Sciences (India)

Transportation problem (TP) is an important network structured linear programming problem that arises in several contexts and has deservedly received a great deal of attention in the literature. The central concept in this problem is to find the least total transportation cost of a commodity in order to satisfy demands at ...

Some applications of linear difference equations in finance with wolfram|alpha and maple

Directory of Open Access Journals (Sweden)

Dana Rıhová

2014-12-01

Full Text Available The principle objective of this paper is to show how linear difference equations can be applied to solve some issues of financial mathematics. We focus on the area of compound interest and annuities. In both cases we determine appropriate recursive rules, which constitute the first order linear difference equations with constant coefficients, and derive formulas required for calculating examples. Finally, we present possibilities of application of two selected computer algebra systems Wolfram|Alpha and Maple in this mathematical area.

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)

Unsteady Solution of Non-Linear Differential Equations Using Walsh Function Series

Science.gov (United States)

Gnoffo, Peter A.

2015-01-01

Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.

Identification of a time-varying point source in a system of two coupled linear diffusion-advection- reaction equations: application to surface water pollution

International Nuclear Information System (INIS)

Hamdi, Adel

2009-01-01

This paper deals with the identification of a point source (localization of its position and recovering the history of its time-varying intensity function) that constitutes the right-hand side of the first equation in a system of two coupled 1D linear transport equations. Assuming that the source intensity function vanishes before reaching the final control time, we prove the identifiability of the sought point source from recording the state relative to the second coupled transport equation at two observation points framing the source region. Note that at least one of the two observation points should be strategic. We establish an identification method that uses these records to identify the source position as the root of a continuous and strictly monotonic function. Whereas the source intensity function is recovered using a recursive formula without any need of an iterative process. Some numerical experiments on a variant of the surface water pollution BOD–OD coupled model are presented

The transport equation in general geometry

International Nuclear Information System (INIS)

Pomraning, G.C.

1990-01-01

As stated in the introduction to the paper, the motivation for this work was to obtain an explicit form for the streaming operator in the transport equation, which could be used to compute curvature effects in an asymptotic analysis leading to diffusion theory. This sign error was discovered while performing this analysis

Differential equation of exospheric lateral transport and its application to terrestrial hydrogen

Science.gov (United States)

Hodges, R. R., Jr.

1973-01-01

The differential equation description of exospheric lateral transport of Hodges and Johnson is reformulated to extend its utility to light gases. Accuracy of the revised equation is established by applying it to terrestrial hydrogen. The resulting global distributions for several static exobase models are shown to be essentially the same as those that have been computed by Quessette using an integral equation approach. The present theory is subsequently used to elucidate the effects of nonzero lateral flow, exobase rotation, and diurnal tidal winds on the hydrogen distribution. Finally it is shown that the differential equation of exospheric transport is analogous to a diffusion equation. Hence it is practical to consider exospheric transport as a continuation of thermospheric diffusion, a concept that alleviates the need for an artificial exobase dividing thermosphere and exosphere.

Nonlinear and linear wave equations for propagation in media with frequency power law losses

Science.gov (United States)

Szabo, Thomas L.

2003-10-01

The Burgers, KZK, and Westervelt wave equations used for simulating wave propagation in nonlinear media are based on absorption that has a quadratic dependence on frequency. Unfortunately, most lossy media, such as tissue, follow a more general frequency power law. The authors first research involved measurements of loss and dispersion associated with a modification to Blackstock's solution to the linear thermoviscous wave equation [J. Acoust. Soc. Am. 41, 1312 (1967)]. A second paper by Blackstock [J. Acoust. Soc. Am. 77, 2050 (1985)] showed the loss term in the Burgers equation for plane waves could be modified for other known instances of loss. The authors' work eventually led to comprehensive time-domain convolutional operators that accounted for both dispersion and general frequency power law absorption [Szabo, J. Acoust. Soc. Am. 96, 491 (1994)]. Versions of appropriate loss terms were developed to extend the standard three nonlinear wave equations to these more general losses. Extensive experimental data has verified the predicted phase velocity dispersion for different power exponents for the linear case. Other groups are now working on methods suitable for solving wave equations numerically for these types of loss directly in the time domain for both linear and nonlinear media.

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)

Mathematical problems in the one-velocity theory of particle transport

Energy Technology Data Exchange (ETDEWEB)

Vladimirov, V S

1963-01-15

This paper describes kinetic (transport) equations which describe the process of neutron transport in a substance. These equations are linear, integro-differential equations in partial derivatives of first order.

Mathematical problems in the one-velocity theory of particle transport

International Nuclear Information System (INIS)

Vladimirov, V.S.

1963-01-01

This paper describes kinetic (transport) equations which describe the process of neutron transport in a substance. These equations are linear, integro-differential equations in partial derivatives of first order.

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.

Analysis of a block Gauss-Seidel iterative method for a finite element discretization of the neutron transport equation

International Nuclear Information System (INIS)

Lorence, L.J. Jr.; Martin, W.R.; Luskin, M.

1985-01-01

We prove the convergence of a finite element discretization of the neutron transport equation. The iterative solution of the resulting linear system by a block Gauss-Seidel method is also analyzed. This procedure is shown to require less storage than the direct solution by Gaussian elimination, and an estimate for the rate of convergence is used to show that fewer arithmetic operations are required

Non-monotone positive solutions of second-order linear differential equations: existence, nonexistence and criteria

Directory of Open Access Journals (Sweden)

Mervan Pašić

2016-10-01

Full Text Available We study non-monotone positive solutions of the second-order linear differential equations: $(p(tx'' + q(t x = e(t$, with positive $p(t$ and $q(t$. For the first time, some criteria as well as the existence and nonexistence of non-monotone positive solutions are proved in the framework of some properties of solutions $\\theta (t$ of the corresponding integrable linear equation: $(p(t\\theta''=e(t$. The main results are illustrated by many examples dealing with equations which allow exact non-monotone positive solutions not necessarily periodic. Finally, we pose some open questions.

Quasi-linear equation for magnetoplasma oscillations in the weakly relativistic approximation

International Nuclear Information System (INIS)

Rizzato, F.B.

1985-01-01

Some limitations which are present in the dynamical equations for collisionless plasmas are discussed. Some elementary corrections to the linear theories are obtained in a heuristic form, which directly lead to the so-called quasi-linear theories in its non-relativistic and relativistic forms. The effect of the relativistic variation of the gyrofrequency on the diffusion coefficient is examined in a typically perturbative approximation. (author)

Hyers-Ulam stability of linear second-order differential equations in complex Banach spaces

Directory of Open Access Journals (Sweden)

Yongjin Li

2013-08-01

Full Text Available We prove the Hyers-Ulam stability of linear second-order differential equations in complex Banach spaces. That is, if y is an approximate solution of the differential equation $y''+ alpha y'(t +eta y = 0$ or $y''+ alpha y'(t +eta y = f(t$, then there exists an exact solution of the differential equation near to y.

A three operator split-step method covering a larger set of non-linear partial differential equations

Science.gov (United States)

Zia, Haider

2017-06-01

This paper describes an updated exponential Fourier based split-step method that can be applied to a greater class of partial differential equations than previous methods would allow. These equations arise in physics and engineering, a notable example being the generalized derivative non-linear Schrödinger equation that arises in non-linear optics with self-steepening terms. These differential equations feature terms that were previously inaccessible to model accurately with low computational resources. The new method maintains a 3rd order error even with these additional terms and models the equation in all three spatial dimensions and time. The class of non-linear differential equations that this method applies to is shown. The method is fully derived and implementation of the method in the split-step architecture is shown. This paper lays the mathematical ground work for an upcoming paper employing this method in white-light generation simulations in bulk material.

Solution of second order linear fuzzy difference equation by Lagrange's multiplier method

Directory of Open Access Journals (Sweden)

Sankar Prasad Mondal

2016-06-01

Full Text Available In this paper we execute the solution procedure for second order linear fuzzy difference equation by Lagrange's multiplier method. In crisp sense the difference equation are easy to solve, but when we take in fuzzy sense it forms a system of difference equation which is not so easy to solve. By the help of Lagrange's multiplier we can solved it easily. The results are illustrated by two different numerical examples and followed by two applications.

Self-Adjoint Angular Flux Equation for Coupled Electron-Photon Transport

International Nuclear Information System (INIS)

Liscum-Powell, J.L.; Lorence, L.J. Jr.; Morel, J.E.; Prinja, A.K.

1999-01-01

Recently, Morel and McGhee described an alternate second-order form of the transport equation called the self adjoint angular flux (SAAF) equation that has the angular flux as its unknown. The SAAF formulation has all the advantages of the traditional even- and odd-parity self-adjoint equations, with the added advantages that it yields the full angular flux when it is numerically solved, it is significantly easier to implement reflective and reflective-like boundary conditions, and in the appropriate form it can be solved in void regions. The SAAF equation has the disadvantage that the angular domain is the full unit sphere and, like the even- and odd- parity form, S n source iteration cannot be implemented using the standard sweeping algorithm. Also, problems arise in pure scattering media. Morel and McGhee demonstrated the efficacy of the SAAF formulation for neutral particle transport. Here we apply the SAAF formulation to coupled electron-photon transport problems using multigroup cross-sections from the CEPXS code and S n discretization

Self-adjoint angular flux equation for coupled electron-photon transport

International Nuclear Information System (INIS)

Liscum-Powell, J.L.; Prinja, A.K.; Morel, J.E.; Lorence, L.J. Jr.

1999-01-01

Recently, Morel and McGhee described an alternate second-order form of the transport equation called the self-adjoint angular flux (SAAF) equation that has the angular flux as its unknown. The SAAF formulation has all the advantages of the traditional even- and odd-parity self-adjoint equations, with the added advantages that it yields the full angular flux when it is numerically solved, it is significantly easier to implement reflective and reflective-like boundary conditions, and in the appropriate form it can be solved in void regions. The SAAF equation has the disadvantage that the angular domain is the full unit sphere, and, like the even- and odd-parity form, S n source iteration cannot be implemented using the standard sweeping algorithm. Also, problems arise in pure scattering media. Morel and McGhee demonstrated the efficacy of the SAAF formulation for neutral particle transport. Here, the authors apply the SAAF formulation to coupled electron-photon transport problems using multigroup cross sections from the CEPXS code and S n discretization

On a class of strongly degenerate and singular linear elliptic equation

International Nuclear Information System (INIS)

Duong Minh Duc, D.M.; Le Dung.

1992-11-01

We consider a class of strongly degenerate linear elliptic equation. The boundedness and the Holder regularity of the weak solutions in the weighted Sobolev-Hardy spaces will be studied. (author). 9 refs

A study on linear and nonlinear Schrodinger equations by the variational iteration method

International Nuclear Information System (INIS)

Wazwaz, Abdul-Majid

2008-01-01

In this work, we introduce a framework to obtain exact solutions to linear and nonlinear Schrodinger equations. The He's variational iteration method (VIM) is used for analytic treatment of these equations. Numerical examples are tested to show the pertinent features of this method

Bounded solutions of self-adjoint second order linear difference equations with periodic coeffients

Directory of Open Access Journals (Sweden)

Encinas A.M.

2018-02-01

Full Text Available In this work we obtain easy characterizations for the boundedness of the solutions of the discrete, self–adjoint, second order and linear unidimensional equations with periodic coefficients, including the analysis of the so-called discrete Mathieu equations as particular cases.

Coupled energy-drift and force-balance equations for high-field hot-carrier transport

International Nuclear Information System (INIS)

Huang, Danhong; Alsing, P.M.; Apostolova, T.; Cardimona, D.A.

2005-01-01

Coupled energy-drift and force-balance equations that contain a frictional force for the center-of-mass motion of electrons are derived for hot-electron transport under a strong dc electric field. The frictional force is found to be related to the net rate of phonon emission, which takes away the momentum of a phonon from an electron during each phonon-emission event. The net rate of phonon emission is determined by the Boltzmann scattering equation, which depends on the distribution of electrons interacting with phonons. The work done by the frictional force is included into the energy-drift equation for the electron-relative scattering motion and is found to increase the thermal energy of the electrons. The importance of the hot-electron effect in the energy-drift term under a strong dc field is demonstrated in reducing the field-dependent drift velocity and mobility. The Doppler shift in the energy conservation of scattering electrons interacting with impurities and phonons is found to lead to an anisotropic distribution of electrons in the momentum space along the field direction. The importance of this anisotropic distribution is demonstrated through a comparison with the isotropic energy-balance equation, from which we find that defining a state-independent electron temperature becomes impossible. To the leading order, the energy-drift equation is linearized with a distribution function by expanding it into a Fokker-Planck-type equation, along with the expansions of both the force-balance equation and the Boltzmann scattering equation for hot phonons

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.

Linear triangle finite element formulation for multigroup neutron transport analysis with anisotropic scattering

Energy Technology Data Exchange (ETDEWEB)

Lillie, R.A.; Robinson, J.C.

1976-05-01

The discrete ordinates method is the most powerful and generally used deterministic method to obtain approximate solutions of the Boltzmann transport equation. A finite element formulation, utilizing a canonical form of the transport equation, is here developed to obtain both integral and pointwise solutions to neutron transport problems. The formulation is based on the use of linear triangles. A general treatment of anisotropic scattering is included by employing discrete ordinates-like approximations. In addition, multigroup source outer iteration techniques are employed to perform group-dependent calculations. The ability of the formulation to reduce substantially ray effects and its ability to perform streaming calculations are demonstrated by analyzing a series of test problems. The anisotropic scattering and multigroup treatments used in the development of the formulation are verified by a number of one-dimensional comparisons. These comparisons also demonstrate the relative accuracy of the formulation in predicting integral parameters. The applicability of the formulation to nonorthogonal planar geometries is demonstrated by analyzing a hexagonal-type lattice. A small, high-leakage reactor model is analyzed to investigate the effects of varying both the spatial mesh and order of angular quadrature. This analysis reveals that these effects are more pronounced in the present formulation than in other conventional formulations. However, the insignificance of these effects is demonstrated by analyzing a realistic reactor configuration. In addition, this final analysis illustrates the importance of incorporating anisotropic scattering into the finite element formulation. 8 tables, 29 figures.

Linear triangle finite element formulation for multigroup neutron transport analysis with anisotropic scattering

International Nuclear Information System (INIS)

Lillie, R.A.; Robinson, J.C.

1976-05-01

The discrete ordinates method is the most powerful and generally used deterministic method to obtain approximate solutions of the Boltzmann transport equation. A finite element formulation, utilizing a canonical form of the transport equation, is here developed to obtain both integral and pointwise solutions to neutron transport problems. The formulation is based on the use of linear triangles. A general treatment of anisotropic scattering is included by employing discrete ordinates-like approximations. In addition, multigroup source outer iteration techniques are employed to perform group-dependent calculations. The ability of the formulation to reduce substantially ray effects and its ability to perform streaming calculations are demonstrated by analyzing a series of test problems. The anisotropic scattering and multigroup treatments used in the development of the formulation are verified by a number of one-dimensional comparisons. These comparisons also demonstrate the relative accuracy of the formulation in predicting integral parameters. The applicability of the formulation to nonorthogonal planar geometries is demonstrated by analyzing a hexagonal-type lattice. A small, high-leakage reactor model is analyzed to investigate the effects of varying both the spatial mesh and order of angular quadrature. This analysis reveals that these effects are more pronounced in the present formulation than in other conventional formulations. However, the insignificance of these effects is demonstrated by analyzing a realistic reactor configuration. In addition, this final analysis illustrates the importance of incorporating anisotropic scattering into the finite element formulation. 8 tables, 29 figures

Thermophoresis of a spherical particle: Modeling through moment-based, macroscopic transport equations

Science.gov (United States)

Padrino, Juan C.; Sprittles, James; Lockerby, Duncan

2017-11-01

Thermophoresis refers to the forces on and motions of objects caused by temperature gradients when these objects are exposed to rarefied gases. This phenomenon can occur when the ratio of the gas mean free path to the characteristic physical length scale (Knudsen number) is not negligible. In this work, we obtain the thermophoretic force on a rigid, heat-conducting spherical particle immersed in a rarefied gas resulting from a uniform temperature gradient imposed far from the sphere. To this end, we model the gas dynamics using the steady, linearized version of the so-called regularized 13-moment equations (R13). This set of equations, derived from the Boltzmann equation using the moment method, provides closures to the mass, momentum, and energy conservation laws in the form of constitutive, transport equations for the stress and heat flux that extends the Navier-Stokes-Fourier model to include rarefaction effects. Integration of the pressure and stress on the surface of the sphere leads to the net force as a function of the Knudsen number, dimensionless temperature gradient, and particle-to-gas thermal conductivity ratio. Results from this expression are compared with predictions from other moment-based models as well as from kinetic models. Supported in the UK by the Engineering and Physical Sciences Research Council (EP/N016602/1).

On the classical theory of ordinary linear differential equations of the second order and the Schroedinger equation for power law potentials

International Nuclear Information System (INIS)

Lima, M.L.; Mignaco, J.A.

1983-01-01

The power law potentials in the Schroedinger equation solved recently are shown to come from the classical treatment of the singularities of a linear, second order differential equation. This allows to enlarge the class of solvable power law potentials. (Author) [pt

Radial solutions to semilinear elliptic equations via linearized operators

Directory of Open Access Journals (Sweden)

Phuong Le

2017-04-01

Full Text Available Let $u$ be a classical solution of semilinear elliptic equations in a ball or an annulus in $\\mathbb{R}^N$ with zero Dirichlet boundary condition where the nonlinearity has a convex first derivative. In this note, we prove that if the $N$-th eigenvalue of the linearized operator at $u$ is positive, then $u$ must be radially symmetric.

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

Non-linear mixed-effects pharmacokinetic/pharmacodynamic modelling in NLME using differential equations

DEFF Research Database (Denmark)

Tornøe, Christoffer Wenzel; Agersø, Henrik; Madsen, Henrik

2004-01-01

The standard software for non-linear mixed-effect analysis of pharmacokinetic/phar-macodynamic (PK/PD) data is NONMEM while the non-linear mixed-effects package NLME is an alternative as tong as the models are fairly simple. We present the nlmeODE package which combines the ordinary differential...... equation (ODE) solver package odesolve and the non-Linear mixed effects package NLME thereby enabling the analysis of complicated systems of ODEs by non-linear mixed-effects modelling. The pharmacokinetics of the anti-asthmatic drug theophylline is used to illustrate the applicability of the nlme...

The intelligence of dual simplex method to solve linear fractional fuzzy transportation problem.

Science.gov (United States)

Narayanamoorthy, S; Kalyani, S

2015-01-01

An approach is presented to solve a fuzzy transportation problem with linear fractional fuzzy objective function. In this proposed approach the fractional fuzzy transportation problem is decomposed into two linear fuzzy transportation problems. The optimal solution of the two linear fuzzy transportations is solved by dual simplex method and the optimal solution of the fractional fuzzy transportation problem is obtained. The proposed method is explained in detail with an example.

The Intelligence of Dual Simplex Method to Solve Linear Fractional Fuzzy Transportation Problem

Directory of Open Access Journals (Sweden)

S. Narayanamoorthy

2015-01-01

Full Text Available An approach is presented to solve a fuzzy transportation problem with linear fractional fuzzy objective function. In this proposed approach the fractional fuzzy transportation problem is decomposed into two linear fuzzy transportation problems. The optimal solution of the two linear fuzzy transportations is solved by dual simplex method and the optimal solution of the fractional fuzzy transportation problem is obtained. The proposed method is explained in detail with an example.

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

International Nuclear Information System (INIS)

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

1998-01-01

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

Improvement of the symbolic Monte-Carlo method for the transport equation: P1 extension and coupling with diffusion

International Nuclear Information System (INIS)

Clouet, J.F.; Samba, G.

2005-01-01

We use asymptotic analysis to study the diffusion limit of the Symbolic Implicit Monte-Carlo (SIMC) method for the transport equation. For standard SIMC with piecewise constant basis functions, we demonstrate mathematically that the solution converges to the solution of a wrong diffusion equation. Nevertheless a simple extension to piecewise linear basis functions enables to obtain the correct solution. This improvement allows the calculation in opaque medium on a mesh resolving the diffusion scale much larger than the transport scale. Anyway, the huge number of particles which is necessary to get a correct answer makes this computation time consuming. Thus, we have derived from this asymptotic study an hybrid method coupling deterministic calculation in the opaque medium and Monte-Carlo calculation in the transparent medium. This method gives exactly the same results as the previous one but at a much lower price. We present numerical examples which illustrate the analysis. (authors)

Variance estimates for transport in stochastic media by means of the master equation

International Nuclear Information System (INIS)

Pautz, S. D.; Franke, B. C.; Prinja, A. K.

2013-01-01

The master equation has been used to examine properties of transport in stochastic media. It has been shown previously that not only may the Levermore-Pomraning (LP) model be derived from the master equation for a description of ensemble-averaged transport quantities, but also that equations describing higher-order statistical moments may be obtained. We examine in greater detail the equations governing the second moments of the distribution of the angular fluxes, from which variances may be computed. We introduce a simple closure for these equations, as well as several models for estimating the variances of derived transport quantities. We revisit previous benchmarks for transport in stochastic media in order to examine the error of these new variance models. We find, not surprisingly, that the errors in these variance estimates are at least as large as the corresponding estimates of the average, and sometimes much larger. We also identify patterns in these variance estimates that may help guide the construction of more accurate models. (authors)

Symmetrized neutron transport equation and the fast Fourier transform method

International Nuclear Information System (INIS)

Sinh, N.Q.; Kisynski, J.; Mika, J.

1978-01-01

The differential equation obtained from the neutron transport equation by the application of the source iteration method in two-dimensional rectangular geometry is transformed into a symmetrized form with respect to one of the angular variables. The discretization of the symmetrized equation leads to finite difference equations based on the five-point scheme and solved by use of the fast Fourier transform method. Possible advantages of the approach are shown on test calculations

Advances in the solution of three-dimensional nodal neutron transport equation

International Nuclear Information System (INIS)

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

2003-01-01

In this paper we study the three-dimensional nodal discrete-ordinates approximations of neutron transport equation in a convex domain with piecewise smooth boundaries. 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 obtaining 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. We give numerical results obtained with an algebraic computer system (for N up to 8) and with a code for higher values of N. We compare our results for the geometry of a box with a source in a vertex and a leakage zone in the opposite with others techniques used in this problem. (author)

Piecewise linear emulator of the nonlinear Schroedinger equation and the resulting analytic solutions for Bose-Einstein condensates

International Nuclear Information System (INIS)

Theodorakis, Stavros

2003-01-01

We emulate the cubic term Ψ 3 in the nonlinear Schroedinger equation by a piecewise linear term, thus reducing the problem to a set of uncoupled linear inhomogeneous differential equations. The resulting analytic expressions constitute an excellent approximation to the exact solutions, as is explicitly shown in the case of the kink, the vortex, and a δ function trap. Such a piecewise linear emulation can be used for any differential equation where the only nonlinearity is a Ψ 3 one. In particular, it can be used for the nonlinear Schroedinger equation in the presence of harmonic traps, giving analytic Bose-Einstein condensate solutions that reproduce very accurately the numerically calculated ones in one, two, and three dimensions

An implicit iterative scheme for solving large systems of linear equations

International Nuclear Information System (INIS)

Barry, J.M.; Pollard, J.P.

1986-12-01

An implicit iterative scheme for the solution of large systems of linear equations arising from neutron diffusion studies is presented. The method is applied to three-dimensional reactor studies and its performance is compared with alternative iterative approaches

The Schroedinger equation for central power law potentials and the classical theory of ordinary linear differential equations of the second order

International Nuclear Information System (INIS)

Lima, M.L.; Mignaco, J.A.

1985-01-01

It is shown that the rational power law potentials in the two-body radial Schoedinger equation admit a systematic treatment available from the classical theory of ordinary linear differential equations of the second order. The admissible potentials come into families evolved from equations having a fixed number of elementary singularities. As a consequence, relations are found and discussed among the several potentials in a family. (Author) [pt

Contact symmetries of general linear second-order ordinary differential equations: letter to the editor

NARCIS (Netherlands)

Martini, Ruud; Kersten, P.H.M.

1983-01-01

Using 1-1 mappings, the complete symmetry groups of contact transformations of general linear second-order ordinary differential equations are determined from two independent solutions of those equations, and applied to the harmonic oscillator with and without damping.

Stable solutions of nonlocal electron heat transport equations

International Nuclear Information System (INIS)

Prasad, M.K.; Kershaw, D.S.

1991-01-01

Electron heat transport equations with a nonlocal heat flux are in general ill-posed and intrinsically unstable, as proved by the present authors [Phys. Fluids B 1, 2430 (1989)]. A straightforward numerical solution of these equations will therefore lead to absurd results. It is shown here that by imposing a minimal set of constraints on the problem it is possible to arrive at a globally stable, consistent, and energy conserving numerical solution

On the calculation of linear stability with the aid of asymptotic solutions of Orr-Sommerfeld equation, 1

International Nuclear Information System (INIS)

Fujimura, Kaoru

1980-11-01

The numerical treatment of Orr-Sommerfeld equation which is the fundamental equation of linear hydrodynamic stability theory is described. Present calculation procedure is applied to the two-dimensional quasi-parallel flow for which linearized disturbance equation (Orr-Sommerfeld equation) contains one simple turning point and αR >> 1. The numerical procedure for this problem and one numerical example for Jeffery-Hamel flow (J-H III 1 ) are presented. These treatment can be extended to the other velocity profiles by slight midifications. (author)

About simple nonlinear and linear superpositions of special exact solutions of Veselov-Novikov equation

International Nuclear Information System (INIS)

Dubrovsky, V. G.; Topovsky, A. V.

2013-01-01

New exact solutions, nonstationary and stationary, of Veselov-Novikov (VN) equation in the forms of simple nonlinear and linear superpositions of arbitrary number N of exact special solutions u (n) , n= 1, …, N are constructed via Zakharov and Manakov ∂-dressing method. Simple nonlinear superpositions are represented up to a constant by the sums of solutions u (n) and calculated by ∂-dressing on nonzero energy level of the first auxiliary linear problem, i.e., 2D stationary Schrödinger equation. It is remarkable that in the zero energy limit simple nonlinear superpositions convert to linear ones in the form of the sums of special solutions u (n) . It is shown that the sums u=u (k 1 ) +...+u (k m ) , 1 ⩽k 1 2 m ⩽N of arbitrary subsets of these solutions are also exact solutions of VN equation. The presented exact solutions include as superpositions of special line solitons and also superpositions of plane wave type singular periodic solutions. By construction these exact solutions represent also new exact transparent potentials of 2D stationary Schrödinger equation and can serve as model potentials for electrons in planar structures of modern electronics.

almaBTE : A solver of the space-time dependent Boltzmann transport equation for phonons in structured materials

Science.gov (United States)

Carrete, Jesús; Vermeersch, Bjorn; Katre, Ankita; van Roekeghem, Ambroise; Wang, Tao; Madsen, Georg K. H.; Mingo, Natalio

2017-11-01

almaBTE is a software package that solves the space- and time-dependent Boltzmann transport equation for phonons, using only ab-initio calculated quantities as inputs. The program can predictively tackle phonon transport in bulk crystals and alloys, thin films, superlattices, and multiscale structures with size features in the nm- μm range. Among many other quantities, the program can output thermal conductances and effective thermal conductivities, space-resolved average temperature profiles, and heat-current distributions resolved in frequency and space. Its first-principles character makes almaBTE especially well suited to investigate novel materials and structures. This article gives an overview of the program structure and presents illustrative examples for some of its uses. PROGRAM SUMMARY Program Title:almaBTE Program Files doi:http://dx.doi.org/10.17632/8tfzwgtp73.1 Licensing provisions: Apache License, version 2.0 Programming language: C++ External routines/libraries: BOOST, MPI, Eigen, HDF5, spglib Nature of problem: Calculation of temperature profiles, thermal flux distributions and effective thermal conductivities in structured systems where heat is carried by phonons Solution method: Solution of linearized phonon Boltzmann transport equation, Variance-reduced Monte Carlo

Non self-similar collapses described by the non-linear Schroedinger equation

International Nuclear Information System (INIS)

Berge, L.; Pesme, D.

1992-01-01

We develop a rapid method in order to find the contraction rates of the radially symmetric collapsing solutions of the nonlinear Schroedinger equation defined for space dimensions exceeding a threshold value. We explicitly determine the asymptotic behaviour of these latter solutions by solving the non stationary linear problem relative to the nonlinear Schroedinger equation. We show that the self-similar states associated with the collapsing solutions are characterized by a spatial extent which is bounded from the top by a cut-off radius

Computer programs for the solution of systems of linear algebraic equations

Science.gov (United States)

Sequi, W. T.

1973-01-01

FORTRAN subprograms for the solution of systems of linear algebraic equations are described, listed, and evaluated in this report. Procedures considered are direct solution, iteration, and matrix inversion. Both incore methods and those which utilize auxiliary data storage devices are considered. Some of the subroutines evaluated require the entire coefficient matrix to be in core, whereas others account for banding or sparceness of the system. General recommendations relative to equation solving are made, and on the basis of tests, specific subprograms are recommended.

Geon-type solutions of the non-linear Heisenberg-Klein-Gordon equation

International Nuclear Information System (INIS)

Mielke, E.W.; Scherzer, R.

1980-10-01

As a model for a ''unitary'' field theory of extended particles we consider the non-linear Klein-Gordon equation - associated with a ''squared'' Heisenberg-Pauli-Weyl non-linear spinor equation - coupled to strong gravity. Using a stationary spherical ansatz for the complex scalar field as well as for the background metric generated via Einstein's field equation, we are able to study the effects of the scalar self-interaction as well as of the classical tensor forces. By numerical integration we obtain a continuous spectrum of localized, gravitational solitons resembling the geons previously constructed for the Einstein-Maxwell system by Wheeler. A self-generated curvature potential originating from the curved background partially confines the Schroedinger type wave functions within the ''scalar geon''. For zero angular momentum states and normalized scalar charge the spectrum for the total gravitational energy of these solitons exhibits a branching with respect to the number of nodes appearing in the radial part of the scalar field. Preliminary studies for higher values of the corresponding ''principal quantum number'' reveal that a kind of fine splitting of the energy levels occurs, which may indicate a rich, particle-like structure of these ''quantized geons''. (author)

An integral equation arising in two group neutron transport theory

International Nuclear Information System (INIS)

Cassell, J S; Williams, M M R

2003-01-01

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

A fresh look at linear ordinary differential equations with constant coefficients. Revisiting the impulsive response method using factorization

Science.gov (United States)

Camporesi, Roberto

2016-01-01

We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations of any order based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as well as of the other more advanced approaches: Laplace transform, linear systems, the general theory of linear equations with variable coefficients and variation of parameters. The approach presented here can be used in a first course on differential equations for science and engineering majors.

Mathematical foundations of transport theory

International Nuclear Information System (INIS)

Ershov, Yu.I.; Shikhov, S.B.

1985-01-01

Foundations of mathematical transport theory are presented. Definitions and theorems of functional analysis are given. Linear kinetic equation of neutron transport in multiplication media is derived. A model of neutron interaction with nuclei of medium determining completely the coefficient properties in transport equation is described. Non-stationary problems regarding and without regard of d=e layed neutrons are analyzed. Results of solving Cauchy problem are discussed

Aspheric surface testing by irradiance transport equation

Science.gov (United States)

Shomali, Ramin; Darudi, Ahmad; Nasiri, Sadollah; Asgharsharghi Bonab, Armir

2010-10-01

In this paper a method for aspheric surface testing is presented. The method is based on solving the Irradiance Transport Equation (ITE).The accuracy of ITE normally depends on the amount of the pick to valley of the phase distribution. This subject is investigated by a simulation procedure.

Solving Linear Equations by Classical Jacobi-SR Based Hybrid Evolutionary Algorithm with Uniform Adaptation Technique

OpenAIRE

Jamali, R. M. Jalal Uddin; Hashem, M. M. A.; Hasan, M. Mahfuz; Rahman, Md. Bazlar

2013-01-01

Solving a set of simultaneous linear equations is probably the most important topic in numerical methods. For solving linear equations, iterative methods are preferred over the direct methods especially when the coefficient matrix is sparse. The rate of convergence of iteration method is increased by using Successive Relaxation (SR) technique. But SR technique is very much sensitive to relaxation factor, {\\omega}. Recently, hybridization of classical Gauss-Seidel based successive relaxation t...

Optimal Homotopy Asymptotic Method for Solving the Linear Fredholm Integral Equations of the First Kind

Directory of Open Access Journals (Sweden)

Mohammad Almousa

2013-01-01

Full Text Available The aim of this study is to present the use of a semi analytical method called the optimal homotopy asymptotic method (OHAM for solving the linear Fredholm integral equations of the first kind. Three examples are discussed to show the ability of the method to solve the linear Fredholm integral equations of the first kind. The results indicated that the method is very effective and simple.

The Schroedinger equation for central power law potentials and the classical theory of ordinary linear differential equations of the second order

International Nuclear Information System (INIS)

Lima, M.L.; Mignaco, J.A.

1985-01-01

It is shown that the rational power law potentials in the two-body radial Schrodinger equations admit a systematic treatment available from the classical theory of ordinary linear differential equations of the second order. The resulting potentials come into families evolved from equations having a fixed number of elementary regular singularities. As a consequence, relations are found and discussed among the several potentials in a family. (Author) [pt

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.

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)

Insights into the School Mathematics Tradition from Solving Linear Equations

Science.gov (United States)

Buchbinder, Orly; Chazan, Daniel; Fleming, Elizabeth

2015-01-01

In this article, we explore how the solving of linear equations is represented in English­-language algebra text books from the early nineteenth century when schooling was becoming institutionalized, and then survey contemporary teachers. In the text books, we identify the increasing presence of a prescribed order of steps (a canonical method) for…

Oscillatory solutions of the Cauchy problem for linear differential equations

Directory of Open Access Journals (Sweden)

Gro Hovhannisyan

2015-06-01

Full Text Available We consider the Cauchy problem for second and third order linear differential equations with constant complex coefficients. We describe necessary and sufficient conditions on the data for the existence of oscillatory solutions. It is known that in the case of real coefficients the oscillatory behavior of solutions does not depend on initial values, but we show that this is no longer true in the complex case: hence in practice it is possible to control oscillatory behavior by varying the initial conditions. Our Proofs are based on asymptotic analysis of the zeros of solutions, represented as linear combinations of exponential functions.

Solution of the transport equation with account for inelastic collisions

International Nuclear Information System (INIS)

Kalashnikov, N.P.; Remizovich, V.S.; Ryazanov, M.I.

1980-01-01

The theory of charged particle scattering in a matter with account for inelastic collisions is considered. In ''directly-forward'' approximation the transport equation at the absence of elastic collisions is obtained. The solution of the transport equation is made without and with account for fluctuation of energy losses. Formulas for path-energy relation are given. Energy spectrum and distribution of fast charged particles with respect to paths are studied. The problem of quantum mechanical approach to the theory of multiple scattering of fast charged particles in a matter is discussed briefly

Role of statistical linearization in the solution of nonlinear stochastic equations

International Nuclear Information System (INIS)

Budgor, A.B.

1977-01-01

The solution of a generalized Langevin equation is referred to as a stochastic process. If the external forcing function is Gaussian white noise, the forward Kolmogarov equation yields the transition probability density function. Nonlinear problems must be handled by approximation procedures e.g., perturbation theories, eigenfunction expansions, and nonlinear optimization procedures. After some comments on the first two of these, attention is directed to the third, and the method of statistical linearization is used to demonstrate a relation to the former two. Nonlinear stochastic systems exhibiting sustained or forced oscillations and the centered nonlinear Schroedinger equation in the presence of Gaussian white noise excitation are considered as examples. 5 figures, 2 tables

Analytical approach to linear fractional partial differential equations arising in fluid mechanics

International Nuclear Information System (INIS)

Momani, Shaher; Odibat, Zaid

2006-01-01

In this Letter, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving linear fractional partial differential equations arising in fluid mechanics. The fractional derivatives are described in the Caputo sense. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these methods, the solution takes the form of a convergent series with easily computable components. The corresponding solutions of the integer order equations are found to follow as special cases of those of fractional order equations. Some numerical examples are presented to illustrate the efficiency and reliability of the two methods

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)

Local Ray-Based Traveltime Computation Using the Linearized Eikonal Equation

KAUST Repository

Almubarak, Mohammed S.

2013-05-01

The computation of traveltimes plays a critical role in the conventional implementations of Kirchhoff migration. Finite-difference-based methods are considered one of the most effective approaches for traveltime calculations and are therefore widely used. However, these eikonal solvers are mainly used to obtain early-arrival traveltime. Ray tracing can be used to pick later traveltime branches, besides the early arrivals, which may lead to an improvement in velocity estimation or in seismic imaging. In this thesis, I improved the accuracy of the solution of the linearized eikonal equation by constructing a linear system of equations (LSE) based on finite-difference approximation, which is of second-order accuracy. The ill-conditioned LSE is initially regularized and subsequently solved to calculate the traveltime update. Numerical tests proved that this method is as accurate as the second-order eikonal solver. Later arrivals are picked using ray tracing. These traveltimes are binned to the nearest node on a regular grid and empty nodes are estimated by interpolating the known values. The resulting traveltime field is used as an input to the linearized eikonal algorithm, which improves the accuracy of the interpolated nodes and yields a local ray-based traveltime. This is a preliminary study and further investigation is required to test the efficiency and the convergence of the solutions.

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)

Linear stochastic differential equations with anticipating initial conditions

DEFF Research Database (Denmark)

Khalifa, Narjess; Kuo, Hui-Hsiung; Ouerdiane, Habib

In this paper we use the new stochastic integral introduced by Ayed and Kuo (2008) and the results obtained by Kuo et al. (2012b) to find a solution to a drift-free linear stochastic differential equation with anticipating initial condition. Our solution is based on well-known results from...... classical Itô theory and anticipative Itô formula results from Kue et al. (2012b). We also show that the solution obtained by our method is consistent with the solution obtained by the methods of Malliavin calculus, e.g. Buckdahn and Nualart (1994)....

A general method for enclosing solutions of interval linear equations

Czech Academy of Sciences Publication Activity Database

Rohn, Jiří

2012-01-01

Roč. 6, č. 4 (2012), s. 709-717 ISSN 1862-4472 R&D Projects: GA ČR GA201/09/1957; GA ČR GC201/08/J020 Institutional research plan: CEZ:AV0Z10300504 Keywords : interval linear equations * solution set * enclosure * absolute value inequality Subject RIV: BA - General Mathematics Impact factor: 1.654, year: 2012

An efficient parallel algorithm for the solution of a tridiagonal linear system of equations

Science.gov (United States)

Stone, H. S.

1971-01-01

Tridiagonal linear systems of equations are solved on conventional serial machines in a time proportional to N, where N is the number of equations. The conventional algorithms do not lend themselves directly to parallel computations on computers of the ILLIAC IV class, in the sense that they appear to be inherently serial. An efficient parallel algorithm is presented in which computation time grows as log sub 2 N. The algorithm is based on recursive doubling solutions of linear recurrence relations, and can be used to solve recurrence relations of all orders.

Least-squares finite element discretizations of neutron transport equations in 3 dimensions

Energy Technology Data Exchange (ETDEWEB)

Manteuffel, T.A [Univ. of Colorado, Boulder, CO (United States); Ressel, K.J. [Interdisciplinary Project Center for Supercomputing, Zurich (Switzerland); Starkes, G. [Universtaet Karlsruhe (Germany)

1996-12-31

The least-squares finite element framework to the neutron transport equation introduced in is based on the minimization of a least-squares functional applied to the properly scaled neutron transport equation. Here we report on some practical aspects of this approach for neutron transport calculations in three space dimensions. The systems of partial differential equations resulting from a P{sub 1} and P{sub 2} approximation of the angular dependence are derived. In the diffusive limit, the system is essentially a Poisson equation for zeroth moment and has a divergence structure for the set of moments of order 1. One of the key features of the least-squares approach is that it produces a posteriori error bounds. We report on the numerical results obtained for the minimum of the least-squares functional augmented by an additional boundary term using trilinear finite elements on a uniform tesselation into cubes.

New methods in linear transport theory. Part of a coordinated programme on methods in neutron transport theory

International Nuclear Information System (INIS)

Mika, J.

1975-09-01

Originally the work was oriented towards two main topics: a) difference and integral methods in neutron transport theory. Two computers were used for numerical calculations GIER and CYBER-72. During the first year the main effort was shifted towards basic theoretical investigations. At the first step the ANIS code was adopted and later modified to check various finite difference approaches against each other. Then the general finite element method and the singular perturbation method were developed. The analysis of singularities of the one-dimensional neutron transport equation in spherical geometry has been done and presented. Later the same analysis for the case of cylindrical symmetry has been carried out. The second and the third year programme included the following topics: 1) finite difference methods in stationary neutron transport theory; 2)mathematical fundamentals of approximate methods for solving the transport equation; 3) singular perturbation method for the time-dependent transport equation; 4) investigation of various iterative procedures in reactor calculations. This investigation will help to better understanding of the mathematical basis for existing and developed numerical methods resulting in more effective algorithms for reactor computer codes

Some oscillation criteria for the second-order linear delay differential equation

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2011-01-01

Roč. 136, č. 2 (2011), s. 195-204 ISSN 0862-7959 Institutional research plan: CEZ:AV0Z10190503 Keywords : second-order linear differential equation with a delay * oscillatory solution Subject RIV: BA - General Mathematics http://www.dml.cz/handle/10338.dmlcz/141582

Dependence of the fundamental time eigenvalue of linear transport operator on the system size and other parameters - An application of the Perron-Frobenius theorem

International Nuclear Information System (INIS)

Sahni, D.C.

1991-01-01

Many papers have been devoted to the study of the spectral properties of the linear (neutron) transport equation. Most of the theoretical investigations have concentrated on the existence (or otherwise) of a continuous spectrum, point spectrum, a leading/dominant eigenvalue, and a corresponding positive eigenvector. It is shown that the fundamental time eigenvalue of the linear transport operator increases with the size of the system. This follows from the increase in the largest eigenvalue of a non-negative irreducible matrix whenever any matrix element his increased. This result of matrix analysis is generalized to more general Krein-Rutman operators that leave a cone of vectors invariant

Fourier analysis of a new P1 synthetic acceleration for Sn transport equations

International Nuclear Information System (INIS)

Turcksin, B.; Ragusa, J. C.

2010-10-01

In this work, is derived a new P1 synthetic acceleration scheme (P1SA) for the S N transport equation and analyze its convergence properties through the means of a Fourier analysis. The Fourier analysis is carried out for both continuous (i.e., not spatially discretized) S N equations and linear discontinuous Fem discretization. We show, thanks to the continuous analysis, that the scheme is unstable when the anisotropy is important (μ - >0.5). However, the discrete analysis shows that when cells are large in comparison to the mean free path, the spectral radius decreases and the acceleration scheme becomes effective, even for highly anisotropic scattering. In charged particles transport, scattering is highly anisotropic and mean free paths are very small and, thus, this scheme could be of interest. To use the P1SA when cells are small and anisotropy is important, the scheme is modified by altering the update of the accelerated flux or by using either K transport sweeps before the application of P1SA. The update scheme performs well as long as μ - - ≥0.9, the modified update scheme is unstable. The multiple transport sweeps scheme is convergent with an arbitrary μ - but the spectral radius increases when scattering is isotropic. When anisotropic increases, the frequency of use of the acceleration scheme needs to be decreased. Even if the P1SA is used less often, the spectral radius is significantly smaller when compared with a method that does not use it for high anisotropy (μ - ≥0.5). It is interesting to notice that using P1SA every two iterations gives the same spectral radius than the update method when μ - ≥0.5 but it is much less efficient when μ - <0.5. (Author)

Refined Fuchs inequalities for systems of linear differential equations

International Nuclear Information System (INIS)

Gontsov, R R

2004-01-01

We refine the Fuchs inequalities obtained by Corel for systems of linear meromorphic differential equations given on the Riemann sphere. Fuchs inequalities enable one to estimate the sum of exponents of the system over all its singular points. We refine these well-known inequalities by considering the Jordan structure of the leading coefficient of the Laurent series for the matrix of the right-hand side of the system in the neighbourhood of a singular point

A non linear half space problem for radiative transfer equations. Application to the Rosseland approximation

International Nuclear Information System (INIS)

Sentis, R.

1984-07-01

The radiative transfer equations may be approximated by a non linear diffusion equation (called Rosseland equation) when the mean free paths of the photons are small with respect to the size of the medium. Some technical assomptions are made, namely about the initial conditions, to avoid any problem of initial layer terms

Discrete Ordinates Approximations to the First- and Second-Order Radiation Transport Equations

International Nuclear Information System (INIS)

FAN, WESLEY C.; DRUMM, CLIFTON R.; POWELL, JENNIFER L. email wcfan@sandia.gov

2002-01-01

The conventional discrete ordinates approximation to the Boltzmann transport equation can be described in a matrix form. Specifically, the within-group scattering integral can be represented by three components: a moment-to-discrete matrix, a scattering cross-section matrix and a discrete-to-moment matrix. Using and extending these entities, we derive and summarize the matrix representations of the second-order transport equations

Discrete Ordinates Approximations to the First- and Second-Order Radiation Transport Equations

CERN Document Server

Fan, W C; Powell, J L

2002-01-01

The conventional discrete ordinates approximation to the Boltzmann transport equation can be described in a matrix form. Specifically, the within-group scattering integral can be represented by three components: a moment-to-discrete matrix, a scattering cross-section matrix and a discrete-to-moment matrix. Using and extending these entities, we derive and summarize the matrix representations of the second-order transport equations.

Separation-induced boundary layer transition: Modeling with a non-linear eddy-viscosity model coupled with the laminar kinetic energy equation

International Nuclear Information System (INIS)

Vlahostergios, Z.; Yakinthos, K.; Goulas, A.

2009-01-01

We present an effort to model the separation-induced transition on a flat plate with a semi-circular leading edge, using a cubic non-linear eddy-viscosity model combined with the laminar kinetic energy. A non-linear model, compared to a linear one, has the advantage to resolve the anisotropic behavior of the Reynolds-stresses in the near-wall region and it provides a more accurate expression for the generation of turbulence in the transport equation of the turbulence kinetic energy. Although in its original formulation the model is not able to accurately predict the separation-induced transition, the inclusion of the laminar kinetic energy increases its accuracy. The adoption of the laminar kinetic energy by the non-linear model is presented in detail, together with some additional modifications required for the adaption of the laminar kinetic energy into the basic concepts of the non-linear eddy-viscosity model. The computational results using the proposed combined model are shown together with the ones obtained using an isotropic linear eddy-viscosity model, which adopts also the laminar kinetic energy concept and in comparison with the existing experimental data.

Dynamical and geometric aspects of Hamilton-Jacobi and linearized Monge-Ampère equations VIASM 2016

CERN Document Server

Tran, Hung

2017-01-01

Consisting of two parts, the first part of this volume is an essentially self-contained exposition of the geometric aspects of local and global regularity theory for the Monge–Ampère and linearized Monge–Ampère equations. As an application, we solve the second boundary value problem of the prescribed affine mean curvature equation, which can be viewed as a coupling of the latter two equations. Of interest in its own right, the linearized Monge–Ampère equation also has deep connections and applications in analysis, fluid mechanics and geometry, including the semi-geostrophic equations in atmospheric flows, the affine maximal surface equation in affine geometry and the problem of finding Kahler metrics of constant scalar curvature in complex geometry. Among other topics, the second part provides a thorough exposition of the large time behavior and discounted approximation of Hamilton–Jacobi equations, which have received much attention in the last two decades, and a new approach to the subject, the n...

Scattering theory of the linear Boltzmann operator

International Nuclear Information System (INIS)

Hejtmanek, J.

1975-01-01

In time dependent scattering theory we know three important examples: the wave equation around an obstacle, the Schroedinger and the Dirac equation with a scattering potential. In this paper another example from time dependent linear transport theory is added and considered in full detail. First the linear Boltzmann operator in certain Banach spaces is rigorously defined, and then the existence of the Moeller operators is proved by use of the theorem of Cook-Jauch-Kuroda, that is generalized to the case of a Banach space. (orig.) [de

Solution of linear ordinary differential equations by means of the method of variation of arbitrary constants

DEFF Research Database (Denmark)

Mejlbro, Leif

1997-01-01

An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians.......An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians....

Linear and nonlinear analogues of the Schroedinger equation in the contextual approach in quantum mechanics

International Nuclear Information System (INIS)

Khrennikov, A.Yu.

2005-01-01

One derived the general evolutionary differential equation within the Hilbert space describing dynamics of the wave function. The derived contextual model is more comprehensive in contrast to a quantum one. The contextual equation may be a nonlinear one. Paper presents the conditions ensuring linearity of the evolution and derivation of the Schroedinger equation [ru

Multi-point boundary value problems for linear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Domoshnitsky, A.; Hakl, Robert; Půža, Bedřich

2017-01-01

Roč. 24, č. 2 (2017), s. 193-206 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : boundary value problems * linear functional-differential equations * functional-differential inequalities Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0076/gmj-2016-0076. xml

Multi-point boundary value problems for linear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Domoshnitsky, A.; Hakl, Robert; Půža, Bedřich

2017-01-01

Roč. 24, č. 2 (2017), s. 193-206 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : boundary value problems * linear functional- differential equations * functional- differential inequalities Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0076/gmj-2016-0076.xml

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

Matrix form of Legendre polynomials for solving linear integro-differential equations of high order

Science.gov (United States)

Kammuji, M.; Eshkuvatov, Z. K.; Yunus, Arif A. M.

2017-04-01

This paper presents an effective approximate solution of high order of Fredholm-Volterra integro-differential equations (FVIDEs) with boundary condition. Legendre truncated series is used as a basis functions to estimate the unknown function. Matrix operation of Legendre polynomials is used to transform FVIDEs with boundary conditions into matrix equation of Fredholm-Volterra type. Gauss Legendre quadrature formula and collocation method are applied to transfer the matrix equation into system of linear algebraic equations. The latter equation is solved by Gauss elimination method. The accuracy and validity of this method are discussed by solving two numerical examples and comparisons with wavelet and methods.

First order linear ordinary differential equations in associative algebras

Directory of Open Access Journals (Sweden)

Gordon Erlebacher

2004-01-01

Full Text Available In this paper, we study the linear differential equation $$ frac{dx}{dt}=sum_{i=1}^n a_i(t x b_i(t + f(t $$ in an associative but non-commutative algebra $mathcal{A}$, where the $b_i(t$ form a set of commuting $mathcal{A}$-valued functions expressed in a time-independent spectral basis consisting of mutually annihilating idempotents and nilpotents. Explicit new closed solutions are derived, and examples are presented to illustrate the theory.

Myshkis type oscillation criteria for second-order linear delay differential equations

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2015-01-01

Roč. 178, č. 1 (2015), s. 143-161 ISSN 0026-9255 Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillation criteria Subject RIV: BA - General Mathematics Impact factor: 0.664, year: 2015 http://link.springer.com/article/10.1007%2Fs00605-014-0719-y

Comparison of nonlinearities in oscillation theory of half-linear differential equations

Czech Academy of Sciences Publication Activity Database

Řehák, Pavel

2008-01-01

Roč. 121, č. 2 (2008), s. 93-105 ISSN 0236-5294 R&D Projects: GA AV ČR KJB100190701 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear differential equation * comparison theorem * Riccati technique Subject RIV: BA - General Mathematics Impact factor: 0.317, year: 2008

Linear hyperbolic functional-differential equations with essentially bounded right-hand side

Czech Academy of Sciences Publication Activity Database

Domoshnitsky, A.; Lomtatidze, Alexander; Maghakyan, A.; Šremr, Jiří

2011-01-01

Roč. 2011, - (2011), s. 242965 ISSN 1085-3375 Institutional research plan: CEZ:AV0Z10190503 Keywords : linear functional-differential equation of hyperbolic type * Darboux problem * unique solvability Subject RIV: BA - General Mathematics Impact factor: 1.318, year: 2011 http://www.hindawi.com/journals/ aaa /2011/242965/

Transport equations in an enzymatic glucose fuel cell

Science.gov (United States)

Jariwala, Soham; Krishnamurthy, Balaji

2018-01-01

A mathematical model is developed to study the effects of convective flux and operating temperature on the performance of an enzymatic glucose fuel cell with a membrane. The model assumes isothermal operating conditions and constant feed rate of glucose. The glucose fuel cell domain is divided into five sections, with governing equations describing transport characteristics in each region, namely - anode diffusion layer, anode catalyst layer (enzyme layer), membrane, cathode catalyst layer and cathode diffusion layer. The mass transport is assumed to be one-dimensional and the governing equations are solved numerically. The effects flow rate of glucose feed on the performance of the fuel cell are studied as it contributes significantly to the convective flux. The effects of operating temperature on the performance of a glucose fuel cell are also modeled. The cell performances are compared using cell polarization curves, which were found compliant with experimental observations.

A stochastic solution of the advective transport equation with uncertainty

International Nuclear Information System (INIS)

Williams, M.M.R.

1991-01-01

A model has been developed for calculating the transport of water-borne radionuclides through layers of porous materials, such as rock or clay. The model is based upon a purely advective transport equation, in which the fluid velocity is a random variable, thereby simulating dispersion in a more realistic manner than the ad hoc introduction of a dispersivity. In addition to a random velocity field, which is an observable physical phenomenon, allowance is made for uncertainty in our knowledge of the parameters which enter the equation, e.g. the retardation coefficient. This too, is assumed to be a random variable and contributes to the stochasticity of the resulting partial differential equation of transport. The stochastic differential equation can be solved analytically and then ensemble averages taken over the associated probability distribution of velocity and retardation coefficient. A method based upon a novel form of the central limit theorem of statistics is employed to obtain tractable solutions of a system consisting of many serial legs of varying properties. One interesting conclusion is that the total flux out of a medium is significantly underestimated by using the deterministic solution with an average transit time compared with that from the stochastically averaged solution. The theory is illustrated numerically for a number of physically relevant cases. (author) 8 figs., 4 tabs., 7 refs

Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations

Directory of Open Access Journals (Sweden)

Matt Challacombe

2014-03-01

Full Text Available A new approach to solving the Time-Dependent Self-Consistent-Field equations is developed based on the double quotient formulation of Tsiper 2001 (J. Phys. B. Dual channel, quasi-independent non-linear optimization of these quotients is found to yield convergence rates approaching those of the best case (single channel Tamm-Dancoff approximation. This formulation is variational with respect to matrix truncation, admitting linear scaling solution of the matrix-eigenvalue problem, which is demonstrated for bulk excitons in the polyphenylene vinylene oligomer and the (4,3 carbon nanotube segment.

A novel algebraic procedure for solving non-linear evolution equations of higher order

International Nuclear Information System (INIS)

Huber, Alfred

2007-01-01

We report here a systematic approach that can easily be used for solving non-linear partial differential equations (nPDE), especially of higher order. We restrict the analysis to the so called evolution equations describing any wave propagation. The proposed new algebraic approach leads us to traveling wave solutions and moreover, new class of solution can be obtained. The crucial step of our method is the basic assumption that the solutions satisfy an ordinary differential equation (ODE) of first order that can be easily integrated. The validity and reliability of the method is tested by its application to some non-linear evolution equations. The important aspect of this paper however is the fact that we are able to calculate distinctive class of solutions which cannot be found in the current literature. In other words, using this new algebraic method the solution manifold is augmented to new class of solution functions. Simultaneously we would like to stress the necessity of such sophisticated methods since a general theory of nPDE does not exist. Otherwise, for practical use the algebraic construction of new class of solutions is of fundamental interest

OSCILLATION OF A SECOND-ORDER HALF-LINEAR NEUTRAL DAMPED DIFFERENTIAL EQUATION WITH TIME-DELAY

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

In this paper,the oscillation for a class of second-order half-linear neutral damped differential equation with time-delay is studied.By means of Yang-inequality,the generalized Riccati transformation and a certain function,some new sufficient conditions for the oscillation are given for all solutions to the equation.

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)

Mathematics Literacy of Secondary Students in Solving Simultanenous Linear Equations

Science.gov (United States)

Sitompul, R. S. I.; Budayasa, I. K.; Masriyah

2018-01-01

This study examines the profile of secondary students’ mathematical literacy in solving simultanenous linear equations problems in terms of cognitive style of visualizer and verbalizer. This research is a descriptive research with qualitative approach. The subjects in this research consist of one student with cognitive style of visualizer and one student with cognitive style of verbalizer. The main instrument in this research is the researcher herself and supporting instruments are cognitive style tests, mathematics skills tests, problem-solving tests and interview guidelines. Research was begun by determining the cognitive style test and mathematics skill test. The subjects chosen were given problem-solving test about simultaneous linear equations and continued with interview. To ensure the validity of the data, the researcher conducted data triangulation; the steps of data reduction, data presentation, data interpretation, and conclusion drawing. The results show that there is a similarity of visualizer and verbalizer-cognitive style in identifying and understanding the mathematical structure in the process of formulating. There are differences in how to represent problems in the process of implementing, there are differences in designing strategies and in the process of interpreting, and there are differences in explaining the logical reasons.

Hadronic equation of state in the statistical bootstrap model and linear graph theory

International Nuclear Information System (INIS)

Fre, P.; Page, R.

1976-01-01

Taking a statistical mechanical point og view, the statistical bootstrap model is discussed and, from a critical analysis of the bootstrap volume comcept, it is reached a physical ipothesis, which leads immediately to the hadronic equation of state provided by the bootstrap integral equation. In this context also the connection between the statistical bootstrap and the linear graph theory approach to interacting gases is analyzed

On the Cauchy problem for a Sobolev-type equation with quadratic non-linearity

International Nuclear Information System (INIS)

Aristov, Anatoly I

2011-01-01

We investigate the asymptotic behaviour as t→∞ of the solution of the Cauchy problem for a Sobolev-type equation with quadratic non-linearity and develop ideas used by I. A. Shishmarev and other authors in the study of classical and Sobolev-type equations. Conditions are found under which it is possible to consider the case of an arbitrary dimension of the spatial variable.

A generalized variational algebra and conserved densities for linear evolution equations

International Nuclear Information System (INIS)

Abellanas, L.; Galindo, A.

1978-01-01

The symbolic algebra of Gel'fand and Dikii is generalized to the case of n variables. Using this algebraic approach a rigorous characterization of the polynomial kernel of the variational derivative is given. This is applied to classify all the conservation laws for linear polynomial evolution equations of arbitrary order. (Auth.)

Peculiarities in power type comparison results for half-linear dynamic equations

Czech Academy of Sciences Publication Activity Database

Řehák, Pavel

2012-01-01

Roč. 42, č. 6 (2012), s. 1995-2013 ISSN 0035-7596 R&D Projects: GA AV ČR KJB100190701 Institutional support: RVO:67985840 Keywords : half-linear dynamic equation * time scale * comparison theorem Subject RIV: BA - General Mathematics Impact factor: 0.389, year: 2012 http://projecteuclid.org/euclid.rmjm/1361800616

About simple nonlinear and linear superpositions of special exact solutions of Veselov-Novikov equation

Energy Technology Data Exchange (ETDEWEB)

Dubrovsky, V. G.; Topovsky, A. V. [Novosibirsk State Technical University, Karl Marx prosp. 20, Novosibirsk 630092 (Russian Federation)

2013-03-15

New exact solutions, nonstationary and stationary, of Veselov-Novikov (VN) equation in the forms of simple nonlinear and linear superpositions of arbitrary number N of exact special solutions u{sup (n)}, n= 1, Horizontal-Ellipsis , N are constructed via Zakharov and Manakov {partial_derivative}-dressing method. Simple nonlinear superpositions are represented up to a constant by the sums of solutions u{sup (n)} and calculated by {partial_derivative}-dressing on nonzero energy level of the first auxiliary linear problem, i.e., 2D stationary Schroedinger equation. It is remarkable that in the zero energy limit simple nonlinear superpositions convert to linear ones in the form of the sums of special solutions u{sup (n)}. It is shown that the sums u=u{sup (k{sub 1})}+...+u{sup (k{sub m})}, 1 Less-Than-Or-Slanted-Equal-To k{sub 1} < k{sub 2} < Horizontal-Ellipsis < k{sub m} Less-Than-Or-Slanted-Equal-To N of arbitrary subsets of these solutions are also exact solutions of VN equation. The presented exact solutions include as superpositions of special line solitons and also superpositions of plane wave type singular periodic solutions. By construction these exact solutions represent also new exact transparent potentials of 2D stationary Schroedinger equation and can serve as model potentials for electrons in planar structures of modern electronics.

Chaotic dynamics and diffusion in a piecewise linear equation

International Nuclear Information System (INIS)

Shahrear, Pabel; Glass, Leon; Edwards, Rod

2015-01-01

Genetic interactions are often modeled by logical networks in which time is discrete and all gene activity states update simultaneously. However, there is no synchronizing clock in organisms. An alternative model assumes that the logical network is preserved and plays a key role in driving the dynamics in piecewise nonlinear differential equations. We examine dynamics in a particular 4-dimensional equation of this class. In the equation, two of the variables form a negative feedback loop that drives a second negative feedback loop. By modifying the original equations by eliminating exponential decay, we generate a modified system that is amenable to detailed analysis. In the modified system, we can determine in detail the Poincaré (return) map on a cross section to the flow. By analyzing the eigenvalues of the map for the different trajectories, we are able to show that except for a set of measure 0, the flow must necessarily have an eigenvalue greater than 1 and hence there is sensitive dependence on initial conditions. Further, there is an irregular oscillation whose amplitude is described by a diffusive process that is well-modeled by the Irwin-Hall distribution. There is a large class of other piecewise-linear networks that might be analyzed using similar methods. The analysis gives insight into possible origins of chaotic dynamics in periodically forced dynamical systems

Chaotic dynamics and diffusion in a piecewise linear equation

Science.gov (United States)

Shahrear, Pabel; Glass, Leon; Edwards, Rod

2015-03-01

Genetic interactions are often modeled by logical networks in which time is discrete and all gene activity states update simultaneously. However, there is no synchronizing clock in organisms. An alternative model assumes that the logical network is preserved and plays a key role in driving the dynamics in piecewise nonlinear differential equations. We examine dynamics in a particular 4-dimensional equation of this class. In the equation, two of the variables form a negative feedback loop that drives a second negative feedback loop. By modifying the original equations by eliminating exponential decay, we generate a modified system that is amenable to detailed analysis. In the modified system, we can determine in detail the Poincaré (return) map on a cross section to the flow. By analyzing the eigenvalues of the map for the different trajectories, we are able to show that except for a set of measure 0, the flow must necessarily have an eigenvalue greater than 1 and hence there is sensitive dependence on initial conditions. Further, there is an irregular oscillation whose amplitude is described by a diffusive process that is well-modeled by the Irwin-Hall distribution. There is a large class of other piecewise-linear networks that might be analyzed using similar methods. The analysis gives insight into possible origins of chaotic dynamics in periodically forced dynamical systems.

Dyson-Schwinger equations for the non-linear σ-model

International Nuclear Information System (INIS)

Drouffe, J.M.; Flyvbjerg, H.

1989-08-01

Dyson-Schwinger equations for the O(N)-symmetric non-linear σ-model are derived. They are polynomials in N, hence 1/N-expanded ab initio. A finite, closed set of equations is obtained by keeping only the leading term and the first correction term in this 1/N-series. These equations are solved numerically in two dimensions on square lattices measuring 50x50, 100x100, 200x200, and 400x400. They are also solved analytically at strong coupling and at weak coupling in a finite volume. In these two limits the solution is asymptotically identical to the exact strong- and weak-coupling series through the first three terms. Between these two limits, results for the magnetic susceptibility and the mass gap are identical to the Monte Carlo results available for N=3 and N=4 within a uniform systematic error of O(1/N 3 ), i.e. the results seem good to O(1/N 2 ), though obtained from equations that are exact only to O(1/N). This is understood by seeing the results as summed infinite subseries of the 1/N-series for the exact susceptibility and mass gap. We conclude that the kind of 1/N-expansion presented here converges as well as one might ever hope for, even for N as small as 3. (orig.)

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.

Time dependence linear transport III convergence of the discrete ordinate method

International Nuclear Information System (INIS)

Wilson, D.G.

1983-01-01

In this paper the uniform pointwise convergence of the discrete ordinate method for weak and strong solutions of the time dependent, linear transport equation posed in a multidimensional, rectangular parallelepiped with partially reflecting walls is established. The first result is that a sequence of discrete ordinate solutions converges uniformly on the quadrature points to a solution of the continuous problem provided that the corresponding sequence of truncation errors for the solution of the continuous problem converges to zero in the same manner. The second result is that continuity of the solution with respect to the velocity variables guarantees that the truncation erros in the quadrature formula go the zero and hence that the discrete ordinate approximations converge to the solution of the continuous problem as the discrete ordinate become dense. An existence theory for strong solutions of the the continuous problem follows as a result

The accuracy of time dependent transport equation ergodic approximation

International Nuclear Information System (INIS)

Stancic, V.

1995-01-01

In order to predict the accuracy of the ergodic approximation for solving the time dependent transport equation, a comparison with respect to multiple collision and time finite difference methods, has been considered. (author)

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.

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)

A Comparison of Angular Difference Schemes for One-Dimensional Spherical Geometry SN Equations

International Nuclear Information System (INIS)

Lathrop, K.D.

2000-01-01

To investigate errors caused by angular differencing in approximating the streaming terms of the transport equation, five different approximations are evaluated for three test problems in one-dimensional spherical geometry. The following schemes are compared: diamond, special truncation error minimizing weighted diamond, linear continuous (the original S N scheme), linear discontinuous, and new quadratic continuous. To isolate errors caused by angular differencing, the approximations are derived from the transport equation without spatial differencing, and the resulting coupled ordinary differential equations (ODEs) are solved with an ODE solver. Results from the approximations are compared with analytic solutions derived for two-region purely absorbing spheres. Most of the approximations are derived by taking moments of the conservation form of the transport equation. The quadratic continuous approximation is derived taking the zeroth moment of both the transport equation and the first angular derivative of the transport equation. The advantages of this approach are described. In all of the approximations, the desirability is shown of using an initializing computation of the μ = -1 angular flux to correctly compute the central flux and of having a difference approximation that ensures this central flux is the same for all directions. The behavior of the standard discrete ordinates equations in the diffusion limit is reviewed, and the linear and quadratic continuous approximations are shown to have the correct diffusion limit if an equal interval discrete quadrature is used.In all three test problems, the weighted diamond difference approximation has smaller maximum and average relative flux errors than the diamond or the linear continuous difference approximations. The quadratic continuous approximation and the linear discontinuous approximation are both more accurate than the other approximations, and the quadratic continuous approximation has a decided edge

A comparison of angular difference schemes for one-dimensional spherical geometry SN equations

International Nuclear Information System (INIS)

Lathrop, K.D.

2000-01-01

To investigate errors caused by angular differencing in approximating the streaming terms of the transport equation, five different approximations are evaluated for three test problems in one-dimensional spherical geometry. The following schemes are compared: diamond, special truncation error minimizing weighted diamond, linear continuous (the original S N scheme), linear discontinuous, and new quadratic continuous. To isolate errors caused by angular differencing, the approximations are derived from the transport equation without spatial differencing, and the resulting coupled ordinary differential equations (ODEs) are solved with an ODE solver. Results from the approximations are compared with analytic solutions derived for two-region purely absorbing spheres. Most of the approximations are derived by taking moments of the conservation form of the transport equation. The quadratic continuous approximation is derived taking the zeroth moment of both the transport equation and the first angular derivative of the transport equation. The advantages of this approach are described, In all of the approximations, the desirability is shown of using an initializing computation of the μ = -1 angular flux to correctly compute the central flux and of having a difference approximation that ensures this central flux is the same for all directions. The behavior of the standard discrete ordinates equations in the diffusion limit is reviewed, and the linear and quadratic continuous approximations are shown to have the correct diffusion limit if an equal interval discrete quadrature is used. In all three test problems, the weighted diamond difference approximation has smaller maximum and average relative flux errors than the diamond or the linear continuous difference approximations. The quadratic continuous approximation and the linear discontinuous approximation are both more accurate than the other approximations, and the quadratic continuous approximation has a decided edge

An Improved Recurrent Neural Network for Complex-Valued Systems of Linear Equation and Its Application to Robotic Motion Tracking.

Science.gov (United States)

Ding, Lei; Xiao, Lin; Liao, Bolin; Lu, Rongbo; Peng, Hua

2017-01-01

To obtain the online solution of complex-valued systems of linear equation in complex domain with higher precision and higher convergence rate, a new neural network based on Zhang neural network (ZNN) is investigated in this paper. First, this new neural network for complex-valued systems of linear equation in complex domain is proposed and theoretically proved to be convergent within finite time. Then, the illustrative results show that the new neural network model has the higher precision and the higher convergence rate, as compared with the gradient neural network (GNN) model and the ZNN model. Finally, the application for controlling the robot using the proposed method for the complex-valued systems of linear equation is realized, and the simulation results verify the effectiveness and superiorness of the new neural network for the complex-valued systems of linear equation.

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 piecewise linear finite element discretization of the diffusion equation for arbitrary polyhedral grids

Energy Technology Data Exchange (ETDEWEB)

Bailey, T S; Adams, M L [Texas A M Univ., Dept. of Nuclear Engineering, College Station, TX (United States); Yang, B; Zika, M R [Lawrence Livermore National Lab., Livermore, CA (United States)

2005-07-01

We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses piecewise linear weight and basis functions in the finite element approximation, and it can be applied on arbitrary polygonal (2-dimensional) or polyhedral (3-dimensional) grids. We show that this new PWL method gives solutions comparable to those from Palmer's finite-volume method. However, since the PWL method produces a symmetric positive definite coefficient matrix, it should be substantially more computationally efficient than Palmer's method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids. (authors)

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

Science.gov (United States)

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

2018-05-01

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

New diffusion-like solutions of one-speed transport equations in spherical geometry

International Nuclear Information System (INIS)

Sahni, D.C.

1988-01-01

Stationary, one-speed, spherically symmetric transport equations are considered in a conservative medium. Closed-form expressions are obtained for the angular flux ψ(r, μ) that yield a total flux varying as 1/r by using Sonine transforms. Properties of this solution are studied and it is shown that the solution can not be identified as a diffusion mode solution of the transport equation. Limitations of the Sonine transform technique are noted. (author)

Solutions of half-linear differential equations in the classes Gamma and Pi

Czech Academy of Sciences Publication Activity Database

Řehák, Pavel; Taddei, V.

2016-01-01

Roč. 29, 7-8 (2016), s. 683-714 ISSN 0893-4983 Institutional support: RVO:67985840 Keywords : half-linear differential equation * positive solution * asymptotic formula Subject RIV: BA - General Mathematics Impact factor: 0.565, year: 2016 http://projecteuclid.org/euclid.die/1462298681

Linear Ordinary Differential Equations with Constant Coefficients. Revisiting the Impulsive Response Method Using Factorization

Science.gov (United States)

Camporesi, Roberto

2011-01-01

We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as well as of…

Stochastic differential equation model for linear growth birth and death processes with immigration and emigration

International Nuclear Information System (INIS)

Granita; Bahar, A.

2015-01-01

This paper discusses on linear birth and death with immigration and emigration (BIDE) process to stochastic differential equation (SDE) model. Forward Kolmogorov equation in continuous time Markov chain (CTMC) with a central-difference approximation was used to find Fokker-Planckequation corresponding to a diffusion process having the stochastic differential equation of BIDE process. The exact solution, mean and variance function of BIDE process was found

Stochastic differential equation model for linear growth birth and death processes with immigration and emigration