Master equation approach to the intra-urban passenger flow and application to the Metropolitan Seoul Subway system

Energy Technology Data Exchange (ETDEWEB)

Lee, Keumsook [Department of Geography, Sungshin University, Seoul 136-742 (Korea, Republic of); Goh, Segun; Choi, M Y [Department of Physics and Astronomy and Center for Theoretical Physics, Seoul National University, Seoul 151-747 (Korea, Republic of); Park, Jong Soo [School of Information Technology, Sungshin University, Seoul 136-742 (Korea, Republic of); Jung, Woo-Sung, E-mail: kslee@sungshin.ac.kr, E-mail: mychoi@snu.ac.kr [Department of Physics and Basic Science Research Institute, Pohang University of Science and Technology, Pohang 790-784 (Korea, Republic of)

2011-03-18

The master equation approach is proposed to describe the evolution of passengers in a subway system. With the transition rate constructed from simple geographical consideration, the evolution equation for the distribution of subway passengers is found to bear skew distributions including log-normal, Weibull, and power-law distributions. This approach is then applied to the Metropolitan Seoul Subway system: analysis of the trip data of all passengers in a day reveals that the data in most cases fit well to the log-normal distributions. Implications of the results are also discussed.

Master equation approach to the intra-urban passenger flow and application to the Metropolitan Seoul Subway system

International Nuclear Information System (INIS)

Lee, Keumsook; Goh, Segun; Choi, M Y; Park, Jong Soo; Jung, Woo-Sung

2011-01-01

The master equation approach is proposed to describe the evolution of passengers in a subway system. With the transition rate constructed from simple geographical consideration, the evolution equation for the distribution of subway passengers is found to bear skew distributions including log-normal, Weibull, and power-law distributions. This approach is then applied to the Metropolitan Seoul Subway system: analysis of the trip data of all passengers in a day reveals that the data in most cases fit well to the log-normal distributions. Implications of the results are also discussed.

Statistical mechanics of normal grain growth in one dimension: A partial integro-differential equation model

International Nuclear Information System (INIS)

Ng, Felix S.L.

2016-01-01

We develop a statistical-mechanical model of one-dimensional normal grain growth that does not require any drift-velocity parameterization for grain size, such as used in the continuity equation of traditional mean-field theories. The model tracks the population by considering grain sizes in neighbour pairs; the probability of a pair having neighbours of certain sizes is determined by the size-frequency distribution of all pairs. Accordingly, the evolution obeys a partial integro-differential equation (PIDE) over ‘grain size versus neighbour grain size’ space, so that the grain-size distribution is a projection of the PIDE's solution. This model, which is applicable before as well as after statistically self-similar grain growth has been reached, shows that the traditional continuity equation is invalid outside this state. During statistically self-similar growth, the PIDE correctly predicts the coarsening rate, invariant grain-size distribution and spatial grain size correlations observed in direct simulations. The PIDE is then reducible to the standard continuity equation, and we derive an explicit expression for the drift velocity. It should be possible to formulate similar parameterization-free models of normal grain growth in two and three dimensions.

A simple algebraic cancer equation: calculating how cancers may arise with normal mutation rates

Directory of Open Access Journals (Sweden)

Shibata Darryl

2010-01-01

Full Text Available Abstract Background The purpose of this article is to present a relatively easy to understand cancer model where transformation occurs when the first cell, among many at risk within a colon, accumulates a set of driver mutations. The analysis of this model yields a simple algebraic equation, which takes as inputs the number of stem cells, mutation and division rates, and the number of driver mutations, and makes predictions about cancer epidemiology. Methods The equation [p = 1 - (1 - (1 - (1 - udkNm ] calculates the probability of cancer (p and contains five parameters: the number of divisions (d, the number of stem cells (N × m, the number of critical rate-limiting pathway driver mutations (k, and the mutation rate (u. In this model progression to cancer "starts" at conception and mutations accumulate with cell division. Transformation occurs when a critical number of rate-limiting pathway mutations first accumulates within a single stem cell. Results When applied to several colorectal cancer data sets, parameter values consistent with crypt stem cell biology and normal mutation rates were able to match the increase in cancer with aging, and the mutation frequencies found in cancer genomes. The equation can help explain how cancer risks may vary with age, height, germline mutations, and aspirin use. APC mutations may shorten pathways to cancer by effectively increasing the numbers of stem cells at risk. Conclusions The equation illustrates that age-related increases in cancer frequencies may result from relatively normal division and mutation rates. Although this equation does not encompass all of the known complexity of cancer, it may be useful, especially in a teaching setting, to help illustrate relationships between small and large cancer features.

Integrable systems of partial differential equations determined by structure equations and Lax pair

International Nuclear Information System (INIS)

Bracken, Paul

2010-01-01

It is shown how a system of evolution equations can be developed both from the structure equations of a submanifold embedded in three-space as well as from a matrix SO(6) Lax pair. The two systems obtained this way correspond exactly when a constraint equation is selected and imposed on the system of equations. This allows for the possibility of selecting the coefficients in the second fundamental form in a general way.

Multiple normalized solutions for a planar gauged nonlinear Schrödinger equation

Science.gov (United States)

Luo, Xiao

2018-06-01

We study the existence, multiplicity, quantitative property and asymptotic behavior of normalized solutions for a gauged nonlinear Schrödinger equation arising from the Chern-Simons theory Δ u + ω u +|x|^2u+ λ ( {{h^2}(| x | )}/{{{| x | ^2}}} + \\int \\limits _{| x | }^{ + ∞} {{h(s)}/s} {u^2}(s)ds) u = {| u | ^{p - 2}}u,\\quad x\\in R^2, where ω \\in R, λ >0, p>4 and h(s) = 1/2\\int \\limits _0^s {r{u^2}(r)dr} . Combining constraint minimization method and minimax principle, we prove that the problem possesses at least two normalized solutions: One is a ground state and the other is an excited state. Furthermore, the asymptotic behavior and quantitative property of the ground state are analyzed.

Parallels between control PDE's (Partial Differential Equations) and systems of ODE's (Ordinary Differential Equations)

Science.gov (United States)

Hunt, L. R.; Villarreal, Ramiro

1987-01-01

System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differentail equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralleled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.

Nonlinear integrodifferential equations as discrete systems

Science.gov (United States)

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

1999-06-01

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

Solutions of system of P1 equations without use of auxiliary differential equations coupled

International Nuclear Information System (INIS)

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

2000-01-01

The system of P1 equations is composed by two equations coupled itself one for the neutron flux and other for the current. Usually this system is solved by definitions of two integrals parameters, which are named slowing down densities of the flux and the current. Hence, the system P1 can be change from integral to only two differential equations. However, there are two new differentials equations that may be solved with the initial system. The present work analyzes this procedure and studies a method, which solve the P1 equations directly, without definitions of slowing down densities. (author)

Nonlinear dynamics exploration through normal forms

CERN Document Server

Kahn, Peter B

2014-01-01

Geared toward advanced undergraduates and graduate students, this exposition covers the method of normal forms and its application to ordinary differential equations through perturbation analysis. In addition to its emphasis on the freedom inherent in the normal form expansion, the text features numerous examples of equations, the kind of which are encountered in many areas of science and engineering. The treatment begins with an introduction to the basic concepts underlying the normal forms. Coverage then shifts to an investigation of systems with one degree of freedom that model oscillations

A normalized PID controller in networked control systems with varying time delays.

Science.gov (United States)

Tran, Hoang-Dung; Guan, Zhi-Hong; Dang, Xuan-Kien; Cheng, Xin-Ming; Yuan, Fu-Shun

2013-09-01

It requires not only simplicity and flexibility but also high specified stability and robustness of system to design a PI/PID controller in such complicated networked control systems (NCSs) with delays. By gain and phase margins approach, this paper proposes a novel normalized PI/PID controller for NCSs based on analyzing the stability and robustness of system under the effect of network-induced delays. Specifically, We take into account the total measured network delays to formulate the gain and phase margins of the closed-loop system in the form of a set of equations. With pre-specified values of gain and phase margins, this set of equations is then solved for calculating the closed forms of control parameters which enable us to propose the normalized PI/PID controller simultaneously satisfying the following two requirements: (1) simplicity without re-solving the optimization problem for a new process, (2) high flexibility to cope with large scale of random delays and deal with many different processes in different conditions of network. Furthermore, in our method, the upper bound of random delay can be estimated to indicate the operating domain of proposed PI/PID controller. Finally, simulation results are shown to demonstrate the advantages of our proposed controller in many situations of network-induced delays. Copyright © 2013 ISA. Published by Elsevier Ltd. All rights reserved.

Normal scheme for solving the transport equation independently of spatial discretization

International Nuclear Information System (INIS)

Zamonsky, O.M.

1993-01-01

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

Differential equations a dynamical systems approach ordinary differential equations

CERN Document Server

Hubbard, John H

1991-01-01

This is a corrected third printing of the first part of the text Differential Equations: A Dynamical Systems Approach written by John Hubbard and Beverly West. The authors' main emphasis in this book is on ordinary differential equations. The book is most appropriate for upper level undergraduate and graduate students in the fields of mathematics, engineering, and applied mathematics, as well as the life sciences, physics and economics. Traditional courses on differential equations focus on techniques leading to solutions. Yet most differential equations do not admit solutions which can be written in elementary terms. The authors have taken the view that a differential equations defines functions; the object of the theory is to understand the behavior of these functions. The tools the authors use include qualitative and numerical methods besides the traditional analytic methods. The companion software, MacMath, is designed to bring these notions to life.

The Dirac equation and the normalization of its solutions in a closed Friedmann- Robertson-Walker universe

Energy Technology Data Exchange (ETDEWEB)

Finster, Felix [NWF I-Mathematik, Universitaet Regensburg, D-93040 Regensburg (Germany); Reintjes, Moritz, E-mail: Felix.Finster@mathematik.uni-regensburg.d, E-mail: moritz@math.ucdavis.ed [Mathematics Department, University of California, Davis, CA 95616 (United States)

2009-05-21

We set up the Dirac equation in a Friedmann-Robertson-Walker geometry and separate the spatial and time variables. In the case of a closed universe, the spatial dependence is solved explicitly, giving rise to a discrete set of solutions. We compute the probability integral and analyze a spacetime normalization integral. This analysis allows us to introduce the fermionic projector in a closed Friedmann-Robertson-Walker geometry and to specify its global normalization as well as its local form.

The Dirac equation and the normalization of its solutions in a closed Friedmann- Robertson-Walker universe

International Nuclear Information System (INIS)

Finster, Felix; Reintjes, Moritz

2009-01-01

We set up the Dirac equation in a Friedmann-Robertson-Walker geometry and separate the spatial and time variables. In the case of a closed universe, the spatial dependence is solved explicitly, giving rise to a discrete set of solutions. We compute the probability integral and analyze a spacetime normalization integral. This analysis allows us to introduce the fermionic projector in a closed Friedmann-Robertson-Walker geometry and to specify its global normalization as well as its local form.

Linear integral equations and soliton systems

International Nuclear Information System (INIS)

Quispel, G.R.W.

1983-01-01

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

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.

External validation of equations to estimate resting energy expenditure in 14952 adults with overweight and obesity and 1948 adults with normal weight from Italy.

Science.gov (United States)

Bedogni, Giorgio; Bertoli, Simona; Leone, Alessandro; De Amicis, Ramona; Lucchetti, Elisa; Agosti, Fiorenza; Marazzi, Nicoletta; Battezzati, Alberto; Sartorio, Alessandro

2017-11-24

We cross-validated 28 equations to estimate resting energy expenditure (REE) in a very large sample of adults with overweight or obesity. 14952 Caucasian men and women with overweight or obesity and 1498 with normal weight were studied. REE was measured using indirect calorimetry and estimated using two meta-regression equations and 26 other equations. The correct classification fraction (CCF) was defined as the fraction of subjects whose estimated REE was within 10% of measured REE. The highest CCF was 79%, 80%, 72%, 64%, and 63% in subjects with normal weight, overweight, class 1 obesity, class 2 obesity, and class 3 obesity, respectively. The Henry weight and height and Mifflin equations performed equally well with CCFs of 77% vs. 77% for subjects with normal weight, 80% vs. 80% for those with overweight, 72% vs. 72% for those with class 1 obesity, 64% vs. 63% for those with class 2 obesity, and 61% vs. 60% for those with class 3 obesity. The Sabounchi meta-regression equations offered an improvement over the above equations only for class 3 obesity (63%). The accuracy of REE equations decreases with increasing values of body mass index. The Henry weight & height and Mifflin equations are similarly accurate and the Sabounchi equations offer an improvement only in subjects with class 3 obesity. Copyright © 2017 Elsevier Ltd and European Society for Clinical Nutrition and Metabolism. All rights reserved.

Lorentz-force equations as Heisenberg equations for a quantum system in the euclidean space

International Nuclear Information System (INIS)

Rodriguez D, R.

2007-01-01

In an earlier work, the dynamic equations for a relativistic charged particle under the action of electromagnetic fields were formulated by R. Yamaleev in terms of external, as well as internal momenta. Evolution equations for external momenta, the Lorentz-force equations, were derived from the evolution equations for internal momenta. The mapping between the observables of external and internal momenta are related by Viete formulae for a quadratic polynomial, the characteristic polynomial of the relativistic dynamics. In this paper we show that the system of dynamic equations, can be cast into the Heisenberg scheme for a four-dimensional quantum system. Within this scheme the equations in terms of internal momenta play the role of evolution equations for a state vector, whereas the external momenta obey the Heisenberg equation for an operator evolution. The solutions of the Lorentz-force equation for the motion inside constant electromagnetic fields are presented via pentagonometric functions. (Author)

The Dirac equation and the normalization of its solutions in a closed Friedmann- Robertson-Walker universe

Science.gov (United States)

Finster, Felix; Reintjes, Moritz

2009-05-01

We set up the Dirac equation in a Friedmann-Robertson-Walker geometry and separate the spatial and time variables. In the case of a closed universe, the spatial dependence is solved explicitly, giving rise to a discrete set of solutions. We compute the probability integral and analyze a spacetime normalization integral. This analysis allows us to introduce the fermionic projector in a closed Friedmann-Robertson-Walker geometry and to specify its global normalization as well as its local form. First author supported in part by the Deutsche Forschungsgemeinschaft.

Structural Equation and Mei Conserved Quantity of Mei Symmetry for Appell Equations in Holonomic Systems with Unilateral Constraints

International Nuclear Information System (INIS)

Jia Liqun; Cui Jinchao; Zhang Yaoyu; Luo Shaokai

2009-01-01

Structural equation and Mei conserved quantity of Mei symmetry for Appell equations in holonomic systems with unilateral constraints are investigated. Appell equations and differential equations of motion for holonomic mechanic systems with unilateral constraints are established. The definition and the criterion of Mei symmetry for Appell equations in holonomic systems with unilateral constraints under the infinitesimal transformations of groups are also given. The expressions of the structural equation and Mei conserved quantity of Mei symmetry for Appell equations in holonomic systems with unilateral constraints expressed by Appell functions are obtained. An example is given to illustrate the application of the results. (general)

Soliton equations and Hamiltonian systems

CERN Document Server

Dickey, L A

2002-01-01

The theory of soliton equations and integrable systems has developed rapidly during the last 30 years with numerous applications in mechanics and physics. For a long time, books in this field have not been written but the flood of papers was overwhelming: many hundreds, maybe thousands of them. All this output followed one single work by Gardner, Green, Kruskal, and Mizura on the Korteweg-de Vries equation (KdV), which had seemed to be merely an unassuming equation of mathematical physics describing waves in shallow water. Besides its obvious practical use, this theory is attractive also becau

Planck constant as spectral parameter in integrable systems and KZB equations

Science.gov (United States)

Levin, A.; Olshanetsky, M.; Zotov, A.

2014-10-01

We construct special rational gl N Knizhnik-Zamolodchikov-Bernard (KZB) equations with Ñ punctures by deformation of the corresponding quantum gl N rational R-matrix. They have two parameters. The limit of the first one brings the model to the ordinary rational KZ equation. Another one is τ. At the level of classical mechanics the deformation parameter τ allows to extend the previously obtained modified Gaudin models to the modified Schlesinger systems. Next, we notice that the identities underlying generic (elliptic) KZB equations follow from some additional relations for the properly normalized R-matrices. The relations are noncommutative analogues of identities for (scalar) elliptic functions. The simplest one is the unitarity condition. The quadratic (in R matrices) relations are generated by noncommutative Fay identities. In particular, one can derive the quantum Yang-Baxter equations from the Fay identities. The cubic relations provide identities for the KZB equations as well as quadratic relations for the classical r-matrices which can be treated as halves of the classical Yang-Baxter equation. At last we discuss the R-matrix valued linear problems which provide gl Ñ CM models and Painlevé equations via the above mentioned identities. The role of the spectral parameter plays the Planck constant of the quantum R-matrix. When the quantum gl N R-matrix is scalar ( N = 1) the linear problem reproduces the Krichever's ansatz for the Lax matrices with spectral parameter for the gl Ñ CM models. The linear problems for the quantum CM models generalize the KZ equations in the same way as the Lax pairs with spectral parameter generalize those without it.

The flow equation approach to many-particle systems

CERN Document Server

Kehrein, Stefan; Fujimori, A; Varma, C; Steiner, F

2006-01-01

This self-contained monograph addresses the flow equation approach to many-particle systems. The flow equation approach consists of a sequence of infinitesimal unitary transformations and is conceptually similar to renormalization and scaling methods. Flow equations provide a framework for analyzing Hamiltonian systems where these conventional many-body techniques fail. The text first discusses the general ideas and concepts of the flow equation method. In a second part these concepts are illustrated with various applications in condensed matter theory including strong-coupling problems and non-equilibrium systems. The monograph is accessible to readers familiar with graduate- level solid-state theory.

Discrete ordinate solution of the radiative transfer equation in the 'polarization normal wave representation'

Science.gov (United States)

Kylling, A.

1991-01-01

The transfer equations for normal waves in finite, inhomogeneous and plane-parallel magnetoactive media are solved using the discrete ordinate method. The physical process of absorption, emission, and multiple scattering are accounted for, and the medium may be forced both at the top and bottom boundary by anisotropic radiation as well as by internal anisotropic sources. The computational procedure is numerically stable for arbitrarily large optical depths, and the computer time is independent of optical thickness.

Dirac equations for generalised Yang-Mills systems

International Nuclear Information System (INIS)

Lechtenfeld, O.; Nahm, W.; Tchrakian, D.H.

1985-06-01

We present Dirac equations in 4p dimensions for the generalised Yang-Mills (GYM) theories introduced earlier. These Dirac equations are related to the self-duality equations of the GYM and are checked to be elliptic in a 'BPST' background. In this background these Dirac equations are integrated exactly. The possibility of imposing supersymmetry in the GYM-Dirac system is investigated, with negative results. (orig.)

Algebraic limit cycles in polynomial systems of differential equations

International Nuclear Information System (INIS)

Llibre, Jaume; Zhao Yulin

2007-01-01

Using elementary tools we construct cubic polynomial systems of differential equations with algebraic limit cycles of degrees 4, 5 and 6. We also construct a cubic polynomial system of differential equations having an algebraic homoclinic loop of degree 3. Moreover, we show that there are polynomial systems of differential equations of arbitrary degree that have algebraic limit cycles of degree 3, as well as give an example of a cubic polynomial system of differential equations with two algebraic limit cycles of degree 4

NIMROD: a program for inference via a normal approximation of the posterior in models with random effects based on ordinary differential equations.

Science.gov (United States)

Prague, Mélanie; Commenges, Daniel; Guedj, Jérémie; Drylewicz, Julia; Thiébaut, Rodolphe

2013-08-01

Models based on ordinary differential equations (ODE) are widespread tools for describing dynamical systems. In biomedical sciences, data from each subject can be sparse making difficult to precisely estimate individual parameters by standard non-linear regression but information can often be gained from between-subjects variability. This makes natural the use of mixed-effects models to estimate population parameters. Although the maximum likelihood approach is a valuable option, identifiability issues favour Bayesian approaches which can incorporate prior knowledge in a flexible way. However, the combination of difficulties coming from the ODE system and from the presence of random effects raises a major numerical challenge. Computations can be simplified by making a normal approximation of the posterior to find the maximum of the posterior distribution (MAP). Here we present the NIMROD program (normal approximation inference in models with random effects based on ordinary differential equations) devoted to the MAP estimation in ODE models. We describe the specific implemented features such as convergence criteria and an approximation of the leave-one-out cross-validation to assess the model quality of fit. In pharmacokinetics models, first, we evaluate the properties of this algorithm and compare it with FOCE and MCMC algorithms in simulations. Then, we illustrate NIMROD use on Amprenavir pharmacokinetics data from the PUZZLE clinical trial in HIV infected patients. Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.

Introduction to differential equations with dynamical systems

CERN Document Server

Campbell, Stephen L

2011-01-01

Many textbooks on differential equations are written to be interesting to the teacher rather than the student. Introduction to Differential Equations with Dynamical Systems is directed toward students. This concise and up-to-date textbook addresses the challenges that undergraduate mathematics, engineering, and science students experience during a first course on differential equations. And, while covering all the standard parts of the subject, the book emphasizes linear constant coefficient equations and applications, including the topics essential to engineering students. Stephen Campbell and Richard Haberman--using carefully worded derivations, elementary explanations, and examples, exercises, and figures rather than theorems and proofs--have written a book that makes learning and teaching differential equations easier and more relevant. The book also presents elementary dynamical systems in a unique and flexible way that is suitable for all courses, regardless of length.

Nonlinear von Neumann equations for quantum dissipative systems

International Nuclear Information System (INIS)

Messer, J.; Baumgartner, B.

1978-01-01

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

Nonlinear von Neumann equations for quantum dissipative systems

International Nuclear Information System (INIS)

Messer, J.; Baumgartner, B.

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

ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY

OpenAIRE

Enrique Gonzalo Reyes Garcia

2004-01-01

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

Lie symmetries for systems of evolution equations

Science.gov (United States)

Paliathanasis, Andronikos; Tsamparlis, Michael

2018-01-01

The Lie symmetries for a class of systems of evolution equations are studied. The evolution equations are defined in a bimetric space with two Riemannian metrics corresponding to the space of the independent and dependent variables of the differential equations. The exact relation of the Lie symmetries with the collineations of the bimetric space is determined.

Comparison of SUVs normalized by lean body mass determined by CT with those normalized by lean body mass estimated by predictive equations in normal tissues

Energy Technology Data Exchange (ETDEWEB)

Kim, Woo Hyoung; Kim, Chang Guhn; Kim, Dae Weung [Wonkwang Univ. School of Medicine, Iksan (Korea, Republic of)

2012-09-15

Standardized uptake values (SUVs)normalized by lean body mass (LBM)determined by CT were compared with those normalized by LBM estimated using predictive equations (PEs)in normal liver, spleen, and aorta using {sup 18}F FDG PET/CT. Fluorine 18 fluorodeoxyglucose (F FDG)positron emission tomography/computed tomography (PET/CT)was conducted on 453 patients. LBM determined by CT was defined in 3 ways (LBM{sup CT1}-3). Five PEs were used for comparison (LBM{sup PE1}-5). Tissue SUV normalized by LBM (SUL) was calculated using LBM from each method (SUL{sup CT1}-3, SUL{sup PE1}-5). Agreement between methods was assessed by Bland Altman analysis. Percentage difference and percentage error were also calculated. For all liver SUL{sup CTS} vs. liver SUL{sup PES} except liver SUL{sup PE3}, the range of biases, SDs of percentage difference and percentage errors were -0.17-0.24 SUL, 6.15-10.17%, and 25.07-38.91%, respectively. For liver SUL{sup CTs} vs. liver SUL{sup PE3}, the corresponding figures were 0.47-0.69 SUL, 10.90-11.25%, and 50.85-51.55%, respectively, showing the largest percentage errors and positive biases. Irrespective of magnitudes of the biases, large percentage errors of 25.07-51.55% were observed between liver SUL{sup CT1}-3 and liver SUL{sup PE1}-5. The results of spleen and aorta SUL{sup CTs} and SUL{sup PEs} comparison were almost identical to those for liver. The present study demonstrated substantial errors in individual SUL{sup PEs} compared with SUL{sup CTs} as a reference value. Normalization of SUV by LBM determined by CT rather than PEs may be a useful approach to reduce errors in individual SUL{sup PEs}.

Comparison of SUVs normalized by lean body mass determined by CT with those normalized by lean body mass estimated by predictive equations in normal tissues

International Nuclear Information System (INIS)

Kim, Woo Hyoung; Kim, Chang Guhn; Kim, Dae Weung

2012-01-01

Standardized uptake values (SUVs)normalized by lean body mass (LBM)determined by CT were compared with those normalized by LBM estimated using predictive equations (PEs)in normal liver, spleen, and aorta using 18 F FDG PET/CT. Fluorine 18 fluorodeoxyglucose (F FDG)positron emission tomography/computed tomography (PET/CT)was conducted on 453 patients. LBM determined by CT was defined in 3 ways (LBM CT1 -3). Five PEs were used for comparison (LBM PE1 -5). Tissue SUV normalized by LBM (SUL) was calculated using LBM from each method (SUL CT1 -3, SUL PE1 -5). Agreement between methods was assessed by Bland Altman analysis. Percentage difference and percentage error were also calculated. For all liver SUL CTS vs. liver SUL PES except liver SUL PE3 , the range of biases, SDs of percentage difference and percentage errors were -0.17-0.24 SUL, 6.15-10.17%, and 25.07-38.91%, respectively. For liver SUL CTs vs. liver SUL PE3 , the corresponding figures were 0.47-0.69 SUL, 10.90-11.25%, and 50.85-51.55%, respectively, showing the largest percentage errors and positive biases. Irrespective of magnitudes of the biases, large percentage errors of 25.07-51.55% were observed between liver SUL CT1 -3 and liver SUL PE1 -5. The results of spleen and aorta SUL CTs and SUL PEs comparison were almost identical to those for liver. The present study demonstrated substantial errors in individual SUL PEs compared with SUL CTs as a reference value. Normalization of SUV by LBM determined by CT rather than PEs may be a useful approach to reduce errors in individual SUL PEs

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.

Solution of Large Systems of Linear Equations with Quadratic or Non-Quadratic Matrices and Deconvoiution of Spectra

Energy Technology Data Exchange (ETDEWEB)

Nygaard, K

1967-12-15

The numerical deconvolution of spectra is equivalent to the solution of a (large) system of linear equations with a matrix which is not necessarily a square matrix. The demand that the square sum of the residual errors shall be minimum is not in general sufficient to ensure a unique or 'sound' solution. Therefore other demands which may include the demand for minimum square errors are introduced which lead to 'sound' and 'non-oscillatory' solutions irrespective of the shape of the original matrix and of the determinant of the matrix of the normal equations.

Solving differential–algebraic equation systems by means of index reduction methodology

DEFF Research Database (Denmark)

Sørensen, Kim; Houbak, Niels; Condra, Thomas

2006-01-01

of a number of differential equations and algebraic equations — a so called DAE system. Two of the DAE systems are of index 1 and they can be solved by means of standard DAE-solvers. For the actual application, the equation systems are integrated by means of MATLAB’s solver: ode23t, that solves moderately...... stiff ODEs and index 1 DAEs by means of the trapezoidal rule. The last sub-model that models the boilers steam drum consist of two differential and three algebraic equations. The index of this model is greater than 1, which means that ode23t cannot integrate this equation system. In this paper......, it is shown how the equation system, by means of an index reduction methodology, can be reduced to a system of ordinary differential equations — ODEs....

The 'strength' of a system of differential equations

International Nuclear Information System (INIS)

Hoenselaers, C.

1977-01-01

A review of Einstein's concept of ''strength'' of a system of differential equations is given. As an example the strength of the Einstein-Maxwell equations for non-null Maxwell field is calculated and shown to be the same as for the pure vacuum Einstein equations. (auth.)

Calculation of normal modes of the closed waveguides in general vector case

Science.gov (United States)

Malykh, M. D.; Sevastianov, L. A.; Tiutiunnik, A. A.

2018-04-01

The article is devoted to the calculation of normal modes of the closed waveguides with an arbitrary filling ɛ, μ in the system of computer algebra Sage. Maxwell equations in the cylinder are reduced to the system of two bounded Helmholtz equations, the notion of weak solution of this system is given and then this system is investigated as a system of ordinary differential equations. The normal modes of this system are an eigenvectors of a matrix pencil. We suggest to calculate the matrix elements approximately and to truncate the matrix by usual way but further to solve the truncated eigenvalue problem exactly in the field of algebraic numbers. This approach allows to keep the symmetry of the initial problem and in particular the multiplicity of the eigenvalues. In the work would be presented some results of calculations.

Two-fluid equations for a nuclear system with arbitrary motions

Energy Technology Data Exchange (ETDEWEB)

Kim, Byoung Jae [Chungnam National University, Daejeon (Korea, Republic of); Kim, Kyung Doo [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)

2016-10-15

Ocean nuclear systems include a seabed-type plant, a floating-type plant, and a nuclear-propulsion ship. We asked ourselves, 'What governing equations should be used for ocean nuclear systems?' Since ocean nuclear systems are apt to move arbitrarily, the two-fluid model must be formulated in the non-inertial frame of reference that is undergoing acceleration with respect to an inertial frame. Two-phase flow systems with arbitrary motions are barely reported. Kim et al. (1996) added the centripetal and Euler acceleration forces to the homogeneous equilibrium momentum equation embedded in the RETRAN code. However, they did not look into the mass and energy equations. The purpose of this study is to derive general two-fluid equations in the non-inertial frame of reference, which can be used for safety analysis of ocean nuclear systems. The two-fluid equation forms for scalar properties such as mass, internal energy, and enthalpy equation in the moving frame are the same as those in the absolute frame. On the other hand, the fictitious effect must be included in the momentum equation.

Zero singularities of codimension two and three in delay differential equations

International Nuclear Information System (INIS)

Campbell, Sue Ann; Yuan Yuan

2008-01-01

We give conditions under which a general class of delay differential equations has a point of Bogdanov–Takens or a triple zero bifurcation. We show how a centre manifold projection of the delay equations reduces the dynamics to two- or three-dimensional systems of ordinary differential equations. We put these equations in normal form and determine how the coefficients of the normal forms depend on the original parameters in the model. Finally we apply our results to two neural models and compare the predictions of the theory with numerical bifurcation analysis of the full equations. One model involves a transcritical bifurcation, hence we derive and analyse the appropriate unfoldings for this case

a Recursive Approach to Compute Normal Forms

Science.gov (United States)

HSU, L.; MIN, L. J.; FAVRETTO, L.

2001-06-01

Normal forms are instrumental in the analysis of dynamical systems described by ordinary differential equations, particularly when singularities close to a bifurcation are to be characterized. However, the computation of a normal form up to an arbitrary order is numerically hard. This paper focuses on the computer programming of some recursive formulas developed earlier to compute higher order normal forms. A computer program to reduce the system to its normal form on a center manifold is developed using the Maple symbolic language. However, it should be stressed that the program relies essentially on recursive numerical computations, while symbolic calculations are used only for minor tasks. Some strategies are proposed to save computation time. Examples are presented to illustrate the application of the program to obtain high order normalization or to handle systems with large dimension.

Solving differential-algebraic equation systems by means of index reduction methodology

DEFF Research Database (Denmark)

Sørensen, Kim; Houbak, Niels; Condra, Thomas Joseph

2006-01-01

of a number of differential equations and algebraic equations - a so called DAE system. Two of the DAE systems are of index 1 and they can be solved by means of standard DAE-solvers. For the actual application, the equation systems are integrated by means of MATLAB’s solver: ode23t, that solves moderately...... stiff ODE’s and index 1 DAE’s by means of the trapezoidal rule. The last sub-model that models the boilers steam drum consist of two differential and three algebraic equations. The index of this model is greater than 1, which means that ode23t cannot integrate this equation system. In this paper......, it is shown how the equation system, by means of an index reduction methodology, can be reduced to a system of Ordinary- Differential-Equations - ODE’s....

The respiratory system in equations

CERN Document Server

Maury, Bertrand

2013-01-01

The book proposes an introduction to the mathematical modeling of the respiratory system. A detailed introduction on the physiological aspects makes it accessible to a large audience without any prior knowledge on the lung. Different levels of description are proposed, from the lumped models with a small number of parameters (Ordinary Differential Equations), up to infinite dimensional models based on Partial Differential Equations. Besides these two types of differential equations, two chapters are dedicated to resistive networks, and to the way they can be used to investigate the dependence of the resistance of the lung upon geometrical characteristics. The theoretical analysis of the various models is provided, together with state-of-the-art techniques to compute approximate solutions, allowing comparisons with experimental measurements. The book contains several exercises, most of which are accessible to advanced undergraduate students.

The method of normal forms for singularly perturbed systems of Fredholm integro-differential equations with rapidly varying kernels

Energy Technology Data Exchange (ETDEWEB)

Bobodzhanov, A A; Safonov, V F [National Research University " Moscow Power Engineering Institute" , Moscow (Russian Federation)

2013-07-31

The paper deals with extending the Lomov regularization method to classes of singularly perturbed Fredholm-type integro-differential systems, which have not so far been studied. In these the limiting operator is discretely noninvertible. Such systems are commonly known as problems with unstable spectrum. Separating out the essential singularities in the solutions to these problems presents great difficulties. The principal one is to give an adequate description of the singularities induced by 'instability points' of the spectrum. A methodology for separating singularities by using normal forms is developed. It is applied to the above type of systems and is substantiated in these systems. Bibliography: 10 titles.

Particle Systems and Partial Differential Equations I

CERN Document Server

Gonçalves, Patricia

2014-01-01

This book presents the proceedings of the international conference Particle Systems and Partial Differential Equations I, which took place at the Centre of Mathematics of the University of Minho, Braga, Portugal, from the 5th to the 7th of December, 2012.  The purpose of the conference was to bring together world leaders to discuss their topics of expertise and to present some of their latest research developments in those fields. Among the participants were researchers in probability, partial differential equations and kinetics theory. The aim of the meeting was to present to a varied public the subject of interacting particle systems, its motivation from the viewpoint of physics and its relation with partial differential equations or kinetics theory, and to stimulate discussions and possibly new collaborations among researchers with different backgrounds.  The book contains lecture notes written by François Golse on the derivation of hydrodynamic equations (compressible and incompressible Euler and Navie...

Integrable coupling system of fractional soliton equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

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

2009-10-05

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

Quantum-statistical kinetic equations

International Nuclear Information System (INIS)

Loss, D.; Schoeller, H.

1989-01-01

Considering a homogeneous normal quantum fluid consisting of identical interacting fermions or bosons, the authors derive an exact quantum-statistical generalized kinetic equation with a collision operator given as explicit cluster series where exchange effects are included through renormalized Liouville operators. This new result is obtained by applying a recently developed superoperator formalism (Liouville operators, cluster expansions, symmetrized projectors, P q -rule, etc.) to nonequilibrium systems described by a density operator ρ(t) which obeys the von Neumann equation. By means of this formalism a factorization theorem is proven (being essential for obtaining closed equations), and partial resummations (leading to renormalized quantities) are performed. As an illustrative application, the quantum-statistical versions (including exchange effects due to Fermi-Dirac or Bose-Einstein statistics) of the homogeneous Boltzmann (binary collisions) and Choh-Uhlenbeck (triple collisions) equations are derived

Cross-Validation of a Recently Published Equation Predicting Energy Expenditure to Run or Walk a Mile in Normal-Weight and Overweight Adults

Science.gov (United States)

Morris, Cody E.; Owens, Scott G.; Waddell, Dwight E.; Bass, Martha A.; Bentley, John P.; Loftin, Mark

2014-01-01

An equation published by Loftin, Waddell, Robinson, and Owens (2010) was cross-validated using ten normal-weight walkers, ten overweight walkers, and ten distance runners. Energy expenditure was measured at preferred walking (normal-weight walker and overweight walkers) or running pace (distance runners) for 5 min and corrected to a mile. Energy…

How to derive biological information from the value of the normalization constant in allometric equations.

Science.gov (United States)

Kaitaniemi, Pekka

2008-04-09

Allometric equations are widely used in many branches of biological science. The potential information content of the normalization constant b in allometric equations of the form Y = bX(a) has, however, remained largely neglected. To demonstrate the potential for utilizing this information, I generated a large number of artificial datasets that resembled those that are frequently encountered in biological studies, i.e., relatively small samples including measurement error or uncontrolled variation. The value of X was allowed to vary randomly within the limits describing different data ranges, and a was set to a fixed theoretical value. The constant b was set to a range of values describing the effect of a continuous environmental variable. In addition, a normally distributed random error was added to the values of both X and Y. Two different approaches were then used to model the data. The traditional approach estimated both a and b using a regression model, whereas an alternative approach set the exponent a at its theoretical value and only estimated the value of b. Both approaches produced virtually the same model fit with less than 0.3% difference in the coefficient of determination. Only the alternative approach was able to precisely reproduce the effect of the environmental variable, which was largely lost among noise variation when using the traditional approach. The results show how the value of b can be used as a source of valuable biological information if an appropriate regression model is selected.

A generalized estimating equations approach to quantitative trait locus detection of non-normal traits

Directory of Open Access Journals (Sweden)

Thomson Peter C

2003-05-01

Full Text Available Abstract To date, most statistical developments in QTL detection methodology have been directed at continuous traits with an underlying normal distribution. This paper presents a method for QTL analysis of non-normal traits using a generalized linear mixed model approach. Development of this method has been motivated by a backcross experiment involving two inbred lines of mice that was conducted in order to locate a QTL for litter size. A Poisson regression form is used to model litter size, with allowances made for under- as well as over-dispersion, as suggested by the experimental data. In addition to fixed parity effects, random animal effects have also been included in the model. However, the method is not fully parametric as the model is specified only in terms of means, variances and covariances, and not as a full probability model. Consequently, a generalized estimating equations (GEE approach is used to fit the model. For statistical inferences, permutation tests and bootstrap procedures are used. This method is illustrated with simulated as well as experimental mouse data. Overall, the method is found to be quite reliable, and with modification, can be used for QTL detection for a range of other non-normally distributed traits.

FORSIM-6, Automatic Solution of Coupled Differential Equation System

International Nuclear Information System (INIS)

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

1983-01-01

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

Local bifurcations in differential equations with state-dependent delay.

Science.gov (United States)

Sieber, Jan

2017-11-01

A common task when analysing dynamical systems is the determination of normal forms near local bifurcations of equilibria. As most of these normal forms have been classified and analysed, finding which particular class of normal form one encounters in a numerical bifurcation study guides follow-up computations. This paper builds on normal form algorithms for equilibria of delay differential equations with constant delay that were developed and implemented in DDE-Biftool recently. We show how one can extend these methods to delay-differential equations with state-dependent delay (sd-DDEs). Since higher degrees of regularity of local center manifolds are still open for sd-DDEs, we give an independent (still only partial) argument which phenomena from the truncated normal must persist in the full sd-DDE. In particular, we show that all invariant manifolds with a sufficient degree of normal hyperbolicity predicted by the normal form exist also in the full sd-DDE.

Local bifurcations in differential equations with state-dependent delay

Science.gov (United States)

Sieber, Jan

2017-11-01

A common task when analysing dynamical systems is the determination of normal forms near local bifurcations of equilibria. As most of these normal forms have been classified and analysed, finding which particular class of normal form one encounters in a numerical bifurcation study guides follow-up computations. This paper builds on normal form algorithms for equilibria of delay differential equations with constant delay that were developed and implemented in DDE-Biftool recently. We show how one can extend these methods to delay-differential equations with state-dependent delay (sd-DDEs). Since higher degrees of regularity of local center manifolds are still open for sd-DDEs, we give an independent (still only partial) argument which phenomena from the truncated normal must persist in the full sd-DDE. In particular, we show that all invariant manifolds with a sufficient degree of normal hyperbolicity predicted by the normal form exist also in the full sd-DDE.

Method for the determination of the equation of state of advanced fuels based on the properties of normal fluids

International Nuclear Information System (INIS)

Hecht, M.J.; Catton, I.; Kastenberg, W.E.

1976-12-01

An equation of state based on the properties of normal fluids, the law of rectilinear averages, and the second law of thermodynamics can be derived for advanced LMFBR fuels on the basis of the vapor pressure, enthalpy of vaporization, change in heat capacity upon vaporization, and liquid density at the melting point. The method consists of estimating an equation of state by means of the law of rectilinear averages and the second law of thermodynamics, integrating by means of the second law until an instability is reached, and then extrapolating by means of a self-consistent estimation of the enthalpy of vaporization

Structural equation modeling and natural systems

Science.gov (United States)

Grace, James B.

2006-01-01

This book, first published in 2006, presents an introduction to the methodology of structural equation modeling, illustrates its use, and goes on to argue that it has revolutionary implications for the study of natural systems. A major theme of this book is that we have, up to this point, attempted to study systems primarily using methods (such as the univariate model) that were designed only for considering individual processes. Understanding systems requires the capacity to examine simultaneous influences and responses. Structural equation modeling (SEM) has such capabilities. It also possesses many other traits that add strength to its utility as a means of making scientific progress. In light of the capabilities of SEM, it can be argued that much of ecological theory is currently locked in an immature state that impairs its relevance. It is further argued that the principles of SEM are capable of leading to the development and evaluation of multivariate theories of the sort vitally needed for the conservation of natural systems.

A New Algorithm for System of Integral Equations

Directory of Open Access Journals (Sweden)

Abdujabar Rasulov

2014-01-01

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

Microscopic coefficients for the quantum master equation of a Fermi system

International Nuclear Information System (INIS)

Stefanescu, E.; Sandulescu, A.

2002-01-01

In a previous paper, we derived a master equation for fermions, of Lindblad's form, with coefficients depending on microscopic quantities. In this paper, we study the properties of the dissipative coefficients taking into account the explicit expressions of: (a) the matrix elements of the dissipative potential, evaluated from the condition that, essentially, this potential induces transitions among the system eigenstates without significantly modifying these states, (b) the densities of the environment states according to the Thomas-Fermi model, and (c) the occupation probabilities of these states taken as a Fermi-Dirac distribution. The matrix of these coefficients correctly describes the system dynamics: (a) for a normal, Fermi-Dirac distribution of the environment population, the decays dominate the excitation processes; (b) for an inverted (exotic) distribution of this population, specific to a clustering state, the excitation processes are dominant. (author)

Stochastic equations for complex systems theoretical and computational topics

CERN Document Server

Bessaih, Hakima

2015-01-01

Mathematical analyses and computational predictions of the behavior of complex systems are needed to effectively deal with weather and climate predictions, for example, and the optimal design of technical processes. Given the random nature of such systems and the recognized relevance of randomness, the equations used to describe such systems usually need to involve stochastics.  The basic goal of this book is to introduce the mathematics and application of stochastic equations used for the modeling of complex systems. A first focus is on the introduction to different topics in mathematical analysis. A second focus is on the application of mathematical tools to the analysis of stochastic equations. A third focus is on the development and application of stochastic methods to simulate turbulent flows as seen in reality.  This book is primarily oriented towards mathematics and engineering PhD students, young and experienced researchers, and professionals working in the area of stochastic differential equations ...

Non-Archimedean reaction-ultradiffusion equations and complex hierarchic systems

Science.gov (United States)

Zúñiga-Galindo, W. A.

2018-06-01

We initiate the study of non-Archimedean reaction-ultradiffusion equations and their connections with models of complex hierarchic systems. From a mathematical perspective, the equations studied here are the p-adic counterpart of the integro-differential models for phase separation introduced by Bates and Chmaj. Our equations are also generalizations of the ultradiffusion equations on trees studied in the 1980s by Ogielski, Stein, Bachas, Huberman, among others, and also generalizations of the master equations of the Avetisov et al models, which describe certain complex hierarchic systems. From a physical perspective, our equations are gradient flows of non-Archimedean free energy functionals and their solutions describe the macroscopic density profile of a bistable material whose space of states has an ultrametric structure. Some of our results are p-adic analogs of some well-known results in the Archimedean setting, however, the mechanism of diffusion is completely different due to the fact that it occurs in an ultrametric space.

Methods of mathematical modelling continuous systems and differential equations

CERN Document Server

Witelski, Thomas

2015-01-01

This book presents mathematical modelling and the integrated process of formulating sets of equations to describe real-world problems. It describes methods for obtaining solutions of challenging differential equations stemming from problems in areas such as chemical reactions, population dynamics, mechanical systems, and fluid mechanics. Chapters 1 to 4 cover essential topics in ordinary differential equations, transport equations and the calculus of variations that are important for formulating models. Chapters 5 to 11 then develop more advanced techniques including similarity solutions, matched asymptotic expansions, multiple scale analysis, long-wave models, and fast/slow dynamical systems. Methods of Mathematical Modelling will be useful for advanced undergraduate or beginning graduate students in applied mathematics, engineering and other applied sciences.

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

International Nuclear Information System (INIS)

Chau Wang, L.L.

1981-01-01

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

MHD stability properties of a system of reduced toroidal MHD equations

International Nuclear Information System (INIS)

Maschke, E.K.; Morros Tosas, J.; Urquijo, G.

1993-01-01

A system of reduced toroidal magneto-hydrodynamic (MHD) equations is derived from a general scalar representation of the complete MHD system, using an ordering in terms of the inverse aspect ratio ε of a toroidal plasma. It is shown that the energy principle for the reduced equations is identical with the usual energy principle of the complete MHD system, to the appropriate order in ε. Thus, the reduced equations have the same ideal MHD stability limits as the full MHD equations. (authors). 6 refs

Numerical analysis of systems of ordinary and stochastic differential equations

CERN Document Server

Artemiev, S S

1997-01-01

This text deals with numerical analysis of systems of both ordinary and stochastic differential equations. It covers numerical solution problems of the Cauchy problem for stiff ordinary differential equations (ODE) systems by Rosenbrock-type methods (RTMs).

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

International Nuclear Information System (INIS)

Chen Jinbing; Qiao Zhijun

2011-01-01

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

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

Using Difference Equation to Model Discrete-time Behavior in System Dynamics Modeling

NARCIS (Netherlands)

Hesan, R.; Ghorbani, A.; Dignum, M.V.

2014-01-01

In system dynamics modeling, differential equations have been used as the basic mathematical operator. Using difference equation to build system dynamics models instead of differential equation, can be insightful for studying small organizations or systems with micro behavior. In this paper we

Equations of motion for free-flight systems of rotating-translating bodies

International Nuclear Information System (INIS)

Hodapp, A.E. Jr.

1976-09-01

General vector differential equations of motion are developed for a system of rotating-translating, unbalanced, constant mass bodies. Complete flexibility is provided in placement of the reference coordinates within the system of bodies and in placement of body fixed axes within each body. Example cases are presented to demonstrate the ease in reduction of these equations to the equations of motion for systems of interest

Normal people working in normal organizations with normal equipment: system safety and cognition in a mid-air collision.

Science.gov (United States)

de Carvalho, Paulo Victor Rodrigues; Gomes, José Orlando; Huber, Gilbert Jacob; Vidal, Mario Cesar

2009-05-01

A fundamental challenge in improving the safety of complex systems is to understand how accidents emerge in normal working situations, with equipment functioning normally in normally structured organizations. We present a field study of the en route mid-air collision between a commercial carrier and an executive jet, in the clear afternoon Amazon sky in which 154 people lost their lives, that illustrates one response to this challenge. Our focus was on how and why the several safety barriers of a well structured air traffic system melted down enabling the occurrence of this tragedy, without any catastrophic component failure, and in a situation where everything was functioning normally. We identify strong consistencies and feedbacks regarding factors of system day-to-day functioning that made monitoring and awareness difficult, and the cognitive strategies that operators have developed to deal with overall system behavior. These findings emphasize the active problem-solving behavior needed in air traffic control work, and highlight how the day-to-day functioning of the system can jeopardize such behavior. An immediate consequence is that safety managers and engineers should review their traditional safety approach and accident models based on equipment failure probability, linear combinations of failures, rules and procedures, and human errors, to deal with complex patterns of coincidence possibilities, unexpected links, resonance among system functions and activities, and system cognition.

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

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

Exact non-Markovian master equations for multiple qubit systems: Quantum-trajectory approach

Science.gov (United States)

Chen, Yusui; You, J. Q.; Yu, Ting

2014-11-01

A wide class of exact master equations for a multiple qubit system can be explicitly constructed by using the corresponding exact non-Markovian quantum-state diffusion equations. These exact master equations arise naturally from the quantum decoherence dynamics of qubit system as a quantum memory coupled to a collective colored noisy source. The exact master equations are also important in optimal quantum control, quantum dissipation, and quantum thermodynamics. In this paper, we show that the exact non-Markovian master equation for a dissipative N -qubit system can be derived explicitly from the statistical average of the corresponding non-Markovian quantum trajectories. We illustrated our general formulation by an explicit construction of a three-qubit system coupled to a non-Markovian bosonic environment. This multiple qubit master equation offers an accurate time evolution of quantum systems in various domains, and paves the way to investigate the memory effect of an open system in a non-Markovian regime without any approximation.

Normal-zone detectors for the MFTF-B coils. Revision 1

International Nuclear Information System (INIS)

Owen, E.W.; Shimer, D.W.

1983-01-01

In order to protect a set of inductively coupled superconducting magnets, it is necessary to locate and measure normal zone voltages that are small compared with the mutual and self-induced voltages. The method described in this report uses two sets of voltage measurements to locate and measure one or more normal zones in any number of coupled coils. One set of voltages is the outputs of bridges that balance out the self-induced voltages. The other set of voltages can be the voltages across the coils, although alternatives are possible. The two sets of equations form a single combined set of equations. Each normal zone location or combination of normal zones has a set of these combined equations associated with it. It is demonstrated that the normal zone can be located and the correct set chosen, allowing determination of the size of the normal zone. Only a few operations take place in a working detector: multiplication of a constant, addition, and simple decision-making. In many cases the detector for each coil, although weakly linked to the other detectors, can be considered to be independent. An example of the detector design is given for four coils with realistic parameters. The effect on accuracy of changes in the system parameters is discussed

Nonlinear Schroedinger Approximations for Partial Differential Equations with Quadratic and Quasilinear Terms

Science.gov (United States)

Cummings, Patrick

We consider the approximation of solutions of two complicated, physical systems via the nonlinear Schrodinger equation (NLS). In particular, we discuss the evolution of wave packets and long waves in two physical models. Due to the complicated nature of the equations governing many physical systems and the in-depth knowledge we have for solutions of the nonlinear Schrodinger equation, it is advantageous to use approximation results of this kind to model these physical systems. The approximations are simple enough that we can use them to understand the qualitative and quantitative behavior of the solutions, and by justifying them we can show that the behavior of the approximation captures the behavior of solutions to the original equation, at least for long, but finite time. We first consider a model of the water wave equations which can be approximated by wave packets using the NLS equation. We discuss a new proof that both simplifies and strengthens previous justification results of Schneider and Wayne. Rather than using analytic norms, as was done by Schneider and Wayne, we construct a modified energy functional so that the approximation holds for the full interval of existence of the approximate NLS solution as opposed to a subinterval (as is seen in the analytic case). Furthermore, the proof avoids problems associated with inverting the normal form transform by working with a modified energy functional motivated by Craig and Hunter et al. We then consider the Klein-Gordon-Zakharov system and prove a long wave approximation result. In this case there is a non-trivial resonance that cannot be eliminated via a normal form transform. By combining the normal form transform for small Fourier modes and using analytic norms elsewhere, we can get a justification result on the order 1 over epsilon squared time scale.

A hierarchy of systems of nonlinear equations

International Nuclear Information System (INIS)

Falkensteiner, P.; Grosse, H.

1985-01-01

Imposing isospectral invariance for the one-dimensional Dirac operator yields an infinite hierarchy of systems of chiral invariant nonlinear partial differential equations. The same system is obtained through a Lax pair construction and finally a formulation in terms of Kac-Moody generators is given. (Author)

The Schroedinger-Newton equation as model of self-gravitating quantum systems

International Nuclear Information System (INIS)

Grossardt, Andre

2013-01-01

The Schroedinger-Newton equation (SN equation) describes a quantummechanical one-particle-system with gravitational self-interaction and might play a role answering the question if gravity must be quantised. As non-relativistic limit of semi-classical gravity, it provides testable predictions of the effects that classical gravity has on genuinely quantum mechanical systems in the mass regime between a few thousand proton masses and the Planck mass, which is experimentally unexplored. In this thesis I subsume the mathematical properties of the SN equation and justify it as a physical model. I will give a short outline of the controversial debate around semi-classical gravity as a fundamental theory, along with the idea of the SN equation as a model of quantum state reduction. Subsequently, I will respond to frequent objections against nonlinear Schrodinger equations. I will show how the SN equation can be obtained from Einstein's General Relativity coupled to either a KleinGordon or a Dirac equation, in the same sense as the linear Schroedinger equation can be derived in flat Minkowski space-time. The equation is, to this effect, a non-relativistic approximation of the semi-classical Einstein equations. Additionally, I will discuss, first by means of analytic estimations and later numerically, in which parameter range effects of gravitational selfinteraction - e.g. in molecular-interferometry experiments - should be expected. Besides the one-particle SN equation I will provide justification for a modified equation describing the centre-of-mass wave-function of a many-particle system. Furthermore, for this modified equation, I will examine, numerically, the consequences for experiments. Although one arrives at the conclusion that no effects of the SN equation can be expected for masses up to six or seven orders of magnitude above those considered in contemporary molecular interferometry experiments, tests of the equation, for example in satellite experiments, seem

Modal analysis of inter-area oscillations using the theory of normal modes

Energy Technology Data Exchange (ETDEWEB)

Betancourt, R.J. [School of Electromechanical Engineering, University of Colima, Manzanillo, Col. 28860 (Mexico); Barocio, E. [CUCEI, University of Guadalajara, Guadalajara, Jal. 44480 (Mexico); Messina, A.R. [Graduate Program in Electrical Engineering, Cinvestav, Guadalajara, Jal. 45015 (Mexico); Martinez, I. [State Autonomous University of Mexico, Toluca, Edo. Mex. 50110 (Mexico)

2009-04-15

Based on the notion of normal modes in mechanical systems, a method is proposed for the analysis and characterization of oscillatory processes in power systems. The method is based on the property of invariance of modal subspaces and can be used to represent complex power system modal behavior by a set of decoupled, two-degree-of-freedom nonlinear oscillator equations. Using techniques from nonlinear mechanics, a new approach is outlined, for determining the normal modes (NMs) of motion of a general n-degree-of-freedom nonlinear system. Equations relating the normal modes and the physical velocities and displacements are developed from the linearized system model and numerical issues associated with the application of the technique are discussed. In addition to qualitative insight, this method can be utilized in the study of nonlinear behavior and bifurcation analyses. The application of these procedures is illustrated on a planning model of the Mexican interconnected system using a quadratic nonlinear model. Specifically, the use of normal mode analysis as a basis for identifying modal parameters, including natural frequencies and damping ratios of general, linear systems with n degrees of freedom is discussed. Comparisons to conventional linear analysis techniques demonstrate the ability of the proposed technique to extract the different oscillation modes embedded in the oscillation. (author)

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.

Nonlinear perturbations of systems of partial differential equations with constant coefficients

Directory of Open Access Journals (Sweden)

Carmen J. Vanegas

2000-01-01

Full Text Available In this article, we show the existence of solutions to boundary-value problems, consisting of nonlinear systems of partial differential equations with constant coefficients. For this purpose, we use the right inverse of an associated operator and a fix point argument. As illustrations, we apply this method to Helmholtz equations and to second order systems of elliptic equations.

Normal forms of Hopf-zero singularity

International Nuclear Information System (INIS)

Gazor, Majid; Mokhtari, Fahimeh

2015-01-01

The Lie algebra generated by Hopf-zero classical normal forms is decomposed into two versal Lie subalgebras. Some dynamical properties for each subalgebra are described; one is the set of all volume-preserving conservative systems while the other is the maximal Lie algebra of nonconservative systems. This introduces a unique conservative–nonconservative decomposition for the normal form systems. There exists a Lie-subalgebra that is Lie-isomorphic to a large family of vector fields with Bogdanov–Takens singularity. This gives rise to a conclusion that the local dynamics of formal Hopf-zero singularities is well-understood by the study of Bogdanov–Takens singularities. Despite this, the normal form computations of Bogdanov–Takens and Hopf-zero singularities are independent. Thus, by assuming a quadratic nonzero condition, complete results on the simplest Hopf-zero normal forms are obtained in terms of the conservative–nonconservative decomposition. Some practical formulas are derived and the results implemented using Maple. The method has been applied on the Rössler and Kuramoto–Sivashinsky equations to demonstrate the applicability of our results. (paper)

Normal forms of Hopf-zero singularity

Science.gov (United States)

Gazor, Majid; Mokhtari, Fahimeh

2015-01-01

The Lie algebra generated by Hopf-zero classical normal forms is decomposed into two versal Lie subalgebras. Some dynamical properties for each subalgebra are described; one is the set of all volume-preserving conservative systems while the other is the maximal Lie algebra of nonconservative systems. This introduces a unique conservative-nonconservative decomposition for the normal form systems. There exists a Lie-subalgebra that is Lie-isomorphic to a large family of vector fields with Bogdanov-Takens singularity. This gives rise to a conclusion that the local dynamics of formal Hopf-zero singularities is well-understood by the study of Bogdanov-Takens singularities. Despite this, the normal form computations of Bogdanov-Takens and Hopf-zero singularities are independent. Thus, by assuming a quadratic nonzero condition, complete results on the simplest Hopf-zero normal forms are obtained in terms of the conservative-nonconservative decomposition. Some practical formulas are derived and the results implemented using Maple. The method has been applied on the Rössler and Kuramoto-Sivashinsky equations to demonstrate the applicability of our results.

On Critical Behaviour in Systems of Hamiltonian Partial Differential Equations.

Science.gov (United States)

Dubrovin, Boris; Grava, Tamara; Klein, Christian; Moro, Antonio

2015-01-01

We study the critical behaviour of solutions to weakly dispersive Hamiltonian systems considered as perturbations of elliptic and hyperbolic systems of hydrodynamic type with two components. We argue that near the critical point of gradient catastrophe of the dispersionless system, the solutions to a suitable initial value problem for the perturbed equations are approximately described by particular solutions to the Painlevé-I (P[Formula: see text]) equation or its fourth-order analogue P[Formula: see text]. As concrete examples, we discuss nonlinear Schrödinger equations in the semiclassical limit. A numerical study of these cases provides strong evidence in support of the conjecture.

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

Exact solutions for a system of nonlinear plasma fluid equations

International Nuclear Information System (INIS)

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

1991-04-01

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

Systems of quasilinear equations and their applications to gas dynamics

CERN Document Server

Roždestvenskiĭ, B L; Schulenberger, J R

1983-01-01

This book is essentially a new edition, revised and augmented by results of the last decade, of the work of the same title published in 1968 by "Nauka." It is devoted to mathematical questions of gas dynamics. Topics covered include Foundations of the Theory of Systems of Quasilinear Equations of Hyperbolic Type in Two Independent Variables; Classical and Generalized Solutions of One-Dimensional Gas Dynamics; Difference Methods for Solving the Equations of Gas Dynamics; and Generalized Solutions of Systems of Quasilinear Equations of Hyperbolic Type.

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

CERN Document Server

Aliyu, MDS

2011-01-01

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

Normal form and synchronization of strict-feedback chaotic systems

International Nuclear Information System (INIS)

Wang, Feng; Chen, Shihua; Yu Minghai; Wang Changping

2004-01-01

This study concerns the normal form and synchronization of strict-feedback chaotic systems. We prove that, any strict-feedback chaotic system can be rendered into a normal form with a invertible transform and then a design procedure to synchronize the normal form of a non-autonomous strict-feedback chaotic system is presented. This approach needs only a scalar driving signal to realize synchronization no matter how many dimensions the chaotic system contains. Furthermore, the Roessler chaotic system is taken as a concrete example to illustrate the procedure of designing without transforming a strict-feedback chaotic system into its normal form. Numerical simulations are also provided to show the effectiveness and feasibility of the developed methods

Normal zone detectors for a large number of inductively coupled coils

International Nuclear Information System (INIS)

Owen, E.W.; Shimer, D.W.

1983-01-01

In order to protect a set of inductively coupled superconducting magnets, it is necessary to locate and measure normal zone voltages that are small compared with the mutual and self-induced voltages. The method described in this report uses two sets of voltage measurements to locate and measure one or more normal zones in any number of coupled coils. One set of voltages is the outputs of bridges that balance out the self-induced voltages The other set of voltages can be the voltages across the coils, although alternatives are possible. The two sets of equations form a single combined set of equations. Each normal zone location or combination of normal zones has a set of these combined equations associated with it. It is demonstrated that the normal zone can be located and the correct set chosen, allowing determination of the size of the normal zone. Only a few operations take plae in a working detector: multiplication of a constant, addition, and simple decision-making. In many cases the detector for each coil, although weakly linked to the other detectors, can be considered to be independent. An example of the detector design is given for four coils with realistic parameters. The effect on accuracy of changes in the system parameters is discussed

Vertical discretizations for compressible Euler equation atmospheric models giving optimal representation of normal modes

International Nuclear Information System (INIS)

Thuburn, J.; Woollings, T.J.

2005-01-01

Accurate representation of different kinds of wave motion is essential for numerical models of the atmosphere, but is sensitive to details of the discretization. In this paper, numerical dispersion relations are computed for different vertical discretizations of the compressible Euler equations and compared with the analytical dispersion relation. A height coordinate, an isentropic coordinate, and a terrain-following mass-based coordinate are considered, and, for each of these, different choices of prognostic variables and grid staggerings are considered. The discretizations are categorized according to whether their dispersion relations are optimal, are near optimal, have a single zero-frequency computational mode, or are problematic in other ways. Some general understanding of the factors that affect the numerical dispersion properties is obtained: heuristic arguments concerning the normal mode structures, and the amount of averaging and coarse differencing in the finite difference scheme, are shown to be useful guides to which configurations will be optimal; the number of degrees of freedom in the discretization is shown to be an accurate guide to the existence of computational modes; there is only minor sensitivity to whether the equations for thermodynamic variables are discretized in advective form or flux form; and an accurate representation of acoustic modes is found to be a prerequisite for accurate representation of inertia-gravity modes, which, in turn, is found to be a prerequisite for accurate representation of Rossby modes

Modulation equations for spatially periodic systems: derivation and solutions

NARCIS (Netherlands)

Schielen, R.; Doelman, A.

1996-01-01

We study a class of partial dierential equations in one spatial dimension, which can be seen as model equations for the analysis of pattern formation in physical systems dened on unbounded, weakly oscillating domains. We perform a linear and weakly nonlinear stability analysis for solutions that

Conservation properties and potential systems of vorticity-type equations

International Nuclear Information System (INIS)

Cheviakov, Alexei F.

2014-01-01

Partial differential equations of the form divN=0, N t +curl M=0 involving two vector functions in R 3 depending on t, x, y, z appear in different physical contexts, including the vorticity formulation of fluid dynamics, magnetohydrodynamics (MHD) equations, and Maxwell's equations. It is shown that these equations possess an infinite family of local divergence-type conservation laws involving arbitrary functions of space and time. Moreover, it is demonstrated that the equations of interest have a rather special structure of a lower-degree (degree two) conservation law in R 4 (t,x,y,z). The corresponding potential system has a clear physical meaning. For the Maxwell's equations, it gives rise to the scalar electric and the vector magnetic potentials; for the vorticity equations of fluid dynamics, the potentialization inverts the curl operator to yield the fluid dynamics equations in primitive variables; for MHD equations, the potential equations yield a generalization of the Galas-Bogoyavlenskij potential that describes magnetic surfaces of ideal MHD equilibria. The lower-degree conservation law is further shown to yield curl-type conservation laws and determined potential equations in certain lower-dimensional settings. Examples of new nonlocal conservation laws, including an infinite family of nonlocal material conservation laws of ideal time-dependent MHD equations in 2+1 dimensions, are presented

Parametric Borel summability for some semilinear system of partial differential equations

Directory of Open Access Journals (Sweden)

Hiroshi Yamazawa

2015-01-01

Full Text Available In this paper we study the Borel summability of formal solutions with a parameter of first order semilinear system of partial differential equations with \\(n\\ independent variables. In [Singular perturbation of linear systems with a regular singularity, J. Dynam. Control. Syst. 8 (2002, 313-322], Balser and Kostov proved the Borel summability of formal solutions with respect to a singular perturbation parameter for a linear equation with one independent variable. We shall extend their results to a semilinear system of equations with general independent variables.

Some exact solutions to the potential Kadomtsev-Petviashvili equation and to a system of shallow water wave equations

International Nuclear Information System (INIS)

Inan, Ibrahim E.; Kaya, Dogan

2006-01-01

In this Letter by considering an improved tanh function method, we found some exact solutions of the potential Kadomtsev-Petviashvili equation. Some exact solutions of the system of the shallow water wave equation were also found

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

MINPACK-1, Subroutine Library for Nonlinear Equation System

International Nuclear Information System (INIS)

Garbow, Burton S.

1984-01-01

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

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

Undergraduate Students' Mental Operations in Systems of Differential Equations

Science.gov (United States)

Whitehead, Karen; Rasmussen, Chris

2003-01-01

This paper reports on research conducted to understand undergraduate students' ways of reasoning about systems of differential equations (SDEs). As part of a semester long classroom teaching experiment in a first course in differential equations, we conducted task-based interviews with six students after their study of first order differential…

Null controllability of a cascade system of Schrodinger equations

Directory of Open Access Journals (Sweden)

Marcos Lopez-Garcia

2016-03-01

Full Text Available This article presents a control problem for a cascade system of two linear N-dimensional Schrodinger equations. We address the problem of null controllability by means of a control supported in a region not satisfying the classical geometrical control condition. The proof is based on the application of a Carleman estimate with degenerate weights to each one of the equations and a careful analysis of the system in order to prove null controllability with only one control force.

Periodic Solutions of a System of Delay Differential Equations for a Small Delay

Directory of Open Access Journals (Sweden)

Adu A.M. Wasike

2002-06-01

Full Text Available We prove the existence of an asymptotically stable periodic solution of a system of delay differential equations with a small time delay t > 0. To achieve this, we transform the system of equations into a system of perturbed ordinary differential equations and then use perturbation results to show the existence of an asymptotically stable periodic solution. This approach is contingent on the fact that the system of equations with t = 0 has a stable limit cycle. We also provide a comparative study of the solutions of the original system and the perturbed system.  This comparison lays the ground for proving the existence of periodic solutions of the original system by Schauder's fixed point theorem.

On realization of nonlinear systems described by higher-order differential equations

NARCIS (Netherlands)

van der Schaft, Arjan

1987-01-01

We consider systems of smooth nonlinear differential and algebraic equations in which some of the variables are distinguished as “external variables.” The realization problem is to replace the higher-order implicit differential equations by first-order explicit differential equations and the

Normal zone detectors for a large number of inductively coupled coils. Revision 1

International Nuclear Information System (INIS)

Owen, E.W.; Shimer, D.W.

1983-01-01

In order to protect a set of inductively coupled superconducting magnets, it is necessary to locate and measure normal zone voltages that are small compared with the mutual and self-induced voltages. The method described in this paper uses two sets of voltage measurements to locate and measure one or more normal zones in any number of coupled coils. One set of voltages is the outputs of bridges that balance out the self-induced voltages. The other set of voltages can be the voltages across the coils, although alternatives are possible. The two sets of equations form a single combined set of equations. Each normal zone location or combination of normal zones has a set of these combined equations associated with it. It is demonstrated that the normal zone can be located and the correct set chosen, allowing determination of the size of the normal zone. Only a few operations take place in a working detector: multiplication of a constant, addition, and simple decision-making. In many cases the detector for each coil, although weakly linked to the other detectors, can be considered to be independent. The effect on accuracy of changes in the system parameters is discussed

Relating systems properties of the wave and the Schrödinger equation

NARCIS (Netherlands)

Zwart, Heiko J.; Le Gorrec, Yann; Maschke, B.M.

In this article we show that systems properties of the systems governed by the second order differential equation d2wdt2=−A0w and the first order differential equation dzdt=iA0z are related. This can be used to show that, for instance, exact observability of the N-dimensional wave equation implies

Topological resilience in non-normal networked systems

Science.gov (United States)

Asllani, Malbor; Carletti, Timoteo

2018-04-01

The network of interactions in complex systems strongly influences their resilience and the system capability to resist external perturbations or structural damages and to promptly recover thereafter. The phenomenon manifests itself in different domains, e.g., parasitic species invasion in ecosystems or cascade failures in human-made networks. Understanding the topological features of the networks that affect the resilience phenomenon remains a challenging goal for the design of robust complex systems. We hereby introduce the concept of non-normal networks, namely networks whose adjacency matrices are non-normal, propose a generating model, and show that such a feature can drastically change the global dynamics through an amplification of the system response to exogenous disturbances and eventually impact the system resilience. This early stage transient period can induce the formation of inhomogeneous patterns, even in systems involving a single diffusing agent, providing thus a new kind of dynamical instability complementary to the Turing one. We provide, first, an illustrative application of this result to ecology by proposing a mechanism to mute the Allee effect and, second, we propose a model of virus spreading in a population of commuters moving using a non-normal transport network, the London Tube.

Bandwidth Optimization of Normal Equation Matrix in Bundle Block Adjustment in Multi-baseline Rotational Photography

Directory of Open Access Journals (Sweden)

WANG Xiang

2016-02-01

Full Text Available A new bandwidth optimization method of normal equation matrix in bundle block adjustment in multi-baseline rotational close range photography by image index re-sorting is proposed. The equivalent exposure station of each image is calculated by its object space coverage and the relationship with other adjacent images. Then, according to the coordinate relations between equivalent exposure stations, new logical indices of all images are computed, based on which, the optimized bandwidth value can be obtained. Experimental results show that the bandwidth determined by our proposed method is significantly better than its original value, thus the operational efficiency, as well as the memory consumption of multi-baseline rotational close range photography in real-data applications, is optimized to a certain extent.

The action principle for a system of differential equations

International Nuclear Information System (INIS)

Gitman, D M; Kupriyanov, V G

2007-01-01

We consider the problem of constructing an action functional for physical systems whose classical equations of motion cannot be directly identified with Euler-Lagrange equations for an action principle. Two ways of constructing the action principle are presented. From simple consideration, we derive the necessary and sufficient conditions for the existence of a multiplier matrix which can endow a prescribed set of second-order differential equations with the structure of the Euler-Lagrange equations. An explicit form of the action is constructed if such a multiplier exists. If a given set of differential equations cannot be derived from an action principle, one can reformulate such a set in an equivalent first-order form which can always be treated as the Euler-Lagrange equations of a certain action. We construct such an action explicitly. There exists an ambiguity (not reduced to a total time derivative) in associating a Lagrange function with a given set of equations. We present a complete description of this ambiguity. The general procedure is illustrated by several examples

The action principle for a system of differential equations

Energy Technology Data Exchange (ETDEWEB)

Gitman, D M [Instituto de FIsica, Universidade de Sao Paulo (Brazil); Kupriyanov, V G [Instituto de FIsica, Universidade de Sao Paulo (Brazil)

2007-08-17

We consider the problem of constructing an action functional for physical systems whose classical equations of motion cannot be directly identified with Euler-Lagrange equations for an action principle. Two ways of constructing the action principle are presented. From simple consideration, we derive the necessary and sufficient conditions for the existence of a multiplier matrix which can endow a prescribed set of second-order differential equations with the structure of the Euler-Lagrange equations. An explicit form of the action is constructed if such a multiplier exists. If a given set of differential equations cannot be derived from an action principle, one can reformulate such a set in an equivalent first-order form which can always be treated as the Euler-Lagrange equations of a certain action. We construct such an action explicitly. There exists an ambiguity (not reduced to a total time derivative) in associating a Lagrange function with a given set of equations. We present a complete description of this ambiguity. The general procedure is illustrated by several examples.

Characteristic Equation of the Modified Smith predictor to MIMO Systems

Directory of Open Access Journals (Sweden)

Jorge A. Herrera-Cuartas

2013-11-01

Full Text Available The delay in control systems is a feature frequently in real systems due to the transport of objects or information, a series connection of multiple systems or own processing and sensors delay, among others. Recently there have been several studies to identify the external delay MIMO systems, these works are focused on identification and on-line control of MIMO systems and use a multimodel structure based on modified Smith predictor using different search method. It is clear that for the implementation of the algorithm, and to obtain the convergence and stability analysis, it is necessary to have closed-loop equations of modified Smith predictor. However, in these works is not presented the analytical procedure, not be the main object, displaying only the closed loop equations without the procedure for obtaining it. Therefore, to respond, in this paper, we present an analytical way to derive the closed-loop equations of a modified Smith predictor.  

Generation of static solutions of self-consistent system of Einstein-Maxwell equations

International Nuclear Information System (INIS)

Anchikov, A.M.; Daishev, R.A.

1988-01-01

The theorem, according to which the static solution of the self-consistent system of the Einstein-Maxwell equations is assigned to energy static solution of the Einstein equations with the arbitrary energy-momentum tensor in the right part, is proved. As a consequence of this theorem, the way of the generation of the static solutions of the self-consistent system of the Einstein-Maxwell equations with charged dust as a source of the vacuum solutions of the Einstein equations is shown

An Explicit Formulation of Singularity-Free Dynamic Equations of Mechanical Systems in Lagrangian Form---Part Two: Multibody Systems

Directory of Open Access Journals (Sweden)

Pål Johan From

2012-04-01

Full Text Available This paper presents the explicit dynamic equations of multibody mechanical systems. This is the second paper on this topic. In the first paper the dynamics of a single rigid body from the Boltzmann--Hamel equations were derived. In this paper these results are extended to also include multibody systems. We show that when quasi-velocities are used, the part of the dynamic equations that appear from the partial derivatives of the system kinematics are identical to the single rigid body case, but in addition we get terms that come from the partial derivatives of the inertia matrix, which are not present in the single rigid body case. We present for the first time the complete and correct derivation of multibody systems based on the Boltzmann--Hamel formulation of the dynamics in Lagrangian form where local position and velocity variables are used in the derivation to obtain the singularity-free dynamic equations. The final equations are written in global variables for both position and velocity. The main motivation of these papers is to allow practitioners not familiar with differential geometry to implement the dynamic equations of rigid bodies without the presence of singularities. Presenting the explicit dynamic equations also allows for more insight into the dynamic structure of the system. Another motivation is to correct some errors commonly found in the literature. Unfortunately, the formulation of the Boltzmann-Hamel equations used here are presented incorrectly. This has been corrected by the authors, but we present here, for the first time, the detailed mathematical details on how to arrive at the correct equations. We also show through examples that using the equations presented here, the dynamics of a single rigid body is reduced to the standard equations on a Lagrangian form, for example Euler's equations for rotational motion and Euler--Lagrange equations for free motion.

Normal zone soliton in large composite superconductors

International Nuclear Information System (INIS)

Kupferman, R.; Mints, R.G.; Ben-Jacob, E.

1992-01-01

The study of normal zone of finite size (normal domains) in superconductors, has been continuously a subject of interest in the field of applied superconductivity. It was shown that in homogeneous superconductors normal domains are always unstable, so that if a normal domain nucleates, it will either expand or shrink. While testing the stability of large cryostable composite superconductors, a new phenomena was found, the existence of stable propagating normal solitons. The formation of these propagating domains was shown to be a result of the high Joule power generated in the superconductor during the relatively long process of current redistribution between the superconductor and the stabilizer. Theoretical studies were performed in investigate the propagation of normal domains in large composite super conductors in the cryostable regime. Huang and Eyssa performed numerical calculations simulating the diffusion of heat and current redistribution in the conductor, and showed the existence of stable propagating normal domains. They compared the velocity of normal domain propagation with the experimental data, obtaining a reasonable agreement. Dresner presented an analytical method to solve this problem if the time dependence of the Joule power is given. He performed explicit calculations of normal domain velocity assuming that the Joule power decays exponentially during the process of current redistribution. In this paper, the authors propose a system of two one-dimensional diffusion equations describing the dynamics of the temperature and the current density distributions along the conductor. Numerical simulations of the equations reconfirm the existence of propagating domains in the cryostable regime, while an analytical investigation supplies an explicit formula for the velocity of the normal domain

Transforming parts of a differential equations system to difference equations as a method for run-time savings in NONMEM.

Science.gov (United States)

Petersson, K J F; Friberg, L E; Karlsson, M O

2010-10-01

Computer models of biological systems grow more complex as computing power increase. Often these models are defined as differential equations and no analytical solutions exist. Numerical integration is used to approximate the solution; this can be computationally intensive, time consuming and be a large proportion of the total computer runtime. The performance of different integration methods depend on the mathematical properties of the differential equations system at hand. In this paper we investigate the possibility of runtime gains by calculating parts of or the whole differential equations system at given time intervals, outside of the differential equations solver. This approach was tested on nine models defined as differential equations with the goal to reduce runtime while maintaining model fit, based on the objective function value. The software used was NONMEM. In four models the computational runtime was successfully reduced (by 59-96%). The differences in parameter estimates, compared to using only the differential equations solver were less than 12% for all fixed effects parameters. For the variance parameters, estimates were within 10% for the majority of the parameters. Population and individual predictions were similar and the differences in OFV were between 1 and -14 units. When computational runtime seriously affects the usefulness of a model we suggest evaluating this approach for repetitive elements of model building and evaluation such as covariate inclusions or bootstraps.

Adiabatically steered open quantum systems: Master equation and optimal phase

International Nuclear Information System (INIS)

Salmilehto, J.; Solinas, P.; Ankerhold, J.; Moettoenen, M.

2010-01-01

We introduce an alternative way to derive the generalized form of the master equation recently presented by J. P. Pekola et al. [Phys. Rev. Lett. 105, 030401 (2010)] for an adiabatically steered two-level quantum system interacting with a Markovian environment. The original derivation employed the effective Hamiltonian in the adiabatic basis with the standard interaction picture approach but without the usual secular approximation. Our approach is based on utilizing a master equation for a nonsteered system in the first superadiabatic basis. It is potentially efficient in obtaining higher-order equations. Furthermore, we show how to select the phases of the adiabatic eigenstates to minimize the local adiabatic parameter and how this selection leads to states which are invariant under a local gauge change. We also discuss the effects of the adiabatic noncyclic geometric phase on the master equation.

ON ENTIRE SOLUTIONS OF TWO TYPES OF SYSTEMS OF COMPLEX DIFFERENTIAL-DIFFERENCE EQUATIONS

Institute of Scientific and Technical Information of China (English)

Lingyun GAO

2017-01-01

In this paper,we will mainly investigate entire solutions with finite order of two types of systems of differential-difference equations,and obtain some interesting results.It extends some results concerning complex differential (difference) equations to the systems of differential-difference equations.

Hybrid inverse problems for a system of Maxwell’s equations

International Nuclear Information System (INIS)

Bal, Guillaume; Zhou, Ting

2014-01-01

This paper concerns the quantitative step of the medical imaging modality thermo-acoustic tomography (TAT). We model the radiation propagation by a system of Maxwell’s equations. We show that the index of refraction of light and the absorption coefficient (conductivity) can be uniquely and stably reconstructed from a sufficiently large number of TAT measurements. Our method is based on verifying that the linearization of the inverse problem forms a redundant elliptic system of equations. We also observe that the reconstructions are qualitatively quite different from the setting where radiation is modeled by a scalar Helmholtz equation as in Bal G et al (2011 Inverse Problems 27 055007). (paper)

Lyapunov equation for infinite-dimensional discrete bilinear systems

International Nuclear Information System (INIS)

Costa, O.L.V.; Kubrusly, C.S.

1991-03-01

Mean-square stability for discrete systems requires that uniform convergence is preserved between input and state correlation sequences. Such a convergence preserving property holds for an infinite-dimensional bilinear system if and only if the associate Lyapunov equation has a unique strictly positive solution. (author)

A structure-preserving approach to normal form analysis of power systems; Una propuesta de preservacion de estructura al analisis de su forma normal en sistemas de potencia

Energy Technology Data Exchange (ETDEWEB)

Martinez Carrillo, Irma

2008-01-15

Power system dynamic behavior is inherently nonlinear and is driven by different processes at different time scales. The size and complexity of these mechanisms has stimulated the search for methods that reduce the original dimension but retain a certain degree of accuracy. In this dissertation, a novel nonlinear dynamical analysis method for the analysis of large amplitude oscillations that embraces ideas from normal form theory and singular perturbation techniques is proposed. This approach allows the full potential of the normal form method to be reached, and is suitably general for application to a wide variety of nonlinear systems. Drawing on the formal theory of dynamical systems, a structure-preserving model of the system is developed that preservers network and load characteristics. By exploiting the separation of fast and slow time scales of the model, an efficient approach based on singular perturbation techniques, is then derived for constructing a nonlinear power system representation that accurately preserves network structure. The method requires no reduction of the constraint equations and gives therefore, information about the effect of network and load characteristics on system behavior. Analytical expressions are then developed that provide approximate solutions to system performance near a singularity and techniques for interpreting these solutions in terms of modal functions are given. New insights into the nature of nonlinear oscillations are also offered and criteria for characterizing network effects on nonlinear system behavior are proposed. Theoretical insight into the behavior of dynamic coupling of differential-algebraic equations and the origin of nonlinearity is given, and implications for analyzing for design and placement of power system controllers in complex nonlinear systems are discussed. The extent of applicability of the proposed procedure is demonstrated by analyzing nonlinear behavior in two realistic test power systems

Multiparameter extrapolation and deflation methods for solving equation systems

Directory of Open Access Journals (Sweden)

A. J. Hughes Hallett

1984-01-01

Full Text Available Most models in economics and the applied sciences are solved by first order iterative techniques, usually those based on the Gauss-Seidel algorithm. This paper examines the convergence of multiparameter extrapolations (accelerations of first order iterations, as an improved approximation to the Newton method for solving arbitrary nonlinear equation systems. It generalises my earlier results on single parameter extrapolations. Richardson's generalised method and the deflation method for detecting successive solutions in nonlinear equation systems are also presented as multiparameter extrapolations of first order iterations. New convergence results are obtained for those methods.

On normal modes in classical Hamiltonian systems

NARCIS (Netherlands)

van Groesen, Embrecht W.C.

1983-01-01

Normal modes of Hamittonian systems that are even and of classical type are characterized as the critical points of a normalized kinetic energy functional on level sets of the potential energy functional. With the aid of this constrained variational formulation the existence of at least one family

Feynman-Kac equations for reaction and diffusion processes

Science.gov (United States)

Hou, Ru; Deng, Weihua

2018-04-01

This paper provides a theoretical framework for deriving the forward and backward Feynman-Kac equations for the distribution of functionals of the path of a particle undergoing both diffusion and reaction processes. Once given the diffusion type and reaction rate, a specific forward or backward Feynman-Kac equation can be obtained. The results in this paper include those for normal/anomalous diffusions and reactions with linear/nonlinear rates. Using the derived equations, we apply our findings to compute some physical (experimentally measurable) statistics, including the occupation time in half-space, the first passage time, and the occupation time in half-interval with an absorbing or reflecting boundary, for the physical system with anomalous diffusion and spontaneous evanescence.

A bifurcation analysis for the Lugiato-Lefever equation

Science.gov (United States)

Godey, Cyril

2017-05-01

The Lugiato-Lefever equation is a cubic nonlinear Schrödinger equation, including damping, detuning and driving, which arises as a model in nonlinear optics. We study the existence of stationary waves which are found as solutions of a four-dimensional reversible dynamical system in which the evolutionary variable is the space variable. Relying upon tools from bifurcation theory and normal forms theory, we discuss the codimension 1 bifurcations. We prove the existence of various types of steady solutions, including spatially localized, periodic, or quasi-periodic solutions. Contribution to the Topical Issue: "Theory and Applications of the Lugiato-Lefever Equation", edited by Yanne K. Chembo, Damia Gomila, Mustapha Tlidi, Curtis R. Menyuk.

Nonlocal Symmetries to Systems of Nonlinear Diffusion Equations

International Nuclear Information System (INIS)

Qu Changzheng; Kang Jing

2008-01-01

In this paper, we study potential symmetries to certain systems of nonlinear diffusion equations. Those systems have physical applications in soil science, mathematical biology, and invariant curve flows in R 3 . Lie point symmetries of the potential system, which cannot be projected to vector fields of the given dependent and independent variables, yield potential symmetries. The class of the system that admits potential symmetries is expanded.

3rd International Conference on Particle Systems and Partial Differential Equations

CERN Document Server

Soares, Ana

2016-01-01

The main focus of this book is on different topics in probability theory, partial differential equations and kinetic theory, presenting some of the latest developments in these fields. It addresses mathematical problems concerning applications in physics, engineering, chemistry and biology that were presented at the Third International Conference on Particle Systems and Partial Differential Equations, held at the University of Minho, Braga, Portugal in December 2014. The purpose of the conference was to bring together prominent researchers working in the fields of particle systems and partial differential equations, providing a venue for them to present their latest findings and discuss their areas of expertise. Further, it was intended to introduce a vast and varied public, including young researchers, to the subject of interacting particle systems, its underlying motivation, and its relation to partial differential equations. This book will appeal to probabilists, analysts and those mathematicians whose wor...

Selected equation of state in the acentric factor system

International Nuclear Information System (INIS)

Schreiber, D.R.; Pitzer, K.S.

1988-06-01

A new equation of state in the acentric factor system is developed on the basis of high-precision data. The region in critical temperature T/sub r/, critical density P/sub r/ space is identified where there is good agreement as well as the regions of significant departures. The equation fits very well in the critical region. 10 refs., 6 figs., 3 tabs

On the stability of some systems of exponential difference equations

Directory of Open Access Journals (Sweden)

N. Psarros

2018-01-01

Full Text Available In this paper we prove the stability of the zero equilibria of two systems of difference equations of exponential type, which are some extensions of an one-dimensional biological model. The stability of these systems is investigated in the special case when one of the eigenvalues is equal to -1 and the other eigenvalue has absolute value less than 1, using centre manifold theory. In addition, we study the existence and uniqueness of positive equilibria, the attractivity and the global asymptotic stability of these equilibria of some related systems of difference equations.

Efficient Instantiation of Parameterised Boolean Equation Systems to Parity Games

NARCIS (Netherlands)

Kant, Gijs; van de Pol, Jan Cornelis; Wijs, A.J.; Bošnački, D.; Edelkamp, S.

Parameterised Boolean Equation Systems (PBESs) are sequences of Boolean fixed point equations with data variables, used for, e.g., verification of modal μ-calculus formulae for process algebraic specifications with data. Solving a PBES is usually done by instantiation to a Parity Game and then

Modeling imperfectly repaired system data via grey differential equations with unequal-gapped times

International Nuclear Information System (INIS)

Guo Renkuan

2007-01-01

In this paper, we argue that grey differential equation models are useful in repairable system modeling. The arguments starts with the review on GM(1,1) model with equal- and unequal-spaced stopping time sequence. In terms of two-stage GM(1,1) filtering, system stopping time can be partitioned into system intrinsic function and repair effect. Furthermore, we propose an approach to use grey differential equation to specify a semi-statistical membership function for system intrinsic function times. Also, we engage an effort to use GM(1,N) model to model system stopping times and the associated operating covariates and propose an unequal-gapped GM(1,N) model for such analysis. Finally, we investigate the GM(1,1)-embed systematic grey equation system modeling of imperfectly repaired system operating data. Practical examples are given in step-by-step manner to illustrate the grey differential equation modeling of repairable system data

Poincare map for some polynomial systems of differential equations

International Nuclear Information System (INIS)

Varin, V P

2004-01-01

One approach to the classical problem of distinguishing between a centre and a focus for a system of differential equations with polynomial right-hand sides in the plane is discussed. For a broad class of such systems necessary and sufficient conditions for a centre are expressed in terms of equations in variations of higher order. By contrast with the existing methods of investigation, attention is concentrated on the explicit calculation of the asymptotic behaviour of the Poincare map rather than on finding sufficient centre conditions as such; this also enables one to study bifurcations of birth of arbitrarily strongly degenerate cycles.

Solving the Coupled System Improves Computational Efficiency of the Bidomain Equations

KAUST Repository

Southern, J.A.

2009-10-01

The bidomain equations are frequently used to model the propagation of cardiac action potentials across cardiac tissue. At the whole organ level, the size of the computational mesh required makes their solution a significant computational challenge. As the accuracy of the numerical solution cannot be compromised, efficiency of the solution technique is important to ensure that the results of the simulation can be obtained in a reasonable time while still encapsulating the complexities of the system. In an attempt to increase efficiency of the solver, the bidomain equations are often decoupled into one parabolic equation that is computationally very cheap to solve and an elliptic equation that is much more expensive to solve. In this study, the performance of this uncoupled solution method is compared with an alternative strategy in which the bidomain equations are solved as a coupled system. This seems counterintuitive as the alternative method requires the solution of a much larger linear system at each time step. However, in tests on two 3-D rabbit ventricle benchmarks, it is shown that the coupled method is up to 80% faster than the conventional uncoupled method-and that parallel performance is better for the larger coupled problem.

Solving the Coupled System Improves Computational Efficiency of the Bidomain Equations

KAUST Repository

Southern, J.A.; Plank, G.; Vigmond, E.J.; Whiteley, J.P.

2009-01-01

The bidomain equations are frequently used to model the propagation of cardiac action potentials across cardiac tissue. At the whole organ level, the size of the computational mesh required makes their solution a significant computational challenge. As the accuracy of the numerical solution cannot be compromised, efficiency of the solution technique is important to ensure that the results of the simulation can be obtained in a reasonable time while still encapsulating the complexities of the system. In an attempt to increase efficiency of the solver, the bidomain equations are often decoupled into one parabolic equation that is computationally very cheap to solve and an elliptic equation that is much more expensive to solve. In this study, the performance of this uncoupled solution method is compared with an alternative strategy in which the bidomain equations are solved as a coupled system. This seems counterintuitive as the alternative method requires the solution of a much larger linear system at each time step. However, in tests on two 3-D rabbit ventricle benchmarks, it is shown that the coupled method is up to 80% faster than the conventional uncoupled method-and that parallel performance is better for the larger coupled problem.

Fractional Diffusion Equations and Anomalous Diffusion

Science.gov (United States)

Evangelista, Luiz Roberto; Kaminski Lenzi, Ervin

2018-01-01

Preface; 1. Mathematical preliminaries; 2. A survey of the fractional calculus; 3. From normal to anomalous diffusion; 4. Fractional diffusion equations: elementary applications; 5. Fractional diffusion equations: surface effects; 6. Fractional nonlinear diffusion equation; 7. Anomalous diffusion: anisotropic case; 8. Fractional Schrödinger equations; 9. Anomalous diffusion and impedance spectroscopy; 10. The Poisson–Nernst–Planck anomalous (PNPA) models; References; Index.

Generation of static solutions of the self-consistent system of Einstein-Maxwell equations

International Nuclear Information System (INIS)

Anchikov, A.M.; Daishev, R.A.

1988-01-01

A theorem is proved, according to which to each solution of the Einstein equations with an arbitrary momentum-energy tensor in the right hand side there corresponds a static solution of the self-consistent system of Einstein-Maxwell equations. As a consequence of this theorem, a method is established of generating static solutions of the self-consistent system of Einstein-Maxwell equations with a charged grain as a source of vacuum solutions of the Einstein equations

Solution of underdetermined systems of equations with gridded a priori constraints.

Science.gov (United States)

Stiros, Stathis C; Saltogianni, Vasso

2014-01-01

The TOPINV, Topological Inversion algorithm (or TGS, Topological Grid Search) initially developed for the inversion of highly non-linear redundant systems of equations, can solve a wide range of underdetermined systems of non-linear equations. This approach is a generalization of a previous conclusion that this algorithm can be used for the solution of certain integer ambiguity problems in Geodesy. The overall approach is based on additional (a priori) information for the unknown variables. In the past, such information was used either to linearize equations around approximate solutions, or to expand systems of observation equations solved on the basis of generalized inverses. In the proposed algorithm, the a priori additional information is used in a third way, as topological constraints to the unknown n variables, leading to an R(n) grid containing an approximation of the real solution. The TOPINV algorithm does not focus on point-solutions, but exploits the structural and topological constraints in each system of underdetermined equations in order to identify an optimal closed space in the R(n) containing the real solution. The centre of gravity of the grid points defining this space corresponds to global, minimum-norm solutions. The rationale and validity of the overall approach are demonstrated on the basis of examples and case studies, including fault modelling, in comparison with SVD solutions and true (reference) values, in an accuracy-oriented approach.

Differential Equations Compatible with KZ Equations

International Nuclear Information System (INIS)

Felder, G.; Markov, Y.; Tarasov, V.; Varchenko, A.

2000-01-01

We define a system of 'dynamical' differential equations compatible with the KZ differential equations. The KZ differential equations are associated to a complex simple Lie algebra g. These are equations on a function of n complex variables z i taking values in the tensor product of n finite dimensional g-modules. The KZ equations depend on the 'dual' variable in the Cartan subalgebra of g. The dynamical differential equations are differential equations with respect to the dual variable. We prove that the standard hypergeometric solutions of the KZ equations also satisfy the dynamical equations. As an application we give a new determinant formula for the coordinates of a basis of hypergeometric solutions

Some overdetermined systems of complex partial differential equations

International Nuclear Information System (INIS)

Le Hung Son.

1990-01-01

In this paper we extend some properties of analytic functions on several complex variables to solutions of overdetermined systems of complex partial differential equations. It is proved that many global properties of analytic functions are true for solutions of the Vekua system in special cases. The relation between analytic functions and solutions of quasi-linear systems is discussed in the paper. (author). 8 refs

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

International Nuclear Information System (INIS)

Konopelchenko, B G; Ortenzi, G

2013-01-01

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

The core spline method for solution of quantum-mechanical systems of differential equations for bound states

International Nuclear Information System (INIS)

Aleksandrov, L.; Drenska, M.; Karadzhov, D.

1986-01-01

A generalization of the core spline method is given in the case of solution of the general bound state problem for a system of M linear differential equations with coefficients depending on the spectral parameter. The recursion scheme for construction of basic splines is described. The wave functions are expressed as linear combinations of basic splines, which are approximate partial solutions of the system. The spectral parameter (the eigenvalue) is determined from the condition for existence of a nontrivial solution of a (MxM) linear algebraic system at the last collocation point. The nontrivial solutions of this system determine (M - 1) coefficients of the linear spans, expressing the wave functions. The last unknown coefficient is determined from a boundary (or normalization) condition for the system. The computational aspects of the method are discussed, in particular, its concrete algorithmic realization used in the RODSOL program. The numerical solution of the Dirac system for the bound states of a hydrogen atom is given is an example

New component-based normalization method to correct PET system models

International Nuclear Information System (INIS)

Kinouchi, Shoko; Miyoshi, Yuji; Suga, Mikio; Yamaya, Taiga; Yoshida, Eiji; Nishikido, Fumihiko; Tashima, Hideaki

2011-01-01

Normalization correction is necessary to obtain high-quality reconstructed images in positron emission tomography (PET). There are two basic types of normalization methods: the direct method and component-based methods. The former method suffers from the problem that a huge count number in the blank scan data is required. Therefore, the latter methods have been proposed to obtain high statistical accuracy normalization coefficients with a small count number in the blank scan data. In iterative image reconstruction methods, on the other hand, the quality of the obtained reconstructed images depends on the system modeling accuracy. Therefore, the normalization weighing approach, in which normalization coefficients are directly applied to the system matrix instead of a sinogram, has been proposed. In this paper, we propose a new component-based normalization method to correct system model accuracy. In the proposed method, two components are defined and are calculated iteratively in such a way as to minimize errors of system modeling. To compare the proposed method and the direct method, we applied both methods to our small OpenPET prototype system. We achieved acceptable statistical accuracy of normalization coefficients while reducing the count number of the blank scan data to one-fortieth that required in the direct method. (author)

Modified Van der Waals equation and law of corresponding states

Science.gov (United States)

Zhong, Wei; Xiao, Changming; Zhu, Yongkai

2017-04-01

It is well known that the Van der Waals equation is a modification of the ideal gas law, yet it can be used to describe both gas and liquid, and some important messages can be obtained from this state equation. However, the Van der Waals equation is not a precise state equation, and it does not give a good description of the law of corresponding states. In this paper, we expand the Van der Waals equation into its Taylor's series form, and then modify the fourth order expansion by changing the constant Virial coefficients into their analogous ones. Via this way, a more precise result about the law of corresponding states has been obtained, and the law of corresponding states can then be expressed as: in terms of the reduced variables, all fluids should obey the same equation with the analogous Virial coefficients. In addition, the system of 3 He with quantum effects has also been taken into consideration with our modified Van der Waals equation, and it is found that, for a normal system without quantum effect, the modification on ideal gas law from the Van der Waals equation is more significant than the real case, however, for a system with quantum effect, this modification is less significant than the real case, thus a factor is introduced in this paper to weaken or strengthen the modification of the Van der Waals equation, respectively.

On extension of solutions of a simultaneous system of iterative functional equations

Directory of Open Access Journals (Sweden)

Janusz Matkowski

2009-01-01

Full Text Available Some sufficient conditions which allow to extend every local solution of a simultaneous system of equations in a single variable of the form \\[ \\varphi(x = h (x, \\varphi[f_1(x],\\ldots,\\varphi[f_m(x],\\] \\[\\varphi(x = H (x, \\varphi[F_1(x],\\ldots,\\varphi[F_m(x],\\] to a global one are presented. Extensions of solutions of functional equations, both in single and in several variables, play important role (cf. for instance [M. Kuczma, Functional equations in a single variable, Monografie Mat. 46, Polish Scientific Publishers, Warsaw, 1968, M. Kuczma, B. Choczewski, R. Ger, Iterative functional equations, Encyclopedia of Mathematics and Its Applications v. 32, Cambridge, 1990, J. Matkowski, Iteration groups, commuting functions and simultaneous systems of linear functional equations, Opuscula Math. 28 (2008 4, 531-541].

Iterative solution of large sparse systems of equations

CERN Document Server

Hackbusch, Wolfgang

2016-01-01

In the second edition of this classic monograph, complete with four new chapters and updated references, readers will now have access to content describing and analysing classical and modern methods with emphasis on the algebraic structure of linear iteration, which is usually ignored in other literature. The necessary amount of work increases dramatically with the size of systems, so one has to search for algorithms that most efficiently and accurately solve systems of, e.g., several million equations. The choice of algorithms depends on the special properties the matrices in practice have. An important class of large systems arises from the discretization of partial differential equations. In this case, the matrices are sparse (i.e., they contain mostly zeroes) and well-suited to iterative algorithms. The first edition of this book grew out of a series of lectures given by the author at the Christian-Albrecht University of Kiel to students of mathematics. The second edition includes quite novel approaches.

Experiments and Recommendations for Partitioning Systems of Equations

Directory of Open Access Journals (Sweden)

Mafteiu-Scai Liviu Octavian

2014-06-01

Full Text Available Partitioning the systems of equations is a very important process when solving it on a parallel computer. This paper presents some criteria which leads to more efficient parallelization, that must be taken into consideration. New criteria added to preconditioning process by reducing average bandwidth are pro- posed in this paper. These new criteria lead to a combination between preconditioning and partitioning of systems equations, so no need two distinct algorithms/processes. In our proposed methods - where the preconditioning is done by reducing the average bandwidth- two directions were followed in terms of partitioning: for a given preconditioned system determining the best partitioning (or one as close and the second consist in achieving an adequate preconditioning, depending on a given/desired partitioning. A mixed method it is also proposed. Experimental results, conclusions and recommendations, obtained after parallel implementation of conjugate gradient on IBM BlueGene /P supercomputer- based on a synchronous model of parallelization- are also presented in this paper.

Prolongation Loop Algebras for a Solitonic System of Equations

Directory of Open Access Journals (Sweden)

Maria A. Agrotis

2006-11-01

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

Lyapunov stability and its application to systems of ordinary differential equations

Science.gov (United States)

Kennedy, E. W.

1979-01-01

An outline and a brief introduction to some of the concepts and implications of Lyapunov stability theory are presented. Various aspects of the theory are illustrated by the inclusion of eight examples, including the Cartesian coordinate equations of the two-body problem, linear and nonlinear (Van der Pol's equation) oscillatory systems, and the linearized Kustaanheimo-Stiefel element equations for the unperturbed two-body problem.

Solution of generalized control system equations at steady state

International Nuclear Information System (INIS)

Vilim, R.B.

1987-01-01

Although a number of reactor systems codes feature generalized control system models, none of the models offer a steady-state solution finder. Indeed, if a transient is to begin from steady-state conditions, the user must provide estimates for the control system initial conditions and run a null transient until the plant converges to steady state. Several such transients may have to be run before values for control system demand signals are found that produce the desired plant steady state. The intent of this paper is (a) to present the control system equations assumed in the SASSYS reactor systems code and to identify the appropriate set of initial conditions, (b) to describe the generalized block diagram approach used to represent these equations, and (c) to describe a solution method and algorithm for computing these initial conditions from the block diagram. The algorithm has been installed in the SASSYS code for use with the code's generalized control system model. The solution finder greatly enhances the effectiveness of the code and the efficiency of the user in running it

Reduced equations of motion for quantum systems driven by diffusive Markov processes.

Science.gov (United States)

Sarovar, Mohan; Grace, Matthew D

2012-09-28

The expansion of a stochastic Liouville equation for the coupled evolution of a quantum system and an Ornstein-Uhlenbeck process into a hierarchy of coupled differential equations is a useful technique that simplifies the simulation of stochastically driven quantum systems. We expand the applicability of this technique by completely characterizing the class of diffusive Markov processes for which a useful hierarchy of equations can be derived. The expansion of this technique enables the examination of quantum systems driven by non-Gaussian stochastic processes with bounded range. We present an application of this extended technique by simulating Stark-tuned Förster resonance transfer in Rydberg atoms with nonperturbative position fluctuations.

The Cauchy problem for the Bogolyubov hierarchy of equations. The BCS model

International Nuclear Information System (INIS)

Vidybida, A.K.

1975-01-01

A chain of Bogolyubov's kinetic equations for an infinite quantum system of particles distributed in space with the mean density 1/V and interacting with the BCS model operator is considered as a single abstract equation in some countable normalized space bsup(v) of sequences of integral operators. In this case an unique solution of the Cauchy problem has been obtained at arbitrary initial conditions from bsup(v), stationary solutions of the equation have been derived, and the class of the initial conditions which approach to stationary ones is indicated

Multi criteria evaluation for universal soil loss equation based on geographic information system

Science.gov (United States)

Purwaamijaya, I. M.

2018-05-01

The purpose of this research were to produce(l) a conceptual, functional model designed and implementation for universal soil loss equation (usle), (2) standard operational procedure for multi criteria evaluation of universal soil loss equation (usle) using geographic information system, (3) overlay land cover, slope, soil and rain fall layers to gain universal soil loss equation (usle) using multi criteria evaluation, (4) thematic map of universal soil loss equation (usle) in watershed, (5) attribute table of universal soil loss equation (usle) in watershed. Descriptive and formal correlation methods are used for this research. Cikapundung Watershed, Bandung, West Java, Indonesia was study location. This research was conducted on January 2016 to May 2016. A spatial analysis is used to superimposed land cover, slope, soil and rain layers become universal soil loss equation (usle). Multi criteria evaluation for universal soil loss equation (usle) using geographic information system could be used for conservation program.

Normal modes and continuous spectra

International Nuclear Information System (INIS)

Balmforth, N.J.; Morrison, P.J.

1994-12-01

The authors consider stability problems arising in fluids, plasmas and stellar systems that contain singularities resulting from wave-mean flow or wave-particle resonances. Such resonances lead to singularities in the differential equations determining the normal modes at the so-called critical points or layers. The locations of the singularities are determined by the eigenvalue of the problem, and as a result, the spectrum of eigenvalues forms a continuum. They outline a method to construct the singular eigenfunctions comprising the continuum for a variety of problems

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.

Asymptotic Analysis of a System of Algebraic Equations Arising in Dislocation Theory

KAUST Repository

Hall, Cameron L.; Chapman, S. Jonathan; Ockendon, John R.

2010-01-01

The system of algebraic equations given by σn j=0, j≠=i sgn(xi-xj )|xi-xj|a = 1, i = 1, 2, ⋯ , n, x0 = 0, appears in dislocation theory in models of dislocation pile-ups. Specifically, the case a = 1 corresponds to the simple situation where n dislocations are piled up against a locked dislocation, while the case a = 3 corresponds to n dislocation dipoles piled up against a locked dipole. We present a general analysis of systems of this type for a > 0 and n large. In the asymptotic limit n→∞, it becomes possible to replace the system of discrete equations with a continuum equation for the particle density. For 0 < a < 2, this takes the form of a singular integral equation, while for a > 2 it is a first-order differential equation. The critical case a = 2 requires special treatment, but, up to corrections of logarithmic order, it also leads to a differential equation. The continuum approximation is valid only for i neither too small nor too close to n. The boundary layers at either end of the pile-up are also analyzed, which requires matching between discrete and continuum approximations to the main problem. © 2010 Society for Industrial and Applied Mathematics.

Investigation of the Dirac Equation by Using the Conformable Fractional Derivative

Science.gov (United States)

Mozaffari, F. S.; Hassanabadi, H.; Sobhani, H.; Chung, W. S.

2018-05-01

In this paper,the Dirac equation is constructed using the conformable fractional derivative so that in its limit for the fractional parameter, the normal version is recovered. Then, the Cornell potential is considered as the interaction of the system. In this case, the wave function and the energy eigenvalue equation are derived with the aim of the bi-confluent Heun functions. use of the conformable fractional derivative is proven to lead to a branching treatment for the energy of the system. Such a treatment is obvious for small values of the fractional parameter, and a united value as the fractional parameter approaches unity.

Convergence criteria for systems of nonlinear elliptic partial differential equations

International Nuclear Information System (INIS)

Sharma, R.K.

1986-01-01

This thesis deals with convergence criteria for a special system of nonlinear elliptic partial differential equations. A fixed-point algorithm is used, which iteratively solves one linearized elliptic partial differential equation at a time. Conditions are established that help foresee the convergence of the algorithm. Under reasonable hypotheses it is proved that the algorithm converges for such nonlinear elliptic systems. Extensive experimental results are reported and they show the algorithm converges in a wide variety of cases and the convergence is well correlated with the theoretical conditions introduced in this thesis

Nonlinear analysis of a reaction-diffusion system: Amplitude equations

Energy Technology Data Exchange (ETDEWEB)

Zemskov, E. P., E-mail: zemskov@ccas.ru [Russian Academy of Sciences, Dorodnicyn Computing Center (Russian Federation)

2012-10-15

A reaction-diffusion system with a nonlinear diffusion term is considered. Based on nonlinear analysis, the amplitude equations are obtained in the cases of the Hopf and Turing instabilities in the system. Turing pattern-forming regions in the parameter space are determined for supercritical and subcritical instabilities in a two-component reaction-diffusion system.

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

International Nuclear Information System (INIS)

Liu Chunliang; Xie Xi; Chen Yinbao

1991-01-01

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

The soliton solution of BBGKY quantum kinetic equations chain for different type particles system

International Nuclear Information System (INIS)

Rasulova, M.Yu.; Avazov, U.; Hassan, T.

2006-12-01

In the present paper on the basis of BBGKY chain of quantum kinetic equations the chain of equations for correlation matrices is derived, describing the evolution of a system of different types particles, which interact by pair potential. The series, which is the solution of this chain of equations for correlation matrices, is suggested. Using this series the solution of the last chain of equations is reduced to a solution of a set of homogeneous and nonhomogeneous von-Neumann's kinetic equations (analogue of Vlasov equations for quantum case). The first and second equations of this set of equations coincide with the first and second kinetic equations of the set, which is used in plasma physics. For an potential in the form of Dirac delta function, the solution of von-Neumann equation is defined through soliton solution of nonlinear Schrodinger equations. Based on von-Neumann equation one can define all terms of series, which is a solution of a chain of equations for correlation matrices. On the basis of these correlation matrices for a system of different types of particles we can define exact solution of BBGKY chain of quantum kinetic equations

Hamiltonian partial differential equations and applications

CERN Document Server

Nicholls, David; Sulem, Catherine

2015-01-01

This book is a unique selection of work by world-class experts exploring the latest developments in Hamiltonian partial differential equations and their applications. Topics covered within are representative of the field’s wide scope, including KAM and normal form theories, perturbation and variational methods, integrable systems, stability of nonlinear solutions as well as applications to cosmology, fluid mechanics and water waves. The volume contains both surveys and original research papers and gives a concise overview of the above topics, with results ranging from mathematical modeling to rigorous analysis and numerical simulation. It will be of particular interest to graduate students as well as researchers in mathematics and physics, who wish to learn more about the powerful and elegant analytical techniques for Hamiltonian partial differential equations.

Numerical Treatment of the Boltzmann Equation for Self-Propelled Particle Systems

Directory of Open Access Journals (Sweden)

Florian Thüroff

2014-11-01

Full Text Available Kinetic theories constitute one of the most promising tools to decipher the characteristic spatiotemporal dynamics in systems of actively propelled particles. In this context, the Boltzmann equation plays a pivotal role, since it provides a natural translation between a particle-level description of the system’s dynamics and the corresponding hydrodynamic fields. Yet, the intricate mathematical structure of the Boltzmann equation substantially limits the progress toward a full understanding of this equation by solely analytical means. Here, we propose a general framework to numerically solve the Boltzmann equation for self-propelled particle systems in two spatial dimensions and with arbitrary boundary conditions. We discuss potential applications of this numerical framework to active matter systems and use the algorithm to give a detailed analysis to a model system of self-propelled particles with polar interactions. In accordance with previous studies, we find that spatially homogeneous isotropic and broken-symmetry states populate two distinct regions in parameter space, which are separated by a narrow region of spatially inhomogeneous, density-segregated moving patterns. We find clear evidence that these three regions in parameter space are connected by first-order phase transitions and that the transition between the spatially homogeneous isotropic and polar ordered phases bears striking similarities to liquid-gas phase transitions in equilibrium systems. Within the density-segregated parameter regime, we find a novel stable limit-cycle solution of the Boltzmann equation, which consists of parallel lanes of polar clusters moving in opposite directions, so as to render the overall symmetry of the system’s ordered state nematic, despite purely polar interactions on the level of single particles.

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.

On Robust Stability of Systems of Differential-Algebraic Equations

Directory of Open Access Journals (Sweden)

A. Shcheglova

2016-06-01

The sufficient conditions of robust stability for index-one and index-two systems are obtained. We use the values of real and complex stability radii obtained for system of ordinary differential equations solved with respect to the derivatives. We consider the example illustrating the obtained results.

Classification of all solutions of the algebraic Riccati equations for infinite-dimensional systems

NARCIS (Netherlands)

Iftime, O; Curtain, R; Zwart, H

2003-01-01

We obtain a complete classification of all self-adjoint solution of the control algebraic Riccati equation for infinite-dimensional systems under the following assumptions: the system is output stabilizable, strongly detectable and the filter Riccati equation has an invertible self-adjoint

Lyapunov Functions and Solutions of the Lyapunov Matrix Equation for Marginally Stable Systems

DEFF Research Database (Denmark)

Kliem, Wolfhard; Pommer, Christian

2000-01-01

We consider linear systems of differential equations $I \\ddot{x}+B \\dot{x}+C{x}={0}$ where $I$ is the identity matrix and $B$ and $C$ are general complex $n$ x $n$ matrices. Our main interest is to determine conditions for complete marginalstability of these systems. To this end we find solutions...... of the Lyapunov matrix equation and characterize the set of matrices $(B, C)$ which guarantees marginal stability. The theory is applied to gyroscopic systems, to indefinite damped systems, and to circulatory systems, showing how to choose certain parameter matrices to get sufficient conditions for marginal...... stability.Comparison is made with some known results for equations with real system matrices.Moreover more general cases are investigated and several examples are given....

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

International Nuclear Information System (INIS)

Maccari, A.

1997-01-01

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

Modelling biochemical reaction systems by stochastic differential equations with reflection.

Science.gov (United States)

Niu, Yuanling; Burrage, Kevin; Chen, Luonan

2016-05-07

In this paper, we gave a new framework for modelling and simulating biochemical reaction systems by stochastic differential equations with reflection not in a heuristic way but in a mathematical way. The model is computationally efficient compared with the discrete-state Markov chain approach, and it ensures that both analytic and numerical solutions remain in a biologically plausible region. Specifically, our model mathematically ensures that species numbers lie in the domain D, which is a physical constraint for biochemical reactions, in contrast to the previous models. The domain D is actually obtained according to the structure of the corresponding chemical Langevin equations, i.e., the boundary is inherent in the biochemical reaction system. A variant of projection method was employed to solve the reflected stochastic differential equation model, and it includes three simple steps, i.e., Euler-Maruyama method was applied to the equations first, and then check whether or not the point lies within the domain D, and if not perform an orthogonal projection. It is found that the projection onto the closure D¯ is the solution to a convex quadratic programming problem. Thus, existing methods for the convex quadratic programming problem can be employed for the orthogonal projection map. Numerical tests on several important problems in biological systems confirmed the efficiency and accuracy of this approach. Copyright © 2016 Elsevier Ltd. All rights reserved.

Wavefield extrapolation in caustic-free normal ray coordinates

KAUST Repository

Ma, Xuxin

2012-11-04

Normal ray coordinates are conventionally constructed from ray tracing, which inherently requires smooth velocity profiles. To use rays as coordinates, the velocities have to be smoothed further to avoid caustics, which is detrimental to the mapping process. Solving the eikonal equation numerically for a line source at the surface provides a platform to map normal rays in complex unsmoothed velocity models and avoid caustics. We implement reverse-time migration (RTM) and downward continuation in the new ray coordinate system, which allows us to obtain efficient images and avoid some of the dip limitations of downward continuation.

Bethe-Salpeter equation for fermion-antifermion system in the ladder approximation

International Nuclear Information System (INIS)

Fukui, Ichio; Seto, Noriaki; Yoshida, Toshihiro.

1977-01-01

The Bethe-Salpeter (B-S) equation is important for studying hadron physics. Especially intensive investigation on the fermion-antifermion B-S equation is indispensable for the phenomenological studies of hardrons. However, many components of the B-S amplitude and the Wick-rotated integral kernel of non-Fredholm type have prevented from knowing details the solutions even in the ladder approximation. Some particular solutions are known in case of the vanishing four-momenta of bound states. The B-S equation for the bound state of fermion-anti-fermion system interacting through vector (axial-vector) particle exchange was studied in the ladder approximation with Feynman gauge. The reduced equations were obtained for suitably decomposed amplitude, and it is shown that, in the S-wave case, the coupled equations separate into two parts. In the nonrelativistic limit, large components of the amplitude satisfy the Wick-Cutkosky equation, and small components are expressed in terms of the large ones. Equations are derived for the equal-time amplitudes. (Kobatake, H.)

Nonlinear evolution equations and solving algebraic systems: the importance of computer algebra

International Nuclear Information System (INIS)

Gerdt, V.P.; Kostov, N.A.

1989-01-01

In the present paper we study the application of computer algebra to solve the nonlinear polynomial systems which arise in investigation of nonlinear evolution equations. We consider several systems which are obtained in classification of integrable nonlinear evolution equations with uniform rank. Other polynomial systems are related with the finding of algebraic curves for finite-gap elliptic potentials of Lame type and generalizations. All systems under consideration are solved using the method based on construction of the Groebner basis for corresponding polynomial ideals. The computations have been carried out using computer algebra systems. 20 refs

Chemical Equation Balancing.

Science.gov (United States)

Blakley, G. R.

1982-01-01

Reviews mathematical techniques for solving systems of homogeneous linear equations and demonstrates that the algebraic method of balancing chemical equations is a matter of solving a system of homogeneous linear equations. FORTRAN programs using this matrix method to chemical equation balancing are available from the author. (JN)

New form of the Euler-Bernoulli rod equation applied to robotic systems

Directory of Open Access Journals (Sweden)

Filipović Mirjana

2008-01-01

Full Text Available This paper presents a theoretical background and an example of extending the Euler-Bernoulli equation from several aspects. Euler-Bernoulli equation (based on the known laws of dynamics should be supplemented with all the forces that are participating in the formation of the bending moment of the considered mode. The stiffness matrix is a full matrix. Damping is an omnipresent elasticity characteristic of real systems, so that it is naturally included in the Euler-Bernoulli equation. It is shown that Daniel Bernoulli's particular integral is just one component of the total elastic deformation of the tip of any mode to which we have to add a component of the elastic deformation of a stationary regime in accordance with the complexity requirements of motion of an elastic robot system. The elastic line equation mode of link of a complex elastic robot system is defined based on the so-called 'Euler-Bernoulli Approach' (EBA. It is shown that the equation of equilibrium of all forces present at mode tip point ('Lumped-mass approach' (LMA follows directly from the elastic line equation for specified boundary conditions. This, in turn, proves the essential relationship between LMA and EBA approaches. In the defined mathematical model of a robotic system with multiple DOF (degree of freedom in the presence of the second mode, the phenomenon of elasticity of both links and joints are considered simultaneously with the presence of the environment dynamics - all based on the previously presented theoretical premises. Simulation results are presented. .

Discovering governing equations from data by sparse identification of nonlinear dynamical systems.

Science.gov (United States)

Brunton, Steven L; Proctor, Joshua L; Kutz, J Nathan

2016-04-12

Extracting governing equations from data is a central challenge in many diverse areas of science and engineering. Data are abundant whereas models often remain elusive, as in climate science, neuroscience, ecology, finance, and epidemiology, to name only a few examples. In this work, we combine sparsity-promoting techniques and machine learning with nonlinear dynamical systems to discover governing equations from noisy measurement data. The only assumption about the structure of the model is that there are only a few important terms that govern the dynamics, so that the equations are sparse in the space of possible functions; this assumption holds for many physical systems in an appropriate basis. In particular, we use sparse regression to determine the fewest terms in the dynamic governing equations required to accurately represent the data. This results in parsimonious models that balance accuracy with model complexity to avoid overfitting. We demonstrate the algorithm on a wide range of problems, from simple canonical systems, including linear and nonlinear oscillators and the chaotic Lorenz system, to the fluid vortex shedding behind an obstacle. The fluid example illustrates the ability of this method to discover the underlying dynamics of a system that took experts in the community nearly 30 years to resolve. We also show that this method generalizes to parameterized systems and systems that are time-varying or have external forcing.

EGSIEM combination service: combination of GRACE monthly K-band solutions on normal equation level

Science.gov (United States)

Meyer, Ulrich; Jean, Yoomin; Arnold, Daniel; Jäggi, Adrian

2017-04-01

The European Gravity Service for Improved Emergency Management (EGSIEM) project offers a scientific combination service, combining for the first time monthly GRACE gravity fields of different analysis centers (ACs) on normal equation (NEQ) level and thus taking all correlations between the gravity field coefficients and pre-eliminated orbit and instrument parameters correctly into account. Optimal weights for the individual NEQs are commonly derived by variance component estimation (VCE), as is the case for the products of the International VLBI Service (IVS) or the DTRF2008 reference frame realisation that are also derived by combination on NEQ-level. But variance factors are based on post-fit residuals and strongly depend on observation sampling and noise modeling, which both are very diverse in case of the individual EGSIEM ACs. These variance factors do not necessarily represent the true error levels of the estimated gravity field parameters that are still governed by analysis noise. We present a combination approach where weights are derived on solution level, thereby taking the analysis noise into account.

Wigner distribution functions for complex dynamical systems: the emergence of the Wigner-Boltzmann equation.

Science.gov (United States)

Sels, Dries; Brosens, Fons

2013-10-01

The equation of motion for the reduced Wigner function of a system coupled to an external quantum system is presented for the specific case when the external quantum system can be modeled as a set of harmonic oscillators. The result is derived from the Wigner function formulation of the Feynman-Vernon influence functional theory. It is shown how the true self-energy for the equation of motion is connected with the influence functional for the path integral. Explicit expressions are derived in terms of the bare Wigner propagator. Finally, we show under which approximations the resulting equation of motion reduces to the Wigner-Boltzmann equation.

Effective methods of solving of model equations of certain class of thermal systems

International Nuclear Information System (INIS)

Lach, J.

1985-01-01

A number of topics connected with solving of model equations of certain class of thermal systems by the method of successive approximations is touched. A system of partial differential equations of the first degree, appearing most frequently in practical applications of heat and mass transfer theory is reduced to an equivalent system of Volterra integral equations of the second kind. Among a few sample applications the thermal processes appearing in the fuel channel of nuclear reactor are solved. The theoretical analysis is illustrated by the results of numerical calculations given in tables and diagrams. 111 refs., 17 figs., 16 tabs. (author)

New Quasi-Newton Method for Solving Systems of Nonlinear Equations

Czech Academy of Sciences Publication Activity Database

Lukšan, Ladislav; Vlček, Jan

2017-01-01

Roč. 62, č. 2 (2017), s. 121-134 ISSN 0862-7940 R&D Projects: GA ČR GA13-06684S Institutional support: RVO:67985807 Keywords : nonlinear equations * systems of equations * trust-region methods * quasi-Newton methods * adjoint Broyden methods * numerical algorithms * numerical experiments Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.618, year: 2016 http://hdl.handle.net/10338.dmlcz/146699

On A System of Rational Difference Equation

OpenAIRE

Din Qamar

2014-01-01

In this paper, we study local asymptotic stability, global character and periodic nature of solutions of the system of rational difference equations given by xn+1= , yn=, n=0, 1,…, where the parameters a; b; c; d; e; f ∊ (0; ∞), and with initial conditions x0; y0 ∊ (0; ∞). Some numerical examples are given to illustrate our results.

Differential equations, dynamical systems, and an introduction to chaos

CERN Document Server

Smale, Stephen; Devaney, Robert L

2003-01-01

Thirty years in the making, this revised text by three of the world''s leading mathematicians covers the dynamical aspects of ordinary differential equations. it explores the relations between dynamical systems and certain fields outside pure mathematics, and has become the standard textbook for graduate courses in this area. The Second Edition now brings students to the brink of contemporary research, starting from a background that includes only calculus and elementary linear algebra.The authors are tops in the field of advanced mathematics, including Steve Smale who is a recipient of the Field''s Medal for his work in dynamical systems.* Developed by award-winning researchers and authors* Provides a rigorous yet accessible introduction to differential equations and dynamical systems* Includes bifurcation theory throughout* Contains numerous explorations for students to embark uponNEW IN THIS EDITION* New contemporary material and updated applications* Revisions throughout the text, including simplification...

PV System Performance Evaluation by Clustering Production Data to Normal and Non-Normal Operation

NARCIS (Netherlands)

Tsafarakis, O.; Sinapis, K.; van Sark, W.G.J.H.M.

2018-01-01

The most common method for assessment of a photovoltaic (PV) system performance is by comparing its energy production to reference data (irradiance or neighboring PV system). Ideally, at normal operation, the compared sets of data tend to show a linear relationship. Deviations from this linearity

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

Science.gov (United States)

Yoshimura, K.

1999-03-01

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

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.

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

International Nuclear Information System (INIS)

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

1988-01-01

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

The use of copulas to practical estimation of multivariate stochastic differential equation mixed effects models

International Nuclear Information System (INIS)

Rupšys, P.

2015-01-01

A system of stochastic differential equations (SDE) with mixed-effects parameters and multivariate normal copula density function were used to develop tree height model for Scots pine trees in Lithuania. A two-step maximum likelihood parameter estimation method is used and computational guidelines are given. After fitting the conditional probability density functions to outside bark diameter at breast height, and total tree height, a bivariate normal copula distribution model was constructed. Predictions from the mixed-effects parameters SDE tree height model calculated during this research were compared to the regression tree height equations. The results are implemented in the symbolic computational language MAPLE

The use of copulas to practical estimation of multivariate stochastic differential equation mixed effects models

Energy Technology Data Exchange (ETDEWEB)

Rupšys, P. [Aleksandras Stulginskis University, Studenų g. 11, Akademija, Kaunas district, LT – 53361 Lithuania (Lithuania)

2015-10-28

A system of stochastic differential equations (SDE) with mixed-effects parameters and multivariate normal copula density function were used to develop tree height model for Scots pine trees in Lithuania. A two-step maximum likelihood parameter estimation method is used and computational guidelines are given. After fitting the conditional probability density functions to outside bark diameter at breast height, and total tree height, a bivariate normal copula distribution model was constructed. Predictions from the mixed-effects parameters SDE tree height model calculated during this research were compared to the regression tree height equations. The results are implemented in the symbolic computational language MAPLE.

Effective Schroedinger equations on submanifolds

Energy Technology Data Exchange (ETDEWEB)

Wachsmuth, Jakob

2010-02-11

In this thesis the time dependent Schroedinger equation is considered on a Riemannian manifold A with a potential that localizes a certain class of states close to a fixed submanifold C, the constraint manifold. When the potential is scaled in the directions normal to C by a small parameter epsilon, the solutions concentrate in an epsilon-neighborhood of the submanifold. An effective Schroedinger equation on the submanifold C is derived and it is shown that its solutions, suitably lifted to A, approximate the solutions of the original equation on A up to errors of order {epsilon}{sup 3} vertical stroke t vertical stroke at time t. Furthermore, it is proved that, under reasonable conditions, the eigenvalues of the corresponding Hamiltonians below a certain energy coincide upto errors of order {epsilon}{sup 3}. These results holds in the situation where tangential and normal energies are of the same order, and where exchange between normal and tangential energies occurs. In earlier results tangential energies were assumed to be small compared to normal energies, and rather restrictive assumptions were needed, to ensure that the separation of energies is maintained during the time evolution. The most important consequence of this thesis is that now constraining potentials that change their shape along the submanifold can be treated, which is the typical situation in applications like molecular dynamics and quantum waveguides.

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.

3j Symbols: To Normalize or Not to Normalize?

Science.gov (United States)

van Veenendaal, Michel

2011-01-01

The systematic use of alternative normalization constants for 3j symbols can lead to a more natural expression of quantities, such as vector products and spherical tensor operators. The redefined coupling constants directly equate tensor products to the inner and outer products without any additional square roots. The approach is extended to…

Riccati and Ermakov Equations in Time-Dependent and Time-Independent Quantum Systems

Directory of Open Access Journals (Sweden)

Dieter Schuch

2008-05-01

Full Text Available The time-evolution of the maximum and the width of exact analytic wave packet (WP solutions of the time-dependent Schrödinger equation (SE represents the particle and wave aspects, respectively, of the quantum system. The dynamics of the maximum, located at the mean value of position, is governed by the Newtonian equation of the corresponding classical problem. The width, which is directly proportional to the position uncertainty, obeys a complex nonlinear Riccati equation which can be transformed into a real nonlinear Ermakov equation. The coupled pair of these equations yields a dynamical invariant which plays a key role in our investigation. It can be expressed in terms of a complex variable that linearizes the Riccati equation. This variable also provides the time-dependent parameters that characterize the Green's function, or Feynman kernel, of the corresponding problem. From there, also the relation between the classical and quantum dynamics of the systems can be obtained. Furthermore, the close connection between the Ermakov invariant and the Wigner function will be shown. Factorization of the dynamical invariant allows for comparison with creation/annihilation operators and supersymmetry where the partner potentials fulfil (real Riccati equations. This provides the link to a nonlinear formulation of time-independent quantum mechanics in terms of an Ermakov equation for the amplitude of the stationary state wave functions combined with a conservation law. Comparison with SUSY and the time-dependent problems concludes our analysis.

Resummed memory kernels in generalized system-bath master equations

International Nuclear Information System (INIS)

Mavros, Michael G.; Van Voorhis, Troy

2014-01-01

Generalized master equations provide a concise formalism for studying reduced population dynamics. Usually, these master equations require a perturbative expansion of the memory kernels governing the dynamics; in order to prevent divergences, these expansions must be resummed. Resummation techniques of perturbation series are ubiquitous in physics, but they have not been readily studied for the time-dependent memory kernels used in generalized master equations. In this paper, we present a comparison of different resummation techniques for such memory kernels up to fourth order. We study specifically the spin-boson Hamiltonian as a model system bath Hamiltonian, treating the diabatic coupling between the two states as a perturbation. A novel derivation of the fourth-order memory kernel for the spin-boson problem is presented; then, the second- and fourth-order kernels are evaluated numerically for a variety of spin-boson parameter regimes. We find that resumming the kernels through fourth order using a Padé approximant results in divergent populations in the strong electronic coupling regime due to a singularity introduced by the nature of the resummation, and thus recommend a non-divergent exponential resummation (the “Landau-Zener resummation” of previous work). The inclusion of fourth-order effects in a Landau-Zener-resummed kernel is shown to improve both the dephasing rate and the obedience of detailed balance over simpler prescriptions like the non-interacting blip approximation, showing a relatively quick convergence on the exact answer. The results suggest that including higher-order contributions to the memory kernel of a generalized master equation and performing an appropriate resummation can provide a numerically-exact solution to system-bath dynamics for a general spectral density, opening the way to a new class of methods for treating system-bath dynamics

Cucker-Smale model with normalized communication weights and time delay

KAUST Repository

Choi, Young-Pil

2017-03-06

We study a Cucker-Smale-type system with time delay in which agents interact with each other through normalized communication weights. We construct a Lyapunov functional for the system and provide sufficient conditions for asymptotic flocking, i.e., convergence to a common velocity vector. We also carry out a rigorous limit passage to the mean-field limit of the particle system as the number of particles tends to infinity. For the resulting Vlasov-type equation we prove the existence, stability and large-time behavior of measure-valued solutions. This is, to our best knowledge, the first such result for a Vlasov-type equation with time delay. We also present numerical simulations of the discrete system with few particles that provide further insights into the flocking and oscillatory behaviors of the particle velocities depending on the size of the time delay.

On A System of Rational Difference Equation

Directory of Open Access Journals (Sweden)

Din Qamar

2014-06-01

Full Text Available In this paper, we study local asymptotic stability, global character and periodic nature of solutions of the system of rational difference equations given by xn+1= , yn=, n=0, 1,…, where the parameters a; b; c; d; e; f ∊ (0; ∞, and with initial conditions x0; y0 ∊ (0; ∞. Some numerical examples are given to illustrate our results.

Geometric methods of global attraction in systems of delay differential equations

Science.gov (United States)

El-Morshedy, Hassan A.; Ruiz-Herrera, Alfonso

2017-11-01

In this paper we deduce criteria of global attraction in systems of delay differential equations. Our methodology is new and consists in "dominating" the nonlinear terms of the system by a scalar function and then studying some dynamical properties of that function. One of the crucial benefits of our approach is that we obtain delay-dependent results of global attraction that cover the best delay-independent conditions. We apply our results in a gene regulatory model and the classical Nicholson's blowfly equation with patch structure.

Darboux and binary Darboux transformations for discrete integrable systems I. Discrete potential KdV equation

International Nuclear Information System (INIS)

Shi, Ying; Zhang, Da-jun; Nimmo, Jonathan J C

2014-01-01

The Hirota–Miwa equation can be written in ‘nonlinear’ form in two ways: the discrete KP equation and, by using a compatible continuous variable, the discrete potential KP equation. For both systems, we consider the Darboux and binary Darboux transformations, expressed in terms of the continuous variable, and obtain exact solutions in Wronskian and Grammian form. We discuss reductions of both systems to the discrete KdV and discrete potential KdV equation, respectively, and exploit this connection to find the Darboux and binary Darboux transformations and exact solutions of these equations. (paper)

Differential equations

CERN Document Server

Barbu, Viorel

2016-01-01

This textbook is a comprehensive treatment of ordinary differential equations, concisely presenting basic and essential results in a rigorous manner. Including various examples from physics, mechanics, natural sciences, engineering and automatic theory, Differential Equations is a bridge between the abstract theory of differential equations and applied systems theory. Particular attention is given to the existence and uniqueness of the Cauchy problem, linear differential systems, stability theory and applications to first-order partial differential equations. Upper undergraduate students and researchers in applied mathematics and systems theory with a background in advanced calculus will find this book particularly useful. Supplementary topics are covered in an appendix enabling the book to be completely self-contained.

Solutions of system of P1 equations without use of auxiliary differential equations coupled; Solucoes do sistema de equacoes P1 sem o uso de equacoes diferenciais auxiliares acopladas

Energy Technology Data Exchange (ETDEWEB)

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

2000-07-01

The system of P1 equations is composed by two equations coupled itself one for the neutron flux and other for the current. Usually this system is solved by definitions of two integrals parameters, which are named slowing down densities of the flux and the current. Hence, the system P1 can be change from integral to only two differential equations. However, there are two new differentials equations that may be solved with the initial system. The present work analyzes this procedure and studies a method, which solve the P1 equations directly, without definitions of slowing down densities. (author)

Maxwell's equations, quantum physics and the quantum graviton

International Nuclear Information System (INIS)

Gersten, Alexander; Moalem, Amnon

2011-01-01

Quantum wave equations for massless particles and arbitrary spin are derived by factorizing the d'Alembertian operator. The procedure is extensively applied to the spin one photon equation which is related to Maxwell's equations via the proportionality of the photon wavefunction Ψ to the sum E + iB of the electric and magnetic fields. Thus Maxwell's equations can be considered as the first quantized one-photon equation. The photon wave equation is written in two forms, one with additional explicit subsidiary conditions and second with the subsidiary conditions implicitly included in the main equation. The second equation was obtained by factorizing the d'Alembertian with 4×4 matrix representation of 'relativistic quaternions'. Furthermore, scalar Lagrangian formalism, consistent with quantization requirements is developed using derived conserved current of probability and normalization condition for the wavefunction. Lessons learned from the derivation of the photon equation are used in the derivation of the spin two quantum equation, which we call the quantum graviton. Quantum wave equation with implicit subsidiary conditions, which factorizes the d'Alembertian with 8×8 matrix representation of relativistic quaternions, is derived. Scalar Lagrangian is formulated and conserved probability current and wavefunction normalization are found, both consistent with the definitions of quantum operators and their expectation values. We are showing that the derived equations are the first quantized equations of the photon and the graviton.

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

Amplitude equations for a sub-diffusive reaction-diffusion system

International Nuclear Information System (INIS)

Nec, Y; Nepomnyashchy, A A

2008-01-01

A sub-diffusive reaction-diffusion system with a positive definite memory operator and a nonlinear reaction term is analysed. Amplitude equations (Ginzburg-Landau type) are derived for short wave (Turing) and long wave (Hopf) bifurcation points

Normal form of linear systems depending on parameters

International Nuclear Information System (INIS)

Nguyen Huynh Phan.

1995-12-01

In this paper we resolve completely the problem to find normal forms of linear systems depending on parameters for the feedback action that we have studied for the special case of controllable linear systems. (author). 24 refs

Poincaré-MacMillan Equations of Motion for a Nonlinear Nonholonomic Dynamical System

Science.gov (United States)

Amjad, Hussain; Syed Tauseef, Mohyud-Din; Ahmet, Yildirim

2012-03-01

MacMillan's equations are extended to Poincaré's formalism, and MacMillan's equations for nonlinear nonholonomic systems are obtained in terms of Poincaré parameters. The equivalence of the results obtained here with other forms of equations of motion is demonstrated. An illustrative example of the theory is provided as well.

Equation-free modeling unravels the behavior of complex ecological systems

Science.gov (United States)

DeAngelis, Donald L.; Yurek, Simeon

2015-01-01

Ye et al. (1) address a critical problem confronting the management of natural ecosystems: How can we make forecasts of possible future changes in populations to help guide management actions? This problem is especially acute for marine and anadromous fisheries, where the large interannual fluctuations of populations, arising from complex nonlinear interactions among species and with varying environmental factors, have defied prediction over even short time scales. The empirical dynamic modeling (EDM) described in Ye et al.’s report, the latest in a series of papers by Sugihara and his colleagues, offers a promising quantitative approach to building models using time series to successfully project dynamics into the future. With the term “equation-free” in the article title, Ye et al. (1) are suggesting broader implications of their approach, considering the centrality of equations in modern science. From the 1700s on, nature has been increasingly described by mathematical equations, with differential or difference equations forming the basic framework for describing dynamics. The use of mathematical equations for ecological systems came much later, pioneered by Lotka and Volterra, who showed that population cycles might be described in terms of simple coupled nonlinear differential equations. It took decades for Lotka–Volterra-type models to become established, but the development of appropriate differential equations is now routine in modeling ecological dynamics. There is no question that the injection of mathematical equations, by forcing “clarity and precision into conjecture” (2), has led to increased understanding of population and community dynamics. As in science in general, in ecology equations are a key method of communication and of framing hypotheses. These equations serve as compact representations of an enormous amount of empirical data and can be analyzed by the powerful methods of mathematics.

Boundary-value problems with free boundaries for elliptic systems of equations

CERN Document Server

Monakhov, V N

1983-01-01

This book is concerned with certain classes of nonlinear problems for elliptic systems of partial differential equations: boundary-value problems with free boundaries. The first part has to do with the general theory of boundary-value problems for analytic functions and its applications to hydrodynamics. The second presents the theory of quasiconformal mappings, along with the theory of boundary-value problems for elliptic systems of equations and applications of it to problems in the mechanics of continuous media with free boundaries: problems in subsonic gas dynamics, filtration theory, and problems in elastico-plasticity.

Advection-diffusion model for normal grain growth and the stagnation of normal grain growth in thin films

International Nuclear Information System (INIS)

Lou, C.

2002-01-01

An advection-diffusion model has been set up to describe normal grain growth. In this model grains are divided into different groups according to their topological classes (number of sides of a grain). Topological transformations are modelled by advective and diffusive flows governed by advective and diffusive coefficients respectively, which are assumed to be proportional to topological classes. The ordinary differential equations governing self-similar time-independent grain size distribution can be derived analytically from continuity equations. It is proved that the time-independent distributions obtained by solving the ordinary differential equations have the same form as the time-dependent distributions obtained by solving the continuity equations. The advection-diffusion model is extended to describe the stagnation of normal grain growth in thin films. Grain boundary grooving prevents grain boundaries from moving, and the correlation between neighbouring grains accelerates the stagnation of normal grain growth. After introducing grain boundary grooving and the correlation between neighbouring grains into the model, the grain size distribution is close to a lognormal distribution, which is usually found in experiments. A vertex computer simulation of normal grain growth has also been carried out to make a cross comparison with the advection-diffusion model. The result from the simulation did not verify the assumption that the advective and diffusive coefficients are proportional to topological classes. Instead, we have observed that topological transformations usually occur on certain topological classes. This suggests that the advection-diffusion model can be improved by making a more realistic assumption on topological transformations. (author)

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

International Nuclear Information System (INIS)

Yu Fajun; Zhang Hongqing

2008-01-01

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

On the use of the autonomous Birkhoff equations in Lie series perturbation theory

Science.gov (United States)

Boronenko, T. S.

2017-02-01

In this article, we present the Lie transformation algorithm for autonomous Birkhoff systems. Here, we are referring to Hamiltonian systems that obey a symplectic structure of the general form. The Birkhoff equations are derived from the linear first-order Pfaff-Birkhoff variational principle, which is more general than the Hamilton principle. The use of 1-form in formulating the equations of motion in dynamics makes the Birkhoff method more universal and flexible. Birkhoff's equations have a tensorial character, so their form is independent of the coordinate system used. Two examples of normalization in the restricted three-body problem are given to illustrate the application of the algorithm in perturbation theory. The efficiency of this algorithm for problems of asymptotic integration in dynamics is discussed for the case where there is a need to use non-canonical variables in phase space.

Volume-preserving normal forms of Hopf-zero singularity

International Nuclear Information System (INIS)

Gazor, Majid; Mokhtari, Fahimeh

2013-01-01

A practical method is described for computing the unique generator of the algebra of first integrals associated with a large class of Hopf-zero singularity. The set of all volume-preserving classical normal forms of this singularity is introduced via a Lie algebra description. This is a maximal vector space of classical normal forms with first integral; this is whence our approach works. Systems with a nonzero condition on their quadratic parts are considered. The algebra of all first integrals for any such system has a unique (modulo scalar multiplication) generator. The infinite level volume-preserving parametric normal forms of any nondegenerate perturbation within the Lie algebra of any such system is computed, where it can have rich dynamics. The associated unique generator of the algebra of first integrals are derived. The symmetry group of the infinite level normal forms are also discussed. Some necessary formulas are derived and applied to appropriately modified Rössler and generalized Kuramoto–Sivashinsky equations to demonstrate the applicability of our theoretical results. An approach (introduced by Iooss and Lombardi) is applied to find an optimal truncation for the first level normal forms of these examples with exponentially small remainders. The numerically suggested radius of convergence (for the first integral) associated with a hypernormalization step is discussed for the truncated first level normal forms of the examples. This is achieved by an efficient implementation of the results using Maple. (paper)

Volume-preserving normal forms of Hopf-zero singularity

Science.gov (United States)

Gazor, Majid; Mokhtari, Fahimeh

2013-10-01

A practical method is described for computing the unique generator of the algebra of first integrals associated with a large class of Hopf-zero singularity. The set of all volume-preserving classical normal forms of this singularity is introduced via a Lie algebra description. This is a maximal vector space of classical normal forms with first integral; this is whence our approach works. Systems with a nonzero condition on their quadratic parts are considered. The algebra of all first integrals for any such system has a unique (modulo scalar multiplication) generator. The infinite level volume-preserving parametric normal forms of any nondegenerate perturbation within the Lie algebra of any such system is computed, where it can have rich dynamics. The associated unique generator of the algebra of first integrals are derived. The symmetry group of the infinite level normal forms are also discussed. Some necessary formulas are derived and applied to appropriately modified Rössler and generalized Kuramoto-Sivashinsky equations to demonstrate the applicability of our theoretical results. An approach (introduced by Iooss and Lombardi) is applied to find an optimal truncation for the first level normal forms of these examples with exponentially small remainders. The numerically suggested radius of convergence (for the first integral) associated with a hypernormalization step is discussed for the truncated first level normal forms of the examples. This is achieved by an efficient implementation of the results using Maple.

Existence of a coupled system of fractional differential equations

International Nuclear Information System (INIS)

Ibrahim, Rabha W.; Siri, Zailan

2015-01-01

We manage the existence and uniqueness of a fractional coupled system containing Schrödinger equations. Such a system appears in quantum mechanics. We confirm that the fractional system under consideration admits a global solution in appropriate functional spaces. The solution is shown to be unique. The method is based on analytic technique of the fixed point theory. The fractional differential operator is considered from the virtue of the Riemann-Liouville differential operator

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

Existence of a coupled system of fractional differential equations

Energy Technology Data Exchange (ETDEWEB)

Ibrahim, Rabha W. [Multimedia unit, Department of Computer System and Technology Faculty of Computer Science & IT, University of Malaya, 50603 Kuala Lumpur (Malaysia); Siri, Zailan [Institute of Mathematical Sciences, University of Malaya, 50603 Kuala Lumpur (Malaysia)

2015-10-22

We manage the existence and uniqueness of a fractional coupled system containing Schrödinger equations. Such a system appears in quantum mechanics. We confirm that the fractional system under consideration admits a global solution in appropriate functional spaces. The solution is shown to be unique. The method is based on analytic technique of the fixed point theory. The fractional differential operator is considered from the virtue of the Riemann-Liouville differential operator.

An equations of motion approach for open shell systems

International Nuclear Information System (INIS)

Yeager, D.L.; McKoy, V.

1975-01-01

A straightforward scheme is developed for extending the equations of motion formalism to systems with simple open shell ground states. Equations for open shell random phase approximation (RPA) are given for the cases of one electron outside of a closed shell in a nondegenerate molecular orbital and for the triplet ground state with two electrons outside of a closed shell in degenerate molecular orbitals. Applications to other open shells and extension of the open shell EOM to higher orders are both straightforward. Results for the open shell RPA for lithium atom and oxygen molecule are given

Systems of Differential Equations with Skew-Symmetric, Orthogonal Matrices

Science.gov (United States)

Glaister, P.

2008-01-01

The solution of a system of linear, inhomogeneous differential equations is discussed. The particular class considered is where the coefficient matrix is skew-symmetric and orthogonal, and where the forcing terms are sinusoidal. More general matrices are also considered.

Convex solutions of systems arising from Monge-Ampere equations

Directory of Open Access Journals (Sweden)

Haiyan Wang

2009-10-01

Full Text Available We establish two criteria for the existence of convex solutions to a boundary value problem for weakly coupled systems arising from the Monge-Ampère equations. We shall use fixed point theorems in a cone.

Fully Digital Chaotic Differential Equation-based Systems And Methods

KAUST Repository

Radwan, Ahmed Gomaa Ahmed; Zidan, Mohammed A.; Salama, Khaled N.

2012-01-01

Various embodiments are provided for fully digital chaotic differential equation-based systems and methods. In one embodiment, among others, a digital circuit includes digital state registers and one or more digital logic modules configured to obtain a first value from two or more of the digital state registers; determine a second value based upon the obtained first values and a chaotic differential equation; and provide the second value to set a state of one of the plurality of digital state registers. In another embodiment, a digital circuit includes digital state registers, digital logic modules configured to obtain outputs from a subset of the digital shift registers and to provide the input based upon a chaotic differential equation for setting a state of at least one of the subset of digital shift registers, and a digital clock configured to provide a clock signal for operating the digital shift registers.

Fully Digital Chaotic Differential Equation-based Systems And Methods

KAUST Repository

Radwan, Ahmed Gomaa Ahmed

2012-09-06

Various embodiments are provided for fully digital chaotic differential equation-based systems and methods. In one embodiment, among others, a digital circuit includes digital state registers and one or more digital logic modules configured to obtain a first value from two or more of the digital state registers; determine a second value based upon the obtained first values and a chaotic differential equation; and provide the second value to set a state of one of the plurality of digital state registers. In another embodiment, a digital circuit includes digital state registers, digital logic modules configured to obtain outputs from a subset of the digital shift registers and to provide the input based upon a chaotic differential equation for setting a state of at least one of the subset of digital shift registers, and a digital clock configured to provide a clock signal for operating the digital shift registers.

Numerical solution of stiff systems of ordinary differential equations with applications to electronic circuits

Science.gov (United States)

Rosenbaum, J. S.

1971-01-01

Systems of ordinary differential equations in which the magnitudes of the eigenvalues (or time constants) vary greatly are commonly called stiff. Such systems of equations arise in nuclear reactor kinetics, the flow of chemically reacting gas, dynamics, control theory, circuit analysis and other fields. The research reported develops an A-stable numerical integration technique for solving stiff systems of ordinary differential equations. The method, which is called the generalized trapezoidal rule, is a modification of the trapezoidal rule. However, the method is computationally more efficient than the trapezoidal rule when the solution of the almost-discontinuous segments is being calculated.

Tensor-GMRES method for large sparse systems of nonlinear equations

Science.gov (United States)

Feng, Dan; Pulliam, Thomas H.

1994-01-01

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

Critical Dynamics : The Expansion of the Master Equation Including a Critical Point

NARCIS (Netherlands)

Dekker, H.

1980-01-01

In this thesis it is shown how to solve the master equation for a Markov process including a critical point by means of successive approximations in terms of a small parameter. A critical point occurs if, by adjusting an externally controlled quantity, the system shows a transition from normal

How the 2SLS/IV estimator can handle equality constraints in structural equation models: a system-of-equations approach.

Science.gov (United States)

Nestler, Steffen

2014-05-01

Parameters in structural equation models are typically estimated using the maximum likelihood (ML) approach. Bollen (1996) proposed an alternative non-iterative, equation-by-equation estimator that uses instrumental variables. Although this two-stage least squares/instrumental variables (2SLS/IV) estimator has good statistical properties, one problem with its application is that parameter equality constraints cannot be imposed. This paper presents a mathematical solution to this problem that is based on an extension of the 2SLS/IV approach to a system of equations. We present an example in which our approach was used to examine strong longitudinal measurement invariance. We also investigated the new approach in a simulation study that compared it with ML in the examination of the equality of two latent regression coefficients and strong measurement invariance. Overall, the results show that the suggested approach is a useful extension of the original 2SLS/IV estimator and allows for the effective handling of equality constraints in structural equation models. © 2013 The British Psychological Society.

Extended rate equations

International Nuclear Information System (INIS)

Shore, B.W.

1981-01-01

The equations of motion are discussed which describe time dependent population flows in an N-level system, reviewing the relationship between incoherent (rate) equations, coherent (Schrodinger) equations, and more general partially coherent (Bloch) equations. Approximations are discussed which replace the elaborate Bloch equations by simpler rate equations whose coefficients incorporate long-time consequences of coherence

Normal mode approach to modelling of feedback stabilization of the resistive wall mode

International Nuclear Information System (INIS)

Chu, M.S.; Chance, M.S.; Okabayashi, M.; Glasser, A.H.

2003-01-01

Feedback stabilization of the resistive wall mode (RWM) of a plasma in a general feedback configuration is formulated in terms of the normal modes of the plasma-resistive wall system. The growth/damping rates and the eigenfunctions of the normal modes are determined by an extended energy principle for the plasma during its open (feedback) loop operation. A set of equations are derived for the time evolution of these normal modes with currents in the feedback coils. The dynamics of the feedback system is completed by the prescription of the feedback logic. The feasibility of the feedback is evaluated by using the Nyquist diagram method or by solving the characteristic equations. The elements of the characteristic equations are formed from the growth and damping rates of the normal modes, the sensor matrix of the perturbation fluxes detected by the sensor loops, the excitation matrix of the energy input to the normal modes by the external feedback coils, and the feedback logic. (The RWM is also predicted to be excited by an external error field to a large amplitude when it is close to marginal stability.) This formulation has been implemented numerically and applied to the DIII-D tokamak. It is found that feedback with poloidal sensors is much more effective than feedback with radial sensors. Using radial sensors, increasing the number of feedback coils from a central band on the outboard side to include an upper and a lower band can substantially increase the effectiveness of the feedback system. The strength of the RWM that can be stabilized is increased from γτ w = 1 to 30 (γ is the growth rate of the RWM in the absence of feedback and τ w is the resistive wall time constant) Using poloidal sensors, just one central band of feedback coils is sufficient for the stabilization of the RWM with γτ w = 30. (author)

Gaussian approximations for stochastic systems with delay: Chemical Langevin equation and application to a Brusselator system

International Nuclear Information System (INIS)

Brett, Tobias; Galla, Tobias

2014-01-01

We present a heuristic derivation of Gaussian approximations for stochastic chemical reaction systems with distributed delay. In particular, we derive the corresponding chemical Langevin equation. Due to the non-Markovian character of the underlying dynamics, these equations are integro-differential equations, and the noise in the Gaussian approximation is coloured. Following on from the chemical Langevin equation, a further reduction leads to the linear-noise approximation. We apply the formalism to a delay variant of the celebrated Brusselator model, and show how it can be used to characterise noise-driven quasi-cycles, as well as noise-triggered spiking. We find surprisingly intricate dependence of the typical frequency of quasi-cycles on the delay period

Gaussian approximations for stochastic systems with delay: chemical Langevin equation and application to a Brusselator system.

Science.gov (United States)

Brett, Tobias; Galla, Tobias

2014-03-28

We present a heuristic derivation of Gaussian approximations for stochastic chemical reaction systems with distributed delay. In particular, we derive the corresponding chemical Langevin equation. Due to the non-Markovian character of the underlying dynamics, these equations are integro-differential equations, and the noise in the Gaussian approximation is coloured. Following on from the chemical Langevin equation, a further reduction leads to the linear-noise approximation. We apply the formalism to a delay variant of the celebrated Brusselator model, and show how it can be used to characterise noise-driven quasi-cycles, as well as noise-triggered spiking. We find surprisingly intricate dependence of the typical frequency of quasi-cycles on the delay period.

Systemic sclerosis with normal or nonspecific nailfold capillaroscopy.

Science.gov (United States)

Fichel, Fanny; Baudot, Nathalie; Gaitz, Jean-Pierre; Trad, Salim; Barbe, Coralie; Francès, Camille; Senet, Patricia

2014-01-01

In systemic sclerosis (SSc), a specific nailfold videocapillaroscopy (NVC) pattern is observed in 90% of cases and seems to be associated with severity and progression of the disease. To describe the characteristics of SSc patients with normal or nonspecific (normal/nonspecific) NVC. In a retrospective cohort study, clinical features and visceral involvements of 25 SSc cases with normal/nonspecific NVC were compared to 63 SSc controls with the SSc-specific NVC pattern. Normal/nonspecific NVC versus SSc-specific NVC pattern was significantly associated with absence of skin sclerosis (32 vs. 6.3%, p = 0.004), absence of telangiectasia (47.8 vs. 17.3%, p = 0.006) and absence of sclerodactyly (60 vs. 25.4%, p = 0.002), and less frequent severe pulmonary involvement (26.3 vs. 58.2%, p = 0.017). Normal/nonspecific NVC in SSc patients appears to be associated with less severe skin involvement and less frequent severe pulmonary involvement. © 2014 S. Karger AG, Basel.

Coupled replicator equations for the dynamics of learning in multiagent systems

Science.gov (United States)

Sato, Yuzuru; Crutchfield, James P.

2003-01-01

Starting with a group of reinforcement-learning agents we derive coupled replicator equations that describe the dynamics of collective learning in multiagent systems. We show that, although agents model their environment in a self-interested way without sharing knowledge, a game dynamics emerges naturally through environment-mediated interactions. An application to rock-scissors-paper game interactions shows that the collective learning dynamics exhibits a diversity of competitive and cooperative behaviors. These include quasiperiodicity, stable limit cycles, intermittency, and deterministic chaos—behaviors that should be expected in heterogeneous multiagent systems described by the general replicator equations we derive.

Multilevel solvers of first-order system least-squares for Stokes equations

Energy Technology Data Exchange (ETDEWEB)

Lai, Chen-Yao G. [National Chung Cheng Univ., Chia-Yi (Taiwan, Province of China)

1996-12-31

Recently, The use of first-order system least squares principle for the approximate solution of Stokes problems has been extensively studied by Cai, Manteuffel, and McCormick. In this paper, we study multilevel solvers of first-order system least-squares method for the generalized Stokes equations based on the velocity-vorticity-pressure formulation in three dimensions. The least-squares functionals is defined to be the sum of the L{sup 2}-norms of the residuals, which is weighted appropriately by the Reynolds number. We develop convergence analysis for additive and multiplicative multilevel methods applied to the resulting discrete equations.

Local first integrals for systems of differential equations

International Nuclear Information System (INIS)

Zhang Xiang

2003-01-01

The main purpose of this paper is to provide some sufficient conditions for a system of differential equations to have local first integrals in a certain neighbourhood of a singularity. Our results generalize those given in Kwek et al (2003 Z. Angew. Math. Phys. 54 26) and Li et al (2003 Z. Angew. Math. Phys. 54 235)

Almost periodic solutions to systems of parabolic equations

Directory of Open Access Journals (Sweden)

Janpou Nee

1994-01-01

Full Text Available In this paper we show that the second-order differential solution is 2-almost periodic, provided it is 2-bounded, and the growth of the components of a non-linear function of a system of parabolic equation is bounded by any pair of con-secutive eigenvalues of the associated Dirichlet boundary value problems.

Asymmetric systems described by a pair of local covariant wave equations

Energy Technology Data Exchange (ETDEWEB)

Mallik, S [Bern Univ. (Switzerland). Inst. fuer Theoretische Physik

1979-07-16

A class of asymmetric solutions of the integrability conditions for systems obeying the Leutwyler-Stern pair of covariant wave equations is obtained. The class of unequal-mass systems described by these solutions does not embed the particle-antiparticle system behaving as a relativistic harmonic oscillator.

The coordinate system transformation of a serial kinematic structures and use in the derivation of systems motion equations

Directory of Open Access Journals (Sweden)

Zátopek Jiří

2016-01-01

Full Text Available This text discusses the use of transformation matrices to determine the motion equations of the complex mechanical structure. Use of the transformation matrix does not apply only to motion equations but has the general use in relative positions determine of objects in the 3D space. Analysed model is divided into seven physical objects, the transformation matrix and the corresponding inertia/pseudo-inertia matrix is included in each of them. This matrices are strictly necessary to the system dynamic description using the matrix form of Lagrange Equations of the Second Type. Another possibility to use the transformation matrix is shown in the camera system measurement. Model was designed in 3D CAD system SolidWorks, MATLAB was used for the mathematical calculations.

Semilinear hyperbolic systems and equations with singular initial data

International Nuclear Information System (INIS)

Gramchev, T.

1991-07-01

We study the weak limits of solutions u ε (t, ·) for ε→0 to semilinear strictly hyperbolic systems and wave equations with initial data u ε (0, ·) approximating a distribution κ, 0 ε (t, ·) for ε→0 exists. 13 refs

FMTLxLyLz DIMENSIONAL EQUAT DIMENSIONAL EQUATION ...

African Journals Online (AJOL)

eobe

plant made of 12mm thick steel plate was used in de steel plate ... water treatment plant. ... ameters affecting filtration processes were used to derive an equation usin ..... system. However, in deriving the equation onl terms are incorporated.

A stockability equation for forest land in Siskiyou County, California.

Science.gov (United States)

Neil. McKay

1985-01-01

An equation is presented that estimates the relative stocking capacity of forest land in Siskiyou County, California, from the amount of precipitation and the presence of significant indicator plants. The equation is a toot for identifying sites incapable of supporting normal stocking. Estimated relative stocking capacity may be used to discount normal yields to levels...

On the relation between the fundamental equation of thermodynamics and the energy balance equation in the context of closed and open systems

NARCIS (Netherlands)

Knuiman, J.T.; Barneveld, P.A.

2012-01-01

In this paper, we elaborate on the connection between the fundamental equation of thermodynamics, which is essentially the combination of the First and Second Laws of thermodynamics, and the energy balance equation in the context of closed and open systems. We point out some misinterpretations in

Comprehensive solutions to the Bloch equations and dynamical models for open two-level systems

Science.gov (United States)

Skinner, Thomas E.

2018-01-01

The Bloch equation and its variants constitute the fundamental dynamical model for arbitrary two-level systems. Many important processes, including those in more complicated systems, can be modeled and understood through the two-level approximation. It is therefore of widespread relevance, especially as it relates to understanding dissipative processes in current cutting-edge applications of quantum mechanics. Although the Bloch equation has been the subject of considerable analysis in the 70 years since its inception, there is still, perhaps surprisingly, significant work that can be done. This paper extends the scope of previous analyses. It provides a framework for more fully understanding the dynamics of dissipative two-level systems. A solution is derived that is compact, tractable, and completely general, in contrast to previous results. Any solution of the Bloch equation depends on three roots of a cubic polynomial that are crucial to the time dependence of the system. The roots are typically only sketched out qualitatively, with no indication of their dependence on the physical parameters of the problem. Degenerate roots, which modify the solutions, have been ignored altogether. Here the roots are obtained explicitly in terms of a single real-valued root that is expressed as a simple function of the system parameters. For the conventional Bloch equation, a simple graphical representation of this root is presented that makes evident the explicit time dependence of the system for each point in the parameter space. Several intuitive, visual models of system dynamics are developed. A Euclidean coordinate system is identified in which any generalized Bloch equation is separable, i.e., the sum of commuting rotation and relaxation operators. The time evolution in this frame is simply a rotation followed by relaxation at modified rates that play a role similar to the standard longitudinal and transverse rates. These rates are functions of the applied field, which

From explicit to implicit normal mode initialization of a limited-area model

Energy Technology Data Exchange (ETDEWEB)

Bijlsma, S.J.

2013-02-15

In this note the implicit normal mode initialization of a limited-area model is discussed from a different point of view. To that end it is shown that the equations describing the explicit normal mode initialization applied to the shallow water equations in differentiated form on the sphere can readily be derived in normal mode space if the model equations are separable, but only in the case of stationary Rossby modes can be transformed into the implicit equations in physical space. This is a consequence of the simple relations between the components of the different modes in that case. In addition a simple eigenvalue problem is given for the frequencies of the gravity waves. (orig.)

Conservation properties of numerical integration methods for systems of ordinary differential equations

Science.gov (United States)

Rosenbaum, J. S.

1976-01-01

If a system of ordinary differential equations represents a property conserving system that can be expressed linearly (e.g., conservation of mass), it is then desirable that the numerical integration method used conserve the same quantity. It is shown that both linear multistep methods and Runge-Kutta methods are 'conservative' and that Newton-type methods used to solve the implicit equations preserve the inherent conservation of the numerical method. It is further shown that a method used by several authors is not conservative.

About Navier-Stokes Equation in the Theory of Convective Heat Transfer

Science.gov (United States)

Davidzon, M. Y.

2017-10-01

A system of differential equations (Navier-Stokes, continuity, heat conductivity) is used to solve convective heat transfer problems. While solving Navier-Stokes equation, it is usually assumed that tangent stress is proportional to the velocity gradient. This assumption is valid with a small velocity gradient, for example, near an axis of the channel, but velocity gradient can be very large near the channel wall. Our paper shows that if we accept power law instead of linear law for tangential stress, then the velocity profile for creeping, laminar, and turbulent flow in the channel can be calculated without using Navier-Stokes equation. Also, in this case Navier-Stokes equation itself changes: the coefficient of dynamic viscosity changes its value from normal (in case of the creeping flow) to tending to infinity (in case of the well-developed turbulent flow).

Exact results in the large system size limit for the dynamics of the chemical master equation, a one dimensional chain of equations.

Science.gov (United States)

Martirosyan, A; Saakian, David B

2011-08-01

We apply the Hamilton-Jacobi equation (HJE) formalism to solve the dynamics of the chemical master equation (CME). We found exact analytical expressions (in large system-size limit) for the probability distribution, including explicit expression for the dynamics of variance of distribution. We also give the solution for some simple cases of the model with time-dependent rates. We derived the results of the Van Kampen method from the HJE approach using a special ansatz. Using the Van Kampen method, we give a system of ordinary differential equations (ODEs) to define the variance in a two-dimensional case. We performed numerics for the CME with stationary noise. We give analytical criteria for the disappearance of bistability in the case of stationary noise in one-dimensional CMEs.

Quantum Discord in Two-Qubit System Constructed from the Yang—Baxter Equation

International Nuclear Information System (INIS)

Gou Li-Dan; Wang Xiao-Qian; Sun Yuan-Yuan; Xu Yu-Mei

2014-01-01

Quantum correlations among parts of a composite quantum system are a fundamental resource for several applications in quantum information. In general, quantum discord can measure quantum correlations. In that way, we investigate the quantum discord of the two-qubit system constructed from the Yang—Baxter Equation. The density matrix of this system is generated through the unitary Yang—Baxter matrix R. The analytical expression and numerical result of quantum discord and geometric measure of quantum discord are obtained for the Yang—Baxter system. These results show that quantum discord and geometric measure of quantum discord are only connect with the parameter θ, which is the important spectral parameter in Yang—Baxter equation. (general)

A nonlinear boundary integral equations method for the solving of quasistatic elastic contact problem with Coulomb friction

Directory of Open Access Journals (Sweden)

Yurii M. Streliaiev

2016-06-01

Full Text Available Three-dimensional quasistatic contact problem of two linearly elastic bodies' interaction with Coulomb friction taken into account is considered. The boundary conditions of the problem have been simplified by the modification of the Coulomb's law of friction. This modification is based on the introducing of a delay in normal contact tractions that bound tangent contact tractions in the Coulomb's law of friction expressions. At this statement the problem is reduced to a sequence of similar systems of nonlinear integral equations describing bodies' interaction at each step of loading. A method for an approximate solution of the integral equations system corresponded to each step of loading is applied. This method consists of system regularization, discretization of regularized system and iterative process application for solving the discretized system. A numerical solution of a contact problem of an elastic sphere with an elastic half-space interaction under increasing and subsequently decreasing normal compressive force has been obtained.

Equation-free model reduction for complex dynamical systems

International Nuclear Information System (INIS)

Le Maitre, O. P.; Mathelin, L.; Le Maitre, O. P.

2010-01-01

This paper presents a reduced model strategy for simulation of complex physical systems. A classical reduced basis is first constructed relying on proper orthogonal decomposition of the system. Then, unlike the alternative approaches, such as Galerkin projection schemes for instance, an equation-free reduced model is constructed. It consists in the determination of an explicit transformation, or mapping, for the evolution over a coarse time-step of the projection coefficients of the system state on the reduced basis. The mapping is expressed as an explicit polynomial transformation of the projection coefficients and is computed once and for all in a pre-processing stage using the detailed model equation of the system. The reduced system can then be advanced in time by successive applications of the mapping. The CPU cost of the method lies essentially in the mapping approximation which is performed offline, in a parallel fashion, and only once. Subsequent application of the mapping to perform a time-integration is carried out at a low cost thanks to its explicit character. Application of the method is considered for the 2-D flow around a circular cylinder. We investigate the effectiveness of the reduced model in rendering the dynamics for both asymptotic state and transient stages. It is shown that the method leads to a stable and accurate time-integration for only a fraction of the cost of a detailed simulation, provided that the mapping is properly approximated and the reduced basis remains relevant for the dynamics investigated. (authors)

Forward-backward equations for nonlinear propagation in axially invariant optical systems

International Nuclear Information System (INIS)

Ferrando, Albert; Zacares, Mario; Fernandez de Cordoba, Pedro; Binosi, Daniele; Montero, Alvaro

2005-01-01

We present a general framework to deal with forward and backward components of the electromagnetic field in axially invariant nonlinear optical systems, which include those having any type of linear or nonlinear transverse inhomogeneities. With a minimum amount of approximations, we obtain a system of two first-order equations for forward and backward components, explicitly showing the nonlinear couplings among them. The modal approach used allows for an effective reduction of the dimensionality of the original problem from 3+1 (three spatial dimensions plus one time dimension) to 1+1 (one spatial dimension plus one frequency dimension). The new equations can be written in a spinor Dirac-like form, out of which conserved quantities can be calculated in an elegant manner. Finally, these equations inherently incorporate spatiotemporal couplings, so that they can be easily particularized to deal with purely temporal or purely spatial effects. Nonlinear forward pulse propagation and nonparaxial evolution of spatial structures are analyzed as examples

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

Science.gov (United States)

Abdulwahhab, Muhammad Alim

2016-10-01

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

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

International Nuclear Information System (INIS)

Velarde, G.

1976-01-01

Using some simplifying hypotheses, different expressions of the Boltzmann integrodifferential equation are obtained. Posteriorly, they are applied to some particular cases: slowing down, thermalization, multigroups, critical reactors and virtual critical reactors with k, α and lambda. (author)

On the interpretations of Langevin stochastic equation in different coordinate systems

International Nuclear Information System (INIS)

Martinez, E.; Lopez-Diaz, L.; Torres, L.; Alejos, O.

2004-01-01

The stochastic Langevin Landau-Lifshitz equation is usually utilized in micromagnetics formalism to account for thermal effects. Commonly, two different interpretations of the stochastic integrals can be made: Ito and Stratonovich. In this work, the Langevin-Landau-Lifshitz (LLL) equation is written in both Cartesian and Spherical coordinates. If Spherical coordinates are employed, the noise is additive, and therefore, Ito and Stratonovich solutions are equal. This is not the case when (LLL) equation is written in Cartesian coordinates. In this case, the Langevin equation must be interpreted in the Stratonovich sense in order to reproduce correct statistical results. Nevertheless, the statistics of the numerical results obtained from Euler-Ito and Euler-Stratonovich schemes are equivalent due to the additional numerical constraint imposed in Cartesian system after each time step, which itself assures that the magnitude of the magnetization is preserved

Normal-Mode Splitting in a Weakly Coupled Optomechanical System

Science.gov (United States)

Rossi, Massimiliano; Kralj, Nenad; Zippilli, Stefano; Natali, Riccardo; Borrielli, Antonio; Pandraud, Gregory; Serra, Enrico; Di Giuseppe, Giovanni; Vitali, David

2018-02-01

Normal-mode splitting is the most evident signature of strong coupling between two interacting subsystems. It occurs when two subsystems exchange energy between themselves faster than they dissipate it to the environment. Here we experimentally show that a weakly coupled optomechanical system at room temperature can manifest normal-mode splitting when the pump field fluctuations are antisquashed by a phase-sensitive feedback loop operating close to its instability threshold. Under these conditions the optical cavity exhibits an effectively reduced decay rate, so that the system is effectively promoted to the strong coupling regime.

On Coupled System of Navier-Stokes Equations and Temperature

African Journals Online (AJOL)

Dr. Anthony Peter

ABSTRACT. This paper deals with the coupled system of Navier-Stokes equations and temperature (Thermohydraulics) in a strip in the class of spatially non-decaying (infinite-energy) solutions belonging to the properly chosen uniformly local Sobolev spaces. The global well-posedness and dissipativity of the Navier- ...

Master equations for degenerate systems: electron radiative cascade in a Coulomb potential

International Nuclear Information System (INIS)

Uskov, D B; Pratt, R H

2004-01-01

We examine the effects of degeneracy and its lifting for the problem of electron radiative cascade, described by master equations of the Lindblad form (quantum optical master equations). A weak external field approximation is used to study the resulting gradual transformation of cascade dynamics between degenerate and non-degenerate forms. Exploiting the spherical symmetry properties of the system we demonstrate significant difference between perturbations commuting with angular momentum and perturbations breaking the spherical symmetry, such as a homogeneous external field. We discuss the possibility and the general approach for reduction of the Lindblad master equations in the case of spectral degeneracy to the Pauli balance equations. This determines the appropriate choice of basis as, for example, spherical or parabolic

The equationally-defined commutator a study in equational logic and algebra

CERN Document Server

Czelakowski, Janusz

2015-01-01

This monograph introduces and explores the notions of a commutator equation and the equationally-defined commutator from the perspective of abstract algebraic logic.  An account of the commutator operation associated with equational deductive systems is presented, with an emphasis placed on logical aspects of the commutator for equational systems determined by quasivarieties of algebras.  The author discusses the general properties of the equationally-defined commutator, various centralization relations for relative congruences, the additivity and correspondence properties of the equationally-defined commutator, and its behavior in finitely generated quasivarieties. Presenting new and original research not yet considered in the mathematical literature, The Equationally-Defined Commutator will be of interest to professional algebraists and logicians, as well as graduate students and other researchers interested in problems of modern algebraic logic.

Integrable motion of curves in self-consistent potentials: Relation to spin systems and soliton equations

Energy Technology Data Exchange (ETDEWEB)

Myrzakulov, R.; Mamyrbekova, G.K.; Nugmanova, G.N.; Yesmakhanova, K.R. [Eurasian International Center for Theoretical Physics and Department of General and Theoretical Physics, Eurasian National University, Astana 010008 (Kazakhstan); Lakshmanan, M., E-mail: lakshman@cnld.bdu.ac.in [Centre for Nonlinear Dynamics, School of Physics, Bharathidasan University, Tiruchirapalli 620 024 (India)

2014-06-13

Motion of curves and surfaces in R{sup 3} lead to nonlinear evolution equations which are often integrable. They are also intimately connected to the dynamics of spin chains in the continuum limit and integrable soliton systems through geometric and gauge symmetric connections/equivalence. Here we point out the fact that a more general situation in which the curves evolve in the presence of additional self-consistent vector potentials can lead to interesting generalized spin systems with self-consistent potentials or soliton equations with self-consistent potentials. We obtain the general form of the evolution equations of underlying curves and report specific examples of generalized spin chains and soliton equations. These include principal chiral model and various Myrzakulov spin equations in (1+1) dimensions and their geometrically equivalent generalized nonlinear Schrödinger (NLS) family of equations, including Hirota–Maxwell–Bloch equations, all in the presence of self-consistent potential fields. The associated gauge equivalent Lax pairs are also presented to confirm their integrability. - Highlights: • Geometry of continuum spin chain with self-consistent potentials explored. • Mapping on moving space curves in R{sup 3} in the presence of potential fields carried out. • Equivalent generalized nonlinear Schrödinger (NLS) family of equations identified. • Integrability of identified nonlinear systems proved by deducing appropriate Lax pairs.

Local thermodynamics and the generalized Gibbs-Duhem equation in systems with long-range interactions.

Science.gov (United States)

Latella, Ivan; Pérez-Madrid, Agustín

2013-10-01

The local thermodynamics of a system with long-range interactions in d dimensions is studied using the mean-field approximation. Long-range interactions are introduced through pair interaction potentials that decay as a power law in the interparticle distance. We compute the local entropy, Helmholtz free energy, and grand potential per particle in the microcanonical, canonical, and grand canonical ensembles, respectively. From the local entropy per particle we obtain the local equation of state of the system by using the condition of local thermodynamic equilibrium. This local equation of state has the form of the ideal gas equation of state, but with the density depending on the potential characterizing long-range interactions. By volume integration of the relation between the different thermodynamic potentials at the local level, we find the corresponding equation satisfied by the potentials at the global level. It is shown that the potential energy enters as a thermodynamic variable that modifies the global thermodynamic potentials. As a result, we find a generalized Gibbs-Duhem equation that relates the potential energy to the temperature, pressure, and chemical potential. For the marginal case where the power of the decaying interaction potential is equal to the dimension of the space, the usual Gibbs-Duhem equation is recovered. As examples of the application of this equation, we consider spatially uniform interaction potentials and the self-gravitating gas. We also point out a close relationship with the thermodynamics of small systems.

Biorthogonal Systems Approximating the Solution of the Nonlinear Volterra Integro-Differential Equation

Directory of Open Access Journals (Sweden)

Berenguer MI

2010-01-01

Full Text Available This paper deals with obtaining a numerical method in order to approximate the solution of the nonlinear Volterra integro-differential equation. We define, following a fixed-point approach, a sequence of functions which approximate the solution of this type of equation, due to some properties of certain biorthogonal systems for the Banach spaces and .

Variables and equations in hybrid systems with structural changes

NARCIS (Netherlands)

Beek, van D.A.

2001-01-01

In many models of physical systems, structural changes are common. Such structural changes may cause a variable to change from a differential variable to an algebraic variable, or to a variable that is not defined by an equation at all. Most hybrid modelling languages either restrict the kind of

Multiple positive solutions to a coupled systems of nonlinear fractional differential equations.

Science.gov (United States)

Shah, Kamal; Khan, Rahmat Ali

2016-01-01

In this article, we study existence, uniqueness and nonexistence of positive solution to a highly nonlinear coupled system of fractional order differential equations. Necessary and sufficient conditions for the existence and uniqueness of positive solution are developed by using Perov's fixed point theorem for the considered problem. Further, we also established sufficient conditions for existence of multiplicity results for positive solutions. Also, we developed some conditions under which the considered coupled system of fractional order differential equations has no positive solution. Appropriate examples are also provided which demonstrate our results.

On generalized fractional vibration equation

International Nuclear Information System (INIS)

Dai, Hongzhe; Zheng, Zhibao; Wang, Wei

2017-01-01

Highlights: • The paper presents a generalized fractional vibration equation for arbitrary viscoelastically damped system. • Some classical vibration equations can be derived from the developed equation. • The analytic solution of developed equation is derived under some special cases. • The generalized equation is particularly useful for developing new fractional equivalent linearization method. - Abstract: In this paper, a generalized fractional vibration equation with multi-terms of fractional dissipation is developed to describe the dynamical response of an arbitrary viscoelastically damped system. It is shown that many classical equations of motion, e.g., the Bagley–Torvik equation, can be derived from the developed equation. The Laplace transform is utilized to solve the generalized equation and the analytic solution under some special cases is derived. Example demonstrates the generalized transfer function of an arbitrary viscoelastic system.

A New Numerical Technique for Solving Systems Of Nonlinear Fractional Partial Differential Equations

Directory of Open Access Journals (Sweden)

Mountassir Hamdi Cherif

2017-11-01

Full Text Available In this paper, we apply an efficient method called the Aboodh decomposition method to solve systems of nonlinear fractional partial differential equations. This method is a combined form of Aboodh transform with Adomian decomposition method. The theoretical analysis of this investigated for systems of nonlinear fractional partial differential equations is calculated in the explicit form of a power series with easily computable terms. Some examples are given to shows that this method is very efficient and accurate. This method can be applied to solve others nonlinear systems problems.

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

International Nuclear Information System (INIS)

Frank, T.D.

2003-01-01

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

Research on the optimal dynamical systems of three-dimensional Navier-Stokes equations based on weighted residual

Science.gov (United States)

Peng, NaiFu; Guan, Hui; Wu, ChuiJie

2016-04-01

In this paper, the theory of constructing optimal dynamical systems based on weighted residual presented by Wu & Sha is applied to three-dimensional Navier-Stokes equations, and the optimal dynamical system modeling equations are derived. Then the multiscale global optimization method based on coarse graining analysis is presented, by which a set of approximate global optimal bases is directly obtained from Navier-Stokes equations and the construction of optimal dynamical systems is realized. The optimal bases show good properties, such as showing the physical properties of complex flows and the turbulent vortex structures, being intrinsic to real physical problem and dynamical systems, and having scaling symmetry in mathematics, etc.. In conclusion, using fewer terms of optimal bases will approach the exact solutions of Navier-Stokes equations, and the dynamical systems based on them show the most optimal behavior.

p-Euler equations and p-Navier-Stokes equations

Science.gov (United States)

Li, Lei; Liu, Jian-Guo

2018-04-01

We propose in this work new systems of equations which we call p-Euler equations and p-Navier-Stokes equations. p-Euler equations are derived as the Euler-Lagrange equations for the action represented by the Benamou-Brenier characterization of Wasserstein-p distances, with incompressibility constraint. p-Euler equations have similar structures with the usual Euler equations but the 'momentum' is the signed (p - 1)-th power of the velocity. In the 2D case, the p-Euler equations have streamfunction-vorticity formulation, where the vorticity is given by the p-Laplacian of the streamfunction. By adding diffusion presented by γ-Laplacian of the velocity, we obtain what we call p-Navier-Stokes equations. If γ = p, the a priori energy estimates for the velocity and momentum have dual symmetries. Using these energy estimates and a time-shift estimate, we show the global existence of weak solutions for the p-Navier-Stokes equations in Rd for γ = p and p ≥ d ≥ 2 through a compactness criterion.

Dirac delta representation by exact parametric equations.. Application to impulsive vibration systems

Science.gov (United States)

Chicurel-Uziel, Enrique

2007-08-01

A pair of closed parametric equations are proposed to represent the Heaviside unit step function. Differentiating the step equations results in two additional parametric equations, that are also hereby proposed, to represent the Dirac delta function. These equations are expressed in algebraic terms and are handled by means of elementary algebra and elementary calculus. The proposed delta representation complies exactly with the values of the definition. It complies also with the sifting property and the requisite unit area and its Laplace transform coincides with the most general form given in the tables. Furthermore, it leads to a very simple method of solution of impulsive vibrating systems either linear or belonging to a large class of nonlinear problems. Two example solutions are presented.

Boltzmann’s Six-Moment One-Dimensional Nonlinear System Equations with the Maxwell-Auzhan Boundary Conditions

Directory of Open Access Journals (Sweden)

A. Sakabekov

2016-01-01

Full Text Available We prove existence and uniqueness of the solution of the problem with initial and Maxwell-Auzhan boundary conditions for nonstationary nonlinear one-dimensional Boltzmann’s six-moment system equations in space of functions continuous in time and summable in square by a spatial variable. In order to obtain a priori estimation of the initial and boundary value problem for nonstationary nonlinear one-dimensional Boltzmann’s six-moment system equations we get the integral equality and then use the spherical representation of vector. Then we obtain the initial value problem for Riccati equation. We have managed to obtain a particular solution of this equation in an explicit form.

Interpretation of metabolic memory phenomenon using a physiological systems model: What drives oxidative stress following glucose normalization?

Science.gov (United States)

Voronova, Veronika; Zhudenkov, Kirill; Helmlinger, Gabriel; Peskov, Kirill

2017-01-01

Hyperglycemia is generally associated with oxidative stress, which plays a key role in diabetes-related complications. A complex, quantitative relationship has been established between glucose levels and oxidative stress, both in vitro and in vivo. For example, oxidative stress is known to persist after glucose normalization, a phenomenon described as metabolic memory. Also, uncontrolled glucose levels appear to be more detrimental to patients with diabetes (non-constant glucose levels) vs. patients with high, constant glucose levels. The objective of the current study was to delineate the mechanisms underlying such behaviors, using a mechanistic physiological systems modeling approach that captures and integrates essential underlying pathophysiological processes. The proposed model was based on a system of ordinary differential equations. It describes the interplay between reactive oxygen species production potential (ROS), ROS-induced cell alterations, and subsequent adaptation mechanisms. Model parameters were calibrated using different sources of experimental information, including ROS production in cell cultures exposed to various concentration profiles of constant and oscillating glucose levels. The model adequately reproduced the ROS excess generation after glucose normalization. Such behavior appeared to be driven by positive feedback regulations between ROS and ROS-induced cell alterations. The further oxidative stress-related detrimental effect as induced by unstable glucose levels can be explained by inability of cells to adapt to dynamic environment. Cell adaptation to instable high glucose declines during glucose normalization phases, and further glucose increase promotes similar or higher oxidative stress. In contrast, gradual ROS production potential decrease, driven by adaptation, is observed in cells exposed to constant high glucose.

Planck constant as spectral parameter in integrable systems and KZB equations

OpenAIRE

Levin, A.NRU HSE, Department of Mathematics, Myasnitskaya str. 20, Moscow, 101000, Russia; Olshanetsky, M.(ITEP, B. Cheremushkinskaya str. 25, Moscow, 117218, Russia); Zotov, A.(ITEP, B. Cheremushkinskaya str. 25, Moscow, 117218, Russia)

2014-01-01

We construct special rational ${\\rm gl}_N$ Knizhnik-Zamolodchikov-Bernard (KZB) equations with $\\tilde N$ punctures by deformation of the corresponding quantum ${\\rm gl}_N$ rational $R$-matrix. They have two parameters. The limit of the first one brings the model to the ordinary rational KZ equation. Another one is $\\tau$. At the level of classical mechanics the deformation parameter $\\tau$ allows to extend the previously obtained modified Gaudin models to the modified Schlesinger systems. Ne...

Radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems

International Nuclear Information System (INIS)

Liu, L.H.; Zhang, L.; Tan, H.P.

2006-01-01

In graded index medium, the ray goes along a curved path determined by Fermat principle, and the curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectory, the methods not based on ray-tracing technique need to be developed for the solution of radiative transfer in graded index medium. For this purpose, in this paper the streaming operator along a curved ray trajectory in original radiative transfer equation for graded index medium is transformed and expressed in spatial and angular ordinates and the radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems are derived. The conservative and the non-conservative forms of radiative transfer equation for three-dimensional graded index medium are given, which can be used as base equations to develop the numerical simulation methods, such as finite volume method, discrete ordinates method, and finite element method, for radiative transfer in graded index medium in cylindrical and spherical coordinate systems

Global existence and decay of solutions of a nonlinear system of wave equations

KAUST Repository

Said-Houari, Belkacem

2012-01-01

This work is concerned with a system of two wave equations with nonlinear damping and source terms acting in both equations. Under some restrictions on the nonlinearity of the damping and the source terms, we show that our problem has a unique local solution. Also, we prove that, for some restrictions on the initial data, the rate of decay of the total energy is exponential or polynomial depending on the exponents of the damping terms in both equations.

Global existence and decay of solutions of a nonlinear system of wave equations

KAUST Repository

Said-Houari, Belkacem

2012-03-01

This work is concerned with a system of two wave equations with nonlinear damping and source terms acting in both equations. Under some restrictions on the nonlinearity of the damping and the source terms, we show that our problem has a unique local solution. Also, we prove that, for some restrictions on the initial data, the rate of decay of the total energy is exponential or polynomial depending on the exponents of the damping terms in both equations.

Computing generalized Langevin equations and generalized Fokker-Planck equations.

Science.gov (United States)

Darve, Eric; Solomon, Jose; Kia, Amirali

2009-07-07

The Mori-Zwanzig formalism is an effective tool to derive differential equations describing the evolution of a small number of resolved variables. In this paper we present its application to the derivation of generalized Langevin equations and generalized non-Markovian Fokker-Planck equations. We show how long time scales rates and metastable basins can be extracted from these equations. Numerical algorithms are proposed to discretize these equations. An important aspect is the numerical solution of the orthogonal dynamics equation which is a partial differential equation in a high dimensional space. We propose efficient numerical methods to solve this orthogonal dynamics equation. In addition, we present a projection formalism of the Mori-Zwanzig type that is applicable to discrete maps. Numerical applications are presented from the field of Hamiltonian systems.

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; Sur la stabilite, les solutions periodiques et la resolution de certaines categories d'equations et systemes d'equations differentielles couplees non lineaires apparaissant dans les oscillations betatroniques

Energy Technology Data Exchange (ETDEWEB)

Valat, J [Commissariat a l' Energie Atomique, Saclay (France). Centre d' Etudes Nucleaires

1960-12-15

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) [French] Pour les equations du genre de Hill-Meissner a coefficients creneles, on a calcule des diagrammes universels de stabilite et ceux-ci ont ete verifies experimentalement. L'etude de ces equations dans le plan de phase a permis ensuite d'etendre le calcul des solutions periodiques au cas des equations differentielles non lineaires a coefficients periodiques creneles. Cette theorie a ete verifiee experimentalement. Pour Jes systemes couples non lineaires a coefficients constants, on a d'abord cherche les solutions menant a des mouvements algebriques. Les fonctions elliptiques et fuchsiennes uniformisent de tels mouvements. L'etude de mouvements non algebriques est plus delicate, a part l'etude des mouvements de Lissajous non lineaires. Une analyse fonctionnelle montre qu'il est toutefois possible dans certains cas de decoupler le systeme et de trouver des solutions generales. Pour les

A Nonmonotone Line Search Filter Algorithm for the System of Nonlinear Equations

Directory of Open Access Journals (Sweden)

Zhong Jin

2012-01-01

Full Text Available We present a new iterative method based on the line search filter method with the nonmonotone strategy to solve the system of nonlinear equations. The equations are divided into two groups; some equations are treated as constraints and the others act as the objective function, and the two groups are just updated at the iterations where it is needed indeed. We employ the nonmonotone idea to the sufficient reduction conditions and filter technique which leads to a flexibility and acceptance behavior comparable to monotone methods. The new algorithm is shown to be globally convergent and numerical experiments demonstrate its effectiveness.

Comparison of equations for dosing of medications in renal impairment.

Science.gov (United States)

Khanal, Aarati; Peterson, Gregory M; Jose, Matthew D; Castelino, Ronald L

2017-06-01

The aim of this study is to determine the concordance among the Cockcroft-Gault, the Modification of Diet in Renal Disease and the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equations in hypothetical dosing of renally cleared medications. A total of 2163 patients prescribed at least one of the 31 renally cleared drugs under review were included in the study. Kidney function was estimated using the three equations. We compared actual prescribed dosages of the same drug with recommended dosages based on the kidney function as calculated by each of the equations and applying dosing recommendations in the Australian Medicines Handbook. There was a significant difference in the kidney function values estimated from the three equations (P < 0.001). Despite the good overall agreement in renal drug dosing, we found selected but potentially important discrepancies among the doses rendered from the equations. The CKD-EPI equation non-normalized for body surface area had a greater rate of concordance with the Cockcroft-Gault equation than the Modification of Diet in Renal Disease equation for renal drug dosing. There is need for a long-term multi-centre study in a diverse population to define the clinical effects of the discrepancies among the equations for drug dosing. Given the greater concordance of the non-normalized CKD-EPI equation with the Cockcroft-Gault equation for dosing, the recommendation by Kidney Health Australia and the United States National Kidney Disease Education Program that 'dosing based on either eCrCl or an eGFR with body surface area normalization removed are acceptable' seems suitable and practicable for the purpose of dosing of non-critical drugs in the primary care setting. © 2016 Asian Pacific Society of Nephrology.

Generalized Lorentz-Force equations

International Nuclear Information System (INIS)

Yamaleev, R.M.

2001-01-01

Guided by Nambu (n+1)-dimensional phase space formalism we build a new system of dynamic equations. These equations describe a dynamic state of the corporeal system composed of n subsystems. The dynamic equations are formulated in terms of dynamic variables of the subsystems as well as in terms of dynamic variables of the corporeal system. These two sets of variables are related respectively as roots and coefficients of the n-degree polynomial equation. In the special n=2 case, this formalism reproduces relativistic dynamics for the charged spinning particles

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.

Numerical simulation of stochastic point kinetic equation in the dynamical system of nuclear reactor

International Nuclear Information System (INIS)

Saha Ray, S.

2012-01-01

Highlights: ► In this paper stochastic neutron point kinetic equations have been analyzed. ► Euler–Maruyama method and Strong Taylor 1.5 order method have been discussed. ► These methods are applied for the solution of stochastic point kinetic equations. ► Comparison between the results of these methods and others are presented in tables. ► Graphs for neutron and precursor sample paths are also presented. -- Abstract: In the present paper, the numerical approximation methods, applied to efficiently calculate the solution for stochastic point kinetic equations () in nuclear reactor dynamics, are investigated. A system of Itô stochastic differential equations has been analyzed to model the neutron density and the delayed neutron precursors in a point nuclear reactor. The resulting system of Itô stochastic differential equations are solved over each time-step size. The methods are verified by considering different initial conditions, experimental data and over constant reactivities. The computational results indicate that the methods are simple and suitable for solving stochastic point kinetic equations. In this article, a numerical investigation is made in order to observe the random oscillations in neutron and precursor population dynamics in subcritical and critical reactors.

Normalization in Lie algebras via mould calculus and applications

Science.gov (United States)

Paul, Thierry; Sauzin, David

2017-11-01

We establish Écalle's mould calculus in an abstract Lie-theoretic setting and use it to solve a normalization problem, which covers several formal normal form problems in the theory of dynamical systems. The mould formalism allows us to reduce the Lie-theoretic problem to a mould equation, the solutions of which are remarkably explicit and can be fully described by means of a gauge transformation group. The dynamical applications include the construction of Poincaré-Dulac formal normal forms for a vector field around an equilibrium point, a formal infinite-order multiphase averaging procedure for vector fields with fast angular variables (Hamiltonian or not), or the construction of Birkhoff normal forms both in classical and quantum situations. As a by-product we obtain, in the case of harmonic oscillators, the convergence of the quantum Birkhoff form to the classical one, without any Diophantine hypothesis on the frequencies of the unperturbed Hamiltonians.

Solution of Fractional Order System of Bagley-Torvik Equation Using Evolutionary Computational Intelligence

Directory of Open Access Journals (Sweden)

Muhammad Asif Zahoor Raja

2011-01-01

Full Text Available A stochastic technique has been developed for the solution of fractional order system represented by Bagley-Torvik equation. The mathematical model of the equation was developed with the help of feed-forward artificial neural networks. The training of the networks was made with evolutionary computational intelligence based on genetic algorithm hybrid with pattern search technique. Designed scheme was successfully applied to different forms of the equation. Results are compared with standard approximate analytic, stochastic numerical solvers and exact solutions.

On computation of limiting coderivatives of the normal-cone mapping to inequality systems and their applications

Czech Academy of Sciences Publication Activity Database

Outrata, Jiří; Gfrerer, H.

2016-01-01

Roč. 65, č. 4 (2016), s. 671-700 ISSN 0233-1934 R&D Projects: GA ČR GAP402/12/1309 Institutional support: RVO:67985556 Keywords : limiting normal cone * metric regularity and subregularity * parameterized generalized equations Subject RIV: BA - General Mathematics Impact factor: 0.943, year: 2016 http://library.utia.cas.cz/separaty/2015/MTR/outrata-0447466.pdf

Mellin-Barnes representations of Feynman diagrams, linear systems of differential equations, and polynomial solutions

International Nuclear Information System (INIS)

Kalmykov, Mikhail Yu.; Kniehl, Bernd A.

2012-05-01

We argue that the Mellin-Barnes representations of Feynman diagrams can be used for obtaining linear systems of homogeneous differential equations for the original Feynman diagrams with arbitrary powers of propagators without recourse to the integration-by-parts technique. These systems of differential equation can be used (i) for the differential reductions to sets of basic functions and (ii) for counting the numbers of master-integrals.

Systems-based decomposition schemes for the approximate solution of multi-term fractional differential equations

Science.gov (United States)

Ford, Neville J.; Connolly, Joseph A.

2009-07-01

We give a comparison of the efficiency of three alternative decomposition schemes for the approximate solution of multi-term fractional differential equations using the Caputo form of the fractional derivative. The schemes we compare are based on conversion of the original problem into a system of equations. We review alternative approaches and consider how the most appropriate numerical scheme may be chosen to solve a particular equation.

Quantum dot as a spin-current diode: A master-equation approach

DEFF Research Database (Denmark)

Souza, F.M.; Egues, J.C.; Jauho, Antti-Pekka

2007-01-01

We report a study of spin-dependent transport in a system composed of a quantum dot coupled to a normal metal lead and a ferromagnetic lead NM-QD-FM. We use the master equation approach to calculate the spin-resolved currents in the presence of an external bias and an intradot Coulomb interaction....... We find that for a range of positive external biases current flow from the normal metal to the ferromagnet the current polarization =I↑−I↓ / I↑+I↓ is suppressed to zero, while for the corresponding negative biases current flow from the ferromagnet to the normal metal attains a relative maximum value....... The system thus operates as a rectifier for spin-current polarization. This effect follows from an interplay between Coulomb interaction and nonequilibrium spin accumulation in the dot. In the parameter range considered, we also show that the above results can be obtained via nonequilibrium Green functions...

On separable Pauli equations

International Nuclear Information System (INIS)

Zhalij, Alexander

2002-01-01

We classify (1+3)-dimensional Pauli equations for a spin-(1/2) particle interacting with the electro-magnetic field, that are solvable by the method of separation of variables. As a result, we obtain the 11 classes of vector-potentials of the electro-magnetic field A(t,x(vector sign))=(A 0 (t,x(vector sign)), A(vector sign)(t,x(vector sign))) providing separability of the corresponding Pauli equations. It is established, in particular, that the necessary condition for the Pauli equation to be separable into second-order matrix ordinary differential equations is its equivalence to the system of two uncoupled Schroedinger equations. In addition, the magnetic field has to be independent of spatial variables. We prove that coordinate systems and the vector-potentials of the electro-magnetic field providing the separability of the corresponding Pauli equations coincide with those for the Schroedinger equations. Furthermore, an efficient algorithm for constructing all coordinate systems providing the separability of Pauli equation with a fixed vector-potential of the electro-magnetic field is developed. Finally, we describe all vector-potentials A(t,x(vector sign)) that (a) provide the separability of Pauli equation, (b) satisfy vacuum Maxwell equations without currents, and (c) describe non-zero magnetic field

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

Science.gov (United States)

Gumral, Hasan

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

Results of numerically solving an integral equation for a two-fermion system

International Nuclear Information System (INIS)

Skachkov, N.B.; Solov'eva, T.M.

2003-01-01

A two-particle system is described by integral equations whose kernels are dependent on the total energy of the system. Such equations can be reduced to an eigenvalue problem featuring an eigenvalue-dependent operator. This nonlinear eigenvalue problem is solved by means of an iterative procedure developed by the present authors. The energy spectra of a two-fermion system formed by particles of identical masses are obtained for two cases, that where the total spin of the system is equal to zero and that where the total spin of the system is equal to unity. The splitting of the ground-state levels of positronium and dimuonium, the frequency of the transition from the ground state of orthopositronium to its first excited state, and the probabilities of parapositronium and paradimuonium decays are computed. The results obtained in this way are found to be in good agreement with experimental data

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

International Nuclear Information System (INIS)

Ohtani, Nobuo

1976-01-01

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

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

Science.gov (United States)

Park, K. C.; Belvin, W. Keith

1990-01-01

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

The connection of two-particle relativistic quantum mechanics with the Bethe-Salpeter equation

International Nuclear Information System (INIS)

Sazdjian, H.

1986-02-01

We show the formal equivalence between the wave equations of two-particle relativistic quantum mechanics, based on the manifestly covariant hamiltonian formalism with constraints, and the Bethe-Salpeter equation. This is achieved by algebraically transforming the latter so as to separate it into two independent equations which match the equations of hamiltonian relativistic quantum mechanics. The first equation determines the relative time evolution of the system, while the second one yields a three-dimensional eigenvalue equation. A connection is thus established between the Bethe-Salpeter wave function and its kernel on the one hand and the quantum mechanical wave function and interaction potential on the other. For the sector of solutions of the Bethe-Salpeter equation having non-relativistic limits, this relationship can be evaluated in perturbation theory. We also device a generalized form of the instantaneous approximation which simplifies the various expressions involved in the above relations. It also permits the evaluation of the normalization condition of the quantum mechanical wave function as a three-dimensional integral

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

International Nuclear Information System (INIS)

Cavazos-Cadena, Rolando; Hernandez-Hernandez, Daniel

2006-01-01

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

Model Selection and Risk Estimation with Applications to Nonlinear Ordinary Differential Equation Systems

DEFF Research Database (Denmark)

Mikkelsen, Frederik Vissing

eective computational tools for estimating unknown structures in dynamical systems, such as gene regulatory networks, which may be used to predict downstream eects of interventions in the system. A recommended algorithm based on the computational tools is presented and thoroughly tested in various......Broadly speaking, this thesis is devoted to model selection applied to ordinary dierential equations and risk estimation under model selection. A model selection framework was developed for modelling time course data by ordinary dierential equations. The framework is accompanied by the R software...... package, episode. This package incorporates a collection of sparsity inducing penalties into two types of loss functions: a squared loss function relying on numerically solving the equations and an approximate loss function based on inverse collocation methods. The goal of this framework is to provide...

A computational method for direct integration of motion equations of structural systems

International Nuclear Information System (INIS)

Brusa, L.; Ciacci, R.; Creco, A.; Rossi, F.

1975-01-01

The dynamic analysis of structural systems requires the solution of the matrix equations: Md 2 delta/dt(t) + Cddelta/dt(t) + Kdelta(t) = F(t). Many numerical methods are available for direct integration of this equation and their efficiency is due to the fulfillment of the following requirements: A reasonable order of accuracy must be obtained for the approximation of the response relevant to the first modes: the model contributions relevant to the eigenvalues with large real part must be essentially neglected. This paper presents a step-by-step numerical scheme for the integration of this equation which satisfies the requirements previously mentioned. (Auth.)

Normal zone detectors for a large number of inductively coupled coils

International Nuclear Information System (INIS)

Owen, E.W.; Shimer, D.W.

1983-01-01

In order to protect a set of inductively coupled superconducting magnets, it is necessary to locate and measure normal zone voltages that are small compared with the mutual and self-induced voltages. The method described in this paper uses two sets of voltage measurements to locate and measure one or more normal zones in any number of coupled coils. One set of voltages is the outputs of bridges that balance out the self-induced voltages. The other set of voltages can be the voltages across the coils, although alternatives are possible. The two sets of equations form a single combined set of equations. Each normal zone location or combination of normal zones has a set of these combined equations associated with it. It is demonstrated that the normal zone can be located and the correct set chosen, allowing determination of the size of the normal zone. Only a few operations take place in a working detector: multiplication of a constant, addition, and simple decision-making. In many cases the detector for each coil, although weakly linked to the other detectors, can be considered to be independent

Consistency of a system of equations: What does that mean?

NARCIS (Netherlands)

Still, Georg J.; Kern, Walter; Koelewijn, Jaap; Bomhoff, M.J.

2010-01-01

The concept of (structural) consistency also called structural solvability is an important basic tool for analyzing the structure of systems of equations. Our aim is to provide a sound and practically relevant meaning to this concept. The implications of consistency are expressed in terms of

Scalar potentials and the Dirac equation

International Nuclear Information System (INIS)

Bergerhoff, B.; Soff, G.

1994-01-01

The Dirac equation is solved for various types of scalar potentials. Energy eigenvalues and normalized bound-state wave functions are calculated analytically for a scalar 1/r-potential as well as for a mixed scalar and Coulomb 1/r-potential. Also continuum wave functions for positive and negative energies are derived. Similarly, we investigate the solutions of the Dirac equation for a scalar square-well potential. Relativistic wave functions for scalar Yukawa and exponential potentials are determined numerically. Finally, we also discuss solutions of the Dirac equation for scalar linear and quadratic potentials which are frequently used to simulate quark confinement. (orig.)

Cellular solutions for the Poisson equation in extended systems

International Nuclear Information System (INIS)

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

1994-01-01

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

Evaluating accounting information systems that support multiple GAAP reporting using Normalized Systems Theory

NARCIS (Netherlands)

Vanhoof, E.; Huysmans, P.; Aerts, Walter; Verelst, J.; Aveiro, D.; Tribolet, J.; Gouveia, D.

2014-01-01

This paper uses a mixed methods approach of design science and case study research to evaluate structures of Accounting Information Systems (AIS) that report in multiple Generally Accepted Accounting Principles (GAAP), using Normalized Systems Theory (NST). To comply with regulation, many companies

A New Normal Form for Multidimensional Mode Conversion

International Nuclear Information System (INIS)

Tracy, E. R.; Richardson, A. S.; Kaufman, A. N.; Zobin, N.

2007-01-01

Linear conversion occurs when two wave types, with distinct polarization and dispersion characteristics, are locally resonant in a nonuniform plasma [1]. In recent work, we have shown how to incorporate a ray-based (WKB) approach to mode conversion in numerical algorithms [2,3]. The method uses the ray geometry in the conversion region to guide the reduction of the full NxN-system of wave equations to a 2x2 coupled pair which can be solved and matched to the incoming and outgoing WKB solutions. The algorithm in [2] assumes the ray geometry is hyperbolic and that, in ray phase space, there is an 'avoided crossing', which is the most common type of conversion. Here, we present a new formulation that can deal with more general types of conversion [4]. This formalism is based upon the fact (first proved in [5]) that it is always possible to put the 2x2 wave equation into a 'normal' form, such that the diagonal elements of the dispersion matrix Poisson-commute with the off-diagonals (at leading order). Therefore, if we use the diagonals (rather than the eigenvalues or the determinant) of the dispersion matrix as ray Hamiltonians, the off-diagonals will be conserved quantities. When cast into normal form, the 2x2 dispersion matrix has a very natural physical interpretation: the diagonals are the uncoupled ray hamiltonians and the off-diagonals are the coupling. We discuss how to incorporate the normal form into ray tracing algorithms

Correlations between chaos in a perturbed sine-Gordon equation and a truncated model system

International Nuclear Information System (INIS)

Bishop, A.R.; Flesch, R.; Forests, M.G.; Overman, E.A.

1990-01-01

The purpose of this paper is to present a first step toward providing coordinates and associated dynamics for low-dimensional attractors in nearly integrable partial differential equations (pdes), in particular, where the truncated system reflects salient geometric properties of the pde. This is achieved by correlating: (1) numerical results on the bifurcations to temporal chaos with spatial coherence of the damped, periodically forced sine-Gordon equation with periodic boundary conditions; (2) an interpretation of the spatial and temporal bifurcation structures of this perturbed integrable system with regard to the exact structure of the sine-Gordon phase space; (3) a model dynamical systems problem, which is itself a perturbed integrable Hamiltonian system, derived from the perturbed sine-Gordon equation by a finite mode Fourier truncation in the nonlinear Schroedinger limit; and (4) the bifurcations to chaos in the truncated phase space. In particular, a potential source of chaos in both the pde and the model ordinary differential equation systems is focused on: the existence of homoclinic orbits in the unperturbed integrable phase space and their continuation in the perturbed problem. The evidence presented here supports the thesis that the chaotic attractors of the weakly perturbed periodic sine-Gordon system consists of low-dimensional metastable attacking states together with intermediate states that are O(1) unstable and correspond to homoclinic states in the integrable phase space. It is surmised that the chaotic dynamics on these attractors is due to the perturbation of these homocline integrable configurations

Dirac equation on a curved surface

Energy Technology Data Exchange (ETDEWEB)

Brandt, F.T., E-mail: fbrandt@usp.br; Sánchez-Monroy, J.A., E-mail: antosan@usp.br

2016-09-07

The dynamics of Dirac particles confined to a curved surface is examined employing the thin-layer method. We perform a perturbative expansion to first-order and split the Dirac field into normal and tangential components to the surface. In contrast to the known behavior of second order equations like Schrödinger, Maxwell and Klein–Gordon, we find that there is no geometric potential for the Dirac equation on a surface. This implies that the non-relativistic limit does not commute with the thin-layer method. Although this problem can be overcome when second-order terms are retained in the perturbative expansion, this would preclude the decoupling of the normal and tangential degrees of freedom. Therefore, we propose to introduce a first-order term which rescues the non-relativistic limit and also clarifies the effect of the intrinsic and extrinsic curvatures on the dynamics of the Dirac particles. - Highlights: • The thin-layer method is employed to derive the Dirac equation on a curved surface. • A geometric potential is absent at least to first-order in the perturbative expansion. • The effects of the extrinsic curvature are included to rescue the non-relativistic limit. • The resulting Dirac equation is consistent with the Heisenberg uncertainty principle.

Numerical Solution of Nonlinear Volterra Integral Equations System Using Simpson’s 3/8 Rule

Directory of Open Access Journals (Sweden)

Adem Kılıçman

2012-01-01

Full Text Available The Simpson’s 3/8 rule is used to solve the nonlinear Volterra integral equations system. Using this rule the system is converted to a nonlinear block system and then by solving this nonlinear system we find approximate solution of nonlinear Volterra integral equations system. One of the advantages of the proposed method is its simplicity in application. Further, we investigate the convergence of the proposed method and it is shown that its convergence is of order O(h4. Numerical examples are given to show abilities of the proposed method for solving linear as well as nonlinear systems. Our results show that the proposed method is simple and effective.

The impact of calcium assay change on a local adjusted calcium equation.

Science.gov (United States)

Davies, Sarah L; Hill, Charlotte; Bailey, Lisa M; Davison, Andrew S; Milan, Anna M

2016-03-01

Deriving and validating local adjusted calcium equations is important for ensuring appropriate calcium status classification. We investigated the impact on our local adjusted calcium equation of a change in calcium method by the manufacturer from cresolphthalein complexone to NM-BAPTA. Calcium and albumin results from general practice requests were extracted from the Laboratory Information Management system for a three-month period. Results for which there was evidence of disturbance in calcium homeostasis were excluded leaving 13,482 sets of results for analysis. The adjusted calcium equation was derived following least squares regression analysis of total calcium on albumin and normalized to the mean calcium concentration of the data-set. The revised equation (NM-BAPTA calcium method) was compared with the previous equation (cresolphthalein complexone calcium method). The switch in calcium assay resulted in a small change in the adjusted calcium equation but was not considered to be clinically significant. The calcium reference interval differed from that proposed by Pathology Harmony in the UK. Local adjusted calcium equations should be re-assessed following changes in the calcium method. A locally derived reference interval may differ from the consensus harmonized reference interval. © The Author(s) 2015.

Quarter-Sweep Iteration Concept on Conjugate Gradient Normal Residual Method via Second Order Quadrature - Finite Difference Schemes for Solving Fredholm Integro-Differential Equations

International Nuclear Information System (INIS)

Aruchunan, E.

2015-01-01

In this paper, we have examined the effectiveness of the quarter-sweep iteration concept on conjugate gradient normal residual (CGNR) iterative method by using composite Simpson's (CS) and finite difference (FD) discretization schemes in solving Fredholm integro-differential equations. For comparison purposes, Gauss- Seidel (GS) and the standard or full- and half-sweep CGNR methods namely FSCGNR and HSCGNR are also presented. To validate the efficacy of the proposed method, several analyses were carried out such as computational complexity and percentage reduction on the proposed and existing methods. (author)

On the validity of the Jarzynski equation in quantum systems; Zur Gueltigkeit der Jarzynskigleichung in Quantensystemen

Energy Technology Data Exchange (ETDEWEB)

Nolte, Roman

2009-11-20

Discovered in 1997, the Jarzynski equation is one of several new theorems of nonequilibrium thermodynamics. Not only this equation makes a more severe statement than the second law of thermodynamics, it does also relate process quantities from processes far from equilibrium to equilibrium quantities. In particular during the investigation of very small systems there has been drawn much attention to this equation and the related fluctuation theorems during the last years. Something similar applies for the description of microbiological processes which take place often stationary but rarely in thermodynamical equilibrium. However, especially according to small systems the question of the validity of the equation in the quantum case emerges. Though meanwhile quite comprehensive proofs concerning classical systems have been found, for that case uncertainty and contradictory statements exist, founding on different definitions and interpretations of the quantum analogon of expressions of the equation. Simple examples on which the different approaches can be tested, are so far missing. In this work two such examples are investigated and it is examined, how the results differ from their classical counterparts and which properties of quantum systems influence the result. (orig.)

4th International Conference on Particle Systems and Partial Differential Equations

CERN Document Server

Soares, Ana

2017-01-01

'This book addresses mathematical problems motivated by various applications in physics, engineering, chemistry and biology. It gathers the lecture notes from the mini-course presented by Jean-Christophe Mourrat on the construction of the various stochastic “basic” terms involved in the formulation of the dynamic Ö4  theory in three space dimensions, as well as selected contributions presented at the fourth meeting on Particle Systems and PDEs, which was held at the University of Minho’s Centre of Mathematics in December 2015. The purpose of the conference was to bring together prominent researchers working in the fields of particle systems and partial differential equations, offering them a forum to present their recent results and discuss their topics of expertise. The meeting was also intended to present to a vast and varied public, including young researchers, the area of interacting particle systems, its underlying motivation, and its relation to partial differential equations.  The book w...

Discrete Painlevé equations from Y-systems

International Nuclear Information System (INIS)

Hone, Andrew N W; Inoue, Rei

2014-01-01

We consider T-systems and Y-systems arising from cluster mutations applied to quivers that have the property of being periodic under a sequence of mutations. The corresponding nonlinear recurrences for cluster variables (coefficient-free T-systems) were described in the work of Fordy and Marsh, who completely classified all such quivers in the case of period 1, and characterized them in terms of the skew-symmetric exchange matrix B that defines the quiver. A broader notion of periodicity in general cluster algebras was introduced by Nakanishi, who also described the corresponding Y-systems, and T-systems with coefficients. A result of Fomin and Zelevinsky says that the coefficient-free T-system provides a solution of the Y-system. In this paper, we show that in general there is a discrepancy between these two systems, in the sense that the solution of the former does not correspond to the general solution of the latter. This discrepancy is removed by introducing additional non-autonomous coefficients into the T-system. In particular, we focus on the period 1 case and show that, when the exchange matrix B is degenerate, discrete Painlevé equations can arise from this construction. This article is part of a special issue of Journal of Physics A: Mathematical and Theoretical devoted to ‘Cluster algebras in mathematical physics’. (paper)

Symmetries of th-Order Approximate Stochastic Ordinary Differential Equations

OpenAIRE

Fredericks, E.; Mahomed, F. M.

2012-01-01

Symmetries of $n$ th-order approximate stochastic ordinary differential equations (SODEs) are studied. The determining equations of these SODEs are derived in an Itô calculus context. These determining equations are not stochastic in nature. SODEs are normally used to model nature (e.g., earthquakes) or for testing the safety and reliability of models in construction engineering when looking at the impact of random perturbations.

Perturbation Solutions for Random Linear Structural Systems subject to Random Excitation using Stochastic Differential Equations

DEFF Research Database (Denmark)

Köyluoglu, H.U.; Nielsen, Søren R.K.; Cakmak, A.S.

1994-01-01

perturbation method using stochastic differential equations. The joint statistical moments entering the perturbation solution are determined by considering an augmented dynamic system with state variables made up of the displacement and velocity vector and their first and second derivatives with respect......The paper deals with the first and second order statistical moments of the response of linear systems with random parameters subject to random excitation modelled as white-noise multiplied by an envelope function with random parameters. The method of analysis is basically a second order...... to the random parameters of the problem. Equations for partial derivatives are obtained from the partial differentiation of the equations of motion. The zero time-lag joint statistical moment equations for the augmented state vector are derived from the Itô differential formula. General formulation is given...

Kinetic equations for an unstable plasma; Equations cinetiques d'un plasma instable

Energy Technology Data Exchange (ETDEWEB)

Laval, G; Pellat, R [Commissariat a l' Energie Atomique, Fontenay-aux-Roses (France). Centre d' Etudes Nucleaires

1968-07-01

In this work, we establish the plasma kinetic equations starting from the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy of equations. We demonstrate that relations existing between correlation functions may help to justify the truncation of the hierarchy. Then we obtain the kinetic equations of a stable or unstable plasma. They do not reduce to an equation for the one-body distribution function, but generally involve two coupled equations for the one-body distribution function and the spectral density of the fluctuating electric field. We study limiting cases where the Balescu-Lenard equation, the quasi-linear theory, the Pines-Schrieffer equations and the equations of weak turbulence in the random phase approximation are recovered. At last we generalise the H-theorem for the system of equations and we define conditions for irreversible behaviour. (authors) [French] Dans ce travail nous etablissons les equations cinetiques d'un plasma a partir des equations de la recurrence de Bogoliubov, Born, Green, Kirkwood et Yvon. Nous demontrons qu'entre les fonctions de correlation d'un plasma existent des relations qui permettent de justifier la troncature de la recurrence. Nous obtenons alors les equations cinetiques d'un plasma stable ou instable. En general elles ne se reduisent pas a une equation d'evolution pour la densite simple, mais se composent de deux equations couplees portant sur la densite simple et la densite spectrale du champ electrique fluctuant. Nous etudions le cas limites ou l'on retrouve l'equation de Balescu-Lenard, les equations de la theorie quasi-lineaire, les equations de Pines et Schrieffer et les equations de la turbulence faible dans l'approximation des phases aleatoires. Enfin, nous generalisons le theoreme H pour ce systeme d'equations et nous precisons les conditions d'evolution irreversible. (auteurs)

Some physical applications of fractional Schroedinger equation

International Nuclear Information System (INIS)

Guo Xiaoyi; Xu Mingyu

2006-01-01

The fractional Schroedinger equation is solved for a free particle and for an infinite square potential well. The fundamental solution of the Cauchy problem for a free particle, the energy levels and the normalized wave functions of a particle in a potential well are obtained. In the barrier penetration problem, the reflection coefficient and transmission coefficient of a particle from a rectangular potential wall is determined. In the quantum scattering problem, according to the fractional Schroedinger equation, the Green's function of the Lippmann-Schwinger integral equation is given

Nonchaoticity of Ordinary Differential Equations Describing Autonomous Transcriptional Regulatory Circuits

International Nuclear Information System (INIS)

Li Pengfei; Hu Gang; Chen Runsheng

2008-01-01

Gene transcriptional regulation (TR) processes are often described by coupled nonlinear ordinary differential equations (ODEs). When the dimension of TR circuits is high (e.g. n ≥ 3) the motions of the corresponding ODEs may, very probably, show self-sustained oscillations and chaos. On the other hand, chaoticity may be harmful for the normal biological functions of TR processes. In this letter we numerically study the dynamics of 3-gene TR ODEs in great detail, and investigate many 4-, 5-, and 10-gene TR systems by randomly choosing figures and parameters in the conventionally accepted ranges. And we find that oscillations are very seldom and no chaotic motion is observed, even if the dimension of systems is sufficiently high (n ≥ 3). It is argued that the observation of nonchaoticity of these ODEs agrees with normal functions of actual TR processes

Gravitational field equations on and off a 3-brane world

International Nuclear Information System (INIS)

Aliev, A N; Guemruekcueoglu, A E

2004-01-01

The effective gravitational field equations on and off a 3-brane world possessing a Z 2 mirror symmetry and embedded in a five-dimensional bulk spacetime with cosmological constant were derived by Shiromizu, Maeda and Sasaki (SMS) in the framework of the Gauss-Codazzi projective approach with the subsequent specialization to the Gaussian normal coordinates in the neighbourhood of the brane. However, the Gaussian normal coordinates imply a very special slicing of spacetime and clearly, the consistent analysis of the brane dynamics would benefit from complete freedom in the slicing of spacetime, pushing the layer surfaces in the fifth dimension at any rates of evolution and in arbitrary positions. We rederive the SMS effective gravitational field equations on a 3-brane and generalize the off-brane equations to the case where there is an arbitrary energy-momentum tensor in the bulk. We use a more general setting to allow for acceleration of the normals to the brane surface through the lapse function and the shift vector in the spirit of Arnowitt, Deser and Misner. We show that the gravitational influence of the bulk spacetime on the brane may be described by a traceless second-rank tensor W ij , constructed from the 'electric' part of the bulk Riemann tensor. We also present the evolution equations for the tensor W ij , as well as for the corresponding 'magnetic' part of the bulk curvature. These equations involve terms determined by both the nonvanishing acceleration of normals in the nongeodesic slicing of spacetime and the presence of other fields in the bulk

Randomly transitional phenomena in the system governed by Duffing's equation

International Nuclear Information System (INIS)

Ueda, Yoshisuke.

1978-06-01

This paper deals with turbulent or chaotic phenomena which occur in the system governed by Duffing's equation, a special type of 2-dimensional periodic systems. By using analog and digital computers, experiments are undertaken with special reference to the changes of attractors and of average power spectra of the random processes under the variation of the system parameters. On the basis of the experimental results, an outline of the random process is made clear. The results obtained in this paper will be applied to the phenomena of the same kind which occur in 3-dimensional autonomous systems. (author)

Advanced-Retarded Differential Equations in Quantum Photonic Systems

Science.gov (United States)

Alvarez-Rodriguez, Unai; Perez-Leija, Armando; Egusquiza, Iñigo L.; Gräfe, Markus; Sanz, Mikel; Lamata, Lucas; Szameit, Alexander; Solano, Enrique

2017-01-01

We propose the realization of photonic circuits whose dynamics is governed by advanced-retarded differential equations. Beyond their mathematical interest, these photonic configurations enable the implementation of quantum feedback and feedforward without requiring any intermediate measurement. We show how this protocol can be applied to implement interesting delay effects in the quantum regime, as well as in the classical limit. Our results elucidate the potential of the protocol as a promising route towards integrated quantum control systems on a chip. PMID:28230090

Hybrid quantum-classical master equations

International Nuclear Information System (INIS)

Diósi, Lajos

2014-01-01

We discuss hybrid master equations of composite systems, which are hybrids of classical and quantum subsystems. A fairly general form of hybrid master equations is suggested. Its consistency is derived from the consistency of Lindblad quantum master equations. We emphasize that quantum measurement is a natural example of exact hybrid systems. We derive a heuristic hybrid master equation of time-continuous position measurement (monitoring). (paper)

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

Late effects of normal tissues (lent) scoring system: the soma scale

International Nuclear Information System (INIS)

Mornex, F.; Pavy, J.J.; Denekamp, J.

1997-01-01

Radiation tolerance of normal tissues remains the limiting factor for delivering tumoricidal dose. The late toxicity of normal tissues is the most critical element of an irradiation: somatic, functional and structural alterations occur during the actual treatment itself, but late effects manifest months to years after acute effects heal, and may progress with time. The optimal therapeutic ratio ultimately requires not only complete tumor clearance, but also minimal residual injury to surrounding vital normal tissues. The disparity between the intensity of acute and late effects and the inability to predict the eventual manifestation of late normal tissue injury has made radiation oncologists recognize the importance of careful patient follow-up. There is so far no uniform toxicity scoring system to compare several clinical studies in the absence of a 'common toxicity language'. This justifies the need to establish a precise evaluation system for the analysis of late effects of radiation on normal tissues. The SOMA/LENT scoring system results from an international collaboration. European Organization Treatment of Cancer (EORTC) and Radiation Therapy Oncology Group (RTOG) have created subcommittees with the aim of addressing the question of standardized toxic effects criteria. This effort appeared as a necessity to standardize and improve the data recording, to then describe and evaluate uniform toxicity at regular time intervals. The current proposed scale is not yet validated, and should be used cautiously. (authors)

On Generating Discrete Integrable Systems via Lie Algebras and Commutator Equations

International Nuclear Information System (INIS)

Zhang Yu-Feng; Tam, Honwah

2016-01-01

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

Efficient Instantiation of Parameterised Boolean Equation Systems to Parity Games

Directory of Open Access Journals (Sweden)

Gijs Kant

2012-10-01

Full Text Available Parameterised Boolean Equation Systems (PBESs are sequences of Boolean fixed point equations with data variables, used for, e.g., verification of modal mu-calculus formulae for process algebraic specifications with data. Solving a PBES is usually done by instantiation to a Parity Game and then solving the game. Practical game solvers exist, but the instantiation step is the bottleneck. We enhance the instantiation in two steps. First, we transform the PBES to a Parameterised Parity Game (PPG, a PBES with each equation either conjunctive or disjunctive. Then we use LTSmin, that offers transition caching, efficient storage of states and both distributed and symbolic state space generation, for generating the game graph. To that end we define a language module for LTSmin, consisting of an encoding of variables with parameters into state vectors, a grouped transition relation and a dependency matrix to indicate the dependencies between parts of the state vector and transition groups. Benchmarks on some large case studies, show that the method speeds up the instantiation significantly and decreases memory usage drastically.

Relations between nonlinear Riccati equations and other equations in fundamental physics

International Nuclear Information System (INIS)

Schuch, Dieter

2014-01-01

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

About local fractional three-dimensional compressible Navier-Stokes equations in Cantor-type cylindrical co-ordinate system

Directory of Open Access Journals (Sweden)

Gao Guo-Ping

2016-01-01

Full Text Available In this article, we investigate the local fractional 3-D compressible Navier-Stokes equation via local fractional derivative. We use the Cantor-type cylindrical co-ordinate method to transfer 3-D compressible Navier-Stokes equation from the Cantorian co-ordinate system to the Cantor-type cylindrical co-ordinate system.

Fuchsia : A tool for reducing differential equations for Feynman master integrals to epsilon form

Science.gov (United States)

Gituliar, Oleksandr; Magerya, Vitaly

2017-10-01

We present Fuchsia - an implementation of the Lee algorithm, which for a given system of ordinary differential equations with rational coefficients ∂x J(x , ɛ) = A(x , ɛ) J(x , ɛ) finds a basis transformation T(x , ɛ) , i.e., J(x , ɛ) = T(x , ɛ) J‧(x , ɛ) , such that the system turns into the epsilon form : ∂xJ‧(x , ɛ) = ɛ S(x) J‧(x , ɛ) , where S(x) is a Fuchsian matrix. A system of this form can be trivially solved in terms of polylogarithms as a Laurent series in the dimensional regulator ɛ. That makes the construction of the transformation T(x , ɛ) crucial for obtaining solutions of the initial system. In principle, Fuchsia can deal with any regular systems, however its primary task is to reduce differential equations for Feynman master integrals. It ensures that solutions contain only regular singularities due to the properties of Feynman integrals. Program Files doi:http://dx.doi.org/10.17632/zj6zn9vfkh.1 Licensing provisions: MIT Programming language:Python 2.7 Nature of problem: Feynman master integrals may be calculated from solutions of a linear system of differential equations with rational coefficients. Such a system can be easily solved as an ɛ-series when its epsilon form is known. Hence, a tool which is able to find the epsilon form transformations can be used to evaluate Feynman master integrals. Solution method: The solution method is based on the Lee algorithm (Lee, 2015) which consists of three main steps: fuchsification, normalization, and factorization. During the fuchsification step a given system of differential equations is transformed into the Fuchsian form with the help of the Moser method (Moser, 1959). Next, during the normalization step the system is transformed to the form where eigenvalues of all residues are proportional to the dimensional regulator ɛ. Finally, the system is factorized to the epsilon form by finding an unknown transformation which satisfies a system of linear equations. Additional comments

Kelvin Equation for a Non-Ideal Multicomponent Mixture

DEFF Research Database (Denmark)

Shapiro, Alexander; Stenby, Erling Halfdan

1997-01-01

The Kelvin equation is generalized by application to a case of a multicomponent non-ideal mixture. Such a generalization is necessary in order to describe the two-phase equilibrium in a capillary medium with respect to both normal and retrograde condensation. The equation obtained is applied...... to the equilibrium state of a hydrocarbon mixture ina gas-condensate reservoir....

New Numerical Solution of von Karman Equation of Lengthwise Rolling

Directory of Open Access Journals (Sweden)

Rudolf Pernis

2015-01-01

Full Text Available The calculation of average material contact pressure to rolls base on mathematical theory of rolling process given by Karman equation was solved by many authors. The solutions reported by authors are used simplifications for solution of Karman equation. The simplifications are based on two cases for approximation of the circular arch: (a by polygonal curve and (b by parabola. The contribution of the present paper for solution of two-dimensional differential equation of rolling is based on description of the circular arch by equation of a circle. The new term relative stress as nondimensional variable was defined. The result from derived mathematical models can be calculated following variables: normal contact stress distribution, front and back tensions, angle of neutral point, coefficient of the arm of rolling force, rolling force, and rolling torque during rolling process. Laboratory cold rolled experiment of CuZn30 brass material was performed. Work hardening during brass processing was calculated. Comparison of theoretical values of normal contact stress with values of normal contact stress obtained from cold rolling experiment was performed. The calculations were not concluded with roll flattening.

Simultaneous exact controllability for Maxwell equations and for a second-order hyperbolic system

Directory of Open Access Journals (Sweden)

Boris V. Kapitonov

2010-02-01

Full Text Available We present a result on "simultaneous" exact controllability for two models that describe two hyperbolic dynamics. One is the system of Maxwell equations and the other a vector-wave equation with a pressure term. We obtain the main result using modified multipliers in order to generate a necessary observability estimate which allow us to use the Hilbert Uniqueness Method (HUM introduced by Lions.

Methodology for the hybrid solution of systems of differential equations

International Nuclear Information System (INIS)

Larrinaga, E.F.; Lopez, M.A.

1993-01-01

This work shows a general methodology of solution to systems of differential equations in hybrid computers. Taking into account this methodology, a mathematical model was elaborated. It offers wide possibilities of recording and handling the results on the basis of using the hybrid system IBM-VIDAC 1224 which the ISCTN has. It also presents the results gained when simulating a simple model of a nuclear reactor, which was used in the validation of the results of the computational model

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.

Differential equations extended to superspace

International Nuclear Information System (INIS)

Torres, J.; Rosu, H.C.

2003-01-01

We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)

Differential equations extended to superspace

Energy Technology Data Exchange (ETDEWEB)

Torres, J. [Instituto de Fisica, Universidad de Guanajuato, A.P. E-143, Leon, Guanajuato (Mexico); Rosu, H.C. [Instituto Potosino de Investigacion Cientifica y Tecnologica, A.P. 3-74, Tangamanga, San Luis Potosi (Mexico)

2003-07-01

We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)

On the dynamics of a non-local parabolic equation arising from the Gierer-Meinhardt system

Science.gov (United States)

Kavallaris, Nikos I.; Suzuki, Takashi

2017-05-01

The purpose of the current paper is to contribute to the comprehension of the dynamics of the shadow system of an activator-inhibitor system known as a Gierer-Meinhardt model. Shadow systems are intended to work as an intermediate step between single equations and reaction-diffusion systems. In the case where the inhibitor’s response to the activator’s growth is rather weak, then the shadow system of the Gierer-Meinhardt model is reduced to a single though non-local equation whose dynamics will be investigated. We mainly focus on the derivation of blow-up results for this non-local equation which can be seen as instability patterns of the shadow system. In particular, a diffusion driven instability (DDI), or Turing instability, in the neighbourhood of a constant stationary solution, which it is destabilised via diffusion-driven blow-up, is obtained. The latter actually indicates the formation of some unstable patterns, whilst some stability results of global-in-time solutions towards non-constant steady states guarantee the occurrence of some stable patterns.

Augmentation of DAA Staggered – Solution Equations in Underwater Shock Problems for Singular Structural Mass Matrices

Directory of Open Access Journals (Sweden)

John A. DeRuntz Jr.

2005-01-01

Full Text Available The numerical solution of underwater shock fluid – structure interaction problems using boundary element/finite element techniques became tractable through the development of the family of Doubly Asymptotic Approximations (DAA. Practical implementation of the method has relied on the so-called augmentation of the DAA equations. The fluid and structural systems are respectively coupled by the structural acceleration vector in the surface normal direction on the right hand side of the DAA equations, and the total pressure applied to the structural equations on its right hand side. By formally solving for the acceleration vector from the structural system and substituting it into its place in the DAA equations, the augmentation introduces a term involving the inverse of the structural mass matrix. However there exist at least two important classes of problems in which the structural mass matrix is singular. This paper develops a method to carry out the augmentation for such problems using a generalized inverse technique.

A new hierarchy of generalized derivative nonlinear Schroedinger equations, its bi-Hamiltonian structure and finite-dimensional involutive system

International Nuclear Information System (INIS)

Yan, Z.; Zhang, H.

2001-01-01

In this paper, an isospectral problem and one associated with a new hierarchy of nonlinear evolution equations are presented. As a reduction, a representative system of new generalized derivative nonlinear Schroedinger equations in the hierarchy is given. It is shown that the hierarchy possesses bi-Hamiltonian structures by using the trace identity method and is Liouville integrable. The spectral problem is non linearized as a finite-dimensional completely integrable Hamiltonian system under a constraint between the potentials and spectral functions. Finally, the involutive solutions of the hierarchy of equations are obtained. In particular, the involutive solutions of the system of new generalized derivative nonlinear Schroedinger equations are developed

Normal anatomy of the lymphatic system in the CT-image

International Nuclear Information System (INIS)

Steinbrich, W.; Peters, P.E.

1982-01-01

To evaluate a pathologic process of a lymphatic node, detailed knowledge is required of the normal anatomy of the lumphatic system in an axial CT image. The anatomy is demonstrated in a comparative study before and after lymphography with CT-scans of patients with normal lymphadenographs. Hereby it appears that with the high-resolution scanning method and favourable imaging conditions even small lymphatic nodes can be differentiated without a lymphographic contrast technique. However, nerves and vessels cannot be differentiated. The extreme variability in the size of normal lymphatic nodes makes the differentiation of pathologic processes very difficult. (orig.) [de

Empirical Radiometric Normalization of Road Points from Terrestrial Mobile Lidar System

Directory of Open Access Journals (Sweden)

Tee-Ann Teo

2015-05-01

Full Text Available Lidar data provide both geometric and radiometric information. Radiometric information is influenced by sensor and target factors and should be calibrated to obtain consistent energy responses. The radiometric correction of airborne lidar system (ALS converts the amplitude into a backscatter cross-section with physical meaning value by applying a model-driven approach. The radiometric correction of terrestrial mobile lidar system (MLS is a challenging task because it does not completely follow the inverse square range function at near-range. This study proposed a radiometric normalization workflow for MLS using a data-driven approach. The scope of this study is to normalize amplitude of road points for road surface classification, assuming that road points from different scanners or strips should have similar responses in overlapped areas. The normalization parameters for range effect were obtained from crossroads. The experiment showed that the amplitude difference between scanners and strips decreased after radiometric normalization and improved the accuracy of road surface classification.

Trigonometric Solutions of WDVV Equations and Generalized Calogero-Moser-Sutherland Systems

Directory of Open Access Journals (Sweden)

Misha V. Feigin

2009-09-01

Full Text Available We consider trigonometric solutions of WDVV equations and derive geometric conditions when a collection of vectors with multiplicities determines such a solution. We incorporate these conditions into the notion of trigonometric Veselov system (v-system and we determine all trigonometric v-systems with up to five vectors. We show that generalized Calogero-Moser-Sutherland operator admits a factorized eigenfunction if and only if it corresponds to the trigonometric v-system; this inverts a one-way implication observed by Veselov for the rational solutions.

Superstatistical generalised Langevin equation: non-Gaussian viscoelastic anomalous diffusion

Science.gov (United States)

Ślęzak, Jakub; Metzler, Ralf; Magdziarz, Marcin

2018-02-01

Recent advances in single particle tracking and supercomputing techniques demonstrate the emergence of normal or anomalous, viscoelastic diffusion in conjunction with non-Gaussian distributions in soft, biological, and active matter systems. We here formulate a stochastic model based on a generalised Langevin equation in which non-Gaussian shapes of the probability density function and normal or anomalous diffusion have a common origin, namely a random parametrisation of the stochastic force. We perform a detailed analysis demonstrating how various types of parameter distributions for the memory kernel result in exponential, power law, or power-log law tails of the memory functions. The studied system is also shown to exhibit a further unusual property: the velocity has a Gaussian one point probability density but non-Gaussian joint distributions. This behaviour is reflected in the relaxation from a Gaussian to a non-Gaussian distribution observed for the position variable. We show that our theoretical results are in excellent agreement with stochastic simulations.

Continuous limits for an integrable coupling system of Toda equation hierarchy

International Nuclear Information System (INIS)

Li Li; Yu Fajun

2009-01-01

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

Continuous limits for an integrable coupling system of Toda equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

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

2009-09-21

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

Fitting and verification of frequency modulation systems on children with normal hearing.

Science.gov (United States)

Schafer, Erin C; Bryant, Danielle; Sanders, Katie; Baldus, Nicole; Algier, Katherine; Lewis, Audrey; Traber, Jordan; Layden, Paige; Amin, Aneeqa

2014-06-01

Several recent investigations support the use of frequency modulation (FM) systems in children with normal hearing and auditory processing or listening disorders such as those diagnosed with auditory processing disorders, autism spectrum disorders, attention-deficit hyperactivity disorder, Friedreich ataxia, and dyslexia. The American Academy of Audiology (AAA) published suggested procedures, but these guidelines do not cite research evidence to support the validity of the recommended procedures for fitting and verifying nonoccluding open-ear FM systems on children with normal hearing. Documenting the validity of these fitting procedures is critical to maximize the potential FM-system benefit in the above-mentioned populations of children with normal hearing and those with auditory-listening problems. The primary goal of this investigation was to determine the validity of the AAA real-ear approach to fitting FM systems on children with normal hearing. The secondary goal of this study was to examine speech-recognition performance in noise and loudness ratings without and with FM systems in children with normal hearing sensitivity. A two-group, cross-sectional design was used in the present study. Twenty-six typically functioning children, ages 5-12 yr, with normal hearing sensitivity participated in the study. Participants used a nonoccluding open-ear FM receiver during laboratory-based testing. Participants completed three laboratory tests: (1) real-ear measures, (2) speech recognition performance in noise, and (3) loudness ratings. Four real-ear measures were conducted to (1) verify that measured output met prescribed-gain targets across the 1000-4000 Hz frequency range for speech stimuli, (2) confirm that the FM-receiver volume did not exceed predicted uncomfortable loudness levels, and (3 and 4) measure changes to the real-ear unaided response when placing the FM receiver in the child's ear. After completion of the fitting, speech recognition in noise at a -5

The Mathlet Toolkit: Creating Dynamic Applets for Differential Equations and Dynamical Systems

Science.gov (United States)

Decker, Robert

2011-01-01

Dynamic/interactive graphing applets can be used to supplement standard computer algebra systems such as Maple, Mathematica, Derive, or TI calculators, in courses such as Calculus, Differential Equations, and Dynamical Systems. The addition of this type of software can lead to discovery learning, with students developing their own conjectures, and…

Multigrid solution of diffusion equations on distributed memory multiprocessor systems

International Nuclear Information System (INIS)

Finnemann, H.

1988-01-01

The subject is the solution of partial differential equations for simulation of the reactor core on high-performance computers. The parallelization and implementation of nodal multigrid diffusion algorithms on array and ring configurations of the DIRMU multiprocessor system is outlined. The particular iteration scheme employed in the nodal expansion method appears similarly efficient in serial and parallel environments. The combination of modern multi-level techniques with innovative hardware (vector-multiprocessor systems) provides powerful tools needed for real time simulation of physical systems. The parallel efficiencies range from 70 to 90%. The same performance is estimated for large problems on large multiprocessor systems being designed at present. (orig.) [de

A theory of post-stall transients in axial compression systems. I - Development of equations

Science.gov (United States)

Moore, F. K.; Greitzer, E. M.

1985-01-01

An approximate theory is presented for post-stall transients in multistage axial compression systems. The theory leads to a set of three simultaneous nonlinear third-order partial differential equations for pressure rise, and average and disturbed values of flow coefficient, as functions of time and angle around the compressor. By a Galerkin procedure, angular dependence is averaged, and the equations become first order in time. These final equations are capable of describing the growth and possible decay of a rotating-stall cell during a compressor mass-flow transient. It is shown how rotating-stall-like and surgelike motions are coupled through these equations, and also how the instantaneous compressor pumping characteristic changes during the transient stall process.

Normal variants and non-pathologic findings of the skeletal system in childhood

International Nuclear Information System (INIS)

Schaper, J.; Heinen, W.

2006-01-01

The knowledge of normal variants of the skeletal system in childhood protects the child from false radiologic diagnosis and the resulting malpractice. In this article we present a survey of common variants of childhood skeletal system. The awareness of these entities, e. g. the benign extraaxial fluid collections of infancy, allows accurate clinical and radiological diagnosis. The most important impact of radiological expertise in normal variations results in avoidance of unneccessary examinations and in omission of unsubstantiated parental and patients fears. Nearly all normal variations with definite radiologic findings should be treated according to the ''leave-me-alone'' principle. (orig.)

Contact Geometry of Hyperbolic Equations of Generic Type

Directory of Open Access Journals (Sweden)

Dennis The

2008-08-01

Full Text Available We study the contact geometry of scalar second order hyperbolic equations in the plane of generic type. Following a derivation of parametrized contact-invariants to distinguish Monge-Ampère (class 6-6, Goursat (class 6-7 and generic (class 7-7 hyperbolic equations, we use Cartan's equivalence method to study the generic case. An intriguing feature of this class of equations is that every generic hyperbolic equation admits at most a nine-dimensional contact symmetry algebra. The nine-dimensional bound is sharp: normal forms for the contact-equivalence classes of these maximally symmetric generic hyperbolic equations are derived and explicit symmetry algebras are presented. Moreover, these maximally symmetric equations are Darboux integrable. An enumeration of several submaximally symmetric (eight and seven-dimensional generic hyperbolic structures is also given.

Pion production and the nuclear equation of state

International Nuclear Information System (INIS)

Harris, J.W.; Odyniec, G.; Pugh, H.G.

1984-10-01

There has been considerable recent interest in the nuclear equation of state and how it may be determined in relativistic nucleus-nucleus collisions. In these collisions extremely high temperatures are reached and compression to densities several times that of normal nuclear matter are predicted. This affords us the unique opportunity to study, in a somewhat controlled manner, the behavior of nuclear matter under these extreme conditions. If the observables that are measured in experiments can be related in a quantitative way to state variables of the system then the equation of state can be extracted. This relation plays a very important role in understanding the formation and collapse of supernovae and the stability and structure of neutron stars. Furthermore, it can be used to test and constrain field theoretical approaches to nuclear matter and to help to better understand the dynamics of high energy nucleus-nucleus collisions. In this presentation the relationship between the nuclear equation of state and relativistic nucleus-nucleus collisions will be discussed with an emphasis on how to extract the former. That a high density state of the collision should exist will be shown. One observable, namely the pion multiplicity, will be shown to survive the succeeding stages of the collision process to provide information on the equation of state at high densities. The resulting equation of state will be presented and discussed in the light of recent theoretical development. 34 refs., 12 figs

Systematic tools for chemical equation balancing

International Nuclear Information System (INIS)

Filby, E.E.; Idaho National Engineering Lab., Idaho Falls, ID; Idaho Univ., Idaho Falls, ID

1989-01-01

One of the most important skills that chemists and chemical engineers must develop is the ability to balance chemical equations. The normal first method taught is ''balancing by inspection'', which is sometimes explained as simply ''mental algebra.'' Every textbook surveyed for this paper presents equation balancing first as a matter of trial and error; this includes four very recently published books. Very little further guidance is provided until oxidation-reduction reactions must be balanced. The most commonly taught approaches for balancing, redox equations have been the oxidation state change and ion-electron methods. Unfortunately, redox reactions are usually treated as a new topic, and what the student has teamed about ''ordinary'' equations is of little or no help. All too often, these contradictions simply confuse and frustrate students, and equation balancing is relegated to the status of a black art. This is ironic because such,confusion and frustration is not necessary: Chemical equations can, in fact, be balanced by explicitly definable mathematical methods. The purpose of this paper is to outline the algebraic methods involved

A New RF System for the CEBAF Normal Conducting Cavities

International Nuclear Information System (INIS)

Curt Hovater; Hai Dong; Alicia Hofler; George Lahti; John Musson; Tomasz Plawski

2004-01-01

The CEBAF Accelerator at Jefferson Lab is a 6 GeV five pass electron accelerator consisting of two superconducting linacs joined by independent magnetic transport arcs. CEBAF also has numerous normal conducting cavities for beam conditioning in the injector and for RF extraction to the experimental halls. The RF systems that presently control these cavities are becoming expensive to maintain, therefore a replacement RF control system is now being developed. For the new RF system, cavity field control is maintained digitally using an FPGA which contains the feedback algorithm. The system incorporates digital down conversion, using quadrature under-sampling at an IF frequency of 70 MHz. The VXI bus-crate was chosen as the operating platform because of its excellent RFI/EMI properties and its compatibility with the EPICS control system. The normal conducting cavities operate at both the 1497 MHz accelerating frequency and the sub-harmonic frequency of 499 MHz. To accommodate this, the ne w design will use different receiver-transmitter daughter cards for each frequency. This paper discusses the development of the new RF system and reports on initial results

Solutions for confluent and double-confluent Heun equations and some applications

Energy Technology Data Exchange (ETDEWEB)

El-Jaick, Lea Jaccoud; Figueiredo, Bartolomeu D.B. [Centro Brasileiro de Pesquisas Fisicas (CBPF), Rio de Janeiro, RJ (Brazil)

2008-07-01

This paper examines some solutions for confluent and double-confluent Heun equations and their applications to the Schroedinger equation with quasi-exactly solvable potentials. In the first place, we review two Leaver's solutions in series of regular and irregular confluent hypergeometric functions for the confluent equation [E. W. Leaver, J. Math. Phys. 27, 1238 (1986)] and introduce an additional expansion in series of irregular confluent hypergeometric functions. Then, we find the conditions under which one of these solutions can be written as a linear combination of the others. In the second place, by means of limiting procedures we generate solutions for the double-confluent equation as well as for special limits of both the confluent and double-confluent equations. In the third place, solutions of the Heun equations are used to solve the one-dimensional Schroedinger equation for quasi-exactly solvable potentials. We consider a symmetric and an asymmetric double-Morse potentials which appear in the theory of quantum spin systems [O. B. Zaslavskii and V. V. Ulyanov, Sov. Phys. JETP 60, 991 (1984)], a bottomless volcano-type potential which gives degenerate eigenstates [S. Kar and R. R. Parwani, Europhys. Lett., 80, 30004 (2007)], and a potential which leads to quasi normal modes, that is, to solutions presenting complex energies [H. T. Cho and C. L. Ho, J. Phys. A: Math. Theor. 40, 1325 (2007)]. (author)

Solutions for confluent and double-confluent Heun equations and some applications

International Nuclear Information System (INIS)

El-Jaick, Lea Jaccoud; Figueiredo, Bartolomeu D.B.

2008-01-01

This paper examines some solutions for confluent and double-confluent Heun equations and their applications to the Schroedinger equation with quasi-exactly solvable potentials. In the first place, we review two Leaver's solutions in series of regular and irregular confluent hypergeometric functions for the confluent equation [E. W. Leaver, J. Math. Phys. 27, 1238 (1986)] and introduce an additional expansion in series of irregular confluent hypergeometric functions. Then, we find the conditions under which one of these solutions can be written as a linear combination of the others. In the second place, by means of limiting procedures we generate solutions for the double-confluent equation as well as for special limits of both the confluent and double-confluent equations. In the third place, solutions of the Heun equations are used to solve the one-dimensional Schroedinger equation for quasi-exactly solvable potentials. We consider a symmetric and an asymmetric double-Morse potentials which appear in the theory of quantum spin systems [O. B. Zaslavskii and V. V. Ulyanov, Sov. Phys. JETP 60, 991 (1984)], a bottomless volcano-type potential which gives degenerate eigenstates [S. Kar and R. R. Parwani, Europhys. Lett., 80, 30004 (2007)], and a potential which leads to quasi normal modes, that is, to solutions presenting complex energies [H. T. Cho and C. L. Ho, J. Phys. A: Math. Theor. 40, 1325 (2007)]. (author)

Mei Symmetry and New Conserved Quantities of Tzénoff Equations for the Variable Mass Higher-Order Nonholonomic System

Institute of Scientific and Technical Information of China (English)

ZHENG Shi-Wang; WANG Jian-Bo; CHEN Xiang-Wei; XIE Jia-Fang

2012-01-01

Operational systems of spacecraft are general variable mass mechanics systems,and their symmetries and conserved quantities imply profound physical rules of the space system.We study the Mei symmetry of Tzénoff equations for a variable mass nonholonomic system and the new conserved quantities derived.The function expression of the new conserved quantities and the criterion equation which deduces these conserved quantities are presented.This result has some theoretical values in further research of conservation laws obeyed by the variable mass system.%Operational systems of spacecraft are general variable mass mechanics systems, and their symmetries and conserved quantities imply profound physical rules of the space system. We study the Mei symmetry of Tzenoff equations for a variable mass nonholonomic system and the new conserved quantities derived. The function expression of the new conserved quantities and the criterion equation which deduces these conserved quantities are presented. This result has some theoretical values in further research of conservation laws obeyed by the variable mass system.

Kinematic equations for resolved-rate control of an industrial robot arm

Science.gov (United States)

Barker, L. K.

1983-01-01

An operator can use kinematic, resolved-rate equations to dynamically control a robot arm by watching its response to commanded inputs. Known resolved-rate equations for the control of a particular six-degree-of-freedom industrial robot arm and proceeds to simplify the equations for faster computations are derived. Methods for controlling the robot arm in regions which normally cause mathematical singularities in the resolved-rate equations are discussed.

Regular growth of systems of functions and systems of non-homogeneous convolution equations in convex domains of the complex plane

International Nuclear Information System (INIS)

Krivosheev, A S

2000-01-01

In this paper we introduce the notion of regular growth for a system of entire functions of finite order and type. This is a direct and natural generalization of the classical completely regular growth of an entire function. We obtain sufficient and necessary conditions for the solubility of a system of non-homogeneous convolution equations in convex domains of the complex plane. These conditions depend on whether the system of Laplace transforms of the analytic functionals that generate the convolution equations has regular growth. In the case of smooth convex domains, these solubility conditions form a criterion

On matrix fractional differential equations

Directory of Open Access Journals (Sweden)

Adem Kılıçman

2017-01-01

Full Text Available The aim of this article is to study the matrix fractional differential equations and to find the exact solution for system of matrix fractional differential equations in terms of Riemann–Liouville using Laplace transform method and convolution product to the Riemann–Liouville fractional of matrices. Also, we show the theorem of non-homogeneous matrix fractional partial differential equation with some illustrative examples to demonstrate the effectiveness of the new methodology. The main objective of this article is to discuss the Laplace transform method based on operational matrices of fractional derivatives for solving several kinds of linear fractional differential equations. Moreover, we present the operational matrices of fractional derivatives with Laplace transform in many applications of various engineering systems as control system. We present the analytical technique for solving fractional-order, multi-term fractional differential equation. In other words, we propose an efficient algorithm for solving fractional matrix equation.

On the analytical solutions of the system of conformable time-fractional Robertson equations with 1-D diffusion

International Nuclear Information System (INIS)

Iyiola, O.S.; Tasbozan, O.; Kurt, A.; Çenesiz, Y.

2017-01-01

In this paper, we consider the system of conformable time-fractional Robertson equations with one-dimensional diffusion having widely varying diffusion coefficients. Due to the mismatched nature of the initial and boundary conditions associated with Robertson equation, there are spurious oscillations appearing in many computational algorithms. Our goal is to obtain an approximate solutions of this system of equations using the q-homotopy analysis method (q-HAM) and examine the widely varying diffusion coefficients and the fractional order of the derivative.

Cracking chaos-based encryption systems ruled by nonlinear time delay differential equations

International Nuclear Information System (INIS)

Udaltsov, Vladimir S.; Goedgebuer, Jean-Pierre; Larger, Laurent; Cuenot, Jean-Baptiste; Levy, Pascal; Rhodes, William T.

2003-01-01

We report that signal encoding with high-dimensional chaos produced by delayed feedback systems with a strong nonlinearity can be broken. We describe the procedure and illustrate the method with chaotic waveforms obtained from a strongly nonlinear optical system that we used previously to demonstrate signal encryption/decryption with chaos in wavelength. The method can be extended to any systems ruled by nonlinear time-delayed differential equations

Reconstruction of normal forms by learning informed observation geometries from data.

Science.gov (United States)

Yair, Or; Talmon, Ronen; Coifman, Ronald R; Kevrekidis, Ioannis G

2017-09-19

The discovery of physical laws consistent with empirical observations is at the heart of (applied) science and engineering. These laws typically take the form of nonlinear differential equations depending on parameters; dynamical systems theory provides, through the appropriate normal forms, an "intrinsic" prototypical characterization of the types of dynamical regimes accessible to a given model. Using an implementation of data-informed geometry learning, we directly reconstruct the relevant "normal forms": a quantitative mapping from empirical observations to prototypical realizations of the underlying dynamics. Interestingly, the state variables and the parameters of these realizations are inferred from the empirical observations; without prior knowledge or understanding, they parametrize the dynamics intrinsically without explicit reference to fundamental physical quantities.

Approximate solution to the Kolmogorov equation for a fission chain-reacting system

International Nuclear Information System (INIS)

Ruby, L.; McSwine, T.L.

1986-01-01

An approximate solution has been obtained for the Kolmogorov equation describing a fission chain-reacting system. The method considers the population of neutrons, delayed-neutron precursors, and detector counts. The effect of the detector is separated from the statistics of the chain reaction by a weak coupling assumption that predicts that the detector responds to the average rather than to the instantaneous neutron population. An approximate solution to the remaining equation, involving the populations of neutrons and precursors, predicts a negative-binomial behaviour for the neutron probability distribution

Ray equations of a weak shock in a hyperbolic system of ...

Indian Academy of Sciences (India)

differential form of this system of conservation laws is a hyperbolic system of partial differential equations. A(u)ut + B(α)(u)uxα = 0,. (1.3) where. A(u) = 〈∇u,H〉 and B(α)(u) = 〈∇u, F(α)〉,. (1.4) and we use the summation convention that a repeated symbol in subscripts and super- scripts in a term will mean summation over the ...

Tunneling junction as an open system. Normal tunneling

International Nuclear Information System (INIS)

Ono, Y.

1978-01-01

The method of the tunneling Hamiltonian is reformulated in the case of normal tunneling by introducing two independent particle baths. Due to the baths, it becomes possible to realize a final stationary state where the electron numbers of the two electrodes in the tunneling system are maintained constant and where there exists a stationary current. The effect of the bath-system couplings on the current-voltage characteristics of the junction is discussed in relation to the usual expression of the current as a function of voltage. (Auth.)

On the use of the GRACE normal equation of inter-satellite tracking data for estimation of soil moisture and groundwater in Australia

Directory of Open Access Journals (Sweden)

N. Tangdamrongsub

2018-03-01

Full Text Available An accurate estimation of soil moisture and groundwater is essential for monitoring the availability of water supply in domestic and agricultural sectors. In order to improve the water storage estimates, previous studies assimilated terrestrial water storage variation (ΔTWS derived from the Gravity Recovery and Climate Experiment (GRACE into land surface models (LSMs. However, the GRACE-derived ΔTWS was generally computed from the high-level products (e.g. time-variable gravity fields, i.e. level 2, and land grid from the level 3 product. The gridded data products are subjected to several drawbacks such as signal attenuation and/or distortion caused by a posteriori filters and a lack of error covariance information. The post-processing of GRACE data might lead to the undesired alteration of the signal and its statistical property. This study uses the GRACE least-squares normal equation data to exploit the GRACE information rigorously and negate these limitations. Our approach combines GRACE's least-squares normal equation (obtained from ITSG-Grace2016 product with the results from the Community Atmosphere Biosphere Land Exchange (CABLE model to improve soil moisture and groundwater estimates. This study demonstrates, for the first time, an importance of using the GRACE raw data. The GRACE-combined (GC approach is developed for optimal least-squares combination and the approach is applied to estimate the soil moisture and groundwater over 10 Australian river basins. The results are validated against the satellite soil moisture observation and the in situ groundwater data. Comparing to CABLE, we demonstrate the GC approach delivers evident improvement of water storage estimates, consistently from all basins, yielding better agreement on seasonal and inter-annual timescales. Significant improvement is found in groundwater storage while marginal improvement is observed in surface soil moisture estimates.

On the use of the GRACE normal equation of inter-satellite tracking data for estimation of soil moisture and groundwater in Australia

Science.gov (United States)

Tangdamrongsub, Natthachet; Han, Shin-Chan; Decker, Mark; Yeo, In-Young; Kim, Hyungjun

2018-03-01

An accurate estimation of soil moisture and groundwater is essential for monitoring the availability of water supply in domestic and agricultural sectors. In order to improve the water storage estimates, previous studies assimilated terrestrial water storage variation (ΔTWS) derived from the Gravity Recovery and Climate Experiment (GRACE) into land surface models (LSMs). However, the GRACE-derived ΔTWS was generally computed from the high-level products (e.g. time-variable gravity fields, i.e. level 2, and land grid from the level 3 product). The gridded data products are subjected to several drawbacks such as signal attenuation and/or distortion caused by a posteriori filters and a lack of error covariance information. The post-processing of GRACE data might lead to the undesired alteration of the signal and its statistical property. This study uses the GRACE least-squares normal equation data to exploit the GRACE information rigorously and negate these limitations. Our approach combines GRACE's least-squares normal equation (obtained from ITSG-Grace2016 product) with the results from the Community Atmosphere Biosphere Land Exchange (CABLE) model to improve soil moisture and groundwater estimates. This study demonstrates, for the first time, an importance of using the GRACE raw data. The GRACE-combined (GC) approach is developed for optimal least-squares combination and the approach is applied to estimate the soil moisture and groundwater over 10 Australian river basins. The results are validated against the satellite soil moisture observation and the in situ groundwater data. Comparing to CABLE, we demonstrate the GC approach delivers evident improvement of water storage estimates, consistently from all basins, yielding better agreement on seasonal and inter-annual timescales. Significant improvement is found in groundwater storage while marginal improvement is observed in surface soil moisture estimates.

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

The Use of a Code-generating System for the Derivation of the Equations for Wind Turbine Dynamics

Science.gov (United States)

Ganander, Hans

2003-10-01

For many reasons the size of wind turbines on the rapidly growing wind energy market is increasing. Relations between aeroelastic properties of these new large turbines change. Modifications of turbine designs and control concepts are also influenced by growing size. All these trends require development of computer codes for design and certification. Moreover, there is a strong desire for design optimization procedures, which require fast codes. General codes, e.g. finite element codes, normally allow such modifications and improvements of existing wind turbine models. This is done relatively easy. However, the calculation times of such codes are unfavourably long, certainly for optimization use. The use of an automatic code generating system is an alternative for relevance of the two key issues, the code and the design optimization. This technique can be used for rapid generation of codes of particular wind turbine simulation models. These ideas have been followed in the development of new versions of the wind turbine simulation code VIDYN. The equations of the simulation model were derived according to the Lagrange equation and using Mathematica®, which was directed to output the results in Fortran code format. In this way the simulation code is automatically adapted to an actual turbine model, in terms of subroutines containing the equations of motion, definitions of parameters and degrees of freedom. Since the start in 1997, these methods, constituting a systematic way of working, have been used to develop specific efficient calculation codes. The experience with this technique has been very encouraging, inspiring the continued development of new versions of the simulation code as the need has arisen, and the interest for design optimization is growing.

Existence and uniqueness of solution for a system of equations of ...

African Journals Online (AJOL)

The existence and uniqueness of solution for a system of equations of microwave heating of biologic issue is discussed. Using the Green function approach we establish the existence and uniqueness of solution. Journal of the Nigerian Association of Mathematical Physics Vol. 8 2004: pp. 177-180 ...

UNIFIED MODELS OF ELEMENTS OF POWER SUPPLY SYSTEMS BASED ON EQUATIONS IN PHASE COORDINATES

Directory of Open Access Journals (Sweden)

Yu.N. Vepryk

2015-12-01

Full Text Available Purpose. The models of electrical machines in the phase coordinates, the universal algorithm for the simulation of separate elements in a d-q coordinates system and in a phase-coordinates system are proposed. Methodology. Computer methods of investigation of transients in electrical systems are based on a compilation of systems of differential equations and their numerical integration solution methods. To solve differential equations an implicit method of numerical integration was chosen. Because it provides to complete structural simulation possibility: firstly developing models of separate elements and then forming a model of the complex system. For the mathematical simulation of electromagnetic transients in the elements of the electrical systems has been accepted the implicit Euler-Cauchy method, because it provides a higher precision and stability of the computing processes. Results. In developing the model elements identified two groups of elements: - Static elements and electrical machines in the d-q coordinates; - Rotating electrical machines in phase coordinates. As an example, the paper provides a model of synchronous and asynchronous electric machines in the d-q coordinates system and the phase coordinate system. The generalization algorithm and the unified notation form of equations of elements of an electrical system are obtained. It provides the possibility of using structural methods to develop a mathematical model of power systems under transient conditions. Practical value. In addition, the using of a computer model allows to implement multivariant calculations for research and study of factors affecting the quantitative characteristics of the transients.

Dynamic pathways to mediate reactions buried in thermal fluctuations. I. Time-dependent normal form theory for multidimensional Langevin equation.

Science.gov (United States)

Kawai, Shinnosuke; Komatsuzaki, Tamiki

2009-12-14

We present a novel theory which enables us to explore the mechanism of reaction selectivity and robust functions in complex systems persisting under thermal fluctuation. The theory constructs a nonlinear coordinate transformation so that the equation of motion for the new reaction coordinate is independent of the other nonreactive coordinates in the presence of thermal fluctuation. In this article we suppose that reacting systems subject to thermal noise are described by a multidimensional Langevin equation without a priori assumption for the form of potential. The reaction coordinate is composed not only of all the coordinates and velocities associated with the system (solute) but also of the random force exerted by the environment (solvent) with friction constants. The sign of the reaction coordinate at any instantaneous moment in the region of a saddle determines the fate of the reaction, i.e., whether the reaction will proceed through to the products or go back to the reactants. By assuming the statistical properties of the random force, one can know a priori a well-defined boundary of the reaction which separates the full position-velocity space in the saddle region into mainly reactive and mainly nonreactive regions even under thermal fluctuation. The analytical expression of the reaction coordinate provides the firm foundation on the mechanism of how and why reaction proceeds in thermal fluctuating environments.

Cylindrical shell under impact load including transverse shear and normal stress

International Nuclear Information System (INIS)

Shakeri, M.; Eslami, M.R.; Ghassaa, M.; Ohadi, A.R.

1993-01-01

The general governing equations of shell of revolution under shock loads are reduced to equations describing the elastic behavior of cylindrical shell under axisymmetric impact load. The effect of lateral normal stress, transverse shear, and rotary inertia are included, and the equations are solved by Galerkin finite element method. The results are compared with the previous works of authors. (author)

Conservation laws for certain time fractional nonlinear systems of partial differential equations

Science.gov (United States)

Singla, Komal; Gupta, R. K.

2017-12-01

In this study, an extension of the concept of nonlinear self-adjointness and Noether operators is proposed for calculating conserved vectors of the time fractional nonlinear systems of partial differential equations. In our recent work (J Math Phys 2016; 57: 101504), by proposing the symmetry approach for time fractional systems, the Lie symmetries for some fractional nonlinear systems have been derived. In this paper, the obtained infinitesimal generators are used to find conservation laws for the corresponding fractional systems.

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.

Prediction Equations for Spirometry for Children from Northern India.

Science.gov (United States)

Chhabra, Sunil K; Kumar, Rajeev; Mittal, Vikas

2016-09-08

To develop prediction equations for spirometry for children from northern India using current international guidelines for standardization. Re-analysis of cross-sectional data from a single school. 670 normal children (age 6-17 y; 365 boys) of northern Indian parentage. After screening for normal health, we carried out spirometry with recommended quality assurance according to current guidelines. We developed linear and nonlinear prediction equations using multiple regression analysis. We selected the final models on the basis of the highest coefficient of multiple determination (R2) and statistical validity. Spirometry parameters: FVC, FEV1, PEFR, FEF50, FEF75 and FEF25-75. The equations for the main parameters were as follows: Boys, Ln FVC = -1.687+0.016*height +0.022*age; Ln FEV1 = -1.748+0.015*height+0.031*age. Girls, Ln FVC = -9.989 +(2.018*Ln(height)) + (0.324*Ln(age)); Ln FEV1 = -10.055 +(1.990*Ln(height))+(0.358*Ln(age)). Nonlinear regression yielded substantially greater R2 values compared to linear models except for FEF50 for girls. Height and age were found to be the significant explanatory variables for all parameters on multiple regression with weight making no significant contribution. We developed prediction equations for spirometry for children from northern India. Nonlinear equations were superior to linear equations.

Langevin equation in systems with also negative temperatures

Science.gov (United States)

Baldovin, Marco; Puglisi, Andrea; Vulpiani, Angelo

2018-04-01

We discuss how to derive a Langevin equation (LE) in non standard systems, i.e. when the kinetic part of the Hamiltonian is not the usual quadratic function. This generalization allows to consider also cases with negative absolute temperature. We first give some phenomenological arguments suggesting the shape of the viscous drift, replacing the usual linear viscous damping, and its relation with the diffusion coefficient modulating the white noise term. As a second step, we implement a procedure to reconstruct the drift and the diffusion term of the LE from the time-series of the momentum of a heavy particle embedded in a large Hamiltonian system. The results of our reconstruction are in good agreement with the phenomenological arguments. Applying the method to systems with negative temperature, we can observe that also in this case there is a suitable LE, obtained with a precise protocol, able to reproduce in a proper way the statistical features of the slow variables. In other words, even in this context, systems with negative temperature do not show any pathology.

A Two-Species Cooperative Lotka-Volterra System of Degenerate Parabolic Equations

Directory of Open Access Journals (Sweden)

Jiebao Sun

2011-01-01

parabolic equations. We are interested in the coexistence of the species in a bounded domain. We establish the existence of global generalized solutions of the initial boundary value problem by means of parabolic regularization and also consider the existence of the nontrivial time-periodic solution for this system.

TOEPLITZ, Solution of Linear Equation System with Toeplitz or Circulant Matrix

International Nuclear Information System (INIS)

Garbow, B.

1984-01-01

Description of program or function: TOEPLITZ is a collection of FORTRAN subroutines for solving linear systems Ax=b, where A is a Toeplitz matrix, a Circulant matrix, or has one or several block structures based on Toeplitz or Circulant matrices. Such systems arise in problems of electrodynamics, acoustics, mathematical statistics, algebra, in the numerical solution of integral equations with a difference kernel, and in the theory of stationary time series and signals

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

International Nuclear Information System (INIS)

Pabst, M.

2014-01-01

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

Correction of Bowtie-Filter Normalization and Crescent Artifacts for a Clinical CBCT System.

Science.gov (United States)

Zhang, Hong; Kong, Vic; Huang, Ke; Jin, Jian-Yue

2017-02-01

To present our experiences in understanding and minimizing bowtie-filter crescent artifacts and bowtie-filter normalization artifacts in a clinical cone beam computed tomography system. Bowtie-filter position and profile variations during gantry rotation were studied. Two previously proposed strategies (A and B) were applied to the clinical cone beam computed tomography system to correct bowtie-filter crescent artifacts. Physical calibration and analytical approaches were used to minimize the norm phantom misalignment and to correct for bowtie-filter normalization artifacts. A combined procedure to reduce bowtie-filter crescent artifacts and bowtie-filter normalization artifacts was proposed and tested on a norm phantom, CatPhan, and a patient and evaluated using standard deviation of Hounsfield unit along a sampling line. The bowtie-filter exhibited not only a translational shift but also an amplitude variation in its projection profile during gantry rotation. Strategy B was better than strategy A slightly in minimizing bowtie-filter crescent artifacts, possibly because it corrected the amplitude variation, suggesting that the amplitude variation plays a role in bowtie-filter crescent artifacts. The physical calibration largely reduced the misalignment-induced bowtie-filter normalization artifacts, and the analytical approach further reduced bowtie-filter normalization artifacts. The combined procedure minimized both bowtie-filter crescent artifacts and bowtie-filter normalization artifacts, with Hounsfield unit standard deviation being 63.2, 45.0, 35.0, and 18.8 Hounsfield unit for the best correction approaches of none, bowtie-filter crescent artifacts, bowtie-filter normalization artifacts, and bowtie-filter normalization artifacts + bowtie-filter crescent artifacts, respectively. The combined procedure also demonstrated reduction of bowtie-filter crescent artifacts and bowtie-filter normalization artifacts in a CatPhan and a patient. We have developed a step

Equationally Noetherian property of Ershov algebras

OpenAIRE

Dvorzhetskiy, Yuriy

2014-01-01

This article is about equationally Noetherian and weak equationally Noetherian property of Ershov algebras. Here we show two canonical forms of the system of equations over Ershov algebras and two criteria of equationally Noetherian and weak equationally Noetherian properties.

Exploring the Phase Space of a System of Differential Equations: Different Mathematical Registers

Science.gov (United States)

Dana-Picard, Thierry; Kidron, Ivy

2008-01-01

We describe and analyze a situation involving symbolic representation and graphical visualization of the solution of a system of two linear differential equations, using a computer algebra system. Symbolic solution and graphical representation complement each other. Graphical representation helps to understand the behavior of the symbolic…

A New Distribution-Random Limit Normal Distribution

OpenAIRE

Gong, Xiaolin; Yang, Shuzhen

2013-01-01

This paper introduces a new distribution to improve tail risk modeling. Based on the classical normal distribution, we define a new distribution by a series of heat equations. Then, we use market data to verify our model.

On the projective normality of Artin-Schreier curves

Directory of Open Access Journals (Sweden)

Alberto Ravagnani

2013-11-01

Full Text Available In this paper we study the projective normality of certain Artin-Schreier curves Y_f defined over a field F of characteristic p by the equations y^q+y=f(x, q being a power of p and f in F[x] being a polynomial in x of degree m, with (m,p=1. Many Y_f curves are singular and so, to be precise, here we study the projective normality of appropriate projective models of their normalization.

An Explicit Formulation of Singularity-Free Dynamic Equations of Mechanical Systems in Lagrangian Form---Part one: Single Rigid Bodies

Directory of Open Access Journals (Sweden)

Pål Johan From

2012-04-01

Full Text Available This paper presents the explicit dynamic equations of a mechanical system. The equations are presented so that they can easily be implemented in a simulation software or controller environment and are also well suited for system and controller analysis. The dynamics of a general mechanical system consisting of one or more rigid bodies can be derived from the Lagrangian. We can then use several well known properties of Lie groups to guarantee that these equations are well defined. This will, however, often lead to rather abstract formulation of the dynamic equations that cannot be implemented in a simulation software directly. In this paper we close this gap and show what the explicit dynamic equations look like. These equations can then be implemented directly in a simulation software and no background knowledge on Lie theory and differential geometry on the practitioner's side is required. This is the first of two papers on this topic. In this paper we derive the dynamics for single rigid bodies, while in the second part we study multibody systems. In addition to making the equations more accessible to practitioners, a motivation behind the papers is to correct a few errors commonly found in literature. For the first time, we show the detailed derivations and how to arrive at the correct set of equations. We also show through some simple examples that these correspond with the classical formulations found from Lagrange's equations. The dynamics is derived from the Boltzmann--Hamel equations of motion in terms of local position and velocity variables and the mapping to the corresponding quasi-velocities. Finally we present a new theorem which states that the Boltzmann--Hamel formulation of the dynamics is valid for all transformations with a Lie group topology. This has previously only been indicated through examples, but here we also present the formal proof. The main motivation of these papers is to allow practitioners not familiar with

Trajectory attractors of equations of mathematical physics

International Nuclear Information System (INIS)

Vishik, Marko I; Chepyzhov, Vladimir V

2011-01-01

In this survey the method of trajectory dynamical systems and trajectory attractors is described, and is applied in the study of the limiting asymptotic behaviour of solutions of non-linear evolution equations. This method is especially useful in the study of dissipative equations of mathematical physics for which the corresponding Cauchy initial-value problem has a global (weak) solution with respect to the time but the uniqueness of this solution either has not been established or does not hold. An important example of such an equation is the 3D Navier-Stokes system in a bounded domain. In such a situation one cannot use directly the classical scheme of construction of a dynamical system in the phase space of initial conditions of the Cauchy problem of a given equation and find a global attractor of this dynamical system. Nevertheless, for such equations it is possible to construct a trajectory dynamical system and investigate a trajectory attractor of the corresponding translation semigroup. This universal method is applied for various types of equations arising in mathematical physics: for general dissipative reaction-diffusion systems, for the 3D Navier-Stokes system, for dissipative wave equations, for non-linear elliptic equations in cylindrical domains, and for other equations and systems. Special attention is given to using the method of trajectory attractors in approximation and perturbation problems arising in complicated models of mathematical physics. Bibliography: 96 titles.

Exponential model normalization for electrical capacitance tomography with external electrodes under gap permittivity conditions

International Nuclear Information System (INIS)

Baidillah, Marlin R; Takei, Masahiro

2017-01-01

A nonlinear normalization model which is called exponential model for electrical capacitance tomography (ECT) with external electrodes under gap permittivity conditions has been developed. The exponential model normalization is proposed based on the inherently nonlinear relationship characteristic between the mixture permittivity and the measured capacitance due to the gap permittivity of inner wall. The parameters of exponential equation are derived by using an exponential fitting curve based on the simulation and a scaling function is added to adjust the experiment system condition. The exponential model normalization was applied to two dimensional low and high contrast dielectric distribution phantoms by using simulation and experimental studies. The proposed normalization model has been compared with other normalization models i.e. Parallel, Series, Maxwell and Böttcher models. Based on the comparison of image reconstruction results, the exponential model is reliable to predict the nonlinear normalization of measured capacitance in term of low and high contrast dielectric distribution. (paper)

Development and adjustment of programs for solving systems of linear equations

International Nuclear Information System (INIS)

Fujimura, Toichiro

1978-03-01

Programs for solving the systems of linear equations have been adjusted and developed in expanding the scientific subroutine library SSL. The principal programs adjusted are based on the congruent method, method of product form of the inverse, orthogonal method, Crout's method for sparse system, and acceleration of iterative methods. The programs developed are based on the escalator method, direct parallel residue method and block tridiagonal method for band system. Described are usage of the programs developed and their future improvement. FORTRAN lists with simple examples in tests of the programs are also given. (auth.)

Automatic simplification of systems of reaction-diffusion equations by a posteriori analysis.

Science.gov (United States)

Maybank, Philip J; Whiteley, Jonathan P

2014-02-01

Many mathematical models in biology and physiology are represented by systems of nonlinear differential equations. In recent years these models have become increasingly complex in order to explain the enormous volume of data now available. A key role of modellers is to determine which components of the model have the greatest effect on a given observed behaviour. An approach for automatically fulfilling this role, based on a posteriori analysis, has recently been developed for nonlinear initial value ordinary differential equations [J.P. Whiteley, Model reduction using a posteriori analysis, Math. Biosci. 225 (2010) 44-52]. In this paper we extend this model reduction technique for application to both steady-state and time-dependent nonlinear reaction-diffusion systems. Exemplar problems drawn from biology are used to demonstrate the applicability of the technique. Copyright © 2014 Elsevier Inc. All rights reserved.

Solvability of a class of systems of infinite-dimensional integral equations and their application in statistical mechanics

International Nuclear Information System (INIS)

Gonchar, N.S.

1986-01-01

This paper presents a mathematical method developed for investigating a class of systems of infinite-dimensional integral equations which have application in statistical mechanics. Necessary and sufficient conditions are obtained for the uniqueness and bifurcation of the solution of this class of systems of equations. Problems of equilibrium statistical mechanics are considered on the basis of this method

Perturbation theory for continuous stochastic equations

International Nuclear Information System (INIS)

Chechetkin, V.R.; Lutovinov, V.S.

1987-01-01

The various general perturbational schemes for continuous stochastic equations are considered. These schemes have many analogous features with the iterational solution of Schwinger equation for S-matrix. The following problems are discussed: continuous stochastic evolution equations for probability distribution functionals, evolution equations for equal time correlators, perturbation theory for Gaussian and Poissonian additive noise, perturbation theory for birth and death processes, stochastic properties of systems with multiplicative noise. The general results are illustrated by diffusion-controlled reactions, fluctuations in closed systems with chemical processes, propagation of waves in random media in parabolic equation approximation, and non-equilibrium phase transitions in systems with Poissonian breeding centers. The rate of irreversible reaction X + X → A (Smoluchowski process) is calculated with the use of general theory based on continuous stochastic equations for birth and death processes. The threshold criterion and range of fluctuational region for synergetic phase transition in system with Poissonian breeding centers are also considered. (author)

Performance estimates for personnel access control systems

International Nuclear Information System (INIS)

Bradley, R.G.

1980-10-01

Current performance estimates for personnel access control systems use estimates of Type I and Type II verification errors. A system performance equation which addresses normal operation, the insider, and outside adversary attack is developed. Examination of this equation reveals the inadequacy of classical Type I and II error evaluations which require detailed knowledge of the adversary threat scenario for each specific installation. Consequently, new performance measures which are consistent with the performance equation and independent of the threat are developed as an aid in selecting personnel access control systems

Normal mode splitting and ground state cooling in a Fabry—Perot optical cavity and transmission line resonator

International Nuclear Information System (INIS)

Chen Hua-Jun; Mi Xian-Wu

2011-01-01

Optomechanical dynamics in two systems which are a transmission line resonator and Fabrya—Perot optical cavity via radiation—pressure are investigated by linearized quantum Langevin equation. We work in the resolved sideband regime where the oscillator resonance frequency exceeds the cavity linewidth. Normal mode splittings of the mechanical resonator as a pure result of the coupling interaction in the two optomechanical systems is studied, and we make a comparison of normal mode splitting of mechanical resonator between the two systems. In the optical cavity, the normal mode splitting of the movable mirror approaches the latest experiment very well. In addition, an approximation scheme is introduced to demonstrate the ground state cooling, and we make a comparison of cooling between the two systems dominated by two key factors, which are the initial bath temperature and the mechanical quality factor. Since both the normal mode splitting and cooling require working in the resolved sideband regime, whether the normal mode splitting influences the cooling of the mirror is considered. Considering the size of the mechanical resonator and precooling the system, the mechanical resonator in the transmission line resonator system is easier to achieve the ground state cooling than in optical cavity. (electromagnetism, optics, acoustics, heat transfer, classical mechanics, and fluid dynamics)

FEQinput—An editor for the full equations (FEQ) hydraulic modeling system

Science.gov (United States)

Ancalle, David S.; Ancalle, Pablo J.; Domanski, Marian M.

2017-10-30

IntroductionThe Full Equations Model (FEQ) is a computer program that solves the full, dynamic equations of motion for one-dimensional unsteady hydraulic flow in open channels and through control structures. As a result, hydrologists have used FEQ to design and operate flood-control structures, delineate inundation maps, and analyze peak-flow impacts. To aid in fighting floods, hydrologists are using the software to develop a system that uses flood-plain models to simulate real-time streamflow.Input files for FEQ are composed of text files that contain large amounts of parameters, data, and instructions that are written in a format exclusive to FEQ. Although documentation exists that can aid in the creation and editing of these input files, new users face a steep learning curve in order to understand the specific format and language of the files.FEQinput provides a set of tools to help a new user overcome the steep learning curve associated with creating and modifying input files for the FEQ hydraulic model and the related utility tool, Full Equations Utilities (FEQUTL).

Elements of partial differential equations

CERN Document Server

Sneddon, Ian Naismith

1957-01-01

Geared toward students of applied rather than pure mathematics, this volume introduces elements of partial differential equations. Its focus is primarily upon finding solutions to particular equations rather than general theory.Topics include ordinary differential equations in more than two variables, partial differential equations of the first and second orders, Laplace's equation, the wave equation, and the diffusion equation. A helpful Appendix offers information on systems of surfaces, and solutions to the odd-numbered problems appear at the end of the book. Readers pursuing independent st

Utility rate equations of group population dynamics in biological and social systems.

Directory of Open Access Journals (Sweden)

Vyacheslav I Yukalov

Full Text Available We present a novel system of equations to describe the evolution of self-organized structured societies (biological or human composed of several trait groups. The suggested approach is based on the combination of ideas employed in the theory of biological populations, system theory, and utility theory. The evolution equations are defined as utility rate equations, whose parameters are characterized by the utility of each group with respect to the society as a whole and by the mutual utilities of groups with respect to each other. We analyze in detail the cases of two groups (cooperators and defectors and of three groups (cooperators, defectors, and regulators and find that, in a self-organized society, neither defectors nor regulators can overpass the maximal fractions of about [Formula: see text] each. This is in agreement with the data for bee and ant colonies. The classification of societies by their distance from equilibrium is proposed. We apply the formalism to rank the countries according to the introduced metric quantifying their relative stability, which depends on the cost of defectors and regulators as well as their respective population fractions. We find a remarkable concordance with more standard economic ranking based, for instance, on GDP per capita.

Utility Rate Equations of Group Population Dynamics in Biological and Social Systems

Science.gov (United States)

Yukalov, Vyacheslav I.; Yukalova, Elizaveta P.; Sornette, Didier

2013-01-01

We present a novel system of equations to describe the evolution of self-organized structured societies (biological or human) composed of several trait groups. The suggested approach is based on the combination of ideas employed in the theory of biological populations, system theory, and utility theory. The evolution equations are defined as utility rate equations, whose parameters are characterized by the utility of each group with respect to the society as a whole and by the mutual utilities of groups with respect to each other. We analyze in detail the cases of two groups (cooperators and defectors) and of three groups (cooperators, defectors, and regulators) and find that, in a self-organized society, neither defectors nor regulators can overpass the maximal fractions of about each. This is in agreement with the data for bee and ant colonies. The classification of societies by their distance from equilibrium is proposed. We apply the formalism to rank the countries according to the introduced metric quantifying their relative stability, which depends on the cost of defectors and regulators as well as their respective population fractions. We find a remarkable concordance with more standard economic ranking based, for instance, on GDP per capita. PMID:24386163

System of adjoint P1 equations for neutron moderation

International Nuclear Information System (INIS)

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

2000-01-01

In some applications of perturbation theory, it is necessary know the adjoint neutron flux, which is obtained by the solution of adjoint neutron diffusion equation. However, the multigroup constants used for this are weighted in only the direct neutron flux, from the solution of direct P1 equations. In this work, this procedure is questioned and the adjoint P1 equations are derived by the neutron transport equation, the reversion operators rules and analogies between direct and adjoint parameters. (author)

Recent symbolic summation methods to solve coupled systems of differential and difference equations

International Nuclear Information System (INIS)

Schneider, Carsten; Bluemlein, Johannes; Freitas, Abilio de

2014-07-01

We outline a new algorithm to solve coupled systems of differential equations in one continuous variable x (resp. coupled difference equations in one discrete variable N) depending on a small parameter ε: given such a system and given sufficiently many initial values, we can determine the first coefficients of the Laurent-series solutions in ε if they are expressible in terms of indefinite nested sums and products. This systematic approach is based on symbolic summation algorithms in the context of difference rings/fields and uncoupling algorithms. The proposed method gives rise to new interesting applications in connection with integration by parts (IBP) methods. As an illustrative example, we will demonstrate how one can calculate the ε-expansion of a ladder graph with 6 massive fermion lines.

Recent symbolic summation methods to solve coupled systems of differential and difference equations

Energy Technology Data Exchange (ETDEWEB)

Schneider, Carsten [Johannes Kepler Univ., Linz (Austria). Research Inst. for Symbolic Computation (RISC); Bluemlein, Johannes; Freitas, Abilio de [Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany)

2014-07-15

We outline a new algorithm to solve coupled systems of differential equations in one continuous variable x (resp. coupled difference equations in one discrete variable N) depending on a small parameter ε: given such a system and given sufficiently many initial values, we can determine the first coefficients of the Laurent-series solutions in ε if they are expressible in terms of indefinite nested sums and products. This systematic approach is based on symbolic summation algorithms in the context of difference rings/fields and uncoupling algorithms. The proposed method gives rise to new interesting applications in connection with integration by parts (IBP) methods. As an illustrative example, we will demonstrate how one can calculate the ε-expansion of a ladder graph with 6 massive fermion lines.

Kidney Length in Normal Korean Children

International Nuclear Information System (INIS)

Kim, In One; Cheon, Jung Eun; Lee, Young Seok; Lee, Sun Wha; Kim, Ok Hwa; Kim, Ji Hye; Kim, Hong Dae; Sim, Jung Suk

2010-01-01

Renal length offers important information to detect or follow-up various renal diseases. The purpose of this study was to determine the kidney length of normal Korean children in relation to age, height, weight, body surface area (BSA), and body mass index (BMI). Children between 1 month and 15 years of age without urological abnormality were recruited. Children below 3rd percentile and over 97th percentile for height or weight were excluded. Both renal lengths were measured in the prone position three times and then averaged by experienced radiologists. The mean length and standard deviation for each age group was obtained, and regression equation was calculated between renal length and age, weight, height, BSA, and BMI, respectively. Renal length was measured in 550 children. Renal length grows rapidly until 24 month, while the growth rate is reduced thereafter. The regression equation for age is: renal length (mm) = 45.953 + 1.064 x age (month, ≤ 24 months) (R2 = 0.720) or 62.173 + 0.203 x age (months, > 24 months) (R2 = 0.711). The regression equation for height is: renal length (mm) = 24.494 + 0.457 x height (cm) (R2 = 0.894). The regression equation for weight is: renal length (mm) = 38.342 + 2.117 x weight (kg, ≤18 kg) (R2 = 0.852) or 64.498 + 0.646 x weight (kg, > 18 kg) (R2 = 0.651). The regression equation for BSA is: renal length (mm) = 31.622 + 61.363 x BSA (m2, ≤ 0.7) (R2 = 0.857) or 52.717 + 29.959 x BSA (m2, > 0.7) (R2 = 0.715). The regression equation for BMI is: renal length (mm) = 44.474 + 1.163 x BMI (R2 = 0.079). This study provides data on the normal renal length and its association with age, weight, height, BSA and BMI. The results of this study will guide the detection and follow-up of renal diseases in Korean children

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

Science.gov (United States)

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

2007-12-27

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

A perturbed Lennard-Jones chain equation of state for liquid metals

Energy Technology Data Exchange (ETDEWEB)

Mousazadeh, M H; Marageh, M Ghanadi [AEOI, JIH Research Laboratory, 11365/8486, Tehran (Iran, Islamic Republic of)

2006-05-24

The perturbed Lennard-Jones chain (PLJC) equation of state is formulated based on first-order variational perturbation theory. The model uses two parameters for a monatomic system, segment size, {sigma}, and segment energy, {epsilon}/k. In this work, we employed the PLJC equation to calculate the liquid density of 26 metals, including alkali and alkali earth metals, iron, cobalt, nickel, copper, silver, gold, zinc, cadmium, mercury, aluminium, gallium, indium, thallium, tin, lead, antimony, and bismuth, for which accurate experimental data exist in the literature. The calculations cover a broad range of temperatures ranging from the melting point to close to the critical point and pressures ranging from the vapour-pressure curve up to pressures as high as 2000 bar. The average absolute deviation in the liquid density predicted by the PLJC equation of state in the saturation line compared with experimental data is 1.26%. Also, using the normal melting temperature and liquid density at melting point (T{sub m}, {rho}{sub m}) as input data for the estimation of the equation of state parameters provides a good correlation of liquid density at saturated and compressed pressures.

On existence of soliton solutions of arbitrary-order system of nonlinear Schrodinger equations

International Nuclear Information System (INIS)

Zhestkov, S.V.

2003-01-01

The soliton solutions are constructed for the system of arbitrary-order coupled nonlinear Schrodinger equations . The necessary and sufficient conditions of existence of these solutions are obtained. It is shown that the maximum number of solitons in nondegenerate case is 4L, where L is order of the system. (author)

Construction and accuracy of partial differential equation approximations to the chemical master equation.

Science.gov (United States)

Grima, Ramon

2011-11-01