Conservation form of the equations of fluid dynamics in general nonsteady coordinates

International Nuclear Information System (INIS)

Zhang, H.; Camarero, R.; Kahawita, R.

1985-01-01

Many of the differential equations arising in fluid dynamics may be stated in conservation-law form. A number of investigations have been conducted with the aim to derive the conservation-law form of the Navier-Stokes equations in general nonsteady coordinate systems. The present note has the objective to illustrate a mathematical methodology with which such forms of the equations may be derived in an easier and more general fashion. For numerical applications, the scalar form of the equations is eventually provided. Attention is given to the conservation form of equations in curvilinear coordinates and numerical considerations. 6 references

Nonlinear self-adjointness, conservation laws, and the construction of solutions of partial differential equations using conservation laws

International Nuclear Information System (INIS)

Ibragimov, N Kh; Avdonina, E D

2013-01-01

The method of nonlinear self-adjointness, which was recently developed by the first author, gives a generalization of Noether's theorem. This new method significantly extends approaches to constructing conservation laws associated with symmetries, since it does not require the existence of a Lagrangian. In particular, it can be applied to any linear equations and any nonlinear equations that possess at least one local conservation law. The present paper provides a brief survey of results on conservation laws which have been obtained by this method and published mostly in recent preprints of the authors, along with a method for constructing exact solutions of systems of partial differential equations with the use of conservation laws. In most cases the solutions obtained by the method of conservation laws cannot be found as invariant or partially invariant solutions. Bibliography: 23 titles

On conserved densities and asymptotic behaviour for the potential Kadomtsev-Petviashvili equation

International Nuclear Information System (INIS)

Rosenhaus, V

2006-01-01

We study local conservation laws with non-vanishing conserved densities and corresponding boundary conditions for the potential Kadomtsev-Petviashvili equation. We analyse an infinite symmetry group of the equation, and generate a finite number of conserved densities corresponding to infinite symmetries through appropriate boundary conditions

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)

Infinite sets of conservation laws for linear and nonlinear field equations

International Nuclear Information System (INIS)

Mickelsson, J.

1984-01-01

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

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.

Painleve analysis, conservation laws, and symmetry of perturbed nonlinear equations

International Nuclear Information System (INIS)

Basak, S.; Chowdhury, A.R.

1987-01-01

The authors consider the Lie-Backlund symmetries and conservation laws of a perturbed KdV equation and NLS equation. The arbitrary coefficients of the perturbing terms can be related to the condition of existence of nontrivial LB symmetry generators. When the perturbed KdV equation is subjected to Painleve analysis a la Weiss, it is found that the resonance position changes compared to the unperturbed one. They prove the compatibility of the overdetermined set of equations obtained at the different stages of recursion relations, at least for one branch. All other branches are also indicated and difficulties associated them are discussed considering the perturbation parameter epsilon to be small. They determine the Lax pair for the aforesaid branch through the use of Schwarzian derivative. For the perturbed NLS equation they determine the conservation laws following the approach of Chen and Liu. From the recurrence of these conservation laws a Lax pair is constructed. But the Painleve analysis does not produce a positive answer for the perturbed NLS equation. So here they have two contrasting examples of perturbed nonlinear equations: one passes the Painleve test and its Lax pair can be found from the analysis itself, but the other equation does not meet the criterion of the Painleve test, though its Lax pair is found in another way

A depth semi-averaged model for coastal dynamics

Science.gov (United States)

Antuono, M.; Colicchio, G.; Lugni, C.; Greco, M.; Brocchini, M.

2017-05-01

The present work extends the semi-integrated method proposed by Antuono and Brocchini ["Beyond Boussinesq-type equations: Semi-integrated models for coastal dynamics," Phys. Fluids 25(1), 016603 (2013)], which comprises a subset of depth-averaged equations (similar to Boussinesq-like models) and a Poisson equation that accounts for vertical dynamics. Here, the subset of depth-averaged equations has been reshaped in a conservative-like form and both the Poisson equation formulations proposed by Antuono and Brocchini ["Beyond Boussinesq-type equations: Semi-integrated models for coastal dynamics," Phys. Fluids 25(1), 016603 (2013)] are investigated: the former uses the vertical velocity component (formulation A) and the latter a specific depth semi-averaged variable, ϒ (formulation B). Our analyses reveal that formulation A is prone to instabilities as wave nonlinearity increases. On the contrary, formulation B allows an accurate, robust numerical implementation. Test cases derived from the scientific literature on Boussinesq-type models—i.e., solitary and Stokes wave analytical solutions for linear dispersion and nonlinear evolution and experimental data for shoaling properties—are used to assess the proposed solution strategy. It is found that the present method gives reliable predictions of wave propagation in shallow to intermediate waters, in terms of both semi-averaged variables and conservation properties.

General particle transport equation. Final report

International Nuclear Information System (INIS)

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

1994-12-01

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

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

International Nuclear Information System (INIS)

Niederle, J.

1984-01-01

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

Averaging problem in general relativity, macroscopic gravity and using Einstein's equations in cosmology.

Science.gov (United States)

Zalaletdinov, R. M.

1998-04-01

The averaging problem in general relativity is briefly discussed. A new setting of the problem as that of macroscopic description of gravitation is proposed. A covariant space-time averaging procedure is described. The structure of the geometry of macroscopic space-time, which follows from averaging Cartan's structure equations, is described and the correlation tensors present in the theory are discussed. The macroscopic field equations (averaged Einstein's equations) derived in the framework of the approach are presented and their structure is analysed. The correspondence principle for macroscopic gravity is formulated and a definition of the stress-energy tensor for the macroscopic gravitational field is proposed. It is shown that the physical meaning of using Einstein's equations with a hydrodynamic stress-energy tensor in looking for cosmological models means neglecting all gravitational field correlations. The system of macroscopic gravity equations to be solved when the correlations are taken into consideration is given and described.

Analysis of the Symmetries and Conservation Laws of the Nonlinear Jaulent-Miodek Equation

Directory of Open Access Journals (Sweden)

Mehdi Nadjafikhah

2014-01-01

Full Text Available Lie symmetry method is performed for the nonlinear Jaulent-Miodek equation. We will find the symmetry group and optimal systems of Lie subalgebras. The Lie invariants associated with the symmetry generators as well as the corresponding similarity reduced equations are also pointed out. And conservation laws of the J-M equation are presented with two steps: firstly, finding multipliers for computation of conservation laws and, secondly, symbolic computation of conservation laws will be applied.

Comment on ''Boltzmann equation and the conservation of particle number''

International Nuclear Information System (INIS)

Zanette, D.

1990-09-01

In a recent paper (Z. Banggu, Phys. Rev. A 42, 761 (1990)) it is argued that some solutions of the Boltzmann equation do not satisfy particle conservation as a consequence of the independence of velocity on position. In this comment, the arguments and conclusions of that paper are discussed. In particular, it is stressed that the temporal series used for solving the kinetic equation are generally divergent. A discussion about the particle conservation in its solutions is also provided. (author). 4 refs

Direct Construction of Conservation Laws from Field Equations

International Nuclear Information System (INIS)

Anco, S.C.; Bluman, G.

1997-01-01

This Letter presents an algorithm to obtain all local conservation laws for any system of field equations. The algorithm uses a formula which directly generates the conservation laws and does not depend on the system having a Lagrangian formulation, in contrast to Noether close-quote s theorem which requires a Lagrangian. Several examples are considered including dissipative systems inherently having no Lagrangian. copyright 1997 The American Physical Society

Nonlocal symmetries and nonlocal conservation laws of Maxwell's equations

International Nuclear Information System (INIS)

Anco, S.C.; Bluman, G.

1997-01-01

Nonlocal symmetries are obtained for Maxwell's equations in three space-time dimensions through the use of two potential systems involving scalar and vector potentials for the electromagnetic field. Corresponding nonlocal conservation laws are derived from these symmetries. The conservation laws yield nine functionally independent constants of motion which cannot be expressed in terms of the constants of motion arising from local conservation laws for space-time symmetries. These nine constants of motion represent additional conserved quantities for the electromagnetic field in three space endash time dimensions. copyright 1997 American Institute of Physics

Non-Noether conserved quantity for differential equations of motion in the phase space

Institute of Scientific and Technical Information of China (English)

无

2002-01-01

A non-Noether conserved quantity for the differential equations of motion of mechanical systems in the phase space is studied. The differential equations of motion of the systems are established and the determining equations of Lie symmetry are given. An existence theorem of non-Noether conserved quantity is obtained. An example is given to illustrate the application of the result.

Conservation laws and covariant equations of motion for spinning particles

OpenAIRE

Obukhov, Yuri N.; Puetzfeld, Dirk

2015-01-01

We derive the Noether identities and the conservation laws for general gravitational models with arbitrarily interacting matter and gravitational fields. These conservation laws are used for the construction of the covariant equations of motion for test bodies with minimal and nonminimal coupling.

Stationarity-conservation laws for fractional differential equations with variable coefficients

International Nuclear Information System (INIS)

Klimek, Malgorzata

2002-01-01

In this paper, we study linear fractional differential equations with variable coefficients. It is shown that, by assuming some conditions for the coefficients, the stationarity-conservation laws can be derived. The area where these are valid is restricted by the asymptotic properties of solutions of the respective equation. Applications of the proposed procedure include the fractional Fokker-Planck equation in (1+1)- and (d+1)-dimensional space and the fractional Klein-Kramers equation. (author)

Stationarity-conservation laws for fractional differential equations with variable coefficients

Energy Technology Data Exchange (ETDEWEB)

Klimek, Malgorzata [Institute of Mathematics and Computer Science, Technical University of Czestochowa, Czestochowa (Poland)

2002-08-09

In this paper, we study linear fractional differential equations with variable coefficients. It is shown that, by assuming some conditions for the coefficients, the stationarity-conservation laws can be derived. The area where these are valid is restricted by the asymptotic properties of solutions of the respective equation. Applications of the proposed procedure include the fractional Fokker-Planck equation in (1+1)- and (d+1)-dimensional space and the fractional Klein-Kramers equation. (author)

Conservation Laws and Traveling Wave Solutions of a Generalized Nonlinear ZK-BBM Equation

Directory of Open Access Journals (Sweden)

Khadijo Rashid Adem

2014-01-01

Full Text Available We study a generalized two-dimensional nonlinear Zakharov-Kuznetsov-Benjamin-Bona-Mahony (ZK-BBM equation, which is in fact Benjamin-Bona-Mahony equation formulated in the ZK sense. Conservation laws for this equation are constructed by using the new conservation theorem due to Ibragimov and the multiplier method. Furthermore, traveling wave solutions are obtained by employing the (G'/G-expansion method.

Symmetries and conservation laws for a sixth-order Boussinesq equation

International Nuclear Information System (INIS)

Recio, E.; Gandarias, M.L.; Bruzón, M.S.

2016-01-01

This paper considers a generalization depending on an arbitrary function f(u) of a sixth-order Boussinesq equation which arises in shallow water waves theory. Interestingly, this equation admits a Hamiltonian formulation when written as a system. A classification of point symmetries and conservation laws in terms of the function f(u) is presented for both, the generalized Boussinesq equation and the equivalent Hamiltonian system.

On double reductions from symmetries and conservation laws for a damped Boussinesq equation

International Nuclear Information System (INIS)

Gandarias, M.L.; Rosa, M.

2016-01-01

In this work, we study a Boussinesq equation with a strong damping term from the point of view of the Lie theory. We derive the classical Lie symmetries admitted by the equation as well as the reduced ordinary differential equations. Some nontrivial conservation laws are derived by using the multipliers method. Taking into account the relationship between symmetries and conservation laws and applying the double reduction method, we obtain a direct reduction of order of the ordinary differential equations and in particular a kink solution.

Stochastic Optimal Prediction with Application to Averaged Euler Equations

Energy Technology Data Exchange (ETDEWEB)

Bell, John [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Chorin, Alexandre J. [Univ. of California, Berkeley, CA (United States); Crutchfield, William [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)

2017-04-24

Optimal prediction (OP) methods compensate for a lack of resolution in the numerical solution of complex problems through the use of an invariant measure as a prior measure in the Bayesian sense. In first-order OP, unresolved information is approximated by its conditional expectation with respect to the invariant measure. In higher-order OP, unresolved information is approximated by a stochastic estimator, leading to a system of random or stochastic differential equations. We explain the ideas through a simple example, and then apply them to the solution of Averaged Euler equations in two space dimensions.

Symmetries, Traveling Wave Solutions, and Conservation Laws of a (3+1-Dimensional Boussinesq Equation

Directory of Open Access Journals (Sweden)

Letlhogonolo Daddy Moleleki

2014-01-01

Full Text Available We analyze the (3+1-dimensional Boussinesq equation, which has applications in fluid mechanics. We find exact solutions of the (3+1-dimensional Boussinesq equation by utilizing the Lie symmetry method along with the simplest equation method. The solutions obtained are traveling wave solutions. Moreover, we construct the conservation laws of the (3+1-dimensional Boussinesq equation using the new conservation theorem, which is due to Ibragimov.

Multicomponent fluid flow analysis using a new set of conservation equations

International Nuclear Information System (INIS)

Kamali, Reza; Emdad, Homayoon; Alishahi, Mohammad M

2008-01-01

In this work hydrodynamics of multicomponent ideal gas mixtures have been studied. Starting from the kinetic equations, the Eulerian approach is used to derive a new set of conservation equations for the multicomponent system where each component may have different velocity and kinetic temperature. The equations are based on the Grad's method of moment derived from the kinetic model in a relaxation time approximation (RTA). Based on this model which contains separate equation sets for each component of the system, a computer code has been developed for numerical computation of compressible flows of binary gas mixture in generalized curvilinear boundary conforming coordinates. Since these equations are similar to the Navier-Stokes equations for the single fluid systems, the same numerical methods are applied to these new equations. The Roe's numerical scheme is used to discretize the convective terms of governing fluid flow equations. The prepared algorithm and the computer code are capable of computing and presenting flow fields of each component of the system separately as well as the average flow field of the multicomponent gas system as a whole. Comparison of the present code results with those of a more common algorithm based on the mixture theory in a supersonic converging-diverging nozzle provides the validation of the present formulation. Afterwards, a more involved nozzle cooling problem with a binary ideal gas (helium-xenon) is chosen to compare the present results with those of the ordinary mixture theory. The present model provides the details of the flow fields of each component separately which is not available otherwise. It is also shown that the separate fluids treatment, such as the present study, is crucial when considering time scales on the order of (or shorter than) the intercollisions relaxation times.

Derivation of a volume-averaged neutron diffusion equation; Atomos para el desarrollo de Mexico

Energy Technology Data Exchange (ETDEWEB)

Vazquez R, R.; Espinosa P, G. [UAM-Iztapalapa, Av. San Rafael Atlixco 186, Col. Vicentina, Mexico D.F. 09340 (Mexico); Morales S, Jaime B. [UNAM, Laboratorio de Analisis en Ingenieria de Reactores Nucleares, Paseo Cuauhnahuac 8532, Jiutepec, Morelos 62550 (Mexico)]. e-mail: rvr@xanum.uam.mx

2008-07-01

This paper presents a general theoretical analysis of the problem of neutron motion in a nuclear reactor, where large variations on neutron cross sections normally preclude the use of the classical neutron diffusion equation. A volume-averaged neutron diffusion equation is derived which includes correction terms to diffusion and nuclear reaction effects. A method is presented to determine closure-relationships for the volume-averaged neutron diffusion equation (e.g., effective neutron diffusivity). In order to describe the distribution of neutrons in a highly heterogeneous configuration, it was necessary to extend the classical neutron diffusion equation. Thus, the volume averaged diffusion equation include two corrections factor: the first correction is related with the absorption process of the neutron and the second correction is a contribution to the neutron diffusion, both parameters are related to neutron effects on the interface of a heterogeneous configuration. (Author)

Application of an analytical method for solution of thermal hydraulic conservation equations

Energy Technology Data Exchange (ETDEWEB)

Fakory, M.R. [Simulation, Systems & Services Technologies Company (S3 Technologies), Columbia, MD (United States)

1995-09-01

An analytical method has been developed and applied for solution of two-phase flow conservation equations. The test results for application of the model for simulation of BWR transients are presented and compared with the results obtained from application of the explicit method for integration of conservation equations. The test results show that with application of the analytical method for integration of conservation equations, the Courant limitation associated with explicit Euler method of integration was eliminated. The results obtained from application of the analytical method (with large time steps) agreed well with the results obtained from application of explicit method of integration (with time steps smaller than the size imposed by Courant limitation). The results demonstrate that application of the analytical approach significantly improves the numerical stability and computational efficiency.

Conserving relativistic many-body approach: Equation of state, spectral function, and occupation probabilities of nuclear matter

International Nuclear Information System (INIS)

de Jong, F.; Malfliet, R.

1991-01-01

Starting from a relativistic Lagrangian we derive a ''conserving'' approximation for the description of nuclear matter. We show this to be a nontrivial extension over the relativistic Dirac-Brueckner scheme. The saturation point of the equation of state calculated agrees very well with the empirical saturation point. The conserving character of the approach is tested by means of the Hugenholtz--van Hove theorem. We find the theorem fulfilled very well around saturation. A new value for compression modulus is derived, K=310 MeV. Also we calculate the occupation probabilities at normal nuclear matter densities by means of the spectral function. The average depletion κ of the Fermi sea is found to be κ∼0.11

Generalized internal long wave equations: construction, hamiltonian structure and conservation laws

International Nuclear Information System (INIS)

Lebedev, D.R.

1982-01-01

Some aspects of the theory of the internal long-wave equations (ILW) are considered. A general class of the ILW type equations is constructed by means of the Zakharov-Shabat ''dressing'' method. Hamiltonian structure and infinite numbers of conservation laws are introduced. The considered equations are shown to be Hamiltonian in the so-called second Hamiltonian structu

76 FR 13168 - Energy Conservation Program for Consumer Products: Representative Average Unit Costs of Energy

Science.gov (United States)

2011-03-10

... average unit costs of residential energy in a Federal Register notice entitled, ``Energy Conservation... DEPARTMENT OF ENERGY Office of Energy Efficiency and Renewable Energy Energy Conservation Program for Consumer Products: Representative Average Unit Costs of Energy AGENCY: Office of Energy Efficiency...

78 FR 17648 - Energy Conservation Program for Consumer Products: Representative Average Unit Costs of Energy

Science.gov (United States)

2013-03-22

... Conservation Program for Consumer Products: Representative Average Unit Costs of Energy'', dated April 26, 2012... DEPARTMENT OF ENERGY Office of Energy Efficiency and Renewable Energy Energy Conservation Program for Consumer Products: Representative Average Unit Costs of Energy AGENCY: Office of Energy Efficiency...

Averaging Principle for the Higher Order Nonlinear Schrödinger Equation with a Random Fast Oscillation

Science.gov (United States)

Gao, Peng

2018-04-01

This work concerns the problem associated with averaging principle for a higher order nonlinear Schrödinger equation perturbed by a oscillating term arising as the solution of a stochastic reaction-diffusion equation evolving with respect to the fast time. This model can be translated into a multiscale stochastic partial differential equations. Stochastic averaging principle is a powerful tool for studying qualitative analysis of stochastic dynamical systems with different time-scales. To be more precise, under suitable conditions, we prove that there is a limit process in which the fast varying process is averaged out and the limit process which takes the form of the higher order nonlinear Schrödinger equation is an average with respect to the stationary measure of the fast varying process. Finally, by using the Khasminskii technique we can obtain the rate of strong convergence for the slow component towards the solution of the averaged equation, and as a consequence, the system can be reduced to a single higher order nonlinear Schrödinger equation with a modified coefficient.

Averaging Principle for the Higher Order Nonlinear Schrödinger Equation with a Random Fast Oscillation

Science.gov (United States)

Gao, Peng

2018-06-01

This work concerns the problem associated with averaging principle for a higher order nonlinear Schrödinger equation perturbed by a oscillating term arising as the solution of a stochastic reaction-diffusion equation evolving with respect to the fast time. This model can be translated into a multiscale stochastic partial differential equations. Stochastic averaging principle is a powerful tool for studying qualitative analysis of stochastic dynamical systems with different time-scales. To be more precise, under suitable conditions, we prove that there is a limit process in which the fast varying process is averaged out and the limit process which takes the form of the higher order nonlinear Schrödinger equation is an average with respect to the stationary measure of the fast varying process. Finally, by using the Khasminskii technique we can obtain the rate of strong convergence for the slow component towards the solution of the averaged equation, and as a consequence, the system can be reduced to a single higher order nonlinear Schrödinger equation with a modified coefficient.

Lie symmetry analysis and conservation laws for the time fractional fourth-order evolution equation

Directory of Open Access Journals (Sweden)

Wang Li

2017-06-01

Full Text Available In this paper, we study Lie symmetry analysis and conservation laws for the time fractional nonlinear fourth-order evolution equation. Using the method of Lie point symmetry, we provide the associated vector fields, and derive the similarity reductions of the equation, respectively. The method can be applied wisely and efficiently to get the reduced fractional ordinary differential equations based on the similarity reductions. Finally, by using the nonlinear self-adjointness method and Riemann-Liouville time-fractional derivative operator as well as Euler-Lagrange operator, the conservation laws of the equation are obtained.

Coarse-mesh discretized low-order quasi-diffusion equations for subregion averaged scalar fluxes

International Nuclear Information System (INIS)

Anistratov, D. Y.

2004-01-01

In this paper we develop homogenization procedure and discretization for the low-order quasi-diffusion equations on coarse grids for core-level reactor calculations. The system of discretized equations of the proposed method is formulated in terms of the subregion averaged group scalar fluxes. The coarse-mesh solution is consistent with a given fine-mesh discretization of the transport equation in the sense that it preserves a set of average values of the fine-mesh transport scalar flux over subregions of coarse-mesh cells as well as the surface currents, and eigenvalue. The developed method generates numerical solution that mimics the large-scale behavior of the transport solution within assemblies. (authors)

Noether's theorem and Steudel's conserved currents for the sine-Gordon equation

International Nuclear Information System (INIS)

Shadwick, W.F.

1980-01-01

A version of Noether's theorem appropriate for the extended Hamilton-Cartan formalism for regular first-order Lagrangians is proposed. Steudel's derivation of an infinite collection of conserved currents for the sine-Gordon equation is presented in this context and it is demonstrated that, as a consequence of the commutativity of the sine-Gordon Baecklund transformations, the conserved charges corresponding to these currents are in involution with respect to the natural Poisson bracket provided by the formalism. Thus one obtains the formal 'complete integrability' of the sine-Gordon equation as a consequence of the properties of the Baecklund transformation. (orig.)

Variational symmetries, conserved quantities and identities for several equations of mathematical physics

Energy Technology Data Exchange (ETDEWEB)

Donchev, Veliko, E-mail: velikod@ie.bas.bg [Laboratory “Physical Problems of Electron and Ion Technologies,” Institute of Electronics, Bulgarian Academy of Sciences, 72 Tzarigradsko shosse, 1784 Sofia (Bulgaria)

2014-03-15

We find variational symmetries, conserved quantities and identities for several equations: envelope equation, Böcher equation, the propagation of sound waves with losses, flow of a gas with losses, and the nonlinear Schrödinger equation with losses or gains, and an electro-magnetic interaction. Most of these equations do not have a variational description with the classical variational principle and we find such a description with the generalized variational principle of Herglotz.

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

International Nuclear Information System (INIS)

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

2013-01-01

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

A Note on the Invariance Properties and Conservation Laws of the Kadomstev—Petviashvili Equation with Power Law Nonlinearity

International Nuclear Information System (INIS)

Bokhari A H; Zaman F D; Fakhar K; Kara A H

2011-01-01

First, we studied the invariance properties of the Kadomstev—Petviashvili equation with power law nonlinearity. Then, we determined the complete class of conservation laws and stated the corresponding conserved densities which are useful in finding the conserved quantities of the equation. The point symmetry generators were also used to reduce the equation to an exact solution and to verify the invariance properties of the conserved flows. (general)

Reductions and conservation laws for BBM and modified BBM equations

Directory of Open Access Journals (Sweden)

Khorshidi Maryam

2016-01-01

Full Text Available In this paper, the classical Lie theory is applied to study the Benjamin-Bona-Mahony (BBM and modified Benjamin-Bona-Mahony equations (MBBM to obtain their symmetries, invariant solutions, symmetry reductions and differential invariants. By observation of the the adjoint representation of Mentioned symmetry groups on their Lie algebras, we find the primary classification (optimal system of their group-invariant solutions which provides new exact solutions to BBM and MBBM equations. Finally, conservation laws of the BBM and MBBM equations are presented. Some aspects of their symmetry properties are given too.

Averaging of the Equations of the Standard Cosmological Model over Rapid Oscillations

Science.gov (United States)

Ignat'ev, Yu. G.; Samigullina, A. R.

2017-11-01

An averaging of the equations of the standard cosmological model (SCM) is carried out. It is shown that the main contribution to the macroscopic energy density of the scalar field comes from its microscopic oscillations with the Compton period. The effective macroscopic equation of state of the oscillations of the scalar field corresponds to the nonrelativistic limit.

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

Science.gov (United States)

LeMesurier, Brenton

2012-01-01

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

On a quantum version of conservation laws for derivative nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Sen, S.; Chowdhury, A.R.

1988-01-01

The authors derived the quantum mechanical versions of infinite number of conservation laws associated with Derivative Nonlinear Schrodinger equation with the help of a methodology used in string theory. The renormalised version of the conserved quantities are obtained with explicit forms of the counter terms

Lie symmetry analysis, conservation laws and exact solutions of the seventh-order time fractional Sawada–Kotera–Ito equation

Directory of Open Access Journals (Sweden)

Emrullah Yaşar

Full Text Available In this paper Lie symmetry analysis of the seventh-order time fractional Sawada–Kotera–Ito (FSKI equation with Riemann–Liouville derivative is performed. Using the Lie point symmetries of FSKI equation, it is shown that it can be transformed into a nonlinear ordinary differential equation of fractional order with a new dependent variable. In the reduced equation the derivative is in Erdelyi–Kober sense. Furthermore, adapting the Ibragimov’s nonlocal conservation method to time fractional partial differential equations, we obtain conservation laws of the underlying equation. In addition, we construct some exact travelling wave solutions for the FSKI equation using the sub-equation method. Keywords: Fractional Sawada–Kotera–Ito equation, Lie symmetry, Riemann–Liouville fractional derivative, Conservation laws, Exact solutions

Numerical artifacts in the Generalized Porous Medium Equation: Why harmonic averaging itself is not to blame

Science.gov (United States)

Maddix, Danielle C.; Sampaio, Luiz; Gerritsen, Margot

2018-05-01

The degenerate parabolic Generalized Porous Medium Equation (GPME) poses numerical challenges due to self-sharpening and its sharp corner solutions. For these problems, we show results for two subclasses of the GPME with differentiable k (p) with respect to p, namely the Porous Medium Equation (PME) and the superslow diffusion equation. Spurious temporal oscillations, and nonphysical locking and lagging have been reported in the literature. These issues have been attributed to harmonic averaging of the coefficient k (p) for small p, and arithmetic averaging has been suggested as an alternative. We show that harmonic averaging is not solely responsible and that an improved discretization can mitigate these issues. Here, we investigate the causes of these numerical artifacts using modified equation analysis. The modified equation framework can be used for any type of discretization. We show results for the second order finite volume method. The observed problems with harmonic averaging can be traced to two leading error terms in its modified equation. This is also illustrated numerically through a Modified Harmonic Method (MHM) that can locally modify the critical terms to remove the aforementioned numerical artifacts.

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.

The role of energy conservation in the BFKL equation

International Nuclear Information System (INIS)

Forshaw, J.R.; Harriman, P.N.; Sutton, P.J.

1993-01-01

We study a modification to the BFKL equation at zero momentum transfer due to the imposition of energy conservation. The significance of our modification, which enters in the form of an ultraviolet cutoff, is illustrated directly and is discussed within the context of the gluon diffusion in k T . (Author)

Equations of motion and conservation laws in a theory of stably stratified turbulence

Energy Technology Data Exchange (ETDEWEB)

L' vov, Victor S; Rudenko, Oleksii [Department of Chemical Physics, Weizmann Institute of Science, Rehovot 76100 (Israel)], E-mail: oleksii.rudenko@weizmann.ac.il

2008-12-15

This paper is part of an invited talk given at the international conference 'Turbulent Mixing and Beyond'. We consider non-isothermal fluid flows and revise simplifications of basic hydrodynamic equations for such flows, arriving eventually at a generalization of the Oberbeck-Boussinesq approximation valid for arbitrary equation of state including both non-ideal gases as well as liquids. The proposed approach is based on a suggested general definition of potential temperature. Special attention is paid to the energy conservation principle: the proposed approximation exactly preserves the total mechanical energy by approximate equations of motion. It is emphasized explicitly the importance for any turbulent boundary layer model to respect the conservation laws.

Nonlinear MHD-equations: symmetries, solutions and conservation laws

International Nuclear Information System (INIS)

Samokhin, A.V.

1985-01-01

To investigate stability and nonlinear effects in a high-temperature plasma the system of two scalar nonlinear equations is considered. The algebra of classical symmetries of this system and a certain natural part of its conservation laws are described. It is shown that first, with symmetries one can derive invariant (self-similar) solutions, second, acting with symmetry on the known solution the latter can be included into parametric family

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

International Nuclear Information System (INIS)

Koiv, M.; Rosenhaus, V.

1979-01-01

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

Mass and energy-capital conservation equations to study the price evolution of non-renewable energy resources

International Nuclear Information System (INIS)

Gori, F.

2006-01-01

Mass conservation equation of non-renewable resources is employed to study the resources remaining in the reservoir according to the extraction policy. The energy conservation equation is transformed into an energy-capital conservation equation. The Hotelling rule is shown to be a special case of the general energy-capital conservation equation when the mass flow rate of extracted resources is equal to unity. Mass and energy-capital conservation equations are then coupled and solved together. It is investigated the price evolution of extracted resources. The conclusion of the Hotelling rule for non-extracted resources, i.e. an exponential increase of the price of non-renewable resources at the rate of current interest, is then generalized. A new parameter, called 'Price Increase Factor', PIF, is introduced as the difference between the current interest rate of capital and the mass flow rate of extraction of non-renewable resources. The price of extracted resources can increase exponentially only if PIF is greater than zero or if the mass flow rate of extraction is lower than the current interest rate of capital. The price is constant if PIF is zero or if the mass flow rate of extraction is equal to the current interest rate. The price is decreasing with time if PIF is smaller than zero or if the mass flow rate of extraction is higher than the current interest rate. (author)

A discrete model on Sierpinski gasket substrate for a conserved current equation with a conservative noise

International Nuclear Information System (INIS)

Kim, Dae Ho; Kim, Jin Min

2012-01-01

A conserved discrete model on the Sierpinski gasket substrate is studied. The interface width W in the model follows the Family–Vicsek dynamic scaling form with growth exponent β ≈ 0.0542, roughness exponent α ≈ 0.240 and dynamic exponent z ≈ 4.42. They satisfy a scaling relation α + z = 2z rw , where z rw is the random walk exponent of the fractal substrate. Also, they are in a good agreement with the predicted values of a fractional Langevin equation where η c is a conservative noise. (paper)

Traveling waves and conservation laws for highly nonlinear wave equations modeling Hertz chains

Science.gov (United States)

Przedborski, Michelle; Anco, Stephen C.

2017-09-01

A highly nonlinear, fourth-order wave equation that models the continuum theory of long wavelength pulses in weakly compressed, homogeneous, discrete chains with a general power-law contact interaction is studied. For this wave equation, all solitary wave solutions and all nonlinear periodic wave solutions, along with all conservation laws, are derived. The solutions are explicitly parameterized in terms of the asymptotic value of the wave amplitude in the case of solitary waves and the peak of the wave amplitude in the case of nonlinear periodic waves. All cases in which the solution expressions can be stated in an explicit analytic form using elementary functions are worked out. In these cases, explicit expressions for the total energy and total momentum for all solutions are obtained as well. The derivation of the solutions uses the conservation laws combined with an energy analysis argument to reduce the wave equation directly to a separable first-order differential equation that determines the wave amplitude in terms of the traveling wave variable. This method can be applied more generally to other highly nonlinear wave equations.

Symmetry Reductions, Exact Solutions and Conservation Laws of Asymmetric Nizhnik-Novikov-Veselov Equation

International Nuclear Information System (INIS)

Wang Ling; Dong Zhongzhou; Liu Xiqiang

2008-01-01

By applying a direct symmetry method, we get the symmetry of the asymmetric Nizhnik-Novikov-Veselov equation (ANNV). Taking the special case, we have a finite-dimensional symmetry. By using the equivalent vector of the symmetry, we construct an eight-dimensional symmetry algebra and get the optimal system of group-invariant solutions. To every case of the optimal system, we reduce the ANNV equation and obtain some solutions to the reduced equations. Furthermore, we find some new explicit solutions of the ANNV equation. At last, we give the conservation laws of the ANNV equation.

A general nonlinear evolution equation for irreversible conservative approach to stable equilibrium

International Nuclear Information System (INIS)

Beretta, G.P.

1986-01-01

This paper addresses a mathematical problem relevant to the question of nonequilibrium and irreversibility, namely, that of ''designing'' a general evolution equation capable of describing irreversible but conservative relaxtion towards equilibrium. The objective is to present an interesting mathematical solution to this design problem, namely, a new nonlinear evolution equation that satisfies a set of very stringent relevant requirements. Three different frameworks are defined from which the new equation could be adopted, with entirely different interpretations. Some useful well-known mathematics involving Gram determinants are presented and a nonlinear evolution equation is given which meets the stringent design specifications

Lie symmetry analysis, explicit solutions and conservation laws for the space-time fractional nonlinear evolution equations

Science.gov (United States)

Inc, Mustafa; Yusuf, Abdullahi; Aliyu, Aliyu Isa; Baleanu, Dumitru

2018-04-01

This paper studies the symmetry analysis, explicit solutions, convergence analysis, and conservation laws (Cls) for two different space-time fractional nonlinear evolution equations with Riemann-Liouville (RL) derivative. The governing equations are reduced to nonlinear ordinary differential equation (ODE) of fractional order using their Lie point symmetries. In the reduced equations, the derivative is in Erdelyi-Kober (EK) sense, power series technique is applied to derive an explicit solutions for the reduced fractional ODEs. The convergence of the obtained power series solutions is also presented. Moreover, the new conservation theorem and the generalization of the Noether operators are developed to construct the nonlocal Cls for the equations . Some interesting figures for the obtained explicit solutions are presented.

Justification of the averaging method for parabolic equations containing rapidly oscillating terms with large amplitudes

International Nuclear Information System (INIS)

Levenshtam, V B

2006-01-01

We justify the averaging method for abstract parabolic equations with stationary principal part that contain non-linearities (subordinate to the principal part) some of whose terms are rapidly oscillating in time with zero mean and are proportional to the square root of the frequency of oscillation. Our interest in the exponent 1/2 is motivated by the fact that terms proportional to lower powers of the frequency have no influence on the average. For linear equations of the same type, we justify an algorithm for the study of the stability of solutions in the case when the stationary averaged problem has eigenvalues on the imaginary axis (the critical case)

Three-parameter relativistic dynamics. 1. Equation of motion, energy conservation

International Nuclear Information System (INIS)

Rogachevskii, A.G.

1995-01-01

A formally geometric analog of the relativistic dynamics of a point charged particle is constructed. Time as a function of the spatial coordinates is taken as the trajectory equation, i.e., the trajectory is a hypersurface in Minkowski space. The dynamics is presented. The law of open-quotes energyclose quotes conservation is examined

Relativistic dissipative hydrodynamic equations at the second order for multi-component systems with multiple conserved currents

International Nuclear Information System (INIS)

Monnai, Akihiko; Hirano, Tetsufumi

2010-01-01

We derive the second order hydrodynamic equations for the relativistic system of multi-components with multiple conserved currents by generalizing the Israel-Stewart theory and Grad's moment method. We find that, in addition to the conventional moment equations, extra moment equations associated with conserved currents should be introduced to consistently match the number of equations with that of unknowns and to satisfy the Onsager reciprocal relations. Consistent expansion of the entropy current leads to constitutive equations which involve the terms not appearing in the original Israel-Stewart theory even in the single component limit. We also find several terms which exhibit thermal diffusion such as Soret and Dufour effects. We finally compare our results with those of other existing formalisms.

Numerical solution of conservation equations in the transient model for the system thermal - hydraulics in the Korsar computer code

International Nuclear Information System (INIS)

Yudov, Y.V.

2001-01-01

The functional part of the KORSAR computer code is based on the computational unit for the reactor system thermal-hydraulics and other thermal power systems with water cooling. The two-phase flow dynamics of the thermal-hydraulic network is modelled by KORSAR in one-dimensional two-fluid (non-equilibrium and nonhomogeneous) approximation with the same pressure of both phases. Each phase is characterized by parameters averaged over the channel sections, and described by the conservation equations for mass, energy and momentum. The KORSAR computer code relies upon a novel approach to mathematical modelling of two-phase dispersed-annular flows. This approach allows a two-fluid model to differentiate the effects of the liquid film and droplets in the gas core on the flow characteristics. A semi-implicit numerical scheme has been chosen for deriving discrete analogs the conservation equations in KORSAR. In the semi-implicit numerical scheme, solution of finite-difference equations is reduced to the problem of determining the pressure field at a new time level. For the one-channel case, the pressure field is found from the solution of a system of linear algebraic equations by using the tri-diagonal matrix method. In the branched network calculation, the matrix of coefficients in the equations describing the pressure field is no longer tri-diagonal but has a sparseness structure. In this case, the system of linear equations for the pressure field can be solved with any of the known classical methods. Such an approach is implemented in the existing best-estimate thermal-hydraulic computer codes (TRAC, RELAP5, etc.) For the KORSAR computer code, we have developed a new non-iterative method for calculating the pressure field in the network of any topology. This method is based on the tri-diagonal matrix method and performs well when solving the thermal-hydraulic network problems. (author)

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

International Nuclear Information System (INIS)

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

Lie symmetry analysis, exact solutions and conservation laws for the time fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation

Science.gov (United States)

Baleanu, Dumitru; Inc, Mustafa; Yusuf, Abdullahi; Aliyu, Aliyu Isa

2018-06-01

In this work, we investigate the Lie symmetry analysis, exact solutions and conservation laws (Cls) to the time fractional Caudrey-Dodd-Gibbon-Sawada-Kotera (CDGDK) equation with Riemann-Liouville (RL) derivative. The time fractional CDGDK is reduced to nonlinear ordinary differential equation (ODE) of fractional order. New exact traveling wave solutions for the time fractional CDGDK are obtained by fractional sub-equation method. In the reduced equation, the derivative is in Erdelyi-Kober (EK) sense. Ibragimov's nonlocal conservation method is applied to construct Cls for time fractional CDGDK.

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.

A generalized variational algebra and conserved densities for linear evolution equations

International Nuclear Information System (INIS)

Abellanas, L.; Galindo, A.

1978-01-01

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

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

International Nuclear Information System (INIS)

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

1984-08-01

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

AKNS hierarchy, Darboux transformation and conservation laws of the 1D nonautonomous nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Zhao Dun; Zhang Yujuan; Lou Weiwei; Luo Honggang

2011-01-01

By constructing nonisospectral Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy, we investigate the nonautonomous nonlinear Schroedinger (NLS) equations which have been used to describe the Feshbach resonance management in matter-wave solitons in Bose-Einstein condensate and the dispersion and nonlinearity managements for optical solitons. It is found that these equations are some special cases of a new integrable model of nonlocal nonautonomous NLS equations. Based on the Lax pairs, the Darboux transformation and conservation laws are explored. It is shown that the local external potentials would break down the classical infinite number of conservation laws. The result indicates that the integrability of the nonautonomous NLS systems may be nontrivial in comparison to the conventional concept of integrability in the canonical case.

Stability analysis of explicit entropy viscosity methods for non-linear scalar conservation equations

KAUST Repository

Bonito, Andrea; Guermond, Jean-Luc; Popov, Bojan

2013-01-01

We establish the L2-stability of an entropy viscosity technique applied to nonlinear scalar conservation equations. First-and second-order explicit time-stepping techniques using continuous finite elements in space are considered. The method

Averaging in spherically symmetric cosmology

International Nuclear Information System (INIS)

Coley, A. A.; Pelavas, N.

2007-01-01

The averaging problem in cosmology is of fundamental importance. When applied to study cosmological evolution, the theory of macroscopic gravity (MG) can be regarded as a long-distance modification of general relativity. In the MG approach to the averaging problem in cosmology, the Einstein field equations on cosmological scales are modified by appropriate gravitational correlation terms. We study the averaging problem within the class of spherically symmetric cosmological models. That is, we shall take the microscopic equations and effect the averaging procedure to determine the precise form of the correlation tensor in this case. In particular, by working in volume-preserving coordinates, we calculate the form of the correlation tensor under some reasonable assumptions on the form for the inhomogeneous gravitational field and matter distribution. We find that the correlation tensor in a Friedmann-Lemaitre-Robertson-Walker (FLRW) background must be of the form of a spatial curvature. Inhomogeneities and spatial averaging, through this spatial curvature correction term, can have a very significant dynamical effect on the dynamics of the Universe and cosmological observations; in particular, we discuss whether spatial averaging might lead to a more conservative explanation of the observed acceleration of the Universe (without the introduction of exotic dark matter fields). We also find that the correlation tensor for a non-FLRW background can be interpreted as the sum of a spatial curvature and an anisotropic fluid. This may lead to interesting effects of averaging on astrophysical scales. We also discuss the results of averaging an inhomogeneous Lemaitre-Tolman-Bondi solution as well as calculations of linear perturbations (that is, the backreaction) in an FLRW background, which support the main conclusions of the analysis

Conservation laws derived by the Neutral-Action Method. A simple application to the Schroedinger equation

International Nuclear Information System (INIS)

Nordbrock, U.; Kienzler, R.

2007-01-01

Conservation laws are a recognized tool in physical and engineering sciences. The classical procedure to construct conservation laws is to apply Noether's Theorem. It requires the existence of a Lagrange-function for the system under consideration. Two unknown sets of functions have to be found. A broader class of such laws is obtainable, if Noether's Theorem is used together with the Bessel-Hagen extension, raising the number of sets of unknown functions to three. By using the recently developed Neutral-Action Method, the same conservation laws can be obtained by calculating only one unknown set of functions. Moreover the Neutral Action Method can also be applied in the absence of a Lagrangian, since only the governing differential equations are required for this procedure. In the paper, an application of this method to the Schroedinger equation is presented. (authors)

Students' Understanding of Conservation of Matter, Stoichiometry and Balancing Equations in Indonesia

Science.gov (United States)

Agung, Salamah; Schwartz, Marc S.

2007-01-01

This study examines Indonesian students' understanding of conservation of matter, balancing of equations and stoichiometry. Eight hundred and sixty-seven Grade 12 students from 22 schools across four different cities in two developed provinces in Indonesia participated in the study. Nineteen teachers also participated in order to validate the…

Subchannel analysis program for boiling water reactor fuel bundles based on five conservation equations of two-phase flow

International Nuclear Information System (INIS)

Bessho, Y.; Uchikawa, S.

1985-01-01

A subchannel analysis program, MENUETT, is developed for evaluation of thermal-hydraulic characteristics in boiling water reactor fuel bundles. This program is based on five conservation equations of two-phase flow with the drift-flux correlation. The cross flows are calculated separately for liquid and vapor phases from the lateral momentum conservation equation. The effects of turbulent mixing and void drift are accounted for in the program. The conservation equations are implicitly differentiated with the convective terms by the donor-cell method, and are solved iteratively in the axial and lateral directions. Data of the 3 X 3 rod bundle experiments are used for program verification. The lateral distributions of equilibrium quality and mass flow rate at the bundle exit calculated by the program compare satisfactorily with the experimental results

Flow and transport simulation of Madeira River using three depth-averaged two-equation turbulence closure models

Directory of Open Access Journals (Sweden)

Li-ren Yu

2012-03-01

Full Text Available This paper describes a numerical simulation in the Amazon water system, aiming to develop a quasi-three-dimensional numerical tool for refined modeling of turbulent flow and passive transport of mass in natural waters. Three depth-averaged two-equation turbulence closure models, k˜−ε˜,k˜−w˜, and k˜−ω˜ , were used to close the non-simplified quasi-three dimensional hydrodynamic fundamental governing equations. The discretized equations were solved with the advanced multi-grid iterative method using non-orthogonal body-fitted coarse and fine grids with collocated variable arrangement. Except for steady flow computation, the processes of contaminant inpouring and plume development at the beginning of discharge, caused by a side-discharge of a tributary, have also been numerically investigated. The three depth-averaged two-equation closure models are all suitable for modeling strong mixing turbulence. The newly established turbulence models such as the k˜−ω˜ model, with a higher order of magnitude of the turbulence parameter, provide a possibility for improving computational precision.

Reynolds Averaged Navier-Stokes (RANS) equation solutions of wind turbine wakes

Energy Technology Data Exchange (ETDEWEB)

Ludwig, Daniel Evandro; Horn, Diego Anderson; Petry, Adriane Prisco [Thermal and Energy Study Group, Mechanical Engeneering Department, Federal University of Rio Grande do Sul, Porto Alegre (Brazil)], E-mail: adrianep@mecanica.ufrgs.br

2010-07-01

This paper aims to evaluate the influence of three different turbulence models in the study of a wind turbine wake. Numerical Simulation is used as working tool to characterize the flow through the wind turbines, it is used the numeric simulation. The numerical analysis is based on the finite volume method and the Reynolds Averaged Navier-Stokes (RANS) equations. Three turbulence models are used to represent the total effects of turbulence in the flow: the two equations k-classical and the RNG k- models, based on the turbulent viscosity; and the Shear Stress Transport (SST) model, based on the transport of the Reynolds tensor. The results of the 'u' velocity profiles are compared to experimental data from Vermeer (2003) at distances equivalent to 2, 4, 6, 8, 10 and 16 diameters downstream from the turbine. Results shows that the SST model gives better results until 6 diameters, beyond this distance there is no significant differences between the compared models. (author)

Deriving average soliton equations with a perturbative method

International Nuclear Information System (INIS)

Ballantyne, G.J.; Gough, P.T.; Taylor, D.P.

1995-01-01

The method of multiple scales is applied to periodically amplified, lossy media described by either the nonlinear Schroedinger (NLS) equation or the Korteweg--de Vries (KdV) equation. An existing result for the NLS equation, derived in the context of nonlinear optical communications, is confirmed. The method is then applied to the KdV equation and the result is confirmed numerically

Conservation laws and self-consistent sources for a super-CKdV equation hierarchy

International Nuclear Information System (INIS)

Li Li

2011-01-01

From the super-matrix Lie algebras, we consider a super-extension of the CKdV equation hierarchy in the present Letter, and propose the super-CKdV hierarchy with self-consistent sources. Furthermore, we establish the infinitely many conservation laws for the integrable super-CKdV hierarchy.

Conservation laws and self-consistent sources for a super-CKdV equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

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

2011-03-14

From the super-matrix Lie algebras, we consider a super-extension of the CKdV equation hierarchy in the present Letter, and propose the super-CKdV hierarchy with self-consistent sources. Furthermore, we establish the infinitely many conservation laws for the integrable super-CKdV hierarchy.

Lie symmetry and the generalized Hojman conserved quantity of Nielsen equations for a variable mass holonomic system of relative motion

International Nuclear Information System (INIS)

Zhang Mei-Ling; Wang Xiao-Xiao; Xie Yin-Li; Jia Li-Qun; Sun Xian-Ting

2011-01-01

Lie symmetry and the generalized Hojman conserved quantity of Nielsen equations for a variable mass holonomic system of relative motion are studied. The determining equation of Lie symmetry of Nielsen equations for a variable mass holonomic system of relative motion under the infinitesimal transformations of groups is given. The expression of generalized Hojman conserved quantity deduced directly from Lie symmetry for a variable mass holonomic system of relative motion is obtained. An example is given to illustrate the application of the results. (general)

A discrete model on Sierpinski gasket substrate for a conserved current equation with a conservative noise

Science.gov (United States)

Kim, Dae Ho; Kim, Jin Min

2012-09-01

A conserved discrete model on the Sierpinski gasket substrate is studied. The interface width W in the model follows the Family-Vicsek dynamic scaling form with growth exponent β ≈ 0.0542, roughness exponent α ≈ 0.240 and dynamic exponent z ≈ 4.42. They satisfy a scaling relation α + z = 2zrw, where zrw is the random walk exponent of the fractal substrate. Also, they are in a good agreement with the predicted values of a fractional Langevin equation \\frac{\\partial h}{\\partial t}={\

Discrete conservation laws and the convergence of long time simulations of the mkdv equation

Science.gov (United States)

Gorria, C.; Alejo, M. A.; Vega, L.

2013-02-01

Pseudospectral collocation methods and finite difference methods have been used for approximating an important family of soliton like solutions of the mKdV equation. These solutions present a structural instability which make difficult to approximate their evolution in long time intervals with enough accuracy. The standard numerical methods do not guarantee the convergence to the proper solution of the initial value problem and often fail by approaching solutions associated to different initial conditions. In this frame the numerical schemes that preserve the discrete invariants related to some conservation laws of this equation produce better results than the methods which only take care of a high consistency order. Pseudospectral spatial discretization appear as the most robust of the numerical methods, but finite difference schemes are useful in order to analyze the rule played by the conservation of the invariants in the convergence.

Physics-preserving averaging scheme based on Grunwald-Letnikov formula for gas flow in fractured media

KAUST Repository

Amir, Sahar Z.

2018-01-02

The heterogeneous natures of rock fabrics, due to the existence of multi-scale fractures and geological formations, led to the deviations from unity in the flux-equations fractional-exponent magnitudes. In this paper, the resulting non-Newtonian non-Darcy fractional-derivatives flux equations are solved using physics-preserving averaging schemes that incorporates both, original and shifted, Grunwald-Letnikov (GL) approximation formulas preserving the physics, by reducing the shifting effects, while maintaining the stability of the system, by keeping one shifted expansion. The proposed way of using the GL expansions also generate symmetrical coefficient matrices that significantly reduces the discretization complexities appearing with all shifted cases from literature, and help considerably in 2D and 3D systems. Systems equations derivations and discretization details are discussed. Then, the physics-preserving averaging scheme is explained and illustrated. Finally, results are presented and reviewed. Edge-based original GL expansions are unstable as also illustrated in literatures. Shifted GL expansions are stable but add a lot of additional weights to both discretization sides affecting the physical accuracy. In comparison, the physics-preserving averaging scheme balances the physical accuracy and stability requirements leading to a more physically conservative scheme that is more stable than the original GL approximation but might be slightly less stable than the shifted GL approximations. It is a locally conservative Single-Continuum averaging scheme that applies a finite-volume viewpoint.

Physics-preserving averaging scheme based on Grunwald-Letnikov formula for gas flow in fractured media

KAUST Repository

Amir, Sahar Z.; Sun, Shuyu

2018-01-01

The heterogeneous natures of rock fabrics, due to the existence of multi-scale fractures and geological formations, led to the deviations from unity in the flux-equations fractional-exponent magnitudes. In this paper, the resulting non-Newtonian non-Darcy fractional-derivatives flux equations are solved using physics-preserving averaging schemes that incorporates both, original and shifted, Grunwald-Letnikov (GL) approximation formulas preserving the physics, by reducing the shifting effects, while maintaining the stability of the system, by keeping one shifted expansion. The proposed way of using the GL expansions also generate symmetrical coefficient matrices that significantly reduces the discretization complexities appearing with all shifted cases from literature, and help considerably in 2D and 3D systems. Systems equations derivations and discretization details are discussed. Then, the physics-preserving averaging scheme is explained and illustrated. Finally, results are presented and reviewed. Edge-based original GL expansions are unstable as also illustrated in literatures. Shifted GL expansions are stable but add a lot of additional weights to both discretization sides affecting the physical accuracy. In comparison, the physics-preserving averaging scheme balances the physical accuracy and stability requirements leading to a more physically conservative scheme that is more stable than the original GL approximation but might be slightly less stable than the shifted GL approximations. It is a locally conservative Single-Continuum averaging scheme that applies a finite-volume viewpoint.

Self-consistent field theory of collisions: Orbital equations with asymptotic sources and self-averaged potentials

Energy Technology Data Exchange (ETDEWEB)

Hahn, Y.K., E-mail: ykhahn22@verizon.net

2014-12-15

The self-consistent field theory of collisions is formulated, incorporating the unique dynamics generated by the self-averaged potentials. The bound state Hartree–Fock approach is extended for the first time to scattering states, by properly resolving the principal difficulties of non-integrable continuum orbitals and imposing complex asymptotic conditions. The recently developed asymptotic source theory provides the natural theoretical basis, as the asymptotic conditions are completely transferred to the source terms and the new scattering function is made fullyintegrable. The scattering solutions can then be directly expressed in terms of bound state HF configurations, establishing the relationship between the bound and scattering state solutions. Alternatively, the integrable spin orbitals are generated by constructing the individual orbital equations that contain asymptotic sources and self-averaged potentials. However, the orbital energies are not determined by the equations, and a special channel energy fixing procedure is developed to secure the solutions. It is also shown that the variational construction of the orbital equations has intrinsic ambiguities that are generally associated with the self-consistent approach. On the other hand, when a small subset of open channels is included in the source term, the solutions are only partiallyintegrable, but the individual open channels can then be treated more simply by properly selecting the orbital energies. The configuration mixing and channel coupling are then necessary to complete the solution. The new theory improves the earlier continuum HF model. - Highlights: • First extension of HF to scattering states, with proper asymptotic conditions. • Orbital equations with asymptotic sources and integrable orbital solutions. • Construction of self-averaged potentials, and orbital energy fixing. • Channel coupling and configuration mixing, involving the new orbitals. • Critical evaluation of the

77 FR 24940 - Energy Conservation Program for Consumer Products: Representative Average Unit Costs of Energy

Science.gov (United States)

2012-04-26

... 5-year average ratio with heating oil prices published in the Monthly Energy Review, but the propane... DEPARTMENT OF ENERGY Office of Energy Efficiency and Renewable Energy Energy Conservation Program... and Renewable Energy, Department of Energy. ACTION: Notice. SUMMARY: In this notice, the U.S...

Averaging of nonlinearity-managed pulses

International Nuclear Information System (INIS)

Zharnitsky, Vadim; Pelinovsky, Dmitry

2005-01-01

We consider the nonlinear Schroedinger equation with the nonlinearity management which describes Bose-Einstein condensates under Feshbach resonance. By using an averaging theory, we derive the Hamiltonian averaged equation and compare it with other averaging methods developed for this problem. The averaged equation is used for analytical approximations of nonlinearity-managed solitons

Four-level conservative finite-difference schemes for Boussinesq paradigm equation

Science.gov (United States)

Kolkovska, N.

2013-10-01

In this paper a two-parametric family of four level conservative finite difference schemes is constructed for the multidimensional Boussinesq paradigm equation. The schemes are explicit in the sense that no inner iterations are needed for evaluation of the numerical solution. The preservation of the discrete energy with this method is proved. The schemes have been numerically tested on one soliton propagation model and two solitons interaction model. The numerical experiments demonstrate that the proposed family of schemes has second order of convergence in space and time steps in the discrete maximal norm.

Mass, momentum and energy conserving (MaMEC) discretizations on general grids for the compressible Euler and shallow water equations

NARCIS (Netherlands)

Hof, Bas van ’t; Veldman, Arthur E.P.

2012-01-01

The paper explains a method by which discretizations of the continuity and momentum equations can be designed, such that they can be combined with an equation of state into a discrete energy equation. The resulting 'MaMEC' discretizations conserve mass, momentum as well as energy, although no

75 FR 13123 - Energy Conservation Program for Consumer Products: Representative Average Unit Costs of Energy

Science.gov (United States)

2010-03-18

... that of heating oil, based on the 2004-2008 averages for these two fuels. The source for these price... DEPARTMENT OF ENERGY Office of Energy Efficiency and Renewable Energy Energy Conservation Program... and Renewable Energy, Department of Energy. ACTION: Notice. SUMMARY: In this notice, the U.S...

A conservative spectral method for the Boltzmann equation with anisotropic scattering and the grazing collisions limit

International Nuclear Information System (INIS)

Gamba, Irene M.; Haack, Jeffrey R.

2014-01-01

We present the formulation of a conservative spectral method for the Boltzmann collision operator with anisotropic scattering cross-sections. The method is an extension of the conservative spectral method of Gamba and Tharkabhushanam [17,18], which uses the weak form of the collision operator to represent the collisional term as a weighted convolution in Fourier space. The method is tested by computing the collision operator with a suitably cut-off angular cross section and comparing the results with the solution of the Landau equation. We analytically study the convergence rate of the Fourier transformed Boltzmann collision operator in the grazing collisions limit to the Fourier transformed Landau collision operator under the assumption of some regularity and decay conditions of the solution to the Boltzmann equation. Our results show that the angular singularity which corresponds to the Rutherford scattering cross section is the critical singularity for which a grazing collision limit exists for the Boltzmann operator. Additionally, we numerically study the differences between homogeneous solutions of the Boltzmann equation with the Rutherford scattering cross section and an artificial cross section, which give convergence to solutions of the Landau equation at different asymptotic rates. We numerically show the rate of the approximation as well as the consequences for the rate of entropy decay for homogeneous solutions of the Boltzmann equation and Landau equation

A High-Accuracy Linear Conservative Difference Scheme for Rosenau-RLW Equation

Directory of Open Access Journals (Sweden)

Jinsong Hu

2013-01-01

Full Text Available We study the initial-boundary value problem for Rosenau-RLW equation. We propose a three-level linear finite difference scheme, which has the theoretical accuracy of Oτ2+h4. The scheme simulates two conservative properties of original problem well. The existence, uniqueness of difference solution, and a priori estimates in infinite norm are obtained. Furthermore, we analyze the convergence and stability of the scheme by energy method. At last, numerical experiments demonstrate the theoretical results.

A conservative local discontinuous Galerkin method for the solution of nonlinear Schr(o)dinger equation in two dimensions

Institute of Scientific and Technical Information of China (English)

ZHANG RongPei; YU XiJun; LI MingJun; LI XiangGui

2017-01-01

In this study,we present a conservative local discontinuous Galerkin (LDG) method for numerically solving the two-dimensional nonlinear Schr(o)dinger (NLS) equation.The NLS equation is rewritten as a firstorder system and then we construct the LDG formulation with appropriate numerical flux.The mass and energy conserving laws for the semi-discrete formulation can be proved based on different choices of numerical fluxes such as the central,alternative and upwind-based flux.We will propose two kinds of time discretization methods for the semi-discrete formulation.One is based on Crank-Nicolson method and can be proved to preserve the discrete mass and energy conservation.The other one is Krylov implicit integration factor (ⅡF) method which demands much less computational effort.Various numerical experiments are presented to demonstrate the conservation law of mass and energy,the optimal rates of convergence,and the blow-up phenomenon.

Critical review of conservation equations for two-phase flow in the U.S. NRC TRACE code

International Nuclear Information System (INIS)

Wulff, Wolfgang

2011-01-01

Research highlights: → Field equations as implemented in TRACE are incorrect. → Boundary conditions needed for cooling of nuclear fuel elements are wrong. → The two-fluid model in TRACE is not closed. → Three-dimensional flow modeling in TRACE has no basis. - Abstract: The field equations for two-phase flow in the computer code TRAC/RELAP Advanced Computational Engine or TRACE are examined to determine their validity, their capabilities and limitations in resolving nuclear reactor safety issues. TRACE was developed for the NRC to predict thermohydraulic phenomena in nuclear power plants during operational transients and postulated accidents. TRACE is based on the rigorously derived and well-established two-fluid field equations for 1-D and 3-D two-phase flow. It is shown that: (1)The two-fluid field equations for mass conservation as implemented in TRACE are wrong because local mass balances in TRACE are in conflict with mass conservation for the whole reactor system, as shown in Section . (2)Wrong equations of motion are used in TRACE in place of momentum balances, compromising at branch points the prediction of momentum transfer between, and the coupling of, loops in hydraulic networks by impedance (form loss and wall shear) and by inertia and thereby the simulation of reactor component interactions. (3)Most seriously, TRACE calculation of heat transfer from fuel elements is incorrect for single and two-phase flows, because Eq. of the TRACE Manual is wrong (see Section ). (4)Boundary conditions for momentum and energy balances in TRACE are restricted to flow regimes with single-phase wall contact because TRACE lacks constitutive relations for solid-fluid exchange of momentum and heat in prevailing flow regimes. Without a quantified assessment of consequences from (3) to (4), predictions of phasic fluid velocities, fuel temperatures and important safety parameters, e.g., peak clad temperature, are questionable. Moreover, TRACE cannot predict 3-D single- or

Stability analysis of explicit entropy viscosity methods for non-linear scalar conservation equations

KAUST Repository

Bonito, Andrea

2013-10-03

We establish the L2-stability of an entropy viscosity technique applied to nonlinear scalar conservation equations. First-and second-order explicit time-stepping techniques using continuous finite elements in space are considered. The method is shown to be stable independently of the polynomial degree of the space approximation under the standard CFL condition. © 2013 American Mathematical Society.

A maximum-principle preserving finite element method for scalar conservation equations

KAUST Repository

Guermond, Jean-Luc; Nazarov, Murtazo

2014-01-01

This paper introduces a first-order viscosity method for the explicit approximation of scalar conservation equations with Lipschitz fluxes using continuous finite elements on arbitrary grids in any space dimension. Provided the lumped mass matrix is positive definite, the method is shown to satisfy the local maximum principle under a usual CFL condition. The method is independent of the cell type; for instance, the mesh can be a combination of tetrahedra, hexahedra, and prisms in three space dimensions. © 2014 Elsevier B.V.

A maximum-principle preserving finite element method for scalar conservation equations

KAUST Repository

Guermond, Jean-Luc

2014-04-01

This paper introduces a first-order viscosity method for the explicit approximation of scalar conservation equations with Lipschitz fluxes using continuous finite elements on arbitrary grids in any space dimension. Provided the lumped mass matrix is positive definite, the method is shown to satisfy the local maximum principle under a usual CFL condition. The method is independent of the cell type; for instance, the mesh can be a combination of tetrahedra, hexahedra, and prisms in three space dimensions. © 2014 Elsevier B.V.

Mass, momentum and energy conserving (MaMEC) discretizations on general grids for the compressible Euler and shallow water equations

International Nuclear Information System (INIS)

Hof, Bas van’t; Veldman, Arthur E.P.

2012-01-01

The paper explains a method by which discretizations of the continuity and momentum equations can be designed, such that they can be combined with an equation of state into a discrete energy equation. The resulting ‘MaMEC’ discretizations conserve mass, momentum as well as energy, although no explicit conservation law for the total energy is present. Essential ingredients are (i) discrete convection that leaves the discrete energy invariant, and (ii) discrete consistency between the thermodynamic terms. Of particular relevance is the way in which finite volume fluxes are related to nodal values. The method is an extension of existing methods based on skew-symmetry of discrete operators, because it allows arbitrary equations of state and a larger class of grids than earlier methods. The method is first illustrated with a one-dimensional example on a highly stretched staggered grid, in which the MaMEC method calculates qualitatively correct results and a non-skew-symmetric finite volume method becomes unstable. A further example is a two-dimensional shallow water calculation on a rectilinear grid as well as on an unstructured grid. The conservation of mass, momentum and energy is checked, and losses are found negligible up to machine accuracy.

Integrating factors and conservation theorems for Hamilton's canonical equations of motion of variable mass nonholonomic nonconservative dynamical systems

Institute of Scientific and Technical Information of China (English)

李仁杰; 乔永芬; 刘洋

2002-01-01

We present a general approach to the construction of conservation laws for variable mass nonholonomic noncon-servative systems. First, we give the definition of integrating factors, and we study in detail the necessary conditionsfor the existence of the conserved quantities. Then, we establish the conservation theorem and its inverse theorem forHamilton's canonical equations of motion of variable mass nonholonomic nonconservative dynamical systems. Finally,we give an example to illustrate the application of the results.

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

International Nuclear Information System (INIS)

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

1993-01-01

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

Theoretical Maxwell's Equations, Gauge Field and Their Universality Based on One Conservation Law

Institute of Scientific and Technical Information of China (English)

Liu Changmao

2005-01-01

The notion of the inner product of vectors is extended to tensors of different orders, which may replace the vector product usually. The essences of the differential and the codifferential forms are pointed out: they represent the tangent surface and the normal surface fluxes of a tensor, respectively. The definitions of the divergence and the curl of a 2D surface flux of a tensor are obtained.Maxwell's equations, namely, the construction law of field, which were usually established based on two conservation laws of electric charge and imaginary magnetic charge, are derived by the author only by using one conservation law ( mass or fluid flux quantity and so on) and the feature of central field ( or its composition). By the feature of central field ( or its composition), the curl of 2D flux is zero. Both universality of gauge field and the difficulty of magnetic monopole theory ( a magnetic monopole has no effect on electric current just like a couple basing no effect on the sum of forces) are presented: magnetic monopole has no the feature of magnet. Finally it is pointed out that the base of relation of mass and energy is already involved in Maxwell's equations.

Nongeostrophic theory of zonally averaged circulation. I - Formulation

Science.gov (United States)

Tung, Ka Kit

1986-01-01

A nongeostrophic theory of zonally averaged circulation is formulated using the nonlinear primitive equations (mass conservation, thermodynamics, and zonal momentum) on a sphere. The relationship between the mean meridional circulation and diabatic heating rate is studied. Differences between results of nongeostropic theory and the geostrophic formulation concerning the role of eddy forcing of the diabatic circulation and the nonlinear nearly inviscid limit versus the geostrophic limit are discussed. Consideration is given to the Eliassen-Palm flux divergence, the Eliassen-Palm pseudodivergence, the nonacceleration theorem, and the nonlinear nongeostrophic Taylor relationship.

On the XFEL Schrödinger Equation: Highly Oscillatory Magnetic Potentials and Time Averaging

KAUST Repository

Antonelli, Paolo

2014-01-14

We analyse a nonlinear Schrödinger equation for the time-evolution of the wave function of an electron beam, interacting selfconsistently through a Hartree-Fock nonlinearity and through the repulsive Coulomb interaction of an atomic nucleus. The electrons are supposed to move under the action of a time dependent, rapidly periodically oscillating electromagnetic potential. This can be considered a simplified effective single particle model for an X-ray free electron laser. We prove the existence and uniqueness for the Cauchy problem and the convergence of wave-functions to corresponding solutions of a Schrödinger equation with a time-averaged Coulomb potential in the high frequency limit for the oscillations of the electromagnetic potential. © 2014 Springer-Verlag Berlin Heidelberg.

On the coupling of systems of hyperbolic conservation laws with ordinary differential equations

International Nuclear Information System (INIS)

Borsche, Raul; Colombo, Rinaldo M; Garavello, Mauro

2010-01-01

Motivated by applications to the piston problem, to a manhole model, to blood flow and to supply chain dynamics, this paper deals with a system of conservation laws coupled with a system of ordinary differential equations. The former is defined on a domain with boundary and the coupling is provided by the boundary condition. For each of the examples considered, numerical integrations are provided

Effects of collisions on conservation laws in gyrokinetic field theory

Energy Technology Data Exchange (ETDEWEB)

Sugama, H.; Nunami, M. [National Institute for Fusion Science, Toki 509-5292 (Japan); Department of Fusion Science, SOKENDAI (The Graduate University for Advanced Studies), Toki 509-5292 (Japan); Watanabe, T.-H. [Department of Physics, Nagoya University, Nagoya 464-8602 (Japan)

2015-08-15

Effects of collisions on conservation laws for toroidal plasmas are investigated based on the gyrokinetic field theory. Associating the collisional system with a corresponding collisionless system at a given time such that the two systems have the same distribution functions and electromagnetic fields instantaneously, it is shown how the collisionless conservation laws derived from Noether's theorem are modified by the collision term. Effects of the external source term added into the gyrokinetic equation can be formulated similarly with the collisional effects. Particle, energy, and toroidal momentum balance equations including collisional and turbulent transport fluxes are systematically derived using a novel gyrokinetic collision operator, by which the collisional change rates of energy and canonical toroidal angular momentum per unit volume in the gyrocenter space can be given in the conservative forms. The ensemble-averaged transport equations of particles, energy, and toroidal momentum given in the present work are shown to include classical, neoclassical, and turbulent transport fluxes which agree with those derived from conventional recursive formulations.

Regularity and mass conservation for discrete coagulation–fragmentation equations with diffusion

KAUST Repository

Cañizo, J.A.

2010-03-01

We present a new a priori estimate for discrete coagulation-fragmentation systems with size-dependent diffusion within a bounded, regular domain confined by homogeneous Neumann boundary conditions. Following from a duality argument, this a priori estimate provides a global L2 bound on the mass density and was previously used, for instance, in the context of reaction-diffusion equations. In this paper we demonstrate two lines of applications for such an estimate: On the one hand, it enables to simplify parts of the known existence theory and allows to show existence of solutions for generalised models involving collision-induced, quadratic fragmentation terms for which the previous existence theory seems difficult to apply. On the other hand and most prominently, it proves mass conservation (and thus the absence of gelation) for almost all the coagulation coefficients for which mass conservation is known to hold true in the space homogeneous case. © 2009 Elsevier Masson SAS. All rights reserved.

On the multiple zeros of a real analytic function with applications to the averaging theory of differential equations

Science.gov (United States)

García, Isaac A.; Llibre, Jaume; Maza, Susanna

2018-06-01

In this work we consider real analytic functions , where , Ω is a bounded open subset of , is an interval containing the origin, are parameters, and ε is a small parameter. We study the branching of the zero-set of at multiple points when the parameter ε varies. We apply the obtained results to improve the classical averaging theory for computing T-periodic solutions of λ-families of analytic T-periodic ordinary differential equations defined on , using the displacement functions defined by these equations. We call the coefficients in the Taylor expansion of in powers of ε the averaged functions. The main contribution consists in analyzing the role that have the multiple zeros of the first non-zero averaged function. The outcome is that these multiple zeros can be of two different classes depending on whether the zeros belong or not to the analytic set defined by the real variety associated to the ideal generated by the averaged functions in the Noetheriang ring of all the real analytic functions at . We bound the maximum number of branches of isolated zeros that can bifurcate from each multiple zero z 0. Sometimes these bounds depend on the cardinalities of minimal bases of the former ideal. Several examples illustrate our results and they are compared with the classical theory, branching theory and also under the light of singularity theory of smooth maps. The examples range from polynomial vector fields to Abel differential equations and perturbed linear centers.

Uniqueness of Mass-Conserving Self-similar Solutions to Smoluchowski's Coagulation Equation with Inverse Power Law Kernels

Science.gov (United States)

Laurençot, Philippe

2018-03-01

Uniqueness of mass-conserving self-similar solutions to Smoluchowski's coagulation equation is shown when the coagulation kernel K is given by K(x,x_*)=2(x x_*)^{-α } , (x,x_*)\\in (0,∞)^2 , for some α >0.

A direct Primitive Variable Recovery Scheme for hyperbolic conservative equations: The case of relativistic hydrodynamics.

Science.gov (United States)

Aguayo-Ortiz, A; Mendoza, S; Olvera, D

2018-01-01

In this article we develop a Primitive Variable Recovery Scheme (PVRS) to solve any system of coupled differential conservative equations. This method obtains directly the primitive variables applying the chain rule to the time term of the conservative equations. With this, a traditional finite volume method for the flux is applied in order avoid violation of both, the entropy and "Rankine-Hugoniot" jump conditions. The time evolution is then computed using a forward finite difference scheme. This numerical technique evades the recovery of the primitive vector by solving an algebraic system of equations as it is often used and so, it generalises standard techniques to solve these kind of coupled systems. The article is presented bearing in mind special relativistic hydrodynamic numerical schemes with an added pedagogical view in the appendix section in order to easily comprehend the PVRS. We present the convergence of the method for standard shock-tube problems of special relativistic hydrodynamics and a graphical visualisation of the errors using the fluctuations of the numerical values with respect to exact analytic solutions. The PVRS circumvents the sometimes arduous computation that arises from standard numerical methods techniques, which obtain the desired primitive vector solution through an algebraic polynomial of the charges.

A conservative finite difference method for the numerical solution of plasma fluid equations

International Nuclear Information System (INIS)

Colella, P.; Dorr, M.R.; Wake, D.D.

1999-01-01

This paper describes a numerical method for the solution of a system of plasma fluid equations. The fluid model is similar to those employed in the simulation of high-density, low-pressure plasmas used in semiconductor processing. The governing equations consist of a drift-diffusion model of the electrons, together with an internal energy equation, coupled via Poisson's equation to a system of Euler equations for each ion species augmented with electrostatic force, collisional, and source/sink terms. The time integration of the full system is performed using an operator splitting that conserves space charge and avoids dielectric relaxation timestep restrictions. The integration of the individual ion species and electrons within the time-split advancement is achieved using a second-order Godunov discretization of the hyperbolic terms, modified to account for the significant role of the electric field in the propagation of acoustic waves, combined with a backward Euler discretization of the parabolic terms. Discrete boundary conditions are employed to accommodate the plasma sheath boundary layer on underresolved grids. The algorithm is described for the case of a single Cartesian grid as the first step toward an implementation on a locally refined grid hierarchy in which the method presented here may be applied on each refinement level

Scaling Relations and Self-Similarity of 3-Dimensional Reynolds-Averaged Navier-Stokes Equations.

Science.gov (United States)

Ercan, Ali; Kavvas, M Levent

2017-07-25

Scaling conditions to achieve self-similar solutions of 3-Dimensional (3D) Reynolds-Averaged Navier-Stokes Equations, as an initial and boundary value problem, are obtained by utilizing Lie Group of Point Scaling Transformations. By means of an open-source Navier-Stokes solver and the derived self-similarity conditions, we demonstrated self-similarity within the time variation of flow dynamics for a rigid-lid cavity problem under both up-scaled and down-scaled domains. The strength of the proposed approach lies in its ability to consider the underlying flow dynamics through not only from the governing equations under consideration but also from the initial and boundary conditions, hence allowing to obtain perfect self-similarity in different time and space scales. The proposed methodology can be a valuable tool in obtaining self-similar flow dynamics under preferred level of detail, which can be represented by initial and boundary value problems under specific assumptions.

Negatons, positons, rational-like solutions and conservation laws of the Korteweg-de Vries equation with loss and non-uniformity terms

International Nuclear Information System (INIS)

Zhang Dajun; Chen Dengyuan

2004-01-01

Solitons, negatons, positons, rational-like solutions and mixed solutions of a non-isospectral equation, the Korteweg-de Vries equation with loss and non-uniformity terms, are obtained through the Wronskian technique. The non-isospectral characteristics of the motion behaviours of some solutions are described with some figures made by using Mathematica. We also derive an infinite number of conservation laws of the equation

An Analysis of the Invariance and Conservation Laws of Some Classes of Nonlinear Ostrovsky Equations and Related Systems

International Nuclear Information System (INIS)

Fakhar, K.; Kara, A. H.

2011-01-01

A large class of partial differential equations in the modelling of ocean waves are due to Ostrovsky. We determine the invariance properties (through the Lie point symmetry generators) and construct classes of conservation laws for some of the models. In the latter case, the method involves finding the ‘multipliers’ associated with the conservation laws with a stronger emphasis on the ‘higher-order’ ones. The relationship between the symmetries and conservation laws is investigated by considering the invariance properties of the multipliers. (general)

Conservation Properties of the Hamiltonian Particle-Mesh method for the Quasi-Geostrophic Equations on a sphere

NARCIS (Netherlands)

H. Thorsdottir (Halldora)

2011-01-01

htmlabstractThe Hamiltonian particle-mesh (HPM) method is used to solve the Quasi-Geostrophic model generalized to a sphere, using the Spherepack modeling package to solve the Helmholtz equation on a colatitude-longitude grid with spherical harmonics. The predicted energy conservation of a

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.

Numerical solution of the thermalhydraulic conservation equations from fundamental concepts to multidimensional two-fluid analysis

International Nuclear Information System (INIS)

Carver, M.B.

1995-08-01

The discussion briefly establishes some requisite concepts of differential equation theory, and applies these to describe methods for numerical solution of the thermalhydraulic conservation equations in their various forms. The intent is to cover the general methodology without obscuring the principles with details. As a short overview of computational thermalhydraulics, the material provides an introductory foundation, so that those working on the application of thermalhydraulic codes can begin to understand the many intricacies involved without having to locate and read the references given. Those intending to work in code development will need to read and understand all the references. (author). 49 refs

Periodic solutions of Lienard differential equations via averaging theory of order two.

Science.gov (United States)

Llibre, Jaume; Novaes, Douglas D; Teixeira, Marco A

2015-01-01

For ε ≠ 0 sufficiently small we provide sufficient conditions for the existence of periodic solutions for the Lienard differential equations of the form x'' + f ⁢(x)⁢ x' + n2⁢x + g (x) = ε2p1 ⁢(t) + ε3 ⁢p2(t), where n is a positive integer, f : ℝ → ℝ is a C 3 function, g : ℝ → ℝ is a C 4 function, and p i : ℝ → ℝ for i = 1, 2 are continuous 2π-periodic function. The main tool used in this paper is the averaging theory of second order. We also provide one application of the main result obtained.

Periodic solutions of Lienard differential equations via averaging theory of order two

Directory of Open Access Journals (Sweden)

JAUME LLIBRE

2015-12-01

Full Text Available Abstract For ε ≠ 0sufficiently small we provide sufficient conditions for the existence of periodic solutions for the Lienard differential equations of the form x ′′ + f ( x x ′ + n 2 x + g ( x = ε 2 p 1 ( t + ε 3 p 2 ( t , where n is a positive integer, f : ℝ → ℝis a C 3function, g : ℝ → ℝis a C 4function, and p i : ℝ → ℝfor i = 1 , 2are continuous 2 π–periodic function. The main tool used in this paper is the averaging theory of second order. We also provide one application of the main result obtained.

Higher order supersymmetries and fermionic conservation laws of the supersymmetric extension of the KdV equation

NARCIS (Netherlands)

Kersten, P.H.M.

1988-01-01

By the introduction of nonlocal basonic and fermionic variables we construct a recursion symmetry of the super KdV equation, leading to a hierarchy of bosonic symmetries and one of fermionic symmetries. The hierarchies of bosonic and fermionic conservation laws arise in a natural way in the

Tensor formulation of the model equations on strong conservation form for an incompressible flow in general coordinates

DEFF Research Database (Denmark)

Jørgensen, Bo Hoffmann

2003-01-01

This brief report expresses the basic equations of an incompressible flow model in a form which can be translated easily into the form used by a numerical solver. The application of tensor notation makes is possible to effectively address the issue ofnumerical robustness and stating the model...... equations on a general form which accommodate curvilinear coordinates. Strong conservation form is obtained by formulating the equations so that the flow variables, velocity and pressure, are expressed in thephysical coordinate system while the location of evaluation is expressed within the transformed...... form of the equations is included which allows for special solutions to be developed in the transformedcoordinate system. Examples of applications are atmospheric flows over complex terrain, aerodynamically flows, industrial flows and environmental flows....

Averaging for solitons with nonlinearity management

International Nuclear Information System (INIS)

Pelinovsky, D.E.; Kevrekidis, P.G.; Frantzeskakis, D.J.

2003-01-01

We develop an averaging method for solitons of the nonlinear Schroedinger equation with a periodically varying nonlinearity coefficient, which is used to effectively describe solitons in Bose-Einstein condensates, in the context of the recently proposed technique of Feshbach resonance management. Using the derived local averaged equation, we study matter-wave bright and dark solitons and demonstrate a very good agreement between solutions of the averaged and full equations

Reduced Braginskii equations

International Nuclear Information System (INIS)

Yagi, M.; Horton, W.

1994-01-01

A set of reduced Braginskii equations is derived without assuming flute ordering and the Boussinesq approximation. These model equations conserve the physical energy. It is crucial at finite β that the perpendicular component of Ohm's law be solved to ensure ∇·j=0 for energy conservation

Stochastic line motion and stochastic flux conservation for nonideal hydromagnetic models

International Nuclear Information System (INIS)

Eyink, Gregory L.

2009-01-01

We prove that smooth solutions of nonideal (viscous and resistive) incompressible magnetohydrodynamic (MHD) equations satisfy a stochastic law of flux conservation. This property implies that the magnetic flux through a surface is equal to the average of the magnetic fluxes through an ensemble of surfaces advected backward in time by the plasma velocity perturbed with a random white noise. Our result is an analog of the well-known Alfven theorem of ideal MHD and is valid for any value of the magnetic Prandtl number. A second stochastic conservation law is shown to hold at unit Prandtl number, a random version of the generalized Kelvin theorem derived by Bekenstein and Oron for ideal MHD. These stochastic conservation laws are not only shown to be consequences of the nonideal MHD equations but are proved in fact to be equivalent to those equations. We derive similar results for two more refined hydromagnetic models, Hall MHD and the two-fluid plasma model, still assuming incompressible velocities and isotropic transport coefficients. Finally, we use these results to discuss briefly the infinite-Reynolds-number limit of hydromagnetic turbulence and to support the conjecture that flux conservation remains stochastic in that limit.

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

Science.gov (United States)

Motsepa, Tanki; Masood Khalique, Chaudry

2018-05-01

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

Reduced Braginskii equations

International Nuclear Information System (INIS)

Yagi, M.; Horton, W.

1993-11-01

A set of reduced Braginskii equations is derived without assuming flute ordering and the Boussinesq approximation. These model equations conserve the physical energy. It is crucial at finite β that we solve the perpendicular component of Ohm's law to conserve the physical energy while ensuring the relation ∇ · j = 0

Reduced Braginskii equations

Energy Technology Data Exchange (ETDEWEB)

Yagi, M. [Japan Atomic Energy Research Inst., Naka, Ibaraki (Japan). Naka Fusion Research Establishment; Horton, W. [Texas Univ., Austin, TX (United States). Inst. for Fusion Studies

1993-11-01

A set of reduced Braginskii equations is derived without assuming flute ordering and the Boussinesq approximation. These model equations conserve the physical energy. It is crucial at finite {beta} that we solve the perpendicular component of Ohm`s law to conserve the physical energy while ensuring the relation {del} {center_dot} j = 0.

Boltzmann equation for a mixture of gases with non-conservative processes

International Nuclear Information System (INIS)

Martiarena, M.L.

1989-01-01

The nonlinear and non-isotropic Boltzmann equation (NLBE) including several molecular species, non-conservative channels and external forces. The general solution of that equation is obtained for a spatially homogeneous mixture of L gases, consisting of Maxwell particles, as a Generalized Laguerre expansion, within a Hilbert space. Removal and self-generation effects are included in presence of a time-dependent external force. An exact particular solution is studied generalizing the well-known BKW-mode for a mixture of L gases with inelastic processes. An homogeneous gas of test particles, in d dimension, is considered which interacts with a background host medium in the presence of an external space and time dependent force. Scattering, removal and self-generation collisions are included. The inhomogeneous Boltzmann equation for this system to an homogeneous one is reduced without background or external forces, using a generalized Nilkoskii transform. It is shown that a background of field particles can confine the test gas, even in absence of external forces. Furthermore, the solution of NLBE with non-isotropic singular initial conditions, is analyzed. The NLBE is transformed into an integral equation which is solved iteratively. The evolution of delta and step singularities in the distribution function is discussed during the initial layer and compared with the isotropic case. As an application of the methods abovementioned, the collision of a beam of ions or neutral atoms with a carbon-foil is considered. The electron experimental spectra from a transport equation is described. It is supposed that convoy electron may be produced inside the solid by single ion-atom collisions as ELC or ECC. The produced electrons lost energy by collision with the atoms of the material, which are considered at rest. The electron distribution function is numerically calculated. The ratio between the intrinsic convoy electron peak height to the background electron intensity

Conservation Laws for Gyrokinetic Equations for Large Perturbations and Flows

Science.gov (United States)

Dimits, Andris

2017-10-01

Gyrokinetic theory has proved to be very useful for the understanding of magnetized plasmas, both to simplify analytical treatments and as a basis for efficient numerical simulations. Gyrokinetic theories were previously developed in two extended orderings that are applicable to large fluctuations and flows as may arise in the tokamak edge and scrapeoff layer. In the present work, we cast the resulting equations in a field-theoretical variational form, and derive, up to second order in the respective orderings, the associated global and local energy and (linear and toroidal) momentum conservation relations that result from Noether's theorem. The consequences of these for the various possible choices of numerical discretization used in gyrokinetic simulations are considered. Prepared for US DOE by LLNL under Contract DE-AC52-07NA27344 and supported by the U.S. DOE, OFES.

Total energy and potential enstrophy conserving schemes for the shallow water equations using Hamiltonian methods - Part 1: Derivation and properties

Science.gov (United States)

Eldred, Christopher; Randall, David

2017-02-01

The shallow water equations provide a useful analogue of the fully compressible Euler equations since they have similar characteristics: conservation laws, inertia-gravity and Rossby waves, and a (quasi-) balanced state. In order to obtain realistic simulation results, it is desirable that numerical models have discrete analogues of these properties. Two prototypical examples of such schemes are the 1981 Arakawa and Lamb (AL81) C-grid total energy and potential enstrophy conserving scheme, and the 2007 Salmon (S07) Z-grid total energy and potential enstrophy conserving scheme. Unfortunately, the AL81 scheme is restricted to logically square, orthogonal grids, and the S07 scheme is restricted to uniform square grids. The current work extends the AL81 scheme to arbitrary non-orthogonal polygonal grids and the S07 scheme to arbitrary orthogonal spherical polygonal grids in a manner that allows for both total energy and potential enstrophy conservation, by combining Hamiltonian methods (work done by Salmon, Gassmann, Dubos, and others) and discrete exterior calculus (Thuburn, Cotter, Dubos, Ringler, Skamarock, Klemp, and others). Detailed results of the schemes applied to standard test cases are deferred to part 2 of this series of papers.

Integrating factors and conservation theorems for Hamilton‘s canonical equations of motion of variable mass nonholonmic nonconservative dynamical systems

Institute of Scientific and Technical Information of China (English)

李仁杰; 刘洋; 等

2002-01-01

We present a general approach to the construction of conservation laws for variable mass noholonmic nonconservative systems.First,we give the definition of integrating factors,and we study in detail the necessary conditions for the existence of the conserved quantities,Then,we establish the conservatioin theorem and its inverse theorem for Hamilton's canonical equations of motion of variable mass nonholonomic nonocnservative dynamical systems.Finally,we give an example to illustrate the application of the results.

equateIRT: An R Package for IRT Test Equating

Directory of Open Access Journals (Sweden)

Michela Battauz

2015-12-01

Full Text Available The R package equateIRT implements item response theory (IRT methods for equating different forms composed of dichotomous items. In particular, the IRT models included are the three-parameter logistic model, the two-parameter logistic model, the one-parameter logistic model and the Rasch model. Forms can be equated when they present common items (direct equating or when they can be linked through a chain of forms that present common items in pairs (indirect or chain equating. When two forms can be equated through different paths, a single conversion can be obtained by averaging the equating coefficients. The package calculates direct and chain equating coefficients. The averaging of direct and chain coefficients that link the same two forms is performed through the bisector method. Furthermore, the package provides analytic standard errors of direct, chain and average equating coefficients.

Development of equations to predict the influence of floor space on average daily gain, average daily feed intake and gain : feed ratio of finishing pigs.

Science.gov (United States)

Flohr, J R; Dritz, S S; Tokach, M D; Woodworth, J C; DeRouchey, J M; Goodband, R D

2018-05-01

Floor space allowance for pigs has substantial effects on pig growth and welfare. Data from 30 papers examining the influence of floor space allowance on the growth of finishing pigs was used in a meta-analysis to develop alternative prediction equations for average daily gain (ADG), average daily feed intake (ADFI) and gain : feed ratio (G : F). Treatment means were compiled in a database that contained 30 papers for ADG and 28 papers for ADFI and G : F. The predictor variables evaluated were floor space (m2/pig), k (floor space/final BW0.67), Initial BW, Final BW, feed space (pigs per feeder hole), water space (pigs per waterer), group size (pigs per pen), gender, floor type and study length (d). Multivariable general linear mixed model regression equations were used. Floor space treatments within each experiment were the observational and experimental unit. The optimum equations to predict ADG, ADFI and G : F were: ADG, g=337.57+(16 468×k)-(237 350×k 2)-(3.1209×initial BW (kg))+(2.569×final BW (kg))+(71.6918×k×initial BW (kg)); ADFI, g=833.41+(24 785×k)-(388 998×k 2)-(3.0027×initial BW (kg))+(11.246×final BW (kg))+(187.61×k×initial BW (kg)); G : F=predicted ADG/predicted ADFI. Overall, the meta-analysis indicates that BW is an important predictor of ADG and ADFI even after computing the constant coefficient k, which utilizes final BW in its calculation. This suggests including initial and final BW improves the prediction over using k as a predictor alone. In addition, the analysis also indicated that G : F of finishing pigs is influenced by floor space allowance, whereas individual studies have concluded variable results.

Well-balanced schemes for the Euler equations with gravitation: Conservative formulation using global fluxes

Science.gov (United States)

Chertock, Alina; Cui, Shumo; Kurganov, Alexander; Özcan, Şeyma Nur; Tadmor, Eitan

2018-04-01

We develop a second-order well-balanced central-upwind scheme for the compressible Euler equations with gravitational source term. Here, we advocate a new paradigm based on a purely conservative reformulation of the equations using global fluxes. The proposed scheme is capable of exactly preserving steady-state solutions expressed in terms of a nonlocal equilibrium variable. A crucial step in the construction of the second-order scheme is a well-balanced piecewise linear reconstruction of equilibrium variables combined with a well-balanced central-upwind evolution in time, which is adapted to reduce the amount of numerical viscosity when the flow is at (near) steady-state regime. We show the performance of our newly developed central-upwind scheme and demonstrate importance of perfect balance between the fluxes and gravitational forces in a series of one- and two-dimensional examples.

Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations

International Nuclear Information System (INIS)

Rhebergen, S.; Bokhove, O.; Vegt, J.J.W. van der

2008-01-01

We present space- and space-time discontinuous Galerkin finite element (DGFEM) formulations for systems containing nonconservative products, such as occur in dispersed multiphase flow equations. The main criterium we pose on the weak formulation is that if the system of nonconservative partial differential equations can be transformed into conservative form, then the formulation must reduce to that for conservative systems. Standard DGFEM formulations cannot be applied to nonconservative systems of partial differential equations. We therefore introduce the theory of weak solutions for nonconservative products into the DGFEM formulation leading to the new question how to define the path connecting left and right states across a discontinuity. The effect of different paths on the numerical solution is investigated and found to be small. We also introduce a new numerical flux that is able to deal with nonconservative products. Our scheme is applied to two different systems of partial differential equations. First, we consider the shallow water equations, where topography leads to nonconservative products, in which the known, possibly discontinuous, topography is formally taken as an unknown in the system. Second, we consider a simplification of a depth-averaged two-phase flow model which contains more intrinsic nonconservative products

Amplitude equation for under water sand-ripples in one dimension

DEFF Research Database (Denmark)

observed when the amplitude $d$ is suddenly varied. The equation has the form h_t=- ε(h-mean(h))+((h_x)^2-1)h_(xx)- h_(xxxx)+ δ((h_x)^2)_(xx) which, due to the first term, is neither completely local (it has long-range coupling through the average height mean(h)) nor has local sand conservation. We argue...

Consistent three-equation model for thin films

Science.gov (United States)

Richard, Gael; Gisclon, Marguerite; Ruyer-Quil, Christian; Vila, Jean-Paul

2017-11-01

Numerical simulations of thin films of newtonian fluids down an inclined plane use reduced models for computational cost reasons. These models are usually derived by averaging over the fluid depth the physical equations of fluid mechanics with an asymptotic method in the long-wave limit. Two-equation models are based on the mass conservation equation and either on the momentum balance equation or on the work-energy theorem. We show that there is no two-equation model that is both consistent and theoretically coherent and that a third variable and a three-equation model are required to solve all theoretical contradictions. The linear and nonlinear properties of two and three-equation models are tested on various practical problems. We present a new consistent three-equation model with a simple mathematical structure which allows an easy and reliable numerical resolution. The numerical calculations agree fairly well with experimental measurements or with direct numerical resolutions for neutral stability curves, speed of kinematic waves and of solitary waves and depth profiles of wavy films. The model can also predict the flow reversal at the first capillary trough ahead of the main wave hump.

Space, time and conservation laws

International Nuclear Information System (INIS)

Aronov, R.A.; Ugarov, V.A.

1978-01-01

The Neter theorem establishing correspondence between conservation laws and symmetry properties (space and time in particular) is considered. The theorem is based on one of the possible ways of finding equations of motion for a physical system. From a certain expression (action functional) equations of motion for a system can be obtained which do not contain new physical assertions in principal in comparison with the Newtonian laws. Neter suggested a way of deriving conservation laws by transforming space and time coordinates. Neter theorem consequences raise a number of problems: 1). Are conservation laws (energy, momentum) consequences of space and time symmetry properties. 2). Is it possible to obtain conservation laws in theory neglecting equations of motion. 3). What is of the primary importance: equations of motion, conservation laws or properties of space and time symmetry. It is shown that direct Neter theorem does not testify to stipulation of conservation laws by properties of space and time symmetry and symmetry properties of other non-space -time properties of material systems in objective reality. It says nothing of whether there is any subordination between symmetry properties and conservation laws

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

International Nuclear Information System (INIS)

Li Xinyue; Zhao Qiulan

2009-01-01

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

A higher-order conservation element solution element method for solving hyperbolic differential equations on unstructured meshes

Science.gov (United States)

Bilyeu, David

This dissertation presents an extension of the Conservation Element Solution Element (CESE) method from second- to higher-order accuracy. The new method retains the favorable characteristics of the original second-order CESE scheme, including (i) the use of the space-time integral equation for conservation laws, (ii) a compact mesh stencil, (iii) the scheme will remain stable up to a CFL number of unity, (iv) a fully explicit, time-marching integration scheme, (v) true multidimensionality without using directional splitting, and (vi) the ability to handle two- and three-dimensional geometries by using unstructured meshes. This algorithm has been thoroughly tested in one, two and three spatial dimensions and has been shown to obtain the desired order of accuracy for solving both linear and non-linear hyperbolic partial differential equations. The scheme has also shown its ability to accurately resolve discontinuities in the solutions. Higher order unstructured methods such as the Discontinuous Galerkin (DG) method and the Spectral Volume (SV) methods have been developed for one-, two- and three-dimensional application. Although these schemes have seen extensive development and use, certain drawbacks of these methods have been well documented. For example, the explicit versions of these two methods have very stringent stability criteria. This stability criteria requires that the time step be reduced as the order of the solver increases, for a given simulation on a given mesh. The research presented in this dissertation builds upon the work of Chang, who developed a fourth-order CESE scheme to solve a scalar one-dimensional hyperbolic partial differential equation. The completed research has resulted in two key deliverables. The first is a detailed derivation of a high-order CESE methods on unstructured meshes for solving the conservation laws in two- and three-dimensional spaces. The second is the code implementation of these numerical methods in a computer code. For

A new nonlinear turbulence model based on Partially-Averaged Navier-Stokes Equations

International Nuclear Information System (INIS)

Liu, J T; Wu, Y L; Cai, C; Liu, S H; Wang, L Q

2013-01-01

Partially-averaged Navier-Stokes (PANS) Model was recognized as a Reynolds-averaged Navier-Stokes (RANS) to direct numerical simulation (DNS) bridging method. PANS model was purported for any filter width-from RANS to DNS. PANS method also shared some similarities with the currently popular URANS (unsteady RANS) method. In this paper, a new PANS model was proposed, which was based on RNG k-ε turbulence model. The Standard and RNG k-ε turbulence model were both isotropic models, as well as PANS models. The sheer stress in those PANS models was solved by linear equation. The linear hypothesis was not accurate in the simulation of complex flow, such as stall phenomenon. The sheer stress here was solved by nonlinear method proposed by Ehrhard. Then, the nonlinear PANS model was set up. The pressure coefficient of the suction side of the NACA0015 hydrofoil was predicted. The result of pressure coefficient agrees well with experimental result, which proves that the nonlinear PANS model can capture the high pressure gradient flow. A low specific centrifugal pump was used to verify the capacity of the nonlinear PANS model. The comparison between the simulation results of the centrifugal pump and Particle Image Velocimetry (PIV) results proves that the nonlinear PANS model can be used in the prediction of complex flow field

How to obtain the covariant form of Maxwell's equations from the continuity equation

International Nuclear Information System (INIS)

Heras, Jose A

2009-01-01

The covariant Maxwell equations are derived from the continuity equation for the electric charge. This result provides an axiomatic approach to Maxwell's equations in which charge conservation is emphasized as the fundamental axiom underlying these equations

Angular finite volume method for solving the multigroup transport equation with piecewise average scattering cross sections

International Nuclear Information System (INIS)

Calloo, A.; Vidal, J.F.; Le Tellier, R.; Rimpault, G.

2011-01-01

This paper deals with the solving of the multigroup integro-differential form of the transport equation for fine energy group structure. In that case, multigroup transfer cross sections display strongly peaked shape for light scatterers and the current Legendre polynomial expansion is not well-suited to represent them. Furthermore, even if considering an exact scattering cross sections representation, the scattering source in the discrete ordinates method (also known as the Sn method) being calculated by sampling the angular flux at given directions, may be wrongly computed due to lack of angular support for the angular flux. Hence, following the work of Gerts and Matthews, an angular finite volume solver has been developed for 2D Cartesian geometries. It integrates the multigroup transport equation over discrete volume elements obtained by meshing the unit sphere with a product grid over the polar and azimuthal coordinates and by considering the integrated flux per solid angle element. The convergence of this method has been compared to the S_n method for a highly anisotropic benchmark. Besides, piecewise-average scattering cross sections have been produced for non-bound Hydrogen atoms using a free gas model for thermal neutrons. LWR lattice calculations comparing Legendre representations of the Hydrogen scattering multigroup cross section at various orders and piecewise-average cross sections for this same atom are carried out (while keeping a Legendre representation for all other isotopes). (author)

An upscaled two-equation model of transport in porous media through unsteady-state closure of volume averaged formulations

Science.gov (United States)

Chaynikov, S.; Porta, G.; Riva, M.; Guadagnini, A.

2012-04-01

We focus on a theoretical analysis of nonreactive solute transport in porous media through the volume averaging technique. Darcy-scale transport models based on continuum formulations typically include large scale dispersive processes which are embedded in a pore-scale advection diffusion equation through a Fickian analogy. This formulation has been extensively questioned in the literature due to its inability to depict observed solute breakthrough curves in diverse settings, ranging from the laboratory to the field scales. The heterogeneity of the pore-scale velocity field is one of the key sources of uncertainties giving rise to anomalous (non-Fickian) dispersion in macro-scale porous systems. Some of the models which are employed to interpret observed non-Fickian solute behavior make use of a continuum formulation of the porous system which assumes a two-region description and includes a bimodal velocity distribution. A first class of these models comprises the so-called ''mobile-immobile'' conceptualization, where convective and dispersive transport mechanisms are considered to dominate within a high velocity region (mobile zone), while convective effects are neglected in a low velocity region (immobile zone). The mass exchange between these two regions is assumed to be controlled by a diffusive process and is macroscopically described by a first-order kinetic. An extension of these ideas is the two equation ''mobile-mobile'' model, where both transport mechanisms are taken into account in each region and a first-order mass exchange between regions is employed. Here, we provide an analytical derivation of two region "mobile-mobile" meso-scale models through a rigorous upscaling of the pore-scale advection diffusion equation. Among the available upscaling methodologies, we employ the Volume Averaging technique. In this approach, the heterogeneous porous medium is supposed to be pseudo-periodic, and can be represented through a (spatially) periodic unit cell

Weak self-adjoint differential equations

International Nuclear Information System (INIS)

Gandarias, M L

2011-01-01

The concepts of self-adjoint and quasi self-adjoint equations were introduced by Ibragimov (2006 J. Math. Anal. Appl. 318 742-57; 2007 Arch. ALGA 4 55-60). In Ibragimov (2007 J. Math. Anal. Appl. 333 311-28), a general theorem on conservation laws was proved. In this paper, we generalize the concept of self-adjoint and quasi self-adjoint equations by introducing the definition of weak self-adjoint equations. We find a class of weak self-adjoint quasi-linear parabolic equations. The property of a differential equation to be weak self-adjoint is important for constructing conservation laws associated with symmetries of the differential equation. (fast track communication)

Theoretical and pragmatic modelling of governing equations for a two-phase flow in bubbly and annular flow regimes

International Nuclear Information System (INIS)

Bottoni, M.; Sengpiel, W.

1992-01-01

Starting from the rigorous formulation of the conservation equations for mass, momentum and enthalpy, derived for a two-phase flow by volume averaging microscopic balance equations over Eulerian control cells, the article discusses the formulation of the terms describing exchanges between the phases. Two flow regimes are taken into consideration, bubbly flow, applicable for small or medium void fractions, and annular flow, for large void fractions. When lack of knowledge of volume-averaged physical quantities make the rigorously formulated terms useless for computational purposes, modelling of these terms is discussed. 3 figs., 15 refs

Efficient positive, conservative, Maxwellian preserving and implicit difference schemes for the 1-D isotropic Fokker-Planck-Landau equation; Schemas positifs, implicites, conservant l'energie et les etats d'equilibre pour l'equation de Fokker-Planck-Landau isotrope

Energy Technology Data Exchange (ETDEWEB)

Buet, Ch. [CEA Bruyeres-le-Chatel, Dept. Sciences de la Simulation et de l' Information, Service Numerique Environnement et Constantes, 91 (France); Le Thanh, K.C. [CEA Bruyeres-le-Chatel, Service Physique des Plasmas et Electromagnetisme, 91 (France). Dept. de Physique Theorique et Appliquee

2008-07-01

The aim of this paper is to describe the discretization of the Fokker-Planck-Landau (FPL) collision term in the isotropic case, which models the self-collision for the electrons when they are totally isotropized by heavy particles background such as ions. The discussion focuses on schemes, which could preserve positivity, mass, energy and Maxwellian equilibrium. The Chang and Cooper method is widely used by plasma's physicists for the FPL equation (and for Fokker-Planck type equations). We present a new variant that is both positive and conservative contrary to the existing one's. We propose also a non Chang and Cooper 'type scheme on non-uniform grid, which is also both positive, conservative and equilibrium state preserving contrary to existing one's. The case of Coulombian potential is emphasized. We address also the problem of the time discretization. In particular we show how to recast some implicit methods to get band diagonal system and to solve it by direct method with a linear cost. (authors)

Notes on Conservation Laws, Equations of Motion of Matter, and Particle Fields in Lorentzian and Teleparallel de Sitter Space-Time Structures

Directory of Open Access Journals (Sweden)

Waldyr A. Rodrigues

2016-01-01

Full Text Available We discuss the physics of interacting fields and particles living in a de Sitter Lorentzian manifold (dSLM, a submanifold of a 5-dimensional pseudo-Euclidean (5dPE equipped with a metric tensor inherited from the metric of the 5dPE space. The dSLM is naturally oriented and time oriented and is the arena used to study the energy-momentum conservation law and equations of motion for physical systems living there. Two distinct de Sitter space-time structures MdSL and MdSTP are introduced given dSLM, the first equipped with the Levi-Civita connection of its metric field and the second with a metric compatible parallel connection. Both connections are used only as mathematical devices. Thus, for example, MdSL is not supposed to be the model of any gravitational field in the General Relativity Theory (GRT. Misconceptions appearing in the literature concerning the motion of free particles in dSLM are clarified. Komar currents are introduced within Clifford bundle formalism permitting the presentation of Einstein equation as a Maxwell like equation and proving that in GRT there are infinitely many conserved currents. We prove that in GRT even when the appropriate Killing vector fields exist it is not possible to define a conserved energy-momentum covector as in special relativistic theories.

New formulation of Hardin-Pope equations for aeroacoustics

DEFF Research Database (Denmark)

Ekaterinaris, J.A.

1999-01-01

Dynamics, Vol. 6, No. 5-6, 1994, pp. 334-340). This method requires detailed information about the unsteady aerodynamic flowfield, which usually is obtained from a computational fluid dynamics solution. A new, conservative formulation of the equations governing acoustic disturbances is presented....... The conservative form of the governing equations is obtained after application of a transformation of variables that produces a set of inhomogeneous equations similar to the conservation-law form of the compressible Euler equations. The source term of these equations depends only on the derivatives...... of the hydrodynamic variables. Explicit time marching is performed. A high-order accurate, upwind-biased numerical scheme is used for numerical solution of the conservative equations. The convective fluxes are evaluated using upwind-biased formulas and flux-vector splitting. Solutions are obtained for the acoustic...

Equations for the gravitational field and local conserved quantities in the general theory of relativity

International Nuclear Information System (INIS)

Manoff, S.

1979-07-01

By utilization of the method of Lagrangians with covariant derivatives (MLCD) the different energy-momentum tensors (canonical, generalized canonical, symmetrical) and the relations between them are considered. On this basis, Einstein's theory of gravitation is studied as a field theory with a Lagrangian density of the type Lsub(g)=√-g.Lsub(g)(gsub(ij),Rsub(A)), (Rsub(A)=Rsub(ijkl)). It is shown that the energy-momentum tensors of the gravitational field can be defined for this theory. The symmetrical energy-momentum tensor of the gravitational field sub(gs)Tsub(k)sup(i), which in the general case is not a local conserved quantity (sub(gs)Tsub(k)sup(i)sub(;i) unequal 0) (in contrast to the material fields satisfying condition sub(Ms)Tsub(k)sup(i)sub(;i) = 0), is equal to zero for the gravitational field in vacuum (cosmological constant Λ = 0). Equations of the gravitational field of a new type are suggested, leading to equations of motion (sub(Ms)Tsub(k)sup(i) + sub(gs)Tsub(k)sup(i))sub(;i) = 0. The equations corresponding to the Lagrangian density Lsub(g)=(√-g/kappasub(o)) (R - lambda approximately), lambda approximately = const., are considered. The equations of Einstein Rsub(ij) = 0 are obtained in the case of gravitational field in vacuum. Some particular cases are examined as an illustration to material fields and the corresponding gravitational equations. (author)

Multi-dimensional, fully-implicit, spectral method for the Vlasov-Maxwell equations with exact conservation laws in discrete form

Science.gov (United States)

Delzanno, G. L.

2015-11-01

A spectral method for the numerical solution of the multi-dimensional Vlasov-Maxwell equations is presented. The plasma distribution function is expanded in Fourier (for the spatial part) and Hermite (for the velocity part) basis functions, leading to a truncated system of ordinary differential equations for the expansion coefficients (moments) that is discretized with an implicit, second order accurate Crank-Nicolson time discretization. The discrete non-linear system is solved with a preconditioned Jacobian-Free Newton-Krylov method. It is shown analytically that the Fourier-Hermite method features exact conservation laws for total mass, momentum and energy in discrete form. Standard tests involving plasma waves and the whistler instability confirm the validity of the conservation laws numerically. The whistler instability test also shows that we can step over the fastest time scale in the system without incurring in numerical instabilities. Some preconditioning strategies are presented, showing that the number of linear iterations of the Krylov solver can be drastically reduced and a significant gain in performance can be obtained.

Multi-symplectic Runge-Kutta methods for nonlinear Dirac equations

International Nuclear Information System (INIS)

Hong Jialin; Li Chun

2006-01-01

In this paper, we consider the multi-symplectic Runge-Kutta (MSRK) methods applied to the nonlinear Dirac equation in relativistic quantum physics, based on a discovery of the multi-symplecticity of the equation. In particular, the conservation of energy, momentum and charge under MSRK discretizations is investigated by means of numerical experiments and numerical comparisons with non-MSRK methods. Numerical experiments presented reveal that MSRK methods applied to the nonlinear Dirac equation preserve exactly conservation laws of charge and momentum, and conserve the energy conservation in the corresponding numerical accuracy to the method utilized. It is verified numerically that MSRK methods are stable and convergent with respect to the conservation laws of energy, momentum and charge, and MSRK methods preserve not only the inner geometric structure of the equation, but also some crucial conservative properties in quantum physics. A remarkable advantage of MSRK methods applied to the nonlinear Dirac equation is the precise preservation of charge conservation law

Angular finite volume method for solving the multigroup transport equation with piecewise average scattering cross sections

Energy Technology Data Exchange (ETDEWEB)

Calloo, A.; Vidal, J.F.; Le Tellier, R.; Rimpault, G., E-mail: ansar.calloo@cea.fr, E-mail: jean-francois.vidal@cea.fr, E-mail: romain.le-tellier@cea.fr, E-mail: gerald.rimpault@cea.fr [CEA, DEN, DER/SPRC/LEPh, Saint-Paul-lez-Durance (France)

2011-07-01

This paper deals with the solving of the multigroup integro-differential form of the transport equation for fine energy group structure. In that case, multigroup transfer cross sections display strongly peaked shape for light scatterers and the current Legendre polynomial expansion is not well-suited to represent them. Furthermore, even if considering an exact scattering cross sections representation, the scattering source in the discrete ordinates method (also known as the Sn method) being calculated by sampling the angular flux at given directions, may be wrongly computed due to lack of angular support for the angular flux. Hence, following the work of Gerts and Matthews, an angular finite volume solver has been developed for 2D Cartesian geometries. It integrates the multigroup transport equation over discrete volume elements obtained by meshing the unit sphere with a product grid over the polar and azimuthal coordinates and by considering the integrated flux per solid angle element. The convergence of this method has been compared to the S{sub n} method for a highly anisotropic benchmark. Besides, piecewise-average scattering cross sections have been produced for non-bound Hydrogen atoms using a free gas model for thermal neutrons. LWR lattice calculations comparing Legendre representations of the Hydrogen scattering multigroup cross section at various orders and piecewise-average cross sections for this same atom are carried out (while keeping a Legendre representation for all other isotopes). (author)

Developing A New Predictive Dispersion Equation Based on Tidal Average (TA) Condition in Alluvial Estuaries

Science.gov (United States)

Anak Gisen, Jacqueline Isabella; Nijzink, Remko C.; Savenije, Hubert H. G.

2014-05-01

Dispersion mathematical representation of tidal mixing between sea water and fresh water in The definition of dispersion somehow remains unclear as it is not directly measurable. The role of dispersion is only meaningful if it is related to the appropriate temporal and spatial scale of mixing, which are identified as the tidal period, tidal excursion (longitudinal), width of estuary (lateral) and mixing depth (vertical). Moreover, the mixing pattern determines the salt intrusion length in an estuary. If a physically based description of the dispersion is defined, this would allow the analytical solution of the salt intrusion problem. The objective of this study is to develop a predictive equation for estimating the dispersion coefficient at tidal average (TA) condition, which can be applied in the salt intrusion model to predict the salinity profile for any estuary during different events. Utilizing available data of 72 measurements in 27 estuaries (including 6 recently studied estuaries in Malaysia), regressions analysis has been performed with various combinations of dimensionless parameters . The predictive dispersion equations have been developed for two different locations, at the mouth D0TA and at the inflection point D1TA (where the convergence length changes). Regressions have been carried out with two separated datasets: 1) more reliable data for calibration; and 2) less reliable data for validation. The combination of dimensionless ratios that give the best performance is selected as the final outcome which indicates that the dispersion coefficient is depending on the tidal excursion, tidal range, tidal velocity amplitude, friction and the Richardson Number. A limitation of the newly developed equation is that the friction is generally unknown. In order to compensate this problem, further analysis has been performed adopting the hydraulic model of Cai et. al. (2012) to estimate the friction and depth. Keywords: dispersion, alluvial estuaries, mixing, salt

How to obtain the covariant form of Maxwell's equations from the continuity equation

Energy Technology Data Exchange (ETDEWEB)

Heras, Jose A [Departamento de Ciencias Basicas, Universidad Autonoma Metropolitana, Unidad Azcapotzalco, Av. San Pablo No. 180, Col. Reynosa, 02200, Mexico D. F. (Mexico); Departamento de Fisica y Matematicas, Universidad Iberoamericana, Prolongacion Paseo de la Reforma 880, Mexico D. F. 01210 (Mexico)

2009-07-15

The covariant Maxwell equations are derived from the continuity equation for the electric charge. This result provides an axiomatic approach to Maxwell's equations in which charge conservation is emphasized as the fundamental axiom underlying these equations.

Efficient conservative ADER schemes based on WENO reconstruction and space-time predictor in primitive variables

Science.gov (United States)

Zanotti, Olindo; Dumbser, Michael

2016-01-01

We present a new version of conservative ADER-WENO finite volume schemes, in which both the high order spatial reconstruction as well as the time evolution of the reconstruction polynomials in the local space-time predictor stage are performed in primitive variables, rather than in conserved ones. To obtain a conservative method, the underlying finite volume scheme is still written in terms of the cell averages of the conserved quantities. Therefore, our new approach performs the spatial WENO reconstruction twice: the first WENO reconstruction is carried out on the known cell averages of the conservative variables. The WENO polynomials are then used at the cell centers to compute point values of the conserved variables, which are subsequently converted into point values of the primitive variables. This is the only place where the conversion from conservative to primitive variables is needed in the new scheme. Then, a second WENO reconstruction is performed on the point values of the primitive variables to obtain piecewise high order reconstruction polynomials of the primitive variables. The reconstruction polynomials are subsequently evolved in time with a novel space-time finite element predictor that is directly applied to the governing PDE written in primitive form. The resulting space-time polynomials of the primitive variables can then be directly used as input for the numerical fluxes at the cell boundaries in the underlying conservative finite volume scheme. Hence, the number of necessary conversions from the conserved to the primitive variables is reduced to just one single conversion at each cell center. We have verified the validity of the new approach over a wide range of hyperbolic systems, including the classical Euler equations of gas dynamics, the special relativistic hydrodynamics (RHD) and ideal magnetohydrodynamics (RMHD) equations, as well as the Baer-Nunziato model for compressible two-phase flows. In all cases we have noticed that the new ADER

Ensemble averaged coherent state path integral for disordered bosons with a repulsive interaction (Derivation of mean field equations)

International Nuclear Information System (INIS)

Mieck, B.

2007-01-01

We consider bosonic atoms with a repulsive contact interaction in a trap potential for a Bose-Einstein condensation (BEC) and additionally include a random potential. The ensemble averages for two models of static (I) and dynamic (II) disorder are performed and investigated in parallel. The bosonic many body systems of the two disorder models are represented by coherent state path integrals on the Keldysh time contour which allow exact ensemble averages for zero and finite temperatures. These ensemble averages of coherent state path integrals therefore present alternatives to replica field theories or super-symmetric averaging techniques. Hubbard-Stratonovich transformations (HST) lead to two corresponding self-energies for the hermitian repulsive interaction and for the non-hermitian disorder-interaction. The self-energy of the repulsive interaction is absorbed by a shift into the disorder-self-energy which comprises as an element of a larger symplectic Lie algebra sp(4M) the self-energy of the repulsive interaction as a subalgebra (which is equivalent to the direct product of M x sp(2); 'M' is the number of discrete time intervals of the disorder-self-energy in the generating function). After removal of the remaining Gaussian integral for the self-energy of the repulsive interaction, the first order variations of the coherent state path integrals result in the exact mean field or saddle point equations, solely depending on the disorder-self-energy matrix. These equations can be solved by continued fractions and are reminiscent to the 'Nambu-Gorkov' Green function formalism in superconductivity because anomalous terms or pair condensates of the bosonic atoms are also included into the selfenergies. The derived mean field equations of the models with static (I) and dynamic (II) disorder are particularly applicable for BEC in d=3 spatial dimensions because of the singularity of the density of states at vanishing wavevector. However, one usually starts out from

Conservation of energy for the Euler-Korteweg equations

KAUST Repository

Dębiec, Tomasz

2017-12-30

In this article we study the principle of energy conservation for the Euler-Korteweg system. We formulate an Onsager-type sufficient regularity condition for weak solutions of the Euler-Korteweg system to conserve the total energy. The result applies to the system of Quantum Hydrodynamics.

Conservation of energy for the Euler-Korteweg equations

KAUST Repository

Dębiec, Tomasz; Gwiazda, Piotr; Świerczewska-Gwiazda, Agnieszka; Tzavaras, Athanasios

2017-01-01

In this article we study the principle of energy conservation for the Euler-Korteweg system. We formulate an Onsager-type sufficient regularity condition for weak solutions of the Euler-Korteweg system to conserve the total energy. The result applies to the system of Quantum Hydrodynamics.

Averaging hydraulic head, pressure head, and gravitational head in subsurface hydrology, and implications for averaged fluxes, and hydraulic conductivity

Directory of Open Access Journals (Sweden)

G. H. de Rooij

2009-07-01

Full Text Available Current theories for water flow in porous media are valid for scales much smaller than those at which problem of public interest manifest themselves. This provides a drive for upscaled flow equations with their associated upscaled parameters. Upscaling is often achieved through volume averaging, but the solution to the resulting closure problem imposes severe restrictions to the flow conditions that limit the practical applicability. Here, the derivation of a closed expression of the effective hydraulic conductivity is forfeited to circumvent the closure problem. Thus, more limited but practical results can be derived. At the Representative Elementary Volume scale and larger scales, the gravitational potential and fluid pressure are treated as additive potentials. The necessary requirement that the superposition be maintained across scales is combined with conservation of energy during volume integration to establish consistent upscaling equations for the various heads. The power of these upscaling equations is demonstrated by the derivation of upscaled water content-matric head relationships and the resolution of an apparent paradox reported in the literature that is shown to have arisen from a violation of the superposition principle. Applying the upscaling procedure to Darcy's Law leads to the general definition of an upscaled hydraulic conductivity. By examining this definition in detail for porous media with different degrees of heterogeneity, a series of criteria is derived that must be satisfied for Darcy's Law to remain valid at a larger scale.

Empirical solution of Green-Ampt equation using soil conservation service - curve number values

Science.gov (United States)

Grimaldi, S.; Petroselli, A.; Romano, N.

2012-09-01

The Soil Conservation Service - Curve Number (SCS-CN) method is a popular widely used rainfall-runoff model for quantifying the total stream-flow volume generated by storm rainfall, but its application is not appropriate for sub-daily resolutions. In order to overcome this drawback, the Green-Ampt (GA) infiltration equation is considered and an empirical solution is proposed and evaluated. The procedure, named CN4GA (Curve Number for Green-Ampt), aims to calibrate the Green-Ampt model parameters distributing in time the global information provided by the SCS-CN method. The proposed procedure is evaluated by analysing observed rainfall-runoff events; results show that CN4GA seems to provide better agreement with the observed hydrographs respect to the classic SCS-CN method.

Diffusion Processes Satisfying a Conservation Law Constraint

Directory of Open Access Journals (Sweden)

J. Bakosi

2014-01-01

Full Text Available We investigate coupled stochastic differential equations governing N nonnegative continuous random variables that satisfy a conservation principle. In various fields a conservation law requires a set of fluctuating variables to be nonnegative and (if appropriately normalized sum to one. As a result, any stochastic differential equation model to be realizable must not produce events outside of the allowed sample space. We develop a set of constraints on the drift and diffusion terms of such stochastic models to ensure that both the nonnegativity and the unit-sum conservation law constraints are satisfied as the variables evolve in time. We investigate the consequences of the developed constraints on the Fokker-Planck equation, the associated system of stochastic differential equations, and the evolution equations of the first four moments of the probability density function. We show that random variables, satisfying a conservation law constraint, represented by stochastic diffusion processes, must have diffusion terms that are coupled and nonlinear. The set of constraints developed enables the development of statistical representations of fluctuating variables satisfying a conservation law. We exemplify the results with the bivariate beta process and the multivariate Wright-Fisher, Dirichlet, and Lochner’s generalized Dirichlet processes.

Hojman's theorem of the third-order ordinary differential equation

International Nuclear Information System (INIS)

Hong-Sheng, Lü; Hong-Bin, Zhang; Shu-Long, Gu

2009-01-01

This paper extends Hojman's conservation law to the third-order differential equation. A new conserved quantity is constructed based on the Lie group of transformation generators of the equations of motion. The generators contain variations of the time and generalized coordinates. Two independent non-trivial conserved quantities of the third-order ordinary differential equation are obtained. A simple example is presented to illustrate the applications of the results. (general)

Cubication of conservative nonlinear oscillators

International Nuclear Information System (INIS)

Belendez, Augusto; Alvarez, Mariela L; Fernandez, Elena; Pascual, Inmaculada

2009-01-01

A cubication procedure of the nonlinear differential equation for conservative nonlinear oscillators is analysed and discussed. This scheme is based on the Chebyshev series expansion of the restoring force, and this allows us to approximate the original nonlinear differential equation by a Duffing equation in which the coefficients for the linear and cubic terms depend on the initial amplitude, A, while in a Taylor expansion of the restoring force these coefficients are independent of A. The replacement of the original nonlinear equation by an approximate Duffing equation allows us to obtain an approximate frequency-amplitude relation as a function of the complete elliptic integral of the first kind. Some conservative nonlinear oscillators are analysed to illustrate the usefulness and effectiveness of this scheme.

Nonlinear Conservation Laws and Finite Volume Methods

Science.gov (United States)

Leveque, Randall J.

Introduction Software Notation Classification of Differential Equations Derivation of Conservation Laws The Euler Equations of Gas Dynamics Dissipative Fluxes Source Terms Radiative Transfer and Isothermal Equations Multi-dimensional Conservation Laws The Shock Tube Problem Mathematical Theory of Hyperbolic Systems Scalar Equations Linear Hyperbolic Systems Nonlinear Systems The Riemann Problem for the Euler Equations Numerical Methods in One Dimension Finite Difference Theory Finite Volume Methods Importance of Conservation Form - Incorrect Shock Speeds Numerical Flux Functions Godunov's Method Approximate Riemann Solvers High-Resolution Methods Other Approaches Boundary Conditions Source Terms and Fractional Steps Unsplit Methods Fractional Step Methods General Formulation of Fractional Step Methods Stiff Source Terms Quasi-stationary Flow and Gravity Multi-dimensional Problems Dimensional Splitting Multi-dimensional Finite Volume Methods Grids and Adaptive Refinement Computational Difficulties Low-Density Flows Discrete Shocks and Viscous Profiles Start-Up Errors Wall Heating Slow-Moving Shocks Grid Orientation Effects Grid-Aligned Shocks Magnetohydrodynamics The MHD Equations One-Dimensional MHD Solving the Riemann Problem Nonstrict Hyperbolicity Stiffness The Divergence of B Riemann Problems in Multi-dimensional MHD Staggered Grids The 8-Wave Riemann Solver Relativistic Hydrodynamics Conservation Laws in Spacetime The Continuity Equation The 4-Momentum of a Particle The Stress-Energy Tensor Finite Volume Methods Multi-dimensional Relativistic Flow Gravitation and General Relativity References

Large Wilson loop averages from the Schwinger-Dyson equation

International Nuclear Information System (INIS)

Xue Shesheng

1987-01-01

Using Schwinger-Dyson equations for the large Wilson loop in abelian lattice gauge theories, we evaluate the vacuum expectation values of the Wilson loop of sizes 1x2, 2x2, 2x3, and so on, from which the string tension is extracted. (orig.)

Fractional Schroedinger equation

International Nuclear Information System (INIS)

Laskin, Nick

2002-01-01

Some properties of the fractional Schroedinger equation are studied. We prove the Hermiticity of the fractional Hamilton operator and establish the parity conservation law for fractional quantum mechanics. As physical applications of the fractional Schroedinger equation we find the energy spectra of a hydrogenlike atom (fractional 'Bohr atom') and of a fractional oscillator in the semiclassical approximation. An equation for the fractional probability current density is developed and discussed. We also discuss the relationships between the fractional and standard Schroedinger equations

The conservation laws for deformed classical models

International Nuclear Information System (INIS)

Klimek, M.

1994-01-01

The problem of deriving the conservation laws for deformed linear equations of motion is investigated. The conserved currents are obtained in explicit form and used in the construction of constants of motion. The equations for the set of non-interacting oscillators with arbitrary scale-time as well as the κ-Klein-Gordon equation are considered as an example of application of the method. (author) 9 refs

Guide for the estimation of the α and β coefficients in the Average enrichment equation as burnt function by fuel type

International Nuclear Information System (INIS)

Montes T, J.L.; Cortes C, C.C.

1992-08-01

The objective of the report is to determine manually or by means of a calculation sheet, the coefficients α and β of the average enrichment equation as function of the fuel burnt (B) using the Lineal Reactivity Pattern, with information generated by the RECORD code of the FMS package. (Author)

Scaling and scale invariance of conservation laws in Reynolds transport theorem framework

Science.gov (United States)

Haltas, Ismail; Ulusoy, Suleyman

2015-07-01

Scale invariance is the case where the solution of a physical process at a specified time-space scale can be linearly related to the solution of the processes at another time-space scale. Recent studies investigated the scale invariance conditions of hydrodynamic processes by applying the one-parameter Lie scaling transformations to the governing equations of the processes. Scale invariance of a physical process is usually achieved under certain conditions on the scaling ratios of the variables and parameters involved in the process. The foundational axioms of hydrodynamics are the conservation laws, namely, conservation of mass, conservation of linear momentum, and conservation of energy from continuum mechanics. They are formulated using the Reynolds transport theorem. Conventionally, Reynolds transport theorem formulates the conservation equations in integral form. Yet, differential form of the conservation equations can also be derived for an infinitesimal control volume. In the formulation of the governing equation of a process, one or more than one of the conservation laws and, some times, a constitutive relation are combined together. Differential forms of the conservation equations are used in the governing partial differential equation of the processes. Therefore, differential conservation equations constitute the fundamentals of the governing equations of the hydrodynamic processes. Applying the one-parameter Lie scaling transformation to the conservation laws in the Reynolds transport theorem framework instead of applying to the governing partial differential equations may lead to more fundamental conclusions on the scaling and scale invariance of the hydrodynamic processes. This study will investigate the scaling behavior and scale invariance conditions of the hydrodynamic processes by applying the one-parameter Lie scaling transformation to the conservation laws in the Reynolds transport theorem framework.

An alternative scheme of the Bogolyubov's average method

International Nuclear Information System (INIS)

Ortiz Peralta, T.; Ondarza R, R.; Camps C, E.

1990-01-01

In this paper the average energy and the magnetic moment conservation laws in the Drift Theory of charged particle motion are obtained in a simple way. The approach starts from the energy and magnetic moment conservation laws and afterwards the average is performed. This scheme is more economic from the standpoint of time and algebraic calculations than the usual procedure of Bogolyubov's method. (Author)

Hamiltonian structures of some non-linear evolution equations

International Nuclear Information System (INIS)

Tu, G.Z.

1983-06-01

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

A Second-Order Maximum Principle Preserving Lagrange Finite Element Technique for Nonlinear Scalar Conservation Equations

KAUST Repository

Guermond, Jean-Luc; Nazarov, Murtazo; Popov, Bojan; Yang, Yong

2014-01-01

© 2014 Society for Industrial and Applied Mathematics. This paper proposes an explicit, (at least) second-order, maximum principle satisfying, Lagrange finite element method for solving nonlinear scalar conservation equations. The technique is based on a new viscous bilinear form introduced in Guermond and Nazarov [Comput. Methods Appl. Mech. Engrg., 272 (2014), pp. 198-213], a high-order entropy viscosity method, and the Boris-Book-Zalesak flux correction technique. The algorithm works for arbitrary meshes in any space dimension and for all Lipschitz fluxes. The formal second-order accuracy of the method and its convergence properties are tested on a series of linear and nonlinear benchmark problems.

Asymptotic analysis of the average, steady, isotherml flow in coupled, parallel channels

International Nuclear Information System (INIS)

Lund, K.O.

1976-01-01

The conservation equations of mass and momentum are derived for the average flow of gases in coupled, parallel channels, or rod bundles. In the case of gas-cooled rod bundles the pitch of the rods is relatively large so the flows in the channels are strongly coupled. From this observation a perturbation parameter is derived and the descriptive equations are scaled using this parameter, which represents the ratio of the axial flow area to the transverse flow area, and which is of the order of 10 -3 in current gas-cooled fast breeder reactor designs. By expanding the velocities into perturbation series the equations for two channels are solved as an initial value problem, and the results compared to a finite difference solution of the same problem. The N-channel problem is solved to the lowest order as a two-point boundary value problem with the pressures specified at the inlet and the outlet. It is concluded from the study that asymptotic methods are effective in solving the flow problems of rod bundles; however, further work is required to evaluate the possible computational advantages of the methods

Exact Solutions of Generalized Modified Boussinesq, Kuramoto-Sivashinsky, and Camassa-Holm Equations via Double Reduction Theory

Directory of Open Access Journals (Sweden)

Zulfiqar Ali

2013-01-01

Full Text Available We find exact solutions of the Generalized Modified Boussinesq (GMB equation, the Kuromoto-Sivashinsky (KS equation the and, Camassa-Holm (CH equation by utilizing the double reduction theory related to conserved vectors. The fourth order GMB equation involves the arbitrary function and mixed derivative terms in highest derivative. The partial Noether’s approach yields seven conserved vectors for GMB equation and one conserved for vector KS equation. Due to presence of mixed derivative term the conserved vectors for GMB equation derived by the Noether like theorem do not satisfy the divergence relationship. The extra terms that constitute the trivial part of conserved vectors are adjusted and the resulting conserved vectors satisfy the divergence property. The double reduction theory yields two independent solutions and one reduction for GMB equation and one solution for KS equation. For CH equation two independent solutions are obtained elsewhere by double reduction theory with the help of conserved Vectors.

Conservation Laws for Partially Conservative Variable Mass Systems via d'Alembert's Principle

Institute of Scientific and Technical Information of China (English)

AFTAB Ahmed; NASEER Ahmed; QUDRAT Khan

2008-01-01

Conservation laws for partially conservative variable mass dynamical systems under symmetric infinitesimal transformations are determined. A generalization of Lagrange-d'Alembert's principle for a variable mass system in terms of asynchronous virtual variation is presented. The generalized Killing equations are obtained such that their solution yields the transformations and the associated conservation laws. An example illustrative of the theory is furnished at the end as well.

The impact of the form of the Euler equations for radial flow in cylindrical and spherical coordinates on numerical conservation and accuracy

Science.gov (United States)

Crittenden, P. E.; Balachandar, S.

2018-03-01

The radial one-dimensional Euler equations are often rewritten in what is known as the geometric source form. The differential operator is identical to the Cartesian case, but source terms result. Since the theory and numerical methods for the Cartesian case are well-developed, they are often applied without modification to cylindrical and spherical geometries. However, numerical conservation is lost. In this article, AUSM^+ -up is applied to a numerically conservative (discrete) form of the Euler equations labeled the geometric form, a nearly conservative variation termed the geometric flux form, and the geometric source form. The resulting numerical methods are compared analytically and numerically through three types of test problems: subsonic, smooth, steady-state solutions, Sedov's similarity solution for point or line-source explosions, and shock tube problems. Numerical conservation is analyzed for all three forms in both spherical and cylindrical coordinates. All three forms result in constant enthalpy for steady flows. The spatial truncation errors have essentially the same order of convergence, but the rate constants are superior for the geometric and geometric flux forms for the steady-state solutions. Only the geometric form produces the correct shock location for Sedov's solution, and a direct connection between the errors in the shock locations and energy conservation is found. The shock tube problems are evaluated with respect to feature location using an approximation with a very fine discretization as the benchmark. Extensions to second order appropriate for cylindrical and spherical coordinates are also presented and analyzed numerically. Conclusions are drawn, and recommendations are made. A derivation of the steady-state solution is given in the Appendix.

Effective average action for gauge theories and exact evolution equations

International Nuclear Information System (INIS)

Reuter, M.; Wetterich, C.

1993-11-01

We propose a new nonperturbative evolution equation for Yang-Mills theories. It describes the scale dependence of an effective action. The running of the nonabelian gauge coupling in arbitrary dimension is computed. (orig.)

First-order partial differential equations

CERN Document Server

Rhee, Hyun-Ku; Amundson, Neal R

2001-01-01

This first volume of a highly regarded two-volume text is fully usable on its own. After going over some of the preliminaries, the authors discuss mathematical models that yield first-order partial differential equations; motivations, classifications, and some methods of solution; linear and semilinear equations; chromatographic equations with finite rate expressions; homogeneous and nonhomogeneous quasilinear equations; formation and propagation of shocks; conservation equations, weak solutions, and shock layers; nonlinear equations; and variational problems. Exercises appear at the end of mo

Averaging Robertson-Walker cosmologies

International Nuclear Information System (INIS)

Brown, Iain A.; Robbers, Georg; Behrend, Juliane

2009-01-01

The cosmological backreaction arises when one directly averages the Einstein equations to recover an effective Robertson-Walker cosmology, rather than assuming a background a priori. While usually discussed in the context of dark energy, strictly speaking any cosmological model should be recovered from such a procedure. We apply the scalar spatial averaging formalism for the first time to linear Robertson-Walker universes containing matter, radiation and dark energy. The formalism employed is general and incorporates systems of multiple fluids with ease, allowing us to consider quantitatively the universe from deep radiation domination up to the present day in a natural, unified manner. Employing modified Boltzmann codes we evaluate numerically the discrepancies between the assumed and the averaged behaviour arising from the quadratic terms, finding the largest deviations for an Einstein-de Sitter universe, increasing rapidly with Hubble rate to a 0.01% effect for h = 0.701. For the ΛCDM concordance model, the backreaction is of the order of Ω eff 0 ≈ 4 × 10 −6 , with those for dark energy models being within a factor of two or three. The impacts at recombination are of the order of 10 −8 and those in deep radiation domination asymptote to a constant value. While the effective equations of state of the backreactions in Einstein-de Sitter, concordance and quintessence models are generally dust-like, a backreaction with an equation of state w eff < −1/3 can be found for strongly phantom models

Computational partial differential equations using Matlab

CERN Document Server

Li, Jichun

2008-01-01

Brief Overview of Partial Differential Equations The parabolic equations The wave equations The elliptic equations Differential equations in broader areasA quick review of numerical methods for PDEsFinite Difference Methods for Parabolic Equations Introduction Theoretical issues: stability, consistence, and convergence 1-D parabolic equations2-D and 3-D parabolic equationsNumerical examples with MATLAB codesFinite Difference Methods for Hyperbolic Equations IntroductionSome basic difference schemes Dissipation and dispersion errors Extensions to conservation lawsThe second-order hyperbolic PDE

Derivation of simplified basic equations of gas-liquid two-phase dispersed flow based on two-fluid model

International Nuclear Information System (INIS)

Kataoka, Isao; Tomiyama, Akio

2004-01-01

The simplified and physically reasonable basic equations for the gas-liquid dispersed flow were developed based on some appropriate assumptions and the treatment of dispersed phase as isothermal rigid particles. Based on the local instant formulation of mass, momentum and energy conservation of the dispersed flow, time-averaged equations were obtained assuming that physical quantities in the dispersed phase are uniform. These assumptions are approximately valid when phase change rate and/or chemical reaction rate are not so large at gas-liquid interface and there is no heat generation in within the dispersed phase. Detailed discussions were made on the characteristics of obtained basic equations and physical meanings of terms consisting the basic equations. It is shown that, in the derived averaged momentum equation, the terms of pressure gradient and viscous momentum diffusion do not appear and, in the energy equation, the term of molecular thermal diffusion heat flux does not appear. These characteristics of the derived equations were shown to be very consistent concerning the physical interpretation of the gas-liquid dispersed flow. Furthermore, the obtained basic equations are consistent with experiments for the dispersed flow where most of averaged physical quantities are obtained assuming that the distributions of those are uniform within the dispersed phase. Investigation was made on the problem whether the obtained basic equations are well-posed or ill-posed for the initial value problem. The eigenvalues of the simplified mass and momentum equations are calculated for basic equations obtained here and previous two-fluid basic equations with one pressure model. Well-posedness and ill-posedness are judged whether the eigenvalues are real or imaginary. The result indicated the newly developed basic equations always constitute the well-posed initial value problem while the previous two-fluid basic equations based on one pressure model constitutes ill

General Navier–Stokes-like momentum and mass-energy equations

Energy Technology Data Exchange (ETDEWEB)

Monreal, Jorge, E-mail: jmonreal@mail.usf.edu

2015-03-15

A new system of general Navier–Stokes-like equations is proposed to model electromagnetic flow utilizing analogues of hydrodynamic conservation equations. Such equations are intended to provide a different perspective and, potentially, a better understanding of electromagnetic mass, energy and momentum behaviour. Under such a new framework additional insights into electromagnetism could be gained. To that end, we propose a system of momentum and mass-energy conservation equations coupled through both momentum density and velocity vectors.

Invariant relations in Boussinesq-type equations

International Nuclear Information System (INIS)

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

2004-01-01

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

Conservation Laws for Partially Conservative Variable Mass Systems via d'Alembert's Principle

International Nuclear Information System (INIS)

Ahmed, Aftab; Ahmed, Naseer; Khan, Qudrat

2008-01-01

Conservation laws for partially conservative variable mass dynamical systems under symmetric infinitesimal transformations are determined. A generalization of Lagrange-d'Alembert's principle for a variable mass system in terms of asynchronous virtual variation is presented. The generalized Killing equations are obtained such that their solution yields the transformations and the associated conservation laws. An example illustrative of the theory is furnished at the end as well. (the physics of elementary particles and fields)

Variational approaches to conservation laws for a nonlinear ...

African Journals Online (AJOL)

The conservation laws of a nonlinear evolution equation of time dependent variable coefficients of damping and dispersion is studied. The equation under consideration is not derivable from a variational principle which means that one cannot appeal to the Noether theorem to determine the conservation laws. We utilize the ...

Deduction of Einstein equation from homogeneity of Riemann spacetime

Science.gov (United States)

Ni, Jun

2012-03-01

The symmetry of spacetime translation leads to the energy-momentum conservation. However, the Lagrange depends on spacetime coordinates, which makes the symmetry of spacetime translation different with other symmetry invariant explicitly under symmetry transformation. We need an equation to guarantee the symmetry of spacetime translation. In this talk, I will show that the Einstein equation can be deduced purely from the general covariant principle and the homogeneity of spacetime in the frame of quantum field theory. The Einstein equation is shown to be the equation to guarantee the symmetry of spacetime translation. Gravity is an apparent force due to the curvature of spacetime resulted from the conservation of energy-momentum. In the action of quantum field, only electroweak-strong interactions appear with curved spacetime metric determined by the Einstein equation.. The general covariant principle and the homogeneity of spacetime are merged into one basic principle: Any Riemann spacetime metric guaranteeing the energy-momentum conservation are equivalent, which can be called as the conserved general covariant principle. [4pt] [1] Jun Ni, Chin. Phys. Lett. 28, 110401 (2011).

Lagrangian averaging with geodesic mean.

Science.gov (United States)

Oliver, Marcel

2017-11-01

This paper revisits the derivation of the Lagrangian averaged Euler (LAE), or Euler- α equations in the light of an intrinsic definition of the averaged flow map as the geodesic mean on the volume-preserving diffeomorphism group. Under the additional assumption that first-order fluctuations are statistically isotropic and transported by the mean flow as a vector field, averaging of the kinetic energy Lagrangian of an ideal fluid yields the LAE Lagrangian. The derivation presented here assumes a Euclidean spatial domain without boundaries.

Newtonian hydrodynamic equations with relativistic pressure and velocity

Energy Technology Data Exchange (ETDEWEB)

Hwang, Jai-chan [Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Daegu 702-701 (Korea, Republic of); Noh, Hyerim [Korea Astronomy and Space Science Institute, Daejeon 305-348 (Korea, Republic of); Fabris, Júlio; Piattella, Oliver F.; Zimdahl, Winfried, E-mail: jchan@knu.ac.kr, E-mail: hr@kasi.re.kr, E-mail: fabris@pq.cnpq.br, E-mail: oliver.piattella@pq.cnpq.br, E-mail: winfried.zimdahl@pq.cnpq.br [Departamento de Fisica, Universidade Federal do Espirito Santo, Vitória (Brazil)

2016-07-01

We present a new approximation to include fully general relativistic pressure and velocity in Newtonian hydrodynamics. The energy conservation, momentum conservation and two Poisson's equations are consistently derived from Einstein's gravity in the zero-shear gauge assuming weak gravity and action-at-a-distance limit. The equations show proper special relativity limit in the absence of gravity. Our approximation is complementary to the post-Newtonian approximation and the equations are valid in fully nonlinear situations.

Partial differential equations

CERN Document Server

Evans, Lawrence C

2010-01-01

This text gives a comprehensive survey of modern techniques in the theoretical study of partial differential equations (PDEs) with particular emphasis on nonlinear equations. The exposition is divided into three parts: representation formulas for solutions; theory for linear partial differential equations; and theory for nonlinear partial differential equations. Included are complete treatments of the method of characteristics; energy methods within Sobolev spaces; regularity for second-order elliptic, parabolic, and hyperbolic equations; maximum principles; the multidimensional calculus of variations; viscosity solutions of Hamilton-Jacobi equations; shock waves and entropy criteria for conservation laws; and, much more.The author summarizes the relevant mathematics required to understand current research in PDEs, especially nonlinear PDEs. While he has reworked and simplified much of the classical theory (particularly the method of characteristics), he primarily emphasizes the modern interplay between funct...

Nonlinear Dirac Equations

Directory of Open Access Journals (Sweden)

Wei Khim Ng

2009-02-01

Full Text Available We construct nonlinear extensions of Dirac's relativistic electron equation that preserve its other desirable properties such as locality, separability, conservation of probability and Poincaré invariance. We determine the constraints that the nonlinear term must obey and classify the resultant non-polynomial nonlinearities in a double expansion in the degree of nonlinearity and number of derivatives. We give explicit examples of such nonlinear equations, studying their discrete symmetries and other properties. Motivated by some previously suggested applications we then consider nonlinear terms that simultaneously violate Lorentz covariance and again study various explicit examples. We contrast our equations and construction procedure with others in the literature and also show that our equations are not gauge equivalent to the linear Dirac equation. Finally we outline various physical applications for these equations.

Discontinuous Galerkin Method for Hyperbolic Conservation Laws

KAUST Repository

Mousikou, Ioanna

2016-11-11

Hyperbolic conservation laws form a special class of partial differential equations. They describe phenomena that involve conserved quantities and their solutions show discontinuities which reflect the formation of shock waves. We consider one-dimensional systems of hyperbolic conservation laws and produce approximations using finite difference, finite volume and finite element methods. Due to stability issues of classical finite element methods for hyperbolic conservation laws, we study the discontinuous Galerkin method, which was recently introduced. The method involves completely discontinuous basis functions across each element and it can be considered as a combination of finite volume and finite element methods. We illustrate the implementation of discontinuous Galerkin method using Legendre polynomials, in case of scalar equations and in case of quasi-linear systems, and we review important theoretical results about stability and convergence of the method. The applications of finite volume and discontinuous Galerkin methods to linear and non-linear scalar equations, as well as to the system of elastodynamics, are exhibited.

Discontinuous Galerkin Method for Hyperbolic Conservation Laws

KAUST Repository

Mousikou, Ioanna

2016-01-01

Hyperbolic conservation laws form a special class of partial differential equations. They describe phenomena that involve conserved quantities and their solutions show discontinuities which reflect the formation of shock waves. We consider one-dimensional systems of hyperbolic conservation laws and produce approximations using finite difference, finite volume and finite element methods. Due to stability issues of classical finite element methods for hyperbolic conservation laws, we study the discontinuous Galerkin method, which was recently introduced. The method involves completely discontinuous basis functions across each element and it can be considered as a combination of finite volume and finite element methods. We illustrate the implementation of discontinuous Galerkin method using Legendre polynomials, in case of scalar equations and in case of quasi-linear systems, and we review important theoretical results about stability and convergence of the method. The applications of finite volume and discontinuous Galerkin methods to linear and non-linear scalar equations, as well as to the system of elastodynamics, are exhibited.

Geometric approach to soliton equations

International Nuclear Information System (INIS)

Sasaki, R.

1979-09-01

A class of nonlinear equations that can be solved in terms of nxn scattering problem is investigated. A systematic geometric method of exploiting conservation laws and related equations, the so-called prolongation structure, is worked out. The nxn problem is reduced to nsub(n-1)x(n-1) problems and finally to 2x2 problems, which have been comprehensively investigated recently by the author. A general method of deriving the infinite numbers of polynomial conservation laws for an nxn problem is presented. The cases of 3x3 and 2x2 problems are discussed explicitly. (Auth.)

A systematic comparison of two-equation Reynolds-averaged Navier-Stokes turbulence models applied to shock-cloud interactions

Science.gov (United States)

Goodson, Matthew D.; Heitsch, Fabian; Eklund, Karl; Williams, Virginia A.

2017-07-01

Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-averaged Navier-Stokes turbulence models, and we implement six models in the athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixing layer width. We then test the validity of each model into the supersonic regime, showing that compressibility corrections can improve agreement with experiment. For models with buoyancy effects, we also verify our implementation via the growth of the Rayleigh-Taylor instability in a stratified medium. The models are then applied to the ubiquitous astrophysical shock-cloud interaction in three dimensions. We focus on the mixing of shock and cloud material, comparing results from turbulence models to high-resolution simulations (up to 200 cells per cloud radius) and ensemble-averaged simulations. We find that the turbulence models lead to increased spreading and mixing of the cloud, although no two models predict the same result. Increased mixing is also observed in inviscid simulations at resolutions greater than 100 cells per radius, which suggests that the turbulent mixing begins to be resolved.

Two-level schemes for the advection equation

Science.gov (United States)

Vabishchevich, Petr N.

2018-06-01

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

Relativistic extension of a charge-conservative finite element solver for time-dependent Maxwell-Vlasov equations

Science.gov (United States)

Na, D.-Y.; Moon, H.; Omelchenko, Y. A.; Teixeira, F. L.

2018-01-01

Accurate modeling of relativistic particle motion is essential for physical predictions in many problems involving vacuum electronic devices, particle accelerators, and relativistic plasmas. A local, explicit, and charge-conserving finite-element time-domain (FETD) particle-in-cell (PIC) algorithm for time-dependent (non-relativistic) Maxwell-Vlasov equations on irregular (unstructured) meshes was recently developed by Moon et al. [Comput. Phys. Commun. 194, 43 (2015); IEEE Trans. Plasma Sci. 44, 1353 (2016)]. Here, we extend this FETD-PIC algorithm to the relativistic regime by implementing and comparing three relativistic particle-pushers: (relativistic) Boris, Vay, and Higuera-Cary. We illustrate the application of the proposed relativistic FETD-PIC algorithm for the analysis of particle cyclotron motion at relativistic speeds, harmonic particle oscillation in the Lorentz-boosted frame, and relativistic Bernstein modes in magnetized charge-neutral (pair) plasmas.

Stochastic optimal control, forward-backward stochastic differential equations and the Schroedinger equation

Energy Technology Data Exchange (ETDEWEB)

Paul, Wolfgang; Koeppe, Jeanette [Institut fuer Physik, Martin Luther Universitaet, 06099 Halle (Germany); Grecksch, Wilfried [Institut fuer Mathematik, Martin Luther Universitaet, 06099 Halle (Germany)

2016-07-01

The standard approach to solve a non-relativistic quantum problem is through analytical or numerical solution of the Schroedinger equation. We show a way to go around it. This way is based on the derivation of the Schroedinger equation from conservative diffusion processes and the establishment of (several) stochastic variational principles leading to the Schroedinger equation under the assumption of a kinematics described by Nelson's diffusion processes. Mathematically, the variational principle can be considered as a stochastic optimal control problem linked to the forward-backward stochastic differential equations of Nelson's stochastic mechanics. The Hamilton-Jacobi-Bellmann equation of this control problem is the Schroedinger equation. We present the mathematical background and how to turn it into a numerical scheme for analyzing a quantum system without using the Schroedinger equation and exemplify the approach for a simple 1d problem.

Constitutive equations for two-phase flows

International Nuclear Information System (INIS)

Boure, J.A.

1974-12-01

The mathematical model of a system of fluids consists of several kinds of equations complemented by boundary and initial conditions. The first kind equations result from the application to the system, of the fundamental conservation laws (mass, momentum, energy). The second kind equations characterize the fluid itself, i.e. its intrinsic properties and in particular its mechanical and thermodynamical behavior. They are the mathematical model of the particular fluid under consideration, the laws they expressed are so called the constitutive equations of the fluid. In practice the constitutive equations cannot be fully stated without reference to the conservation laws. Two classes of model have been distinguished: mixture model and two-fluid models. In mixture models, the mixture is considered as a single fluid. Besides the usual friction factor and heat transfer correlations, a single constitutive law is necessary. In diffusion models, the mixture equation of state is replaced by the phasic equations of state and by three consitutive laws, for phase change mass transfer, drift velocity and thermal non-equilibrium respectively. In the two-fluid models, the two phases are considered separately; two phasic equations of state, two friction factor correlations, two heat transfer correlations and four constitutive laws are included [fr

Bridging the Knowledge Gaps between Richards' Equation and Budyko Equation

Science.gov (United States)

Wang, D.

2017-12-01

The empirical Budyko equation represents the partitioning of mean annual precipitation into evaporation and runoff. Richards' equation, based on Darcy's law, represents the movement of water in unsaturated soils. The linkage between Richards' equation and Budyko equation is presented by invoking the empirical Soil Conservation Service curve number (SCS-CN) model for computing surface runoff at the event-scale. The basis of the SCS-CN method is the proportionality relationship, i.e., the ratio of continuing abstraction to its potential is equal to the ratio of surface runoff to its potential value. The proportionality relationship can be derived from the Richards' equation for computing infiltration excess and saturation excess models at the catchment scale. Meanwhile, the generalized proportionality relationship is demonstrated as the common basis of SCS-CN method, monthly "abcd" model, and Budyko equation. Therefore, the linkage between Darcy's law and the emergent pattern of mean annual water balance at the catchment scale is presented through the proportionality relationship.

The average action for scalar fields near phase transitions

International Nuclear Information System (INIS)

Wetterich, C.

1991-08-01

We compute the average action for fields in two, three and four dimensions, including the effects of wave function renormalization. A study of the one loop evolution equations for the scale dependence of the average action gives a unified picture of the qualitatively different behaviour in various dimensions for discrete as well as abelian and nonabelian continuous symmetry. The different phases and the phase transitions can be infered from the evolution equation. (orig.)

The Use of Nonlinear Constitutive Equations to Evaluate Draw Resistance and Filter Ventilation

Directory of Open Access Journals (Sweden)

Eitzinger B

2014-12-01

Full Text Available This study investigates by nonlinear constitutive equations the influence of tipping paper, cigarette paper, filter, and tobacco rod on the degree of filter ventilation and draw resistance. Starting from the laws of conservation, the path to the theory of fluid dynamics in porous media and Darcy's law is reviewed and, as an extension to Darcy's law, two different nonlinear pressure drop-flow relations are proposed. It is proven that these relations are valid constitutive equations and the partial differential equations for the stationary flow in an unlit cigarette covering anisotropic, inhomogeneous and nonlinear behaviour are derived. From these equations a system of ordinary differential equations for the one-dimensional flow in the cigarette is derived by averaging pressure and velocity over the cross section of the cigarette. By further integration, the concept of an electrical analog is reached and discussed in the light of nonlinear pressure drop-flow relations. By numerical calculations based on the system of ordinary differential equations, it is shown that the influence of nonlinearities cannot be neglected because variations in the degree of filter ventilation can reach up to 20% of its nominal value.

On energy conservation in extended magnetohydrodynamics

International Nuclear Information System (INIS)

Kimura, Keiji; Morrison, P. J.

2014-01-01

A systematic study of energy conservation for extended magnetohydrodynamic models that include Hall terms and electron inertia is performed. It is observed that commonly used models do not conserve energy in the ideal limit, i.e., when viscosity and resistivity are neglected. In particular, a term in the momentum equation that is often neglected is seen to be needed for conservation of energy

Conservation laws in disordered electron systems: Thermodynamic limit and configurational averaging

Czech Academy of Sciences Publication Activity Database

Janiš, Václav; Kolorenč, Jindřich

2004-01-01

Roč. 241, č. 9 (2004), s. 2032-2042 ISSN 0370-1972 R&D Projects: GA ČR GA202/04/1055 Institutional research plan: CEZ:AV0Z1010914 Keywords : conservation laws * noninteracting disordered electrons * diffusion pole Subject RIV: BE - Theoretical Physics Impact factor: 0.982, year: 2004

Conservative numerical schemes for Euler-Lagrange equations

Energy Technology Data Exchange (ETDEWEB)

Vazquez, L. [Universidad Complutense, Madrid (Spain). Dept. de Matematica Aplicada; Jimenez, S. [Universidad Alfonso X El Sabio, Madrid (Spain). Dept. de Matematica Aplicada

1999-05-01

As a preliminary step to study magnetic field lines, the authors seek numerical schemes that reproduce at discrete level the significant feature of the continuous model, based on an underling Lagrangian structure. The resulting scheme give discrete counterparts of the variation law for the energy as well of as the Euler-Lagrange equations and their symmetries.

Comments on the integrability of the loop-space chiral equations

International Nuclear Information System (INIS)

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

1980-01-01

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

Conservative fourth-order time integration of non-linear dynamic systems

DEFF Research Database (Denmark)

Krenk, Steen

2015-01-01

An energy conserving time integration algorithm with fourth-order accuracy is developed for dynamic systems with nonlinear stiffness. The discrete formulation is derived by integrating the differential state-space equations of motion over the integration time increment, and then evaluating...... the resulting time integrals of the inertia and stiffness terms via integration by parts. This process introduces the time derivatives of the state space variables, and these are then substituted from the original state-space differential equations. The resulting discrete form of the state-space equations...... is a direct fourth-order accurate representation of the original differential equations. This fourth-order form is energy conserving for systems with force potential in the form of a quartic polynomial in the displacement components. Energy conservation for a force potential of general form is obtained...

Estimating glomerular filtration rate (GFR) in children. The average between a cystatin C- and a creatinine-based equation improves estimation of GFR in both children and adults and enables diagnosing Shrunken Pore Syndrome.

Science.gov (United States)

Leion, Felicia; Hegbrant, Josefine; den Bakker, Emil; Jonsson, Magnus; Abrahamson, Magnus; Nyman, Ulf; Björk, Jonas; Lindström, Veronica; Larsson, Anders; Bökenkamp, Arend; Grubb, Anders

2017-09-01

Estimating glomerular filtration rate (GFR) in adults by using the average of values obtained by a cystatin C- (eGFR cystatin C ) and a creatinine-based (eGFR creatinine ) equation shows at least the same diagnostic performance as GFR estimates obtained by equations using only one of these analytes or by complex equations using both analytes. Comparison of eGFR cystatin C and eGFR creatinine plays a pivotal role in the diagnosis of Shrunken Pore Syndrome, where low eGFR cystatin C compared to eGFR creatinine has been associated with higher mortality in adults. The present study was undertaken to elucidate if this concept can also be applied in children. Using iohexol and inulin clearance as gold standard in 702 children, we studied the diagnostic performance of 10 creatinine-based, 5 cystatin C-based and 3 combined cystatin C-creatinine eGFR equations and compared them to the result of the average of 9 pairs of a eGFR cystatin C and a eGFR creatinine estimate. While creatinine-based GFR estimations are unsuitable in children unless calibrated in a pediatric or mixed pediatric-adult population, cystatin C-based estimations in general performed well in children. The average of a suitable creatinine-based and a cystatin C-based equation generally displayed a better diagnostic performance than estimates obtained by equations using only one of these analytes or by complex equations using both analytes. Comparing eGFR cystatin and eGFR creatinine may help identify pediatric patients with Shrunken Pore Syndrome.

Conservation laws with coinciding smooth solutions but different conserved variables

Science.gov (United States)

Colombo, Rinaldo M.; Guerra, Graziano

2018-04-01

Consider two hyperbolic systems of conservation laws in one space dimension with the same eigenvalues and (right) eigenvectors. We prove that solutions to Cauchy problems with the same initial data differ at third order in the total variation of the initial datum. As a first application, relying on the classical Glimm-Lax result (Glimm and Lax in Decay of solutions of systems of nonlinear hyperbolic conservation laws. Memoirs of the American Mathematical Society, No. 101. American Mathematical Society, Providence, 1970), we obtain estimates improving those in Saint-Raymond (Arch Ration Mech Anal 155(3):171-199, 2000) on the distance between solutions to the isentropic and non-isentropic inviscid compressible Euler equations, under general equations of state. Further applications are to the general scalar case, where rather precise estimates are obtained, to an approximation by Di Perna of the p-system and to a traffic model.

Lorentz invariance violation and charge (non)conservation: A general theoretical frame for extensions of the Maxwell equations

International Nuclear Information System (INIS)

Laemmerzahl, Claus; Macias, Alfredo; Mueller, Holger

2005-01-01

All quantum gravity approaches lead to small modifications in the standard laws of physics which in most cases lead to violations of Lorentz invariance. One particular example is the extended standard model (SME). Here, a general phenomenological approach for extensions of the Maxwell equations is presented which turns out to be more general than the SME and which covers charge nonconservation (CNC), too. The new Lorentz invariance violating terms cannot be probed by optical experiments but need, instead, the exploration of the electromagnetic field created by a point charge or a magnetic dipole. Some scalar tensor theories and higher dimensional brane theories predict CNC in four dimensions and some models violating special relativity have been shown to be connected with CNC. Its relation to the Einstein Equivalence Principle has been discussed. Because of this upcoming interest, the experimental status of electric charge conservation is reviewed. Up to now there seem to exist no unique tests of charge conservation. CNC is related to the precession of polarization, to a modification of the 1/r-Coulomb potential, and to a time dependence of the fine structure constant. This gives the opportunity to describe a dedicated search for CNC

Equations of motion for two-phase flow in a pin bundle of a nuclear reactor

International Nuclear Information System (INIS)

Chawla, T.C.; Ishii, M.

1978-01-01

By performing Eulerian area averaging over a channel area of the local continuity, momentum, and energy equations for single phase turbulent flow and assuming each phase in two-phase flows to be continuum but coupled by the appropriate 'jump' conditions at the interface, the corresponding axial macroscopic balances for two-fluid model in a pin bundle are obtained. To determine the crossflow, a momentum equation in transverse (to the gap between the pins) direction is obtained for each phase by carrying out Eulerian segment averaging of the local momentum equation, where the segment is taken parallel to the gap. By considering the mixture as a whole, a diffusion model based on drift-flux velocity is formulated. In the axial direction it is expressed in terms of three mixture conservation equations of mass, momentum, and energy with one additional continuity equation for the vapor phase. For the determination of crossflow, transverse momentum equation for a mixture is obtained. It is considered that the previous formulation of the two-phase flow based on the 'slip' flow model and the integral subchannel balances using finite control volumes is inadequate in that the model is heuristic and, a priori, assumes the order of magnitude of the terms, also the model is incomplete and incorrect when applied to two-phase mixtures in thermal non-equilibrium such as during accidental depressurization of a water cooled reactor. The governing equations presented are shown to be a very formal and sound physical basis and are indispensable for physically correct methods of analyzing two-phase flows in a pin bundle. (author)

Characteristics of phase-averaged equations for modulated wave groups

NARCIS (Netherlands)

Klopman, G.; Petit, H.A.H.; Battjes, J.A.

2000-01-01

The project concerns the influence of long waves on coastal morphology. The modelling of the combined motion of the long waves and short waves in the horizontal plane is done by phase-averaging over the short wave motion and using intra-wave modelling for the long waves, see e.g. Roelvink (1993).

Modeling of frequency-domain scalar wave equation with the average-derivative optimal scheme based on a multigrid-preconditioned iterative solver

Science.gov (United States)

Cao, Jian; Chen, Jing-Bo; Dai, Meng-Xue

2018-01-01

An efficient finite-difference frequency-domain modeling of seismic wave propagation relies on the discrete schemes and appropriate solving methods. The average-derivative optimal scheme for the scalar wave modeling is advantageous in terms of the storage saving for the system of linear equations and the flexibility for arbitrary directional sampling intervals. However, using a LU-decomposition-based direct solver to solve its resulting system of linear equations is very costly for both memory and computational requirements. To address this issue, we consider establishing a multigrid-preconditioned BI-CGSTAB iterative solver fit for the average-derivative optimal scheme. The choice of preconditioning matrix and its corresponding multigrid components is made with the help of Fourier spectral analysis and local mode analysis, respectively, which is important for the convergence. Furthermore, we find that for the computation with unequal directional sampling interval, the anisotropic smoothing in the multigrid precondition may affect the convergence rate of this iterative solver. Successful numerical applications of this iterative solver for the homogenous and heterogeneous models in 2D and 3D are presented where the significant reduction of computer memory and the improvement of computational efficiency are demonstrated by comparison with the direct solver. In the numerical experiments, we also show that the unequal directional sampling interval will weaken the advantage of this multigrid-preconditioned iterative solver in the computing speed or, even worse, could reduce its accuracy in some cases, which implies the need for a reasonable control of directional sampling interval in the discretization.

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

Analysis of Influence of the Thermal Dependence of Air Thermophysical Properties on the Accuracy of Simulation of Heat Transfer in a Turbulent Flow in Case of Applying Different Methods of Averaging Navier-Stokes Equations

Directory of Open Access Journals (Sweden)

A. D. Kliukvin

2014-01-01

Full Text Available There is theoretically investigated the influence of thermal dependence of air thermophysical properties on accuracy of heat transfer problems solution in a turbulent flow when using different methods of averaging the Navier-Stokes equations.There is analyzed the practicability of using particular method of averaging the NavierStokes equations when it’s necessary to clarify the solution of heat transfer problem taking into account the variability of air thermophysical properties.It’s shown that Reynolds and Favre averaging (the most common methods of averaging the Navier-Stokes equations are not effective in this case because these methods inaccurately describe behavior of large scale turbulent structures which strongly depends on geometry of particular flow. Thus it’s necessary to use more universal methods of turbulent flow simulation which are not based on averaging of all turbulent scales.In the article it’s shown that instead of Reynold and Favre averaging it’s possible to use large eddy simulation whereby turbulent structures are divided into small-scale and large-scale ones with subsequent modelling of small-scale ones only. But this approach leads to the necessarity of increasing the computational power by 2-3 orders.For different methods of averaging the form of additional terms of averaged Navier-Stokes equations in case of accounting pulsation of thermophysical properties of the air is obtained.On the example of a submerged heated air jet the errors (which occur when neglecting the thermal dependence of air thermophysical properties on averaged flow temperature in determination of convectional and conductive components of heat flux and viscous stresses are evaluated. It’s shown that the greatest increase of solution accuracy can be obtained in case of the flows with high temperature gradients.Finally using infinite Teylor series it’s found that underestimation of convective and conductive components of heat flux and

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.

18 CFR 301.7 - Average System Cost methodology functionalization.

Science.gov (United States)

2010-04-01

... 18 Conservation of Power and Water Resources 1 2010-04-01 2010-04-01 false Average System Cost... REGULATORY COMMISSION, DEPARTMENT OF ENERGY REGULATIONS FOR FEDERAL POWER MARKETING ADMINISTRATIONS AVERAGE SYSTEM COST METHODOLOGY FOR SALES FROM UTILITIES TO BONNEVILLE POWER ADMINISTRATION UNDER NORTHWEST POWER...

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

International Nuclear Information System (INIS)

Yu Guofu; Tam, H.-W.

2006-01-01

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

Formal derivation of a 6 equation macro scale model for two-phase flows - link with the 4 equation macro scale model implemented in Flica 4; Etablissement formel d'un modele diphasique macroscopique a 6 equations - lien avec le modele macroscopique a 4 equations flica 4

Energy Technology Data Exchange (ETDEWEB)

Gregoire, O

2008-07-01

In order to simulate nuclear reactor cores, we presently use the 4 equation model implemented within FLICA4 code. This model is complemented with 2 algebraic closures for thermal disequilibrium and relative velocity between phases. Using such closures, means an 'a priori' knowledge of flows calculated in order to ensure that modelling assumptions apply. In order to improve the degree of universality to our macroscopic modelling, we propose in the report to derive a more general 6 equation model (balance equations for mass, momentum and enthalpy for each phase) for 2-phase flows. We apply the up-scaling procedure (Whitaker, 1999) classically used in porous media analysis to the statistically averaged equations (Aniel-Buchheit et al., 2003). By doing this, we apply the double-averaging procedure (Pedras and De Lemos, 2001 and Pinson et al. 2006): statistical and spatial averages. Then, using weighted averages (analogous to Favre's average) we extend the spatial averaging concept to variable density and 2-phase flows. This approach allows the global recovering of the structure of the systems of equations implemented in industrial codes. Supplementary contributions, such as dispersion, are also highlighted. Mechanical and thermal exchanges between solids and fluid are formally derived. Then, thanks to realistic simplifying assumptions, we show how it is possible to derive the original 4 equation model from the full 6 equation model. (author)

The conservation laws of nonrelativistic classical and quantum mechanics for a system of interacting particles

International Nuclear Information System (INIS)

Havas, P.

1978-01-01

The various classical or quantum mechanical equations describing a system of N particles with central two-body interactions are invariant under the 10 transformations of the Galilei group, and for interaction potential inversely proportional to the squares of the particle separations also under two further transformations. From the invariance of the corresponding classical and quantum mechanical variation principles under this 12-parameter conformal extension of the Galilei group, the 'Jacobi-Schroedinger group', the 12 well-known conservation laws of Newtonian dynamics as well as 12 local conservation laws implied by the Schroedinger equation are obtained via Noether's theorem. Under appropriate conditions on the wave functions, these local laws yield 12 global conservation laws which are analogous to the Newtonian ones. The Hamiltonian-Jacobi equation implies a classical equation differing from the Schroedinger equation only by a potential-like term involving the Van Vleck determinant, from which 12 local balance equations and the corresponding global equations are obtained, which under certain conditions reduce the true conservation laws. (Auth.)

The average kinetic energy of the heavy quark in Λb in the Bethe-Salpeter equation approach

International Nuclear Information System (INIS)

Guo, X.-H.; Wu, H.-K.

2007-01-01

In the previous paper, based on the SU(2) f xSU(2) s heavy quark symmetries of the QCD Lagrangian in the heavy quark limit, the Bethe-Salpeter equation for the heavy baryon Λ b was established with the picture that Λ b is composed of a heavy quark and a scalar light diquark. In the present work, we apply this model to calculate μ π 2 for Λ b , the average kinetic energy of the heavy quark inside Λ b . This quantity is particularly interesting since it can be measured in experiments and since it contributes to the inclusive semileptonic decays of Λ b when contributions from higher order terms in 1/M b expansions are taken into account and consequently influences the determination of the Cabibbo-Kobayashi-Maskawa matrix elements V ub and V cb . We find that μ π 2 for Λ b is 0.25GeV 2 ∼0.95GeV 2 , depending on the parameters in the model including the light diquark mass and the interaction strength between the heavy quark and the light diquark in the kernel of the BS equation. We also find that this result is consistent with the value of μ π 2 for Λ b which is derived from the experimental value of μ π 2 for the B meson with the aid of the heavy quark effective theory

An interpolation between the wave and diffusion equations through the fractional evolution equations Dirac like

International Nuclear Information System (INIS)

Pierantozzi, T.; Vazquez, L.

2005-01-01

Through fractional calculus and following the method used by Dirac to obtain his well-known equation from the Klein-Gordon equation, we analyze a possible interpolation between the Dirac and the diffusion equations in one space dimension. We study the transition between the hyperbolic and parabolic behaviors by means of the generalization of the D'Alembert formula for the classical wave equation and the invariance under space and time inversions of the interpolating fractional evolution equations Dirac like. Such invariance depends on the values of the fractional index and is related to the nonlocal property of the time fractional differential operator. For this system of fractional evolution equations, we also find an associated conserved quantity analogous to the Hamiltonian for the classical Dirac case

Conservative multi-implicit integral deferred correction methods with adaptive mesh refinement

International Nuclear Information System (INIS)

Layton, A.T.

2004-01-01

In most models of reacting gas dynamics, the characteristic time scales of chemical reactions are much shorter than the hydrodynamic and diffusive time scales, rendering the reaction part of the model equations stiff. Moreover, nonlinear forcings may introduce into the solutions sharp gradients or shocks, the robust behavior and correct propagation of which require the use of specialized spatial discretization procedures. This study presents high-order conservative methods for the temporal integration of model equations of reacting flows. By means of a method of lines discretization on the flux difference form of the equations, these methods compute approximations to the cell-averaged or finite-volume solution. The temporal discretization is based on a multi-implicit generalization of integral deferred correction methods. The advection term is integrated explicitly, and the diffusion and reaction terms are treated implicitly but independently, with the splitting errors present in traditional operator splitting methods reduced via the integral deferred correction procedure. To reduce computational cost, time steps used to integrate processes with widely-differing time scales may differ in size. (author)

Introduction to the method of average magnitude analysis and application to natural convection in cavities

International Nuclear Information System (INIS)

Lykoudis, P.S.

1995-01-01

The method of Average Magnitude Analysis is a mixture of the Integral Method and the Order of Magnitude Analysis. The paper shows how the differential equations of conservation for steady-state, laminar, boundary layer flows are converted to a system of algebraic equations, where the result is a sum of the order of magnitude of each term, multiplied by, a weight coefficient. These coefficients are determined from integrals containing the assumed velocity and temperature profiles. The method is illustrated by applying it to the case of drag and heat transfer over an infinite flat plate. It is then applied to the case of natural convection over an infinite flat plate with and without the presence of a horizontal magnetic field, and subsequently to enclosures of aspect ratios of one or higher. The final correlation in this instance yields the Nusselt number as a function of the aspect ratio and the Rayleigh and Prandtl numbers. This correlation is tested against a wide range of small and large values of these parameters. 19 refs., 4 figs

ADM pseudotensors, conserved quantities and covariant conservation laws in general relativity

International Nuclear Information System (INIS)

Fatibene, L.; Ferraris, M.; Francaviglia, M.; Lusanna, L.

2012-01-01

The ADM formalism is reviewed and techniques for decomposing generic components of metric, connection and curvature are obtained. These techniques will turn out to be enough to decompose not only Einstein equations but also covariant conservation laws. Then a number of independent sets of hypotheses that are sufficient (though not necessary) to obtain standard ADM quantities (and Hamiltonian) from covariant conservation laws are considered. This determines explicitly the range in which standard techniques are equivalent to covariant conserved quantities. The Schwarzschild metric in different coordinates is then considered, showing how the standard ADM quantities fail dramatically in non-Cartesian coordinates or even worse when asymptotically flatness is not manifest; while, in view of their covariance, covariant conservation laws give the correct result in all cases. - Highlights: ► In the paper ADM conserved quantities for GR are obtained from augmented conservation laws. ► Boundary conditions for this to be possible are considered and compared with the literature. ► Some different forms of Schwarzschild solutions are considered as simple examples of different boundary conditions.

Kinetic equations in dirty superconductors

International Nuclear Information System (INIS)

Kraehenbuehl, Y.

1981-01-01

Kinetic equations for superconductors in the dirty limit are derived using a method developed for superfluid systems, which allows a systematic expansion in small parameters; exact charge conservation is obeyed. (orig.)

Control of Stochastic Master Equation Models of Genetic Regulatory Networks by Approximating Their Average Behavior

Science.gov (United States)

Umut Caglar, Mehmet; Pal, Ranadip

2010-10-01

The central dogma of molecular biology states that ``information cannot be transferred back from protein to either protein or nucleic acid.'' However, this assumption is not exactly correct in most of the cases. There are a lot of feedback loops and interactions between different levels of systems. These types of interactions are hard to analyze due to the lack of data in the cellular level and probabilistic nature of interactions. Probabilistic models like Stochastic Master Equation (SME) or deterministic models like differential equations (DE) can be used to analyze these types of interactions. SME models based on chemical master equation (CME) can provide detailed representation of genetic regulatory system, but their use is restricted by the large data requirements and computational costs of calculations. The differential equations models on the other hand, have low calculation costs and much more adequate to generate control procedures on the system; but they are not adequate to investigate the probabilistic nature of interactions. In this work the success of the mapping between SME and DE is analyzed, and the success of a control policy generated by DE model with respect to SME model is examined. Index Terms--- Stochastic Master Equation models, Differential Equation Models, Control Policy Design, Systems biology

Generalized continuity equations from two-field Schrödinger Lagrangians

Science.gov (United States)

Spourdalakis, A. G. B.; Pappas, G.; Morfonios, C. Â. V.; Kalozoumis, P. A.; Diakonos, F. K.; Schmelcher, P.

2016-11-01

A variational scheme for the derivation of generalized, symmetry-induced continuity equations for Hermitian and non-Hermitian quantum mechanical systems is developed. We introduce a Lagrangian which involves two complex wave fields and whose global invariance under dilation and phase variations leads to a mixed continuity equation for the two fields. In combination with discrete spatial symmetries of the underlying Hamiltonian, the mixed continuity equation is shown to produce bilocal conservation laws for a single field. This leads to generalized conserved charges for vanishing boundary currents and to divergenceless bilocal currents for stationary states. The formalism reproduces the bilocal continuity equation obtained in the special case of P T -symmetric quantum mechanics and paraxial optics.

On the regularization in the Callan-Symanzik equation

International Nuclear Information System (INIS)

Fujii, Yasunori; Takahashi, Yasushi

1975-01-01

The conservative approach of canonical theory of broken scale invariance to the Callan-Symanzik equation is pushed further with the Pauli-Villars regulators. The authors confirm that the Callan-Symanzik equation is derived in a completely general manner. (BMS) [de

Comparison of Fully-Compressible Equation Sets for Atmospheric Dynamics

Science.gov (United States)

Ahmad, Nashat N.

2016-01-01

Traditionally, the equation for the conservation of energy used in atmospheric models is based on potential temperature and is used in place of the total energy conservation. This paper compares the application of the two equations sets for both the Euler and the Navier-Stokes solutions using several benchmark test cases. A high-resolution wave-propagation method which accurately takes into account the source term due to gravity is used for computing the non-hydrostatic atmospheric flows. It is demonstrated that there is little to no difference between the results obtained using the two different equation sets for Euler as well as Navier-Stokes solutions.

DKE: a fast numerical solver for the 3-D relativistic bounce-averaged electron drift kinetic equation

Energy Technology Data Exchange (ETDEWEB)

Decker, J.; Peysson, Y

2004-12-01

A new original code for solving the 3-D relativistic and bounce-averaged electron drift kinetic equation is presented. It designed for the current drive problem in tokamak with an arbitrary magnetic equilibrium. This tool allows self-consistent calculations of the bootstrap current in presence of other external current sources. RF current drive for arbitrary type of waves may be used. Several moments of the electron distribution function are determined, like the exact and effective fractions of trapped electrons, the plasma current, absorbed RF power, runaway and magnetic ripple loss rates and non-thermal Bremsstrahlung. Advanced numerical techniques have been used to make it the first fully implicit (reverse time) 3-D solver, particularly well designed for implementation in a chain of code for realistic current drive calculations in high {beta}{sub p} plasmas. All the details of the physics background and the numerical scheme are presented, as well a some examples to illustrate main code capabilities. Several important numerical points are addressed concerning code stability and potential numerical and physical limitations. (authors)

DKE: a fast numerical solver for the 3-D relativistic bounce-averaged electron drift kinetic equation

International Nuclear Information System (INIS)

Decker, J.; Peysson, Y.

2004-12-01

A new original code for solving the 3-D relativistic and bounce-averaged electron drift kinetic equation is presented. It designed for the current drive problem in tokamak with an arbitrary magnetic equilibrium. This tool allows self-consistent calculations of the bootstrap current in presence of other external current sources. RF current drive for arbitrary type of waves may be used. Several moments of the electron distribution function are determined, like the exact and effective fractions of trapped electrons, the plasma current, absorbed RF power, runaway and magnetic ripple loss rates and non-thermal Bremsstrahlung. Advanced numerical techniques have been used to make it the first fully implicit (reverse time) 3-D solver, particularly well designed for implementation in a chain of code for realistic current drive calculations in high β p plasmas. All the details of the physics background and the numerical scheme are presented, as well a some examples to illustrate main code capabilities. Several important numerical points are addressed concerning code stability and potential numerical and physical limitations. (authors)

Invariance analysis and conservation laws of the wave equation on ...

Indian Academy of Sciences (India)

in [7], the more interesting case being the latter since these lead to conservation laws via ... obtained and, hence, more conservation laws are classified. .... −2r2 sin θurt − 2r sin θut + 2r sin θ. (. 1 −. 2t r. ) ur + 2t sin θur. +r2 sin θ. (. 1 −. 2t r. ) urr + cos θuθ + sin θuθθ = 0,. (15) and then ¯X2 = u∂u + t∂t + r∂r leads to dt t. = dr.

Towards the ultimate variance-conserving convection scheme

International Nuclear Information System (INIS)

Os, J.J.A.M. van; Uittenbogaard, R.E.

2004-01-01

In the past various arguments have been used for applying kinetic energy-conserving advection schemes in numerical simulations of incompressible fluid flows. One argument is obeying the programmed dissipation by viscous stresses or by sub-grid stresses in Direct Numerical Simulation and Large Eddy Simulation, see e.g. [Phys. Fluids A 3 (7) (1991) 1766]. Another argument is that, according to e.g. [J. Comput. Phys. 6 (1970) 392; 1 (1966) 119], energy-conserving convection schemes are more stable i.e. by prohibiting a spurious blow-up of volume-integrated energy in a closed volume without external energy sources. In the above-mentioned references it is stated that nonlinear instability is due to spatial truncation rather than to time truncation and therefore these papers are mainly concerned with the spatial integration. In this paper we demonstrate that discretized temporal integration of a spatially variance-conserving convection scheme can induce non-energy conserving solutions. In this paper the conservation of the variance of a scalar property is taken as a simple model for the conservation of kinetic energy. In addition, the derivation and testing of a variance-conserving scheme allows for a clear definition of kinetic energy-conserving advection schemes for solving the Navier-Stokes equations. Consequently, we first derive and test a strictly variance-conserving space-time discretization for the convection term in the convection-diffusion equation. Our starting point is the variance-conserving spatial discretization of the convection operator presented by Piacsek and Williams [J. Comput. Phys. 6 (1970) 392]. In terms of its conservation properties, our variance-conserving scheme is compared to other spatially variance-conserving schemes as well as with the non-variance-conserving schemes applied in our shallow-water solver, see e.g. [Direct and Large-eddy Simulation Workshop IV, ERCOFTAC Series, Kluwer Academic Publishers, 2001, pp. 409-287

Structure preserving transformations for Newtonian Lie-admissible equations

International Nuclear Information System (INIS)

Cantrijn, F.

1979-01-01

Recently, a new formulation of non-conservative mechanics has been presented in terms of Hamilton-admissible equations which constitute a generalization of the conventional Hamilton equations. The algebraic structure entering the Hamilton-admissible description of a non-conservative system is that of a Lie-admissible algebra. The corresponding geometrical treatment is related to the existence of a so-called symplectic-admissible form. The transformation theory for Hamilton-admissible systems is currently investigated. The purpose of this paper is to describe one aspect of this theory by identifying the class of transformations which preserve the structure of Hamilton-admissible equations. Necessary and sufficient conditions are established for a transformation to be structure preserving. Some particular cases are discussed and an example is worked out

Introduction to partial differential equations

CERN Document Server

Borthwick, David

2016-01-01

This modern take on partial differential equations does not require knowledge beyond vector calculus and linear algebra. The author focuses on the most important classical partial differential equations, including conservation equations and their characteristics, the wave equation, the heat equation, function spaces, and Fourier series, drawing on tools from analysis only as they arise.Within each section the author creates a narrative that answers the five questions: (1) What is the scientific problem we are trying to understand? (2) How do we model that with PDE? (3) What techniques can we use to analyze the PDE? (4) How do those techniques apply to this equation? (5) What information or insight did we obtain by developing and analyzing the PDE? The text stresses the interplay between modeling and mathematical analysis, providing a thorough source of problems and an inspiration for the development of methods.

Asymptotic Conservation Laws in Classical Field Theory

International Nuclear Information System (INIS)

Anderson, I.M.; Torre, C.G.

1996-01-01

A new, general, field theoretic approach to the derivation of asymptotic conservation laws is presented. In this approach asymptotic conservation laws are constructed directly from the field equations according to a universal prescription which does not rely upon the existence of Noether identities or any Lagrangian or Hamiltonian formalisms. The resulting general expressions of the conservation laws enjoy important invariance properties and synthesize all known asymptotic conservation laws, such as the Arnowitt-Deser-Misner energy in general relativity. copyright 1996 The American Physical Society

Integrability of an extended (2+1)-dimensional shallow water wave equation with Bell polynomials

International Nuclear Information System (INIS)

Wang Yun-Hu; Chen Yong

2013-01-01

We investigate the extended (2+1)-dimensional shallow water wave equation. The binary Bell polynomials are used to construct bilinear equation, bilinear Bäcklund transformation, Lax pair, and Darboux covariant Lax pair for this equation. Moreover, the infinite conservation laws of this equation are found by using its Lax pair. All conserved densities and fluxes are given with explicit recursion formulas. The N-soliton solutions are also presented by means of the Hirota bilinear method. (general)

High Weak Order Methods for Stochastic Differential Equations Based on Modified Equations

KAUST Repository

Abdulle, Assyr

2012-01-01

© 2012 Society for Industrial and Applied Mathematics. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new methods of weak order two, in particular, semi-implicit integrators well suited for stiff (meansquare stable) stochastic problems, and implicit integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of our methodology.

Bianchi-Baecklund transformations, conservation laws, and linearization of various field theories

International Nuclear Information System (INIS)

Chau Wang, L.L.

1980-01-01

The discussion includes: the Sine-Gordon equation, parametric Bianchi-Baecklund transformations and the derivation of local conservation laws; chiral fields, parametric Bianchi-Baecklund transformations, local and non-local conservation laws, and linearization; super chiral fields, a parallel development similar to the chiral field; and self-dual Yang-Mills fields in 4-dimensional Euclidean space; loop-cpace chiral equations, parallel development but with subtlety

Symmetric truncations of the shallow-water equations

International Nuclear Information System (INIS)

Rouhi, A.; Abarbanel, H.D.I.

1993-01-01

Conservation of potential vorticity in Eulerian fluids reflects particle interchange symmetry in the Lagrangian fluid version of the same theory. The algebra associated with this symmetry in the shallow-water equations is studied here, and we give a method for truncating the degrees of freedom of the theory which preserves a maximal number of invariants associated with this algebra. The finite-dimensional symmetry associated with keeping only N modes of the shallow-water flow is SU(N). In the limit where the number of modes goes to infinity (N→∞) all the conservation laws connected with potential vorticity conservation are recovered. We also present a Hamiltonian which is invariant under this truncated symmetry and which reduces to the familiar shallow-water Hamiltonian when N→∞. All this provides a finite-dimensional framework for numerical work with the shallow-water equations which preserves not only energy and enstrophy but all other known conserved quantities consistent with the finite number of degrees of freedom. The extension of these ideas to other nearly two-dimensional flows is discussed

Making conservation work for everyone

Energy Technology Data Exchange (ETDEWEB)

Wiersma, J. [Veridian Corp., Ajax, ON (Canada)

2004-07-01

This presentation discussed the economic value of conservation, the optimal deployment of energy conservation. A sample load profile was presented to demonstrate how much electricity the average residential customer uses on a summer day. The average customer does not have the tools to understand the financial consequences of conservation for different types of equipment at different times of the day. Smart metering technology could help in this regard. Accurate unsubsidized prices are also considered to be the best incentive to conserve because customers will reduce electricity use when the prices are high. It was also suggested that standards for new appliances should be increased effectively to their economic value. The enablers to energy conservation include solid consumer education programs, real time metering in places where it is cost effective, real time pricing in places where it is practical, and power rates that reflect real costs. Barriers to energy conservation include the residual economic advantage that may be insufficient to justify investment; support from local distribution companies and transmission companies if the lost revenue adjustment mechanism (LRAM) is not sufficient to recover lost revenue and if LDCs are not sufficiently involved in the design of the electricity conservation program. 7 figs.

On the method of inverse scattering problem and Baecklund transformations for supersymmetric equations

International Nuclear Information System (INIS)

Chaichian, M.; Kulish, P. P.

1978-04-01

Supersymmetric Liouville and sine-Gordon equations are studied. We write down for these models the system of linear equations for which the method of inverse scattering problem should be applicable. Expressions for an infinite set of conserved currents are explicitly given. Supersymmetric Baecklund transformations and generalized conservation laws are constructed. (author)

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

International Nuclear Information System (INIS)

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

2003-01-01

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

Combined Sinh-Cosh-Gordon equation: Symmetry reductions, exact ...

African Journals Online (AJOL)

Combined Sinh-Cosh-Gordon equation: Symmetry reductions, exact solutions and conservation laws. ... In this paper we study the combined sinh-cosh-Gordon equation, which arises in mathematical physics and has a wide range of scientific applications that range from chemical reactions to water surface gravity waves.

Complex nonlinear Lagrangian for the Hasegawa-Mima equation

International Nuclear Information System (INIS)

Dewar, R.L.; Abdullatif, R.F.; Sangeetha, G.G.

2005-01-01

The Hasegawa-Mima equation is the simplest nonlinear single-field model equation that captures the essence of drift wave dynamics. Like the Schroedinger equation it is first order in time. However its coefficients are real, so if the potential φ is initially real it remains real. However, by embedding φ in the space of complex functions a simple Lagrangian is found from which the Hasegawa-Mima equation may be derived from Hamilton's Principle. This Lagrangian is used to derive an action conservation equation which agrees with that of Biskamp and Horton. (author)

Symbolic Detection of Permutation and Parity Symmetries of Evolution Equations

KAUST Repository

Alghamdi, Moataz

2017-06-18

We introduce a symbolic computational approach to detecting all permutation and parity symmetries in any general evolution equation, and to generating associated invariant polynomials, from given monomials, under the action of these symmetries. Traditionally, discrete point symmetries of differential equations are systemically found by solving complicated nonlinear systems of partial differential equations; in the presence of Lie symmetries, the process can be simplified further. Here, we show how to find parity- and permutation-type discrete symmetries purely based on algebraic calculations. Furthermore, we show that such symmetries always form groups, thereby allowing for the generation of new group-invariant conserved quantities from known conserved quantities. This work also contains an implementation of the said results in Mathematica. In addition, it includes, as a motivation for this work, an investigation of the connection between variational symmetries, described by local Lie groups, and conserved quantities in Hamiltonian systems.

Rotational averaging of multiphoton absorption cross sections

Energy Technology Data Exchange (ETDEWEB)

Friese, Daniel H., E-mail: daniel.h.friese@uit.no; Beerepoot, Maarten T. P.; Ruud, Kenneth [Centre for Theoretical and Computational Chemistry, University of Tromsø — The Arctic University of Norway, N-9037 Tromsø (Norway)

2014-11-28

Rotational averaging of tensors is a crucial step in the calculation of molecular properties in isotropic media. We present a scheme for the rotational averaging of multiphoton absorption cross sections. We extend existing literature on rotational averaging to even-rank tensors of arbitrary order and derive equations that require only the number of photons as input. In particular, we derive the first explicit expressions for the rotational average of five-, six-, and seven-photon absorption cross sections. This work is one of the required steps in making the calculation of these higher-order absorption properties possible. The results can be applied to any even-rank tensor provided linearly polarized light is used.

Kinks and the Dirac equation

International Nuclear Information System (INIS)

Skyrme, T.H.R.

1994-01-01

In a model quantum theory of interacting mesons, the motion of certain conserved particle-like structures is discussed. It is shown how collective coordinates may be introduced to describe them, leading, in lowest approximation, to a Dirac equation. (author)

On the non-stationary generalized Langevin equation

Science.gov (United States)

Meyer, Hugues; Voigtmann, Thomas; Schilling, Tanja

2017-12-01

In molecular dynamics simulations and single molecule experiments, observables are usually measured along dynamic trajectories and then averaged over an ensemble ("bundle") of trajectories. Under stationary conditions, the time-evolution of such averages is described by the generalized Langevin equation. By contrast, if the dynamics is not stationary, it is not a priori clear which form the equation of motion for an averaged observable has. We employ the formalism of time-dependent projection operator techniques to derive the equation of motion for a non-equilibrium trajectory-averaged observable as well as for its non-stationary auto-correlation function. The equation is similar in structure to the generalized Langevin equation but exhibits a time-dependent memory kernel as well as a fluctuating force that implicitly depends on the initial conditions of the process. We also derive a relation between this memory kernel and the autocorrelation function of the fluctuating force that has a structure similar to a fluctuation-dissipation relation. In addition, we show how the choice of the projection operator allows us to relate the Taylor expansion of the memory kernel to data that are accessible in MD simulations and experiments, thus allowing us to construct the equation of motion. As a numerical example, the procedure is applied to Brownian motion initialized in non-equilibrium conditions and is shown to be consistent with direct measurements from simulations.

Statistical Average of Spin Operators for Calculation of Three-Component Magnetization (II): Solution of Equation

International Nuclear Information System (INIS)

Wang Huaiyu; Long Yao; Chen Nanxian

2006-01-01

In this paper, the solution of Chebyshev equation with its argument being greater than 1 is obtained. The initial value of the derivative of the solution is the expression of magnetization, which is valid for any spin quantum number S. The Chebyshev equation is transformed from an ordinary differential equation obtained when we dealt with Heisenberg model, in order to calculate all three components of magnetization, by many-body Green's function under random phase approximation. The Chebyshev functions with argument being greater than 1 are discussed. This paper shows that the Chebyshev polynomials with their argument being greater than 1 have their physical application.

Layer-Mean Quantities, Local Conservation Laws, and Vorticity

International Nuclear Information System (INIS)

Camassa, R.; Levermore, C.D.

1997-01-01

We derive local conservation laws for layer-mean quantities in two general settings. When applied to Euler flows, the first of these settings yields well-known local conservation laws for quantities averaged between material surfaces. The second, however, leads to new local conservation laws for quantities involving the vorticity that are averaged between arbitrary surfaces. These produce the crucial vorticity conservation laws in shallow water models that admit nonhydrostatic and noncolumnar motion. Moreover, they seem to lie outside the Hamiltonian paradigm of fluid dynamics. The formalism generalizes to skew-symmetric matrix fields; applications to electromagnetism are suggested. copyright 1997 The American Physical Society

Equations of radiation hydrodynamics

International Nuclear Information System (INIS)

Mihalas, D.

1982-01-01

The purpose of this paper is to give an overview of the role of radiation in the transport of energy and momentum in a combined matter-radiation fluid. The transport equation for a moving radiating fluid is presented in both a fully Eulerian and a fully Lagrangian formulation, along with conservation equations describing the dynamics of the fluid. Special attention is paid to the problem of deriving equations that are mutually consistent in each frame, and between frames, to 0(v/c). A detailed analysis is made to show that in situations of broad interest, terms that are formally of 0(v/c) actually dominate the solution, demonstrating that it is esential (1) to pay scrupulous attention to the question of the frame dependence in formulating the equations; and (2) to solve the equations to 0(v/c) in quite general circumstances. These points are illustrated in the context of the nonequilibrium radiation diffusion limit, and a sketch of how the Lagrangian equations are to be solved will be presented

Gauge-invariant flow equation

Science.gov (United States)

Wetterich, C.

2018-06-01

We propose a closed gauge-invariant functional flow equation for Yang-Mills theories and quantum gravity that only involves one macroscopic gauge field or metric. It is based on a projection on physical and gauge fluctuations. Deriving this equation from a functional integral we employ the freedom in the precise choice of the macroscopic field and the effective average action in order to realize a closed and simple form of the flow equation.

Development of Multi-Dimensional RELAP5 with Conservative Momentum Flux

Energy Technology Data Exchange (ETDEWEB)

Jang, Hyung Wook; Lee, Sang Yong [KINGS, Ulsan (Korea, Republic of)

2016-10-15

The non-conservative form of the momentum equations are used in many codes. It tells us that using the non-conservative form in the non-porous or open body problem may not be good. In this paper, two aspects concerning the multi-dimensional codes will be discussed. Once the validity of the modified code is confirmed, it is applied to the analysis of the large break LOCA for APR-1400. One of them is the properness of the type of the momentum equations. The other discussion will be the implementation of the conservative momentum flux term in RELAP5. From the present study and former, it is shown that the RELAP5 Multi-D with conservative convective terms is applicable to LOCA analysis. And the implementation of the conservative convective terms in RELAP5 seems to be successful. Further efforts have to be made on making it more robust.

Diesel conservation: GSRTC'S experience

Energy Technology Data Exchange (ETDEWEB)

Ramesh Kumar, I V

1980-01-01

The Gujarat State Road Transport Corporation (GSRTC) in India has a fleet of about 6000 buses. The increasing cost of fuel and lubricants added to uncertainty in supplies, has necessitated the need for conserving High Speed Diesel Oil (HSD). GSRTC had achieved an overall average Kilometre Per Litre (kmpl) of 4.44 in the year 1976-1977 due to a variety of measures. In the year 1978-1979 the average kmpl was 4.52 and it is expected to be 4.60 for 1979-1980. The case study outlined describes the measures taken by GSRTC in conserving high speed diesel oil by various methods.

The symmetries and conservation laws of some Gordon-type

Indian Academy of Sciences (India)

Conservation laws; Milne space-time; Gordon-type equations. Abstract. In this letter, the Lie point symmetries of a class of Gordon-type wave equations that arise in the Milne space-time are presented ... Pramana – Journal of Physics | News.

Generalized equations of gravitational field

International Nuclear Information System (INIS)

Stanyukovich, K.P.; Borisova, L.B.

1985-01-01

Equations for gravitational fields are obtained on the basis of a generalized Lagrangian Z=f(R) (R is the scalar curvature). Such an approach permits to take into account the evolution of a gravitation ''constant''. An expression for the force Fsub(i) versus the field variability is obtained. Conservation laws are formulated differing from the standard ones by the fact that in the right part of new equations the value Fsub(i) is present that goes to zero at an ultimate passage to the standard Einstein theory. An equation of state is derived for cosmological metrics for a particular case, f=bRsup(1+α) (b=const, α=const)

Multi-symplectic Preissmann methods for generalized Zakharov-Kuznetsov equation

International Nuclear Information System (INIS)

Wang Junjie; Yang Kuande; Wang Liantang

2012-01-01

Generalized Zakharov-Kuznetsov equation, a typical nonlinear wave equation, was studied based on the multi-symplectic theory in Hamilton space. The multi-symplectic formulations of generalized Zakharov-Kuznetsov equation with several conservation laws are presented. The multi-symplectic Preissmann method is used to discretize the formulations. The numerical experiment is given, and the results verify the efficiency of the multi-symplectic scheme. (authors)

An energy-stable convex splitting for the phase-field crystal equation

KAUST Repository

Vignal, P.; Dalcin, L.; Brown, D. L.; Collier, N.; Calo, V. M.

2015-01-01

Abstract The phase-field crystal equation, a parabolic, sixth-order and nonlinear partial differential equation, has generated considerable interest as a possible solution to problems arising in molecular dynamics. Nonetheless, solving this equation is not a trivial task, as energy dissipation and mass conservation need to be verified for the numerical solution to be valid. This work addresses these issues, and proposes a novel algorithm that guarantees mass conservation, unconditional energy stability and second-order accuracy in time. Numerical results validating our proofs are presented, and two and three dimensional simulations involving crystal growth are shown, highlighting the robustness of the method. © 2015 Elsevier Ltd.

An energy-stable convex splitting for the phase-field crystal equation

KAUST Repository

Vignal, P.

2015-10-01

Abstract The phase-field crystal equation, a parabolic, sixth-order and nonlinear partial differential equation, has generated considerable interest as a possible solution to problems arising in molecular dynamics. Nonetheless, solving this equation is not a trivial task, as energy dissipation and mass conservation need to be verified for the numerical solution to be valid. This work addresses these issues, and proposes a novel algorithm that guarantees mass conservation, unconditional energy stability and second-order accuracy in time. Numerical results validating our proofs are presented, and two and three dimensional simulations involving crystal growth are shown, highlighting the robustness of the method. © 2015 Elsevier Ltd.

Gyrofluid potential vorticity equation and turbulent equipartion states

DEFF Research Database (Denmark)

Madsen, Jens; Juul Rasmussen, Jens; Naulin, Volker

2015-01-01

. The equation is relevant for transport barriers in magnetically confined plasmas because particle density, ion temperature and the radial electric field are mutually coupled through the potential vorticity. The potential vorticity equation is derived from an energy conserving, four-field, electrostatic, full......An equation governing potential vorticity in a magnetized plasmas is derived. The equation is analogous to Ertel's theorem. In the long wave-length limit the potential vorticity equals the ratio of the gyro-frequency plus the E × B- and diamagnetic polarization densities to the particle density...

Integrable peakon equations with cubic nonlinearity

International Nuclear Information System (INIS)

Hone, Andrew N W; Wang, J P

2008-01-01

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

Properties of bright solitons in averaged and unaveraged models for SDG fibres

Science.gov (United States)

Kumar, Ajit; Kumar, Atul

1996-04-01

Using the slowly varying envelope approximation and averaging over the fibre cross-section the evolution equation for optical pulses in semiconductor-doped glass (SDG) fibres is derived from the nonlinear wave equation. Bright soliton solutions of this equation are obtained numerically and their properties are studied and compared with those of the bright solitons in the unaveraged model.

Adaptive integral equation methods in transport theory

International Nuclear Information System (INIS)

Kelley, C.T.

1992-01-01

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

Derivation of Grad’s Thirteen Regularized Moment Equations Using a Hermite Polynomial Representation of Velocity Distribution Function (Preprint)

Science.gov (United States)

2010-06-16

B4) Substituting tui  / and tVT  /2 from the momentum and energy conservation law equations, Eqs...B9) Substituting tui  / and tVT  /2 from the momentum and energy conservation law equations, Eqs. (15...Substituting tui  / and tVT  /2 from the momentum and energy conservation law equations, Eqs. (15) and (16), into Eq. (B13) and then dropping all

Discrete conservation properties for shallow water flows using mixed mimetic spectral elements

Science.gov (United States)

Lee, D.; Palha, A.; Gerritsma, M.

2018-03-01

A mixed mimetic spectral element method is applied to solve the rotating shallow water equations. The mixed method uses the recently developed spectral element histopolation functions, which exactly satisfy the fundamental theorem of calculus with respect to the standard Lagrange basis functions in one dimension. These are used to construct tensor product solution spaces which satisfy the generalized Stokes theorem, as well as the annihilation of the gradient operator by the curl and the curl by the divergence. This allows for the exact conservation of first order moments (mass, vorticity), as well as higher moments (energy, potential enstrophy), subject to the truncation error of the time stepping scheme. The continuity equation is solved in the strong form, such that mass conservation holds point wise, while the momentum equation is solved in the weak form such that vorticity is globally conserved. While mass, vorticity and energy conservation hold for any quadrature rule, potential enstrophy conservation is dependent on exact spatial integration. The method possesses a weak form statement of geostrophic balance due to the compatible nature of the solution spaces and arbitrarily high order spatial error convergence.

Conservation laws for steady flow and solitons in a multifluid plasma revisited

International Nuclear Information System (INIS)

Mace, R. L.; McKenzie, J. F.; Webb, G. M.

2007-01-01

The conservation laws used in constructing the governing equations for planar solitons in multifluid plasmas are revisited. In particular, the concept of generalized vorticity facilitates the derivation of some general ''Bernoulli theorems,'' which reduce, in specific instances, to conservation laws previously deduced by other means. These theorems clarify the underlying physical principles that give rise to the conserved quantities. As an example of the usefulness of the techniques, even for relatively simple flows and progressive waves, the equations governing stationary nonlinear whistler waves propagating parallel to an ambient magnetic field are derived using generalized vorticity concepts

Study of nonequilibrium dispersed two phase flow

International Nuclear Information System (INIS)

Reyes, J.N. Jr.

1986-01-01

Understanding the behavior of liquid droplets in a superheated steam environment is essential to the accurate prediction of nuclear fuel rod surface temperatures during the blowdown and reflood phase of a loss-of-coolant-accident (LOCA). In response to this need, this treatise presents several original and significant contributions to the field of thermofluid physics. The research contained herein presents a statistical derivation of the two-phase mass, momentum, and energy-conservation equations using a droplet continuity equation analogous to that used in the Kinetic Theory of Gases. Unlike the Eulerian volume and time-averaged conservation equations generally used to describe dispersed two-phase flow behavior, this statistical averaging approach results in an additional mass momentum or energy term in each of the respective conservation equations. Further, this study demonstrates that current definitions of the volumetric vapor generation rate used in the mass conservation equation are inappropriate results under certain circumstances. The mass conservation equation derived herein is used to obtain a new definition for the volumetric vapor-generation rate. Last, a simple two phase phenomenological model, based on the statistically averaged conservation equations, is presented and solved analytically. It is shown that the actual quality and vapor temperature, under these circumstances, depend on a single dimensionless group

Parameterization of Time-Averaged Suspended Sediment Concentration in the Nearshore

Directory of Open Access Journals (Sweden)

Hyun-Doug Yoon

2015-11-01

Full Text Available To quantify the effect of wave breaking turbulence on sediment transport in the nearshore, the vertical distribution of time-averaged suspended sediment concentration (SSC in the surf zone was parameterized in terms of the turbulent kinetic energy (TKE at different cross-shore locations, including the bar crest, bar trough, and inner surf zone. Using data from a large-scale laboratory experiment, a simple relationship was developed between the time-averaged SSC and the time-averaged TKE. The vertical variation of the time-averaged SSC was fitted to an equation analogous to the turbulent dissipation rate term. At the bar crest, the proposed equation was slightly modified to incorporate the effect of near-bed sediment processes and yielded reasonable agreement. This parameterization yielded the best agreement at the bar trough, with a coefficient of determination R2 ≥ 0.72 above the bottom boundary layer. The time-averaged SSC in the inner surf zone showed good agreement near the bed but poor agreement near the water surface, suggesting that there is a different sedimentation mechanism that controls the SSC in the inner surf zone.

Regularity and energy conservation for the compressible Euler equations

Czech Academy of Sciences Publication Activity Database

Feireisl, Eduard; Gwiazda, P.; Swierczewska-Gwiazda, A.; Wiedemann, E.

2017-01-01

Roč. 223, č. 3 (2017), s. 1375-1395 ISSN 0003-9527 EU Projects: European Commission(XE) 320078 - MATHEF Institutional support: RVO:67985840 Keywords : compressible Euler equations Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 2.392, year: 2016 http://link.springer.com/article/10.1007%2Fs00205-016-1060-5

Variation estimation of the averaged cross sections in the direct and adjoint fluxes

International Nuclear Information System (INIS)

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

1995-01-01

There are several applications of the perturbation theory to specifics problems of reactor physics, such as nonuniform fuel burnup, nonuniform poison accumulation and evaluations of Doppler effects on reactivity. The neutron fluxes obtained from the solutions of direct and adjoint diffusion equations, are used in these applications. In the adjoint diffusion equation has been used the group constants averaged in the energy-dependent direct neutron flux, that it is not theoretically consistent. In this paper it is presented a method to calculate the energy-dependent adjoint neutron flux, to obtain the average group-constant that will be used in the adjoint diffusion equation. The method is based on the solution of the adjoint neutron balance equations, that were derived for a two regions cell. (author). 5 refs, 2 figs, 1 tab

Balance equations for a relativistic plasma. Pt. 1

International Nuclear Information System (INIS)

Hebenstreit, H.

1983-01-01

Relativistic power moments of the four-momentum are decomposed according to a macroscopic four-velocity. The thus obtained quantities are identified as relativistic generalization of the nonrelativistic orthogonal moments, e.g. diffusion flow, heat flow, pressure, etc. From the relativistic Boltzmann equation we then derive balance equations for these quantities. Explicit expressions for the relativistic mass conservation, energy balance, pressure balance, heat flow balance are presented. The weak relativistic limit is discussed. The derivation of higher order balance equations is sketched. (orig.)

Michaelis - Menten equation for degradation of insoluble substrate

DEFF Research Database (Denmark)

Andersen, Morten; Kari, Jeppe; Borch, Kim

2017-01-01

substrate it is difficult to assess whether the requirement of the MM equation is met. In this paper we study a simple kinetic model, where removal of attack sites expose new ones which preserve the total accessible substrate, and denote this approach the substrate conserving model. The kinetic equations...... are solved in closed form, both steady states and progress curves, for any admissible values of initial conditions and rate constants. The model is shown to merge with the MM equation and the reverse MM equation when these are valid. The relation between available molar concentration of attack sites and mass...

Approach in Theory of Nonlinear Evolution Equations: The Vakhnenko-Parkes Equation

Directory of Open Access Journals (Sweden)

V. O. Vakhnenko

2016-01-01

Full Text Available A variety of methods for examining the properties and solutions of nonlinear evolution equations are explored by using the Vakhnenko equation (VE as an example. The VE, which arises in modelling the propagation of high-frequency waves in a relaxing medium, has periodic and solitary traveling wave solutions some of which are loop-like in nature. The VE can be written in an alternative form, known as the Vakhnenko-Parkes equation (VPE, by a change of independent variables. The VPE has an N-soliton solution which is discussed in detail. Individual solitons are hump-like in nature whereas the corresponding solution to the VE comprises N-loop-like solitons. Aspects of the inverse scattering transform (IST method, as applied originally to the KdV equation, are used to find one- and two-soliton solutions to the VPE even though the VPE’s spectral equation is third-order and not second-order. A Bäcklund transformation for the VPE is used to construct conservation laws. The standard IST method for third-order spectral problems is used to investigate solutions corresponding to bound states of the spectrum and to a continuous spectrum. This leads to N-soliton solutions and M-mode periodic solutions, respectively. Interactions between these types of solutions are investigated.

Wave equations on anti self dual (ASD) manifolds

Science.gov (United States)

Bashingwa, Jean-Juste; Kara, A. H.

2017-11-01

In this paper, we study and perform analyses of the wave equation on some manifolds with non diagonal metric g_{ij} which are of neutral signatures. These include the invariance properties, variational symmetries and conservation laws. In the recent past, wave equations on the standard (space time) Lorentzian manifolds have been performed but not on the manifolds from metrics of neutral signatures.

Neoclassical MHD equations for tokamaks

International Nuclear Information System (INIS)

Callen, J.D.; Shaing, K.C.

1986-03-01

The moment equation approach to neoclassical-type processes is used to derive the flows, currents and resistive MHD-like equations for studying equilibria and instabilities in axisymmetric tokamak plasmas operating in the banana-plateau collisionality regime (ν* approx. 1). The resultant ''neoclassical MHD'' equations differ from the usual reduced equations of resistive MHD primarily by the addition of the important viscous relaxation effects within a magnetic flux surface. The primary effects of the parallel (poloidal) viscous relaxation are: (1) Rapid (approx. ν/sub i/) damping of the poloidal ion flow so the residual flow is only toroidal; (2) addition of the bootstrap current contribution to Ohm's laws; and (3) an enhanced (by B 2 /B/sub theta/ 2 ) polarization drift type term and consequent enhancement of the perpendicular dielectric constant due to parallel flow inertia, which causes the equations to depend only on the poloidal magnetic field B/sub theta/. Gyroviscosity (or diamagnetic vfiscosity) effects are included to properly treat the diamagnetic flow effects. The nonlinear form of the neoclassical MHD equations is derived and shown to satisfy an energy conservation equation with dissipation arising from Joule and poloidal viscous heating, and transport due to classical and neoclassical diffusion

Parametrized post-Newtonian approximation and Rastall's gravitational field equations

International Nuclear Information System (INIS)

Smalley, L.L.

1978-01-01

The parametrized post-Newtonian (PPN) approximation is generalized to accomodate Rastall's modification of Einstein's theory of gravity, which allows nonzero divergence of the energy-momentum tensor. Rastall's theory is then shown to have consistent field equations, gauge conditions, and the correct Newtonian limit of the equations of motion. The PPN parameters are obtained and shown to agree experimentally with those for the Einstein theory. In light of the nonzero divergence condition, integral conservation laws are investigated and shown to yield conserved energy-momentum and angular-momentum. We conclude that the above generalization of metric theories, within the PPN framework, is a natural extension of the concept of metric theories

Nonstandard conserved Hamiltonian structures in dissipative/damped systems: Nonlinear generalizations of damped harmonic oscillator

International Nuclear Information System (INIS)

Pradeep, R. Gladwin; Chandrasekar, V. K.; Senthilvelan, M.; Lakshmanan, M.

2009-01-01

In this paper we point out the existence of a remarkable nonlocal transformation between the damped harmonic oscillator and a modified Emden-type nonlinear oscillator equation with linear forcing, xe+αxx+βx 3 +γx=0, which preserves the form of the time independent integral, conservative Hamiltonian, and the equation of motion. Generalizing this transformation we prove the existence of nonstandard conservative Hamiltonian structure for a general class of damped nonlinear oscillators including Lienard-type systems. Further, using the above Hamiltonian structure for a specific example, namely, the generalized modified Emden equation xe+αx q x+βx 2q+1 =0, where α, β, and q are arbitrary parameters, the general solution is obtained through appropriate canonical transformations. We also present the conservative Hamiltonian structure of the damped Mathews-Lakshmanan oscillator equation. The associated Lagrangian description for all the above systems is also briefly discussed.

A new multi-symplectic scheme for the generalized Kadomtsev-Petviashvili equation

Science.gov (United States)

Li, Haochen; Sun, Jianqiang

2012-09-01

We propose a new scheme for the generalized Kadomtsev-Petviashvili (KP) equation. The multi-symplectic conservation property of the new scheme is proved. Back error analysis shows that the new multi-symplectic scheme has second order accuracy in space and time. Numerical application on studying the KPI equation and the KPII equation are presented in detail.

Out-of-equilibrium quantum fields with conserved charge

International Nuclear Information System (INIS)

Bedingham, D.J.

2004-01-01

We study the out-of-equilibrium evolution of an O(2)-invariant scalar field in which a conserved charge is stored. We apply a loop expansion of the 2-particle irreducible effective action to 3-loop order. Equations of motion are derived which conserve both total charge and total energy yet allow for the effects of scattering whereby charge and energy can transfer between modes. Working in 1+1 dimensions we solve the equations of motion numerically for a system knocked out of equilibrium by a sudden temperature quench. We examine the initial stages of the charge and energy redistribution. This provides a basis from which we can understand the formation of Bose-Einstein condensates from first principles

Vapor-droplet flow equations

International Nuclear Information System (INIS)

Crowe, C.T.

1975-01-01

General features of a vapor-droplet flow are discussed and the equations expressing the conservation of mass, momentum, and energy for the vapor, liquid, and mixture using the control volume approach are derived. The phenomenological laws describing the exchange of mass, momentum, and energy between phases are also reviewed. The results have application to development of water-dominated geothermal resources

A stochastic Galerkin method for the Euler equations with Roe variable transformation

KAUST Repository

Pettersson, Per; Iaccarino, Gianluca; Nordströ m, Jan

2014-01-01

The Euler equations subject to uncertainty in the initial and boundary conditions are investigated via the stochastic Galerkin approach. We present a new fully intrusive method based on a variable transformation of the continuous equations. Roe variables are employed to get quadratic dependence in the flux function and a well-defined Roe average matrix that can be determined without matrix inversion.In previous formulations based on generalized polynomial chaos expansion of the physical variables, the need to introduce stochastic expansions of inverse quantities, or square roots of stochastic quantities of interest, adds to the number of possible different ways to approximate the original stochastic problem. We present a method where the square roots occur in the choice of variables, resulting in an unambiguous problem formulation.The Roe formulation saves computational cost compared to the formulation based on expansion of conservative variables. Moreover, the Roe formulation is more robust and can handle cases of supersonic flow, for which the conservative variable formulation fails to produce a bounded solution. For certain stochastic basis functions, the proposed method can be made more effective and well-conditioned. This leads to increased robustness for both choices of variables. We use a multi-wavelet basis that can be chosen to include a large number of resolution levels to handle more extreme cases (e.g. strong discontinuities) in a robust way. For smooth cases, the order of the polynomial representation can be increased for increased accuracy. © 2013 Elsevier Inc.

Preservation of support and positivity for solutions of degenerate evolution equations

International Nuclear Information System (INIS)

Ambrose, David M; Wright, J Douglas

2010-01-01

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

Scale dependence of the average potential around the maximum in Φ4 theories

International Nuclear Information System (INIS)

Tetradis, N.; Wetterich, C.

1992-04-01

The average potential describes the physics at a length scale k - 1 by averaging out the degrees of freedom with characteristic moments larger than k. The dependence on k can be described by differential evolution equations. We solve these equations for the nonconvex part of the potential around the origin in φ 4 theories, in the phase with spontaneous symmetry breaking. The average potential is real and approaches the convex effective potential in the limit k → 0. Our calculation is relevant for processes for which the shape of the potential at a given scale is important, such as tunneling phenomena or inflation. (orig.)

Wave function collapse implies divergence of average displacement

OpenAIRE

Marchewka, A.; Schuss, Z.

2005-01-01

We show that propagating a truncated discontinuous wave function by Schr\\"odinger's equation, as asserted by the collapse axiom, gives rise to non-existence of the average displacement of the particle on the line. It also implies that there is no Zeno effect. On the other hand, if the truncation is done so that the reduced wave function is continuous, the average coordinate is finite and there is a Zeno effect. Therefore the collapse axiom of measurement needs to be revised.

Coupled force-balance and particle-occupation rate equations for high-field electron transport

International Nuclear Information System (INIS)

Lei, X. L.

2008-01-01

It is pointed out that in the framework of balance-equation approach, the coupled force-balance and particle-occupation rate equations can be used as a complete set of equations to determine the high-field transport of semiconductors in both strong and weak electron-electron interaction limits. We call to attention that the occupation rate equation conserves the total particle number and maintains the energy balance of the relative electron system, and there is no need to introduce any other term in it. The addition of an energy-drift term in the particle-occupation rate equation [Phys. Rev. B 71, 195205 (2005)] is physically inadequate for the violation of the total particle-number conservation and the energy balance. It may lead to a substantial unphysical increase of the total particle number by the application of a dc electric field

Conservation laws for multidimensional systems and related linear algebra problems

International Nuclear Information System (INIS)

Igonin, Sergei

2002-01-01

We consider multidimensional systems of PDEs of generalized evolution form with t-derivatives of arbitrary order on the left-hand side and with the right-hand side dependent on lower order t-derivatives and arbitrary space derivatives. For such systems we find an explicit necessary condition for the existence of higher conservation laws in terms of the system's symbol. For systems that violate this condition we give an effective upper bound on the order of conservation laws. Using this result, we completely describe conservation laws for viscous transonic equations, for the Brusselator model and the Belousov-Zhabotinskii system. To achieve this, we solve over an arbitrary field the matrix equations SA=A t S and SA=-A t S for a quadratic matrix A and its transpose A t , which may be of independent interest

Modification of the FEM3 model to ensure mass conservation

International Nuclear Information System (INIS)

Gresho, P.M.

1987-01-01

The problem of global mass conservation (lack thereof) in the current anelastic equations solved by FEM3 is described and its cause explained. The additional equations necessary to solve the problem are presented and methods for their incorporation into the current code are suggested. 14 refs

Optimal bounds and extremal trajectories for time averages in nonlinear dynamical systems

Science.gov (United States)

Tobasco, Ian; Goluskin, David; Doering, Charles R.

2018-02-01

For any quantity of interest in a system governed by ordinary differential equations, it is natural to seek the largest (or smallest) long-time average among solution trajectories, as well as the extremal trajectories themselves. Upper bounds on time averages can be proved a priori using auxiliary functions, the optimal choice of which is a convex optimization problem. We prove that the problems of finding maximal trajectories and minimal auxiliary functions are strongly dual. Thus, auxiliary functions provide arbitrarily sharp upper bounds on time averages. Moreover, any nearly minimal auxiliary function provides phase space volumes in which all nearly maximal trajectories are guaranteed to lie. For polynomial equations, auxiliary functions can be constructed by semidefinite programming, which we illustrate using the Lorenz system.

Average beta measurement in EXTRAP T1

International Nuclear Information System (INIS)

Hedin, E.R.

1988-12-01

Beginning with the ideal MHD pressure balance equation, an expression for the average poloidal beta, Β Θ , is derived. A method for unobtrusively measuring the quantities used to evaluate Β Θ in Extrap T1 is described. The results if a series of measurements yielding Β Θ as a function of externally applied toroidal field are presented. (author)

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

Number-conserving random phase approximation with analytically integrated matrix elements

International Nuclear Information System (INIS)

Kyotoku, M.; Schmid, K.W.; Gruemmer, F.; Faessler, A.

1990-01-01

In the present paper a number conserving random phase approximation is derived as a special case of the recently developed random phase approximation in general symmetry projected quasiparticle mean fields. All the occurring integrals induced by the number projection are performed analytically after writing the various overlap and energy matrices in the random phase approximation equation as polynomials in the gauge angle. In the limit of a large number of particles the well-known pairing vibration matrix elements are recovered. We also present a new analytically number projected variational equation for the number conserving pairing problem

Current conserving theory at the operator level

Science.gov (United States)

Yuan, Jiangtao; Wang, Yin; Wang, Jian

The basic assumption of quantum transport in mesoscopic systems is that the total charge inside the scattering region is zero. This means that the potential deep inside reservoirs is effectively screened and therefore the electric field at interface of scattering region is zero. Thus the current conservation condition can be satisfied automatically which is an important condition in mesoscopic transport. So far the current conserving ac theory is well developed by considering the displacement current which is due to Coulomb interaction if we just focus on the average current. However, the frequency dependent shot noise does not satisfy the conservation condition since we do not consider the current conservation at the operator level. In this work, we formulate a generalized current conserving theory at the operator level using non-equilibrium Green's function theory which could be applied to both average current and frequency dependent shot noise. A displacement operator is derived for the first time so that the frequency dependent correlation of displacement currents could be investigated. Moreover, the equilibrium shot noise is investigated and a generalized fluctuation-dissipation relationship is presented.

Equations for estimating bankfull channel geometry and discharge for streams in Massachusetts

Science.gov (United States)

Bent, Gardner C.; Waite, Andrew M.

2013-01-01

Regression equations were developed for estimating bankfull geometry—width, mean depth, cross-sectional area—and discharge for streams in Massachusetts. The equations provide water-resource and conservation managers with methods for estimating bankfull characteristics at specific stream sites in Massachusetts. This information can be used for the adminstration of the Commonwealth of Massachusetts Rivers Protection Act of 1996, which establishes a protected riverfront area extending from the mean annual high-water line corresponding to the elevation of bankfull discharge along each side of a perennial stream. Additionally, information on bankfull channel geometry and discharge are important to Federal, State, and local government agencies and private organizations involved in stream assessment and restoration projects. Regression equations are based on data from stream surveys at 33 sites (32 streamgages and 1 crest-stage gage operated by the U.S. Geological Survey) in and near Massachusetts. Drainage areas of the 33 sites ranged from 0.60 to 329 square miles (mi2). At 27 of the 33 sites, field data were collected and analyses were done to determine bankfull channel geometry and discharge as part of the present study. For 6 of the 33 sites, data on bankfull channel geometry and discharge were compiled from other studies done by the U.S. Geological Survey, Natural Resources Conservation Service of the U.S. Department of Agriculture, and the Vermont Department of Environmental Conservation. Similar techniques were used for field data collection and analysis for bankfull channel geometry and discharge at all 33 sites. Recurrence intervals of the bankfull discharge, which represent the frequency with which a stream fills its channel, averaged 1.53 years (median value 1.34 years) at the 33 sites. Simple regression equations were developed for bankfull width, mean depth, cross-sectional area, and discharge using drainage area, which is the most significant explanatory

Application of discontinuous Galerkin method for solving a compressible five-equation two-phase flow model

Science.gov (United States)

Saleem, M. Rehan; Ali, Ishtiaq; Qamar, Shamsul

2018-03-01

In this article, a reduced five-equation two-phase flow model is numerically investigated. The formulation of the model is based on the conservation and energy exchange laws. The model is non-conservative and the governing equations contain two equations for the mass conservation, one for the over all momentum and one for the total energy. The fifth equation is the energy equation for one of the two phases that includes a source term on the right hand side for incorporating energy exchange between the two fluids in the form of mechanical and thermodynamical works. A Runge-Kutta discontinuous Galerkin finite element method is applied to solve the model equations. The main attractive features of the proposed method include its formal higher order accuracy, its nonlinear stability, its ability to handle complicated geometries, and its ability to capture sharp discontinuities or strong gradients in the solutions without producing spurious oscillations. The proposed method is robust and well suited for large-scale time-dependent computational problems. Several case studies of two-phase flows are presented. For validation and comparison of the results, the same model equations are also solved by using a staggered central scheme. It was found that discontinuous Galerkin scheme produces better results as compared to the staggered central scheme.

Solitons, Lie Group Analysis and Conservation Laws of a (3+1)-Dimensional Modified Zakharov-Kuznetsov Equation in a Multicomponent Magnetised Plasma

Science.gov (United States)

Du, Xia-Xia; Tian, Bo; Chai, Jun; Sun, Yan; Yuan, Yu-Qiang

2017-11-01

In this paper, we investigate a (3+1)-dimensional modified Zakharov-Kuznetsov equation, which describes the nonlinear plasma-acoustic waves in a multicomponent magnetised plasma. With the aid of the Hirota method and symbolic computation, bilinear forms and one-, two- and three-soliton solutions are derived. The characteristics and interaction of the solitons are discussed graphically. We present the effects on the soliton's amplitude by the nonlinear coefficients which are related to the ratio of the positive-ion mass to negative-ion mass, number densities, initial densities of the lower- and higher-temperature electrons and ratio of the lower temperature to the higher temperature for electrons, as well as by the dispersion coefficient, which is related to the ratio of the positive-ion mass to the negative-ion mass and number densities. Moreover, using the Lie symmetry group theory, we derive the Lie point symmetry generators and the corresponding symmetry reductions, through which certain analytic solutions are obtained via the power series expansion method and the (G'/G) expansion method. We demonstrate that such an equation is strictly self-adjoint, and the conservation laws associated with the Lie point symmetry generators are derived.

Asymptotic behaviour of time averages for non-ergodic Gaussian processes

Science.gov (United States)

Ślęzak, Jakub

2017-08-01

In this work, we study the behaviour of time-averages for stationary (non-ageing), but ergodicity-breaking Gaussian processes using their representation in Fourier space. We provide explicit formulae for various time-averaged quantities, such as mean square displacement, density, and analyse the behaviour of time-averaged characteristic function, which gives insight into rich memory structure of the studied processes. Moreover, we show applications of the ergodic criteria in Fourier space, determining the ergodicity of the generalised Langevin equation's solutions.

On the XFEL Schrödinger Equation: Highly Oscillatory Magnetic Potentials and Time Averaging

KAUST Repository

Antonelli, Paolo; Athanassoulis, Agisillaos; Hajaiej, Hichem; Markowich, Peter A.

2014-01-01

We analyse a nonlinear Schrödinger equation for the time-evolution of the wave function of an electron beam, interacting selfconsistently through a Hartree-Fock nonlinearity and through the repulsive Coulomb interaction of an atomic nucleus

Lie-Mei symmetry and conserved quantities of the Rosenberg problem

International Nuclear Information System (INIS)

Liu Xiao-Wei; Li Yuan-Cheng

2011-01-01

The Rosenberg problem is a typical but not too complicated problem of nonholonomic mechanical systems. The Lie—Mei symmetry and the conserved quantities of the Rosenberg problem are studied. For the Rosenberg problem, the Lie and the Mei symmetries for the equation are obtained, the conserved quantities are deduced from them and then the definition and the criterion for the Lie—Mei symmetry of the Rosenberg problem are derived. Finally, the Hojman conserved quantity and the Mei conserved quantity are deduced from the Lie—Mei symmetry. (general)

Lie-Mei symmetry and conserved quantities of the Rosenberg problem

Science.gov (United States)

Liu, Xiao-Wei; Li, Yuan-Cheng

2011-07-01

The Rosenberg problem is a typical but not too complicated problem of nonholonomic mechanical systems. The Lie—Mei symmetry and the conserved quantities of the Rosenberg problem are studied. For the Rosenberg problem, the Lie and the Mei symmetries for the equation are obtained, the conserved quantities are deduced from them and then the definition and the criterion for the Lie—Mei symmetry of the Rosenberg problem are derived. Finally, the Hojman conserved quantity and the Mei conserved quantity are deduced from the Lie—Mei symmetry.

Nonlinearity, Conservation Law and Shocks

Indian Academy of Sciences (India)

However, genuine nonlinearity is always present in an ideal gas. The conservation form of the equation (25) brings in shocks which cut off the growing part of the amplitUde as shown in. Figure 15. Acknowledgements. The author sincerely thanks the two referees whose valuable comments led to an improvement of the ...

On dynamic equations for interaction of the affinor field with affine connection

International Nuclear Information System (INIS)

Pestov, A.B.

1987-01-01

The Lagrangian of interaction of an affinor field with an affine connection is constructed and the equations of motion and conservation laws are derived. It is shown that there exists a symmetric conserved tensor of the affine-connection energy-momentum

General relativistic continuum mechanics and the post-Newtonian equations of motion

International Nuclear Information System (INIS)

Morrill, T.H.

1991-01-01

Aspects are examined of general relativistic continuum mechanics. Perfectly elastic materials are dealt with but not exclusively. The derivation of their equations of motion is emphasized, in the post-Newtonian approximation. A reformulation is presented based on the tetrad formalism, of Carter and Quintana's theory of general relativistic elastic continua. A field Lagrangian is derived describing perfect material media; show that the usual covariant conservations law for perfectly elastic media is fully equivalent to the Euler-Lagrange equations describing these same media; and further show that the equations of motion for such materials follow directly from Einstein's field equations. In addition, a version of this principle shows that the local mass density in curved space-time partially depends on the amount and distribution of mass energy in the entire universe and is related to the mass density that would occur if space-time were flat. The total Lagrangian was also expanded in an EIH (Einstein, Infeld, Hoffmann) series to obtain a total post-Newtonian Lagrangian. The results agree with those found by solving Einstein's equations for the metric coefficients and by deriving the post-Newtonian equations of motion from the covariant conservation law

A simple mass-conserved level set method for simulation of multiphase flows

Science.gov (United States)

Yuan, H.-Z.; Shu, C.; Wang, Y.; Shu, S.

2018-04-01

In this paper, a modified level set method is proposed for simulation of multiphase flows with large density ratio and high Reynolds number. The present method simply introduces a source or sink term into the level set equation to compensate the mass loss or offset the mass increase. The source or sink term is derived analytically by applying the mass conservation principle with the level set equation and the continuity equation of flow field. Since only a source term is introduced, the application of the present method is as simple as the original level set method, but it can guarantee the overall mass conservation. To validate the present method, the vortex flow problem is first considered. The simulation results are compared with those from the original level set method, which demonstrates that the modified level set method has the capability of accurately capturing the interface and keeping the mass conservation. Then, the proposed method is further validated by simulating the Laplace law, the merging of two bubbles, a bubble rising with high density ratio, and Rayleigh-Taylor instability with high Reynolds number. Numerical results show that the mass is a well-conserved by the present method.

Maximum-principle-satisfying space-time conservation element and solution element scheme applied to compressible multifluids

KAUST Repository

Shen, Hua; Wen, Chih-Yung; Parsani, Matteo; Shu, Chi-Wang

2016-01-01

A maximum-principle-satisfying space-time conservation element and solution element (CE/SE) scheme is constructed to solve a reduced five-equation model coupled with the stiffened equation of state for compressible multifluids. We first derive a sufficient condition for CE/SE schemes to satisfy maximum-principle when solving a general conservation law. And then we introduce a slope limiter to ensure the sufficient condition which is applicative for both central and upwind CE/SE schemes. Finally, we implement the upwind maximum-principle-satisfying CE/SE scheme to solve the volume-fraction-based five-equation model for compressible multifluids. Several numerical examples are carried out to carefully examine the accuracy, efficiency, conservativeness and maximum-principle-satisfying property of the proposed approach.

Maximum-principle-satisfying space-time conservation element and solution element scheme applied to compressible multifluids

KAUST Repository

Shen, Hua

2016-10-19

A maximum-principle-satisfying space-time conservation element and solution element (CE/SE) scheme is constructed to solve a reduced five-equation model coupled with the stiffened equation of state for compressible multifluids. We first derive a sufficient condition for CE/SE schemes to satisfy maximum-principle when solving a general conservation law. And then we introduce a slope limiter to ensure the sufficient condition which is applicative for both central and upwind CE/SE schemes. Finally, we implement the upwind maximum-principle-satisfying CE/SE scheme to solve the volume-fraction-based five-equation model for compressible multifluids. Several numerical examples are carried out to carefully examine the accuracy, efficiency, conservativeness and maximum-principle-satisfying property of the proposed approach.

On the Definition of Energy for a Continuum, Its Conservation Laws, and the Energy-Momentum Tensor

Directory of Open Access Journals (Sweden)

Mayeul Arminjon

2016-01-01

Full Text Available We review the energy concept in the case of a continuum or a system of fields. First, we analyze the emergence of a true local conservation equation for the energy of a continuous medium, taking the example of an isentropic continuum in Newtonian gravity. Next, we consider a continuum or a system of fields in special relativity: we recall that the conservation of the energy-momentum tensor contains two local conservation equations of the same kind as before. We show that both of these equations depend on the reference frame and that, however, they can be given a rigorous meaning. Then, we review the definitions of the canonical and Hilbert energy-momentum tensors from a Lagrangian through the principle of stationary action in general space-time. Using relatively elementary mathematics, we prove precise results regarding the definition of the Hilbert tensor field, its uniqueness, and its tensoriality. We recall the meaning of its covariant conservation equation. We end with a proof of uniqueness of the energy density and flux, when both depend polynomially on the fields.

Deriving the equations of motion of porous isotropic media

International Nuclear Information System (INIS)

Pride, S.R.; Gangi, A.F.; Morgan, F.D.

1992-01-01

The equations of motion and stress/strain relations for the linear dynamics of a two-phase, fluid/solid, isotropic, porous material have been derived by a direct volume averaging of the equations of motion and stress-strain relations known to apply in each phase. The equations thus obtained are shown to be consistent with Biot's equations of motion and stress/strain relations; however, the effective fluid density in the equation of relative flow has an unambiguous definition in terms of the tractions acting on the pore walls. The stress/strain relations of the theory correspond to 'quasistatic' stressing (i.e., inertial effects are ignored). It is demonstrated that using such quasistatic stress/strain relations in the equations of motion is justified whenever the wavelengths are greater than a length characteristic of the averaging volume size. 37 refs., 2 figs

[Perinatal result with conservative treatment in preeclampsia-eclampsia].

Science.gov (United States)

Briones-Garduño, Jesús Carlos; de León-Ponce, Manuel Díaz; González-Vargas, Angel; Briones-Vega, Carlos Gabriel

2003-01-01

Conservative treatment in severe preeclampsia has been documented by several authors citing significant improvement in neonatal outcome lacking a significant increase in maternal complications. Our objective was to inform of our preliminary results using protocolized conservative management in women with preeclampsia-eclampsia, favoring better neonate conditions. We included 34 patients with average age of 28.2 years with documented severe preeclampsia-eclampsia complicating a 36-weeks or less pregnancy, admitted in the obstetric intensive care unit (OICU) between October 2001 and February 2002. Patients received protocolized management consisting of intravascular colume expansion, anti-hypertensive control, target organ protection, monitoring, and clinical observation. We considered conservative management as a 24 or more period offered to patients with satisfactory response to medical treatment and no evidence of binomial compromise. Of our group, 85% corresponded to severe preeclampsia, 9% to eclampsia, 3% to imminence of eclampsia, and 3% to HELLP syndrome. Average stay in OICU was 5.5 days with 3.5 days average management before pregnancy was interrupted. These patients presented mean gestational age of 32.8 weeks during which we observed anemia, low platelets, D dimmer increments, MAP average of 112.8, PCOc 18.6, and BI 0.15. We obtained 36 live newborns of whom 12% four died, two were extremely immatures (510 g and 600 g, respectively); one 980-g newborn presented intraventricular hemorrhage, and a 1,450-g newborn had multiple organ failure. Conservative treatment in patients with severe preeclampsia-eclampsia is a feasible alternative in hospitals with an ICU. Conservative management can improve neonatal survival and prognosis in preterm newborns.

Conservation caring: measuring the influence of zoo visitors' connection to wildlife on pro-conservation behaviors.

Science.gov (United States)

Skibins, Jeffrey C; Powell, Robert B

2013-01-01

Zoos in the 21st century are striving to make effective contributions to conservation. Although zoos are extremely popular and host over 600 million visitors worldwide, one challenge zoos face is how to effectively engage visitors and raise awareness and action for conservation. To this end, zoos commonly rely on charismatic megafauna, which have been shown to elicit a connection with zoo visitors. However, little is known about how to measure a connection to a species or how this connection may influence conservation behaviors. This study had two sequential objectives. The first was to develop a scale to measure visitors' connection to a species (Conservation Caring). The second was to investigate the relationship of Conservation Caring to pro-conservation behaviors, following a zoo experience. Pre- (n = 411) and post-visit (n = 452) responses were collected from three sites in order to assess the reliability and validity of a scale to measure Conservation Caring. Structural equation modeling was used to explore the relationship between Conservation Caring and pro-conservation behaviors. Conservation Caring was deemed a valid and reliable scale and was a strong predictor of species oriented behaviors (β = 0.62), for example, "adopting" an animal, but a weak predictor for biodiversity oriented behaviors (β = 0.07), for example, supporting sustainability policies. Results support the role zoos can play in fostering a connection to wildlife and stimulating pro-conservation behaviors. Additionally, visitors connected to a wide array of animals. On the basis of these results, zoos may recruit a wider assemblage of species as potential flagships. © 2013 Wiley Periodicals, Inc.

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

International Nuclear Information System (INIS)

Chau, L.L.

1987-01-01

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

Comparison of Numerical Approaches to a Steady-State Landscape Equation

Science.gov (United States)

Bachman, S.; Peckham, S.

2008-12-01

A mathematical model of an idealized fluvial landscape has been developed, in which a land surface will evolve to preserve dendritic channel networks as the surface is lowered. The physical basis for this model stems from the equations for conservation of mass for water and sediment. These equations relate the divergence of the 2D vector fields showing the unit-width discharge of water and sediment to the excess rainrate and tectonic uplift on the land surface. The 2D flow direction is taken to be opposite to the water- surface gradient vector. These notions are combined with a generalized Manning-type flow resistance formula and a generalized sediment transport law to give a closed mathematical system that can, in principle, be solved for all variables of interest: discharge of water and sediment, land surface height, vertically- averaged flow velocity, water depth, and shear stress. The hydraulic geometry equations (Leopold et. al, 1964, 1995) are used to incorporate width, depth, velocity, and slope of river channels as powers of the mean-annual river discharge. Combined, they give the unit- width discharge of the stream as a power, γ, of the water surface slope. The simplified steady-state model takes into account three components among those listed above: conservation of mass for water, flow opposite the gradient, and a slope-discharge exponent γ = -1 to reflect mature drainage networks. The mathematical representation of this model appears as a second-order hyperbolic partial differential equation (PDE) where the diffusivity is inversely proportional to the square of the local surface slope. The highly nonlinear nature of this PDE has made it very difficult to solve both analytically and numerically. We present simplistic analytic solutions to this equation which are used to test the validity of the numerical algorithms. We also present three such numerical approaches which have been used in solving the differential equation. The first is based on a

Averaging Bias Correction for Future IPDA Lidar Mission MERLIN

Science.gov (United States)

Tellier, Yoann; Pierangelo, Clémence; Wirth, Martin; Gibert, Fabien

2018-04-01

The CNES/DLR MERLIN satellite mission aims at measuring methane dry-air mixing ratio column (XCH4) and thus improving surface flux estimates. In order to get a 1% precision on XCH4 measurements, MERLIN signal processing assumes an averaging of data over 50 km. The induced biases due to the non-linear IPDA lidar equation are not compliant with accuracy requirements. This paper analyzes averaging biases issues and suggests correction algorithms tested on realistic simulated scenes.

The finite volume element (FVE) and multigrid method for the incompressible Navier-Stokes equations

International Nuclear Information System (INIS)

Gu Lizhen; Bao Weizhu

1992-01-01

The authors apply FVE method to discrete INS equations with the original variable, in which the bilinear square finite element and the square finite volume are chosen. The discrete schemes of INS equations are presented. The FMV multigrid algorithm is applied to solve that discrete system, where DGS iteration is used as smoother, DGS distributive mode for the INS discrete system is also presented. The sample problems for the square cavity flow with Reynolds number Re≤100 are successfully calculated. The numerical solutions show that the results with 1 FMV is satisfactory and when Re is not large, The FVE discrete scheme of the conservative INS equations and that of non-conservative INS equations with linearization both can provide almost same accuracy

Front tracking for hyperbolic conservation laws

CERN Document Server

Holden, Helge

2002-01-01

Hyperbolic conservation laws are central in the theory of nonlinear partial differential equations and in science and technology. The reader is given a self-contained presentation using front tracking, which is also a numerical method. The multidimensional scalar case and the case of systems on the line are treated in detail. A chapter on finite differences is included. "It is already one of the few best digests on this topic. The present book is an excellent compromise between theory and practice. Students will appreciate the lively and accurate style." D. Serre, MathSciNet "I have read the book with great pleasure, and I can recommend it to experts as well as students. It can also be used for reliable and very exciting basis for a one-semester graduate course." S. Noelle, Book review, German Math. Soc. "Making it an ideal first book for the theory of nonlinear partial differential equations...an excellent reference for a graduate course on nonlinear conservation laws." M. Laforest, Comp. Phys. Comm.

Electromagnetic-field equations in the six-dimensional space-time R6

International Nuclear Information System (INIS)

Teli, M.T.; Palaskar, D.

1984-01-01

Maxwell's equations (without monopoles) for electromagnetic fields are obtained in six-dimensional space-time. The equations possess structural symmetry in space and time, field and source densities. Space-time-symmetric conservation laws and field solutions are obtained. The results are successfully correlated with their four-dimensional space-time counterparts

Transfer equations for spectral densities of inhomogeneous MHD turbulence

International Nuclear Information System (INIS)

Tu, C.-Y.; Marsch, E.

1990-01-01

On the basis of the dynamic equations governing the evolution of magnetohydrodynamic fluctuations expressed in terms of Elsaesser variables and of their correlation functions derived by Marsch and Tu, a new set of equations is presented describing the evolutions of the energy spectrum e ± and of the residual energy spectra e R and e S of MHD turbulence in an inhomogeneous magnetofluid. The nonlinearities associated with triple correlations in these equations are analysed in detail and evaluated approximately. The resulting energy-transfer functions across wavenumber space are discussed. For e ± they are shown to be approximately energy-conserving if the gradients of the flow speed and density are weak. New cascading functions are heuristically determined by an appropriate dimensional analysis and plausible physical arguments, following the standard phenomenology of fluid turbulence. However, for e R the triple correlations do not correspond to an 'energy' conserving process, but also represent a nonlinear source term for e R . If this source term can be neglected, the spectrum equations are found to be closed. The problem of dealing with the nonlinear source terms remains to be solved in future investigations. (author)

Nonlocal symmetry generators and explicit solutions of some partial differential equations

International Nuclear Information System (INIS)

Qin Maochang

2007-01-01

The nonlocal symmetry of a partial differential equation is studied in this paper. The partial differential equation written as a conservation law can be transformed into an equivalent system by introducing a suitable potential. The nonlocal symmetry group generators of original partial differential equations can be obtained through their equivalent system. Further, new explicit solutions can be constructed from the newly obtained symmetry generators. The Burgers equation is chosen as an example; many new valuable explicit solutions and nonlocal symmetry generators are presented

Relative Stabilities of Conserved and Non-Conserved Structures in the OB-Fold Superfamily

Directory of Open Access Journals (Sweden)

Andrei T. Alexandrescu

2009-05-01

Full Text Available The OB-fold is a diverse structure superfamily based on a β-barrel motif that is often supplemented with additional non-conserved secondary structures. Previous deletion mutagenesis and NMR hydrogen exchange studies of three OB-fold proteins showed that the structural stabilities of sites within the conserved β-barrels were larger than sites in non-conserved segments. In this work we examined a database of 80 representative domain structures currently classified as OB-folds, to establish the basis of this effect. Residue-specific values were obtained for the number of Cα-Cα distance contacts, sequence hydrophobicities, crystallographic B-factors, and theoretical B-factors calculated from a Gaussian Network Model. All four parameters point to a larger average flexibility for the non-conserved structures compared to the conserved β-barrels. The theoretical B-factors and contact densities show the highest sensitivity.Our results suggest a model of protein structure evolution in which novel structural features develop at the periphery of conserved motifs. Core residues are more resistant to structural changes during evolution since their substitution would disrupt a larger number of interactions. Similar factors are likely to account for the differences in stability to unfolding between conserved and non-conserved structures.

Some Properties of the M3D-C1 Form of the 3D Magnetohydrodynamics Equations

International Nuclear Information System (INIS)

Breslau, J.; Ferraro, N.; Jardin, S.

2009-01-01

We introduce a set of scalar variables and projection operators for the vector momentum and magnetic field evolution equations that have several unique and desirable properties, making them a preferred system for solving the magnetohydrodynamics equations in a torus with a strong toroidal magnetic field. We derive a 'weak form' of these equations that explicitly conserves energy and is suitable for a Galerkin finite element formulation provided the basis elements have C1 continuity. Systems of reduced equations are discussed, along with their energy conservation properties. An implicit time advance is presented that adds diagonally dominant self-adjoint energy terms to the mass matrix to obtain numerical stability.

Rarefaction and shock waves for multi-dimensional hyperbolic conservation laws

International Nuclear Information System (INIS)

Dening, Li

1991-01-01

In this paper, the author wants to show the local existence of a solution of combination of shock and rarefaction waves for the multi-dimensional hyperbolic system of conservation laws. The typical example he has in mind is the Euler equations for compressible fluid. More generally, he studies the hyperbolic system of conservation laws ∂ t F 0 (u) + Σ j=1 n ∂ x j F j (u)=0 where u=(u 1 ....,u m ) and F j (u), j=0,...,n are m-dimensional vector-valued functions. He'll impose some conditions in the following on the systems (1.2). All these conditions are satisfied by the Euler equations

Lie symmetry and their conserved quantities of Tzénoff equations for the vairable mass nonholonomic systems%变质量非完整系统Tzénoff方程的Lie对称性与其导出的守恒量

Institute of Scientific and Technical Information of China (English)

郑世旺; 王建波; 陈向炜; 李彦敏; 解加芳

2012-01-01

航天器运行系统大都属于变质量力学系统,变质量力学系统的对称性和守恒量隐含着航天系统更深刻的物理规律.本文首先导出了变质量非完整力学系统的Tzénoff方程,然后研究了变质量非完整力学系统Tzénoff方程的Lie对称性及其所导出的守恒量,给出了这种守恒量的函数表达式和导出这种守恒量的判据方程.该研究结果对进一步探究变质量系统所遵循的守恒规律具有一定的理论价值.%The operational system of the spacecraft is general a variable mass one,of which the symmetry and the conserved quantity imply physical rules of the space system.In this paper,Tzénoff equations of the variable mass nonholonomic system are derived,from which the Lie symmetries of Tzénoff equations for the variable mass nonholonomic system and conserved quantities are derived and are researched.The function expressions of conserved quantities and the criterion equations which deduce these conserved quantities are presented.This result has some theoretical value for further research of the conservation laws obeyed by the variable mass system.

Discrete Exterior Calculus Discretization of Incompressible Navier-Stokes Equations

KAUST Repository

Mohamed, Mamdouh S.; Hirani, Anil N.; Samtaney, Ravi

2017-01-01

A conservative discretization of incompressible Navier-Stokes equations over surface simplicial meshes is developed using discrete exterior calculus (DEC). Numerical experiments for flows over surfaces reveal a second order accuracy

Charge-conserving FEM-PIC schemes on general grids

International Nuclear Information System (INIS)

Campos Pinto, M.; Jund, S.; Salmon, S.; Sonnendruecker, E.

2014-01-01

Particle-In-Cell (PIC) solvers are a major tool for the understanding of the complex behavior of a plasma or a particle beam in many situations. An important issue for electromagnetic PIC solvers, where the fields are computed using Maxwell's equations, is the problem of discrete charge conservation. In this article, we aim at proposing a general mathematical formulation for charge-conserving finite-element Maxwell solvers coupled with particle schemes. In particular, we identify the finite-element continuity equations that must be satisfied by the discrete current sources for several classes of time-domain Vlasov-Maxwell simulations to preserve the Gauss law at each time step, and propose a generic algorithm for computing such consistent sources. Since our results cover a wide range of schemes (namely curl-conforming finite element methods of arbitrary degree, general meshes in two or three dimensions, several classes of time discretization schemes, particles with arbitrary shape factors and piecewise polynomial trajectories of arbitrary degree), we believe that they provide a useful roadmap in the design of high-order charge-conserving FEM-PIC numerical schemes. (authors)

Propagation of multidimensional nonlinear waves and kinematical conservation laws

CERN Document Server

Prasad, Phoolan

2017-01-01

This book formulates the kinematical conservation laws (KCL), analyses them and presents their applications to various problems in physics. Finally, it addresses one of the most challenging problems in fluid dynamics: finding successive positions of a curved shock front. The topics discussed are the outcome of collaborative work that was carried out mainly at the Indian Institute of Science, Bengaluru, India. The theory presented in the book is supported by referring to extensive numerical results. The book is organised into ten chapters. Chapters 1–4 offer a summary of and briefly discuss the theory of hyperbolic partial differential equations and conservation laws. Formulation of equations of a weakly nonlinear wavefront and those of a shock front are briefly explained in Chapter 5, while Chapter 6 addresses KCL theory in space of arbitrary dimensions. The remaining chapters examine various analyses and applications of KCL equations ending in the ultimate goal-propagation of a three-dimensional curved sho...

The NO-pair equation---Fundamental problems, numerical solutions and applications

International Nuclear Information System (INIS)

Martensson-Pendrill, A.

1989-01-01

The solution of inhomogeneous two-particle differential equations is a powerful method for treating correlation effects in the non-relativistic case and is reviewed briefly. The relativistic generalization, the so-called Dirac-Coulomb equation, is complicated by the need to distinguish between positive and negative energy states. The construction and use of projection operators onto positive energy states are presented and the application of the pair equation to other properties such as parity non-conservation is discussed

Stable solutions of nonlocal electron heat transport equations

International Nuclear Information System (INIS)

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

1991-01-01

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

A generalised groundwater flow equation using the concept of non ...

African Journals Online (AJOL)

The classical Darcy law is generalised by regarding the water flow as a function of a non-integer order derivative of the piezometric head. This generalised law and the law of conservation of mass are then used to derive a new equation for groundwater flow. Numerical solutions of this equation for various fractional orders of ...

A stable penalty method for the compressible Navier-Stokes equations: II: One-dimensional domain decomposition schemes

DEFF Research Database (Denmark)

Hesthaven, Jan

1997-01-01

This paper presents asymptotically stable schemes for patching of nonoverlapping subdomains when approximating the compressible Navier-Stokes equations given on conservation form. The scheme is a natural extension of a previously proposed scheme for enforcing open boundary conditions and as a res......This paper presents asymptotically stable schemes for patching of nonoverlapping subdomains when approximating the compressible Navier-Stokes equations given on conservation form. The scheme is a natural extension of a previously proposed scheme for enforcing open boundary conditions...... and as a result the patching of subdomains is local in space. The scheme is studied in detail for Burgers's equation and developed for the compressible Navier-Stokes equations in general curvilinear coordinates. The versatility of the proposed scheme for the compressible Navier-Stokes equations is illustrated...

HYDRA-II: A hydrothermal analysis computer code: Volume 1, Equations and numerics

International Nuclear Information System (INIS)

McCann, R.A.

1987-04-01

HYDRA-II is a hydrothermal computer code capable of three-dimensional analysis of coupled conduction, convection, and thermal radiation problems. This code is especially appropriate for simulating the steady-state performance of spent fuel storage systems. The code has been evaluated for this application for the US Department of Energy's Commercial Spent Fuel Management Program. HYDRA-II provides a finite difference solution in Cartesian coordinates to the equations governing the conservation of mass, momentum, and energy. A cylindrical coordinate system may also be used to enclose the Cartesian coordinate system. This exterior coordinate system is useful for modeling cylindrical cask bodies. The difference equations for conservation of momentum are enhanced by the incorporation of directional porosities and permeabilities that aid in modeling solid structures whose dimensions may be smaller than the computational mesh. The equation for conservation of energy permits of modeling of orthotropic physical properties and film resistances. Several automated procedures are available to model radiation transfer within enclosures and from fuel rod to fuel rod. The documentation of HYDRA-II is presented in three separate volumes. This volume, Volume I - Equations and Numerics, describes the basic differential equations, illustrates how the difference equations are formulated, and gives the solution procedures employed. Volume II - User's Manual contains code flow charts, discusses the code structure, provides detailed instructions for preparing an input file, and illustrates the operation of the code by means of a model problem. The final volume, Volume III - Verification/Validation Assessments, presents results of numerical simulations of single- and multiassembly storage systems and comparisons with experimental data. 4 refs

Conservation of energy and momentum in nonrelativistic plasmas

International Nuclear Information System (INIS)

Sugama, H.; Watanabe, T.-H.; Nunami, M.

2013-01-01

Conservation laws of energy and momentum for nonrelativistic plasmas are derived from applying Noether's theorem to the action integral for the Vlasov-Poisson-Ampère system [Sugama, Phys. Plasmas 7, 466 (2000)]. The symmetric pressure tensor is obtained from modifying the asymmetric canonical pressure tensor with using the rotational symmetry of the action integral. Differences between the resultant conservation laws and those for the Vlasov-Maxwell system including the Maxwell displacement current are clarified. These results provide a useful basis for gyrokinetic conservation laws because gyrokinetic equations are derived as an approximation of the Vlasov-Poisson-Ampère system.

Stochastic TDHF and the Boltzman-Langevin equation

International Nuclear Information System (INIS)

Suraud, E.; Reinhard, P.G.

1991-01-01

Outgoing from a time-dependent theory of correlations, we present a stochastic differential equation for the propagation of ensembles of Slater determinants, called Stochastic Time-Dependent Hartree-Fock (Stochastic TDHF). These ensembles are allowed to develop large fluctuations in the Hartree-Fock mean fields. An alternative stochastic differential equation, the Boltzmann-Langevin equation, can be derived from Stochastic TDHF by averaging over subensembles with small fluctuations

Report on energy conservation. Potential feasibility study on energy conservation at the national textile factories (Energy-saving equipment)

Energy Technology Data Exchange (ETDEWEB)

NONE

2001-03-01

In consideration of the future CDM project, feasibility study was made on the energy conservation/cost reduction in national textile factories in Egypt. In textile factories in Egypt, the steam by oil fueled boiler is used in dyeing, drying and other production processes. In this study, studied were the energy conservation effect, greenhouse effect gas reduction effect and economical efficiency in the case of using electric power/steam by the natural gas fueled cogeneration system. As a result of the study, the energy conservation effect in 6 national textile factories became 13% on average, which equals to the conservation of 15,000 toe per year. The greenhouse effect gas reduction effect also became 13%, which equals to the reduction of approximately 46,000 t-CO2 per year. As to the economical efficiency, the depreciation period of investment became 4.0-8.6 years (5.6 years on average). In the light of the operation period of cogeneration facilities of 15 years, it was indicated that the system was very advantageous. (NEDO)

A low-dispersion, exactly energy-charge-conserving semi-implicit relativistic particle-in-cell algorithm

Science.gov (United States)

Chen, Guangye; Luis, Chacon; Bird, Robert; Stark, David; Yin, Lin; Albright, Brian

2017-10-01

Leap-frog based explicit algorithms, either ``energy-conserving'' or ``momentum-conserving'', do not conserve energy discretely. Time-centered fully implicit algorithms can conserve discrete energy exactly, but introduce large dispersion errors in the light-wave modes, regardless of timestep sizes. This can lead to intolerable simulation errors where highly accurate light propagation is needed (e.g. laser-plasma interactions, LPI). In this study, we selectively combine the leap-frog and Crank-Nicolson methods to produce a low-dispersion, exactly energy-and-charge-conserving PIC algorithm. Specifically, we employ the leap-frog method for Maxwell equations, and the Crank-Nicolson method for particle equations. Such an algorithm admits exact global energy conservation, exact local charge conservation, and preserves the dispersion properties of the leap-frog method for the light wave. The algorithm has been implemented in a code named iVPIC, based on the VPIC code developed at LANL. We will present numerical results that demonstrate the properties of the scheme with sample test problems (e.g. Weibel instability run for 107 timesteps, and LPI applications.

Energy preserving integration of bi-Hamiltonian partial differential equations

NARCIS (Netherlands)

Karasozen, B.; Simsek, G.

2013-01-01

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

Approximations for transport parameters and self-averaging properties for point-like injections in heterogeneous media

International Nuclear Information System (INIS)

Eberhard, Jens

2004-01-01

We focus on transport parameters in heterogeneous media with a flow modelled by an ensemble of periodic and Gaussian random fields. The parameters are determined by ensemble averages. We study to what extent these averages represent the behaviour in a single realization. We calculate the centre-of-mass velocity and the dispersion coefficient using approximations based on a perturbative expansion for the transport equation, and on the iterative solution of the Langevin equation. Compared with simulations, the perturbation theory reproduces the numerical results only poorly, whereas the iterative solution yields good results. Using these approximations, we investigate the self-averaging properties. The ensemble average of the velocity characterizes the behaviour of a realization for large times in both ensembles. The dispersion coefficient is not self-averaging in the ensemble of periodic fields. For the Gaussian ensemble the asymptotic dispersion coefficient is self-averaging. For finite times, however, the fluctuations are so large that the average does not represent the behaviour in a single realization

Averaging Bias Correction for Future IPDA Lidar Mission MERLIN

Directory of Open Access Journals (Sweden)

Tellier Yoann

2018-01-01

Full Text Available The CNES/DLR MERLIN satellite mission aims at measuring methane dry-air mixing ratio column (XCH4 and thus improving surface flux estimates. In order to get a 1% precision on XCH4 measurements, MERLIN signal processing assumes an averaging of data over 50 km. The induced biases due to the non-linear IPDA lidar equation are not compliant with accuracy requirements. This paper analyzes averaging biases issues and suggests correction algorithms tested on realistic simulated scenes.

Evaluation of abutment scour prediction equations with field data

Science.gov (United States)

Benedict, S.T.; Deshpande, N.; Aziz, N.M.

2007-01-01

The U.S. Geological Survey, in cooperation with FHWA, compared predicted abutment scour depths, computed with selected predictive equations, with field observations collected at 144 bridges in South Carolina and at eight bridges from the National Bridge Scour Database. Predictive equations published in the 4th edition of Evaluating Scour at Bridges (Hydraulic Engineering Circular 18) were used in this comparison, including the original Froehlich, the modified Froehlich, the Sturm, the Maryland, and the HIRE equations. The comparisons showed that most equations tended to provide conservative estimates of scour that at times were excessive (as large as 158 ft). Equations also produced underpredictions of scour, but with less frequency. Although the equations provide an important resource for evaluating abutment scour at bridges, the results of this investigation show the importance of using engineering judgment in conjunction with these equations.

A RED modified weighted moving average for soft real-time application

Directory of Open Access Journals (Sweden)

Domanśka Joanna

2014-09-01

Full Text Available The popularity of TCP/IP has resulted in an increase in usage of best-effort networks for real-time communication. Much effort has been spent to ensure quality of service for soft real-time traffic over IP networks. The Internet Engineering Task Force has proposed some architecture components, such as Active Queue Management (AQM. The paper investigates the influence of the weighted moving average on packet waiting time reduction for an AQM mechanism: the RED algorithm. The proposed method for computing the average queue length is based on a difference equation (a recursive equation. Depending on a particular optimality criterion, proper parameters of the modified weighted moving average function can be chosen. This change will allow reducing the number of violations of timing constraints and better use of this mechanism for soft real-time transmissions. The optimization problem is solved through simulations performed in OMNeT++ and later verified experimentally on a Linux implementation

BHR equations re-derived with immiscible particle effects

Energy Technology Data Exchange (ETDEWEB)

Schwarzkopf, John Dennis [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Horwitz, Jeremy A. [Stanford Univ., CA (United States)

2015-05-01

Compressible and variable density turbulent flows with dispersed phase effects are found in many applications ranging from combustion to cloud formation. These types of flows are among the most challenging to simulate. While the exact equations governing a system of particles and fluid are known, computational resources limit the scale and detail that can be simulated in this type of problem. Therefore, a common method is to simulate averaged versions of the flow equations, which still capture salient physics and is relatively less computationally expensive. Besnard developed such a model for variable density miscible turbulence, where ensemble-averaging was applied to the flow equations to yield a set of filtered equations. Besnard further derived transport equations for the Reynolds stresses, the turbulent mass flux, and the density-specific volume covariance, to help close the filtered momentum and continuity equations. We re-derive the exact BHR closure equations which include integral terms owing to immiscible effects. Physical interpretations of the additional terms are proposed along with simple models. The goal of this work is to extend the BHR model to allow for the simulation of turbulent flows where an immiscible dispersed phase is non-trivially coupled with the carrier phase.

Group-theoretical model of developed turbulence and renormalization of the Navier-Stokes equation.

Science.gov (United States)

Saveliev, V L; Gorokhovski, M A

2005-07-01

On the basis of the Euler equation and its symmetry properties, this paper proposes a model of stationary homogeneous developed turbulence. A regularized averaging formula for the product of two fields is obtained. An equation for the averaged turbulent velocity field is derived from the Navier-Stokes equation by renormalization-group transformation.

Lax pairs and conservation laws for two differential-difference systems

International Nuclear Information System (INIS)

Li Chunxia

2003-01-01

A coupled extended Lotka-Volterra lattice and a special Toda lattice are derived from the existing bilinear equations. Starting from the corresponding bilinear Baecklund transformation, Lax pairs for these two differential-difference systems are obtained. Furthermore, an infinite number of conservation laws for the differential-difference equations are deduced from the Lax pairs in a systematic way

Predictive Trailing-Edge Modulation Average Current Control in DC-DC Converters

Directory of Open Access Journals (Sweden)

LASCU, D.

2013-11-01

Full Text Available The paper investigates predictive digital average current control (PDACC in dc/dc converters using trailing-edge modulation (TEM. The study is focused on the recurrence duty cycle equation and then stability analysis is performed. It is demonstrated that average current control using trailing-edge modulation is stable on the whole range of the duty cycle and thus design problems are highly reduced. The analysis is carried out in a general manner, independent of converter topology and therefore the results can then be easily applied for a certain converter (buck, boost, buck-boost, etc.. The theoretical considerations are confirmed for a boost converter first using the MATLAB program based on state-space equations and finally with the CASPOC circuit simulation package.

On the evolution equations, solvable through the inverse scattering method

International Nuclear Information System (INIS)

Gerdjikov, V.S.; Khristov, E.Kh.

1979-01-01

The nonlinear evolution equations (NLEE), related to the one-parameter family of Dirac operators are considered in a uniform manner. The class of NLEE solvable through the inverse scatterina method and their conservation laws are described. The description of the hierarchy of Hamiltonian structures and the proof of complete integrability of the NLEE is presented. The class of Baecklund transformations for these NLEE is derived. The general formulae are illustrated by two important examples: the nonlinear Schroedinger equation and the sine-Gordon equation

Fokker-Planck equation for the non-Markovian Brownian motion in the presence of a magnetic field

Science.gov (United States)

Das, Joydip; Mondal, Shrabani; Bag, Bidhan Chandra

2017-10-01

In the present study, we have proposed the Fokker-Planck equation in a simple way for a Langevin equation of motion having ordinary derivative (OD), the Gaussian random force and a generalized frictional memory kernel. The equation may be associated with or without conservative force field from harmonic potential. We extend this method for a charged Brownian particle in the presence of a magnetic field. Thus, the present method is applicable for a Langevin equation of motion with OD, the Gaussian colored thermal noise and any kind of linear force field that may be conservative or not. It is also simple to apply this method for the colored Gaussian noise that is not related to the damping strength.

Quantitative Compactness Estimates for Hamilton-Jacobi Equations

Science.gov (United States)

Ancona, Fabio; Cannarsa, Piermarco; Nguyen, Khai T.

2016-02-01

We study quantitative compactness estimates in {W^{1,1}_{loc}} for the map {S_t}, {t > 0} that is associated with the given initial data {u_0in Lip (R^N)} for the corresponding solution {S_t u_0} of a Hamilton-Jacobi equation u_t+Hbig(nabla_{x} ubig)=0, qquad t≥ 0,quad xinR^N, with a uniformly convex Hamiltonian {H=H(p)}. We provide upper and lower estimates of order {1/\\varepsilon^N} on the Kolmogorov {\\varepsilon}-entropy in {W^{1,1}} of the image through the map S t of sets of bounded, compactly supported initial data. Estimates of this type are inspired by a question posed by Lax (Course on Hyperbolic Systems of Conservation Laws. XXVII Scuola Estiva di Fisica Matematica, Ravello, 2002) within the context of conservation laws, and could provide a measure of the order of "resolution" of a numerical method implemented for this equation.

Structural Equation Modeling of Multivariate Time Series

Science.gov (United States)

du Toit, Stephen H. C.; Browne, Michael W.

2007-01-01

The covariance structure of a vector autoregressive process with moving average residuals (VARMA) is derived. It differs from other available expressions for the covariance function of a stationary VARMA process and is compatible with current structural equation methodology. Structural equation modeling programs, such as LISREL, may therefore be…

Covariant field equations, gauge fields and conservation laws from Yang-Mills matrix models

International Nuclear Information System (INIS)

Steinacker, Harold

2009-01-01

The effective geometry and the gravitational coupling of nonabelian gauge and scalar fields on generic NC branes in Yang-Mills matrix models is determined. Covariant field equations are derived from the basic matrix equations of motions, known as Yang-Mills algebra. Remarkably, the equations of motion for the Poisson structure and for the nonabelian gauge fields follow from a matrix Noether theorem, and are therefore protected from quantum corrections. This provides a transparent derivation and generalization of the effective action governing the SU(n) gauge fields obtained in [1], including the would-be topological term. In particular, the IKKT matrix model is capable of describing 4-dimensional NC space-times with a general effective metric. Metric deformations of flat Moyal-Weyl space are briefly discussed.

Discrete Exterior Calculus Discretization of Incompressible Navier-Stokes Equations

KAUST Repository

Mohamed, Mamdouh S.

2017-05-23

A conservative discretization of incompressible Navier-Stokes equations over surface simplicial meshes is developed using discrete exterior calculus (DEC). Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured-triangular meshes, and first order accuracy otherwise. The mimetic character of many of the DEC operators provides exact conservation of both mass and vorticity, in addition to superior kinetic energy conservation. The employment of barycentric Hodge star allows the discretization to admit arbitrary simplicial meshes. The discretization scheme is presented along with various numerical test cases demonstrating its main characteristics.

BRST, generalized Maurer-Cartan equations and CFT

Energy Technology Data Exchange (ETDEWEB)

Zeitlin, Anton M. [Department of Mathematics, Yale University, 442 Dunham Lab, 10 Hillhouse Ave., New Haven, CT 06511 (United States); St. Petersburg Department of Steklov Mathematical Institute, Fontanka, 27, St. Petersburg 191023 (Russian Federation)]. E-mail: zam@math.ipme.ru

2006-12-25

The paper is devoted to the study of BRST charge in perturbed two-dimensional conformal field theory. The main goal is to write the operator equation expressing the conservation law of BRST charge in perturbed theory in terms of purely algebraic operations on the corresponding operator algebra, which are defined via the OPE. The corresponding equations are constructed and their symmetries are studied up to the second order in formal coupling constant. It appears that the obtained equations can be interpreted as generalized Maurer-Cartan ones. We study two concrete examples in detail: the bosonic nonlinear sigma model and perturbed first order theory. In particular, we show that the Einstein equations, which are the conformal invariance conditions for both these perturbed theories, expanded up to the second order, can be rewritten in such generalized Maurer-Cartan form.

On conservation laws for models in discrete, noncommutative and fractional differential calculus

International Nuclear Information System (INIS)

Klimek, M.

2001-01-01

We present the general method of derivation the explicit form of conserved currents for equations built within the framework of discrete, noncommutative or fractional differential calculus. The procedure applies to linear models with variable coefficients including also nonlinear potential part. As an example an equation on quantum plane, nonlinear Toda lattice model and homogeneous equation of fractional diffusion in 1+1 dimensions are studied

Conservative integration of rigid body motion by quaternion parameters with implicit constraints

DEFF Research Database (Denmark)

Nielsen, Martin Bjerre; Krenk, Steen

2012-01-01

An angular momentum and energy‐conserving time integration algorithm for rigid body rotation is formulated in terms of the quaternion parameters and the corresponding four‐component conjugate momentum vector via Hamilton's equations. The introduction of an extended mass matrix leads to a symmetric...... these equations via the set of momentum equations. Initially, the normalization of the quaternion array is introduced via a Lagrange multiplier. However, this Lagrange multiplier can be expressed explicitly in terms of the gradient of the external load potential, and elimination of the Lagrange multiplier from...... the final format leaves only an explicit projection applied to the external load potential gradient. An algorithm is developed by forming a finite increment of the Hamiltonian. This procedure identifies the proper selection of increments and mean values, and leads to an algorithm with conservation...

Random walk and the heat equation

CERN Document Server

Lawler, Gregory F

2010-01-01

The heat equation can be derived by averaging over a very large number of particles. Traditionally, the resulting PDE is studied as a deterministic equation, an approach that has brought many significant results and a deep understanding of the equation and its solutions. By studying the heat equation by considering the individual random particles, however, one gains further intuition into the problem. While this is now standard for many researchers, this approach is generally not presented at the undergraduate level. In this book, Lawler introduces the heat equation and the closely related notion of harmonic functions from a probabilistic perspective. The theme of the first two chapters of the book is the relationship between random walks and the heat equation. The first chapter discusses the discrete case, random walk and the heat equation on the integer lattice; and the second chapter discusses the continuous case, Brownian motion and the usual heat equation. Relationships are shown between the two. For exa...

A method for solving the KDV equation and some numerical experiments

International Nuclear Information System (INIS)

Chang Jinjiang.

1993-01-01

In this paper, by means of difference method for discretization of space partial derivatives of KDV equation, an initial value problem in ordinary differential equations of large dimensions is produced. By using this ordinary differential equations the existence and the uniqueness of the solution of the KDV equation and the conservation of scheme are proved. This ordinary differential equation can be solved by using implicit Runge-Kutta methods, so a new method for finding the numerical solution of the KDV equation is presented. Numerical experiments not only describe in detail the procedure of two solitons collision, soliton reflex and soliton produce, but also show that this method is very effective. (author). 7 refs, 3 figs

Chaotic Universe, Friedmannian on the average 2

Energy Technology Data Exchange (ETDEWEB)

Marochnik, L S [AN SSSR, Moscow. Inst. Kosmicheskikh Issledovanij

1980-11-01

The cosmological solutions are found for the equations for correlators, describing a statistically chaotic Universe, Friedmannian on the average in which delta-correlated fluctuations with amplitudes h >> 1 are excited. For the equation of state of matter p = n epsilon, the kind of solutions depends on the position of maximum of the spectrum of the metric disturbances. The expansion of the Universe, in which long-wave potential and vortical motions and gravitational waves (modes diverging at t ..-->.. 0) had been excited, tends asymptotically to the Friedmannian one at t ..-->.. identity and depends critically on n: at n < 0.26, the solution for the scalefactor is situated higher than the Friedmannian one, and lower at n > 0.26. The influence of finite at t ..-->.. 0 long-wave fluctuation modes leads to an averaged quasiisotropic solution. The contribution of quantum fluctuations and of short-wave parts of the spectrum of classical fluctuations to the expansion law is considered. Their influence is equivalent to the contribution from an ultrarelativistic gas with corresponding energy density and pressure. The restrictions are obtained for the degree of chaos (the spectrum characteristics) compatible with the observed helium abundance, which could have been retained by a completely chaotic Universe during its expansion up to the nucleosynthesis epoch.

Front tracking for hyperbolic conservation laws

CERN Document Server

Holden, Helge

2015-01-01

This is the second edition of a well-received book providing the fundamentals of the theory hyperbolic conservation laws. Several chapters have been rewritten, new material has been added, in particular, a chapter on space dependent flux functions, and the detailed solution of the Riemann problem for the Euler equations. Hyperbolic conservation laws are central in the theory of nonlinear partial differential equations and in science and technology. The reader is given a self-contained presentation using front tracking, which is also a numerical method. The multidimensional scalar case and the case of systems on the line are treated in detail. A chapter on finite differences is included. From the reviews of the first edition: "It is already one of the few best digests on this topic. The present book is an excellent compromise between theory and practice. Students will appreciate the lively and accurate style." D. Serre, MathSciNet  "I have read the book with great pleasure, and I can recommend it to experts ...

A parametrization of two-dimensional turbulence based on a maximum entropy production principle with a local conservation of energy

International Nuclear Information System (INIS)

Chavanis, Pierre-Henri

2014-01-01

In the context of two-dimensional (2D) turbulence, we apply the maximum entropy production principle (MEPP) by enforcing a local conservation of energy. This leads to an equation for the vorticity distribution that conserves all the Casimirs, the energy, and that increases monotonically the mixing entropy (H-theorem). Furthermore, the equation for the coarse-grained vorticity dissipates monotonically all the generalized enstrophies. These equations may provide a parametrization of 2D turbulence. They do not generally relax towards the maximum entropy state. The vorticity current vanishes for any steady state of the 2D Euler equation. Interestingly, the equation for the coarse-grained vorticity obtained from the MEPP turns out to coincide, after some algebraic manipulations, with the one obtained with the anticipated vorticity method. This shows a connection between these two approaches when the conservation of energy is treated locally. Furthermore, the newly derived equation, which incorporates a diffusion term and a drift term, has a nice physical interpretation in terms of a selective decay principle. This sheds new light on both the MEPP and the anticipated vorticity method. (paper)

LETTER TO THE EDITOR: Bicomplexes and conservation laws in non-Abelian Toda models

Science.gov (United States)

Gueuvoghlanian, E. P.

2001-08-01

A bicomplex structure is associated with the Leznov-Saveliev equation of integrable models. The linear problem associated with the zero-curvature condition is derived in terms of the bicomplex linear equation. The explicit example of a non-Abelian conformal affine Toda model is discussed in detail and its conservation laws are derived from the zero-curvature representation of its equation of motion.

Explicit wave action conservation for water waves on vertically sheared flows

Science.gov (United States)

Quinn, Brenda; Toledo, Yaron; Shrira, Victor

2016-04-01

Water waves almost always propagate on currents with a vertical structure such as currents directed towards the beach accompanied by an under-current directed back toward the deep sea or wind-induced currents which change magnitude with depth due to viscosity effects. On larger scales they also change their direction due to the Coriolis force as described by the Ekman spiral. This implies that the existing wave models, which assume vertically-averaged currents, is an approximation which is far from realistic. In recent years, ocean circulation models have significantly improved with the capability to model vertically-sheared current profiles in contrast with the earlier vertically-averaged current profiles. Further advancements have coupled wave action models to circulation models to relate the mutual effects between the two types of motion. Restricting wave models to vertically-averaged non-turbulent current profiles is obviously problematic in these cases and the primary goal of this work is to derive and examine a general wave action equation which accounts for these shortcoming. The formulation of the wave action conservation equation is made explicit by following the work of Voronovich (1976) and using known asymptotic solutions of the boundary value problem which exploit the smallness of the current magnitude compared to the wave phase velocity and/or its vertical shear and curvature. The adopted approximations are shown to be sufficient for most of the conceivable applications. This provides correction terms to the group velocity and wave action definition accounting for the shear effects, which are fitting for application to operational wave models. In the limit of vanishing current shear, the new formulation reduces to the commonly used Bretherton & Garrett (1968) no-shear wave action equation where the invariant is calculated with the current magnitude taken at the free surface. It is shown that in realistic oceanic conditions, the neglect of the vertical

A stable penalty method for the compressible Navier-Stokes equations: I. Open boundary conditions

DEFF Research Database (Denmark)

Hesthaven, Jan; Gottlieb, D.

1996-01-01

The purpose of this paper is to present asymptotically stable open boundary conditions for the numerical approximation of the compressible Navier-Stokes equations in three spatial dimensions. The treatment uses the conservation form of the Navier-Stokes equations and utilizes linearization...

Agriculture’s Soil Conservation Programs Miss Full Potential in the Fight against Soil Erosion.

Science.gov (United States)

1983-11-28

Soil Loss Equation ( USLE ) and Wind Erosion Equation can be used with a reasonable degree of accuracy. It is the intention of ASCS to expand VC/SL to...HD-R37 495 AGRICULTURE’S SOIL CONSERVATION PROGRAMS MISS FULL i/i POTENTIAL IN THE FIGHT.(U) GENERAL ACCOUNTING OFFICE WASHINGTON DC RESOURCES...GENERAL Report To The Congress OF THE UNITED STATES Agriculture’s Soil Conservation Programs Miss Full Potential In The Fight Against Soil Erosion

Nonholonomic deformation of generalized KdV-type equations

International Nuclear Information System (INIS)

Guha, Partha

2009-01-01

Karasu-Kalkani et al (2008 J. Math. Phys. 49 073516) recently derived a new sixth-order wave equation KdV6, which was shown by Kupershmidt (2008 Phys. Lett. 372A 2634) to have an infinite commuting hierarchy with a common infinite set of conserved densities. Incidentally, this equation was written for the first time by Calogero and is included in the book by Calogero and Degasperis (1982 Lecture Notes in Computer Science vol 144 (Amsterdam: North-Holland) p 516). In this paper, we give a geometric insight into the KdV6 equation. Using Kirillov's theory of coadjoint representation of the Virasoro algebra, we show how to obtain a large class of KdV6-type equations equivalent to the original equation. Using a semidirect product extension of the Virasoro algebra, we propose the nonholonomic deformation of the Ito equation. We also show that the Adler-Kostant-Symes scheme provides a geometrical method for constructing nonholonomic deformed integrable systems. Applying the Adler-Kostant-Symes scheme to loop algebra, we construct a new nonholonomic deformation of the coupled KdV equation.

Measuring the equations of state in a relaxed magnetohydrodynamic plasma

Science.gov (United States)

Kaur, M.; Barbano, L. J.; Suen-Lewis, E. M.; Shrock, J. E.; Light, A. D.; Brown, M. R.; Schaffner, D. A.

2018-01-01

We report measurements of the equations of state of a fully relaxed magnetohydrodynamic (MHD) laboratory plasma. Parcels of magnetized plasma, called Taylor states, are formed in a coaxial magnetized plasma gun, and are allowed to relax and drift into a closed flux conserving volume. Density, ion temperature, and magnetic field are measured as a function of time as the Taylor states compress and heat. The theoretically predicted MHD and double adiabatic equations of state are compared to experimental measurements. We find that the MHD equation of state is inconsistent with our data.

Quasiclassical deformation in KP hierarchy and Benney's long wave equations

International Nuclear Information System (INIS)

Kolokol'tsov, V.N.; Lebedev, D.R.

1987-01-01

In the paper by means of the formal variant of Zakharov-Shabat ''dressing'' method various formulas are obtained for the generating functions of the conservation laws of Kadomtsev-Petvias hierarchy which turn into analogous formulas for Benney hierarchy in the quasiclassical limit. The generating fucntion of the conservation laws of Miura type is constructed for higher Benney equations and the simple proof of the related identities is given

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

Reproducing multi-model ensemble average with Ensemble-averaged Reconstructed Forcings (ERF) in regional climate modeling

Science.gov (United States)

Erfanian, A.; Fomenko, L.; Wang, G.

2016-12-01

Multi-model ensemble (MME) average is considered the most reliable for simulating both present-day and future climates. It has been a primary reference for making conclusions in major coordinated studies i.e. IPCC Assessment Reports and CORDEX. The biases of individual models cancel out each other in MME average, enabling the ensemble mean to outperform individual members in simulating the mean climate. This enhancement however comes with tremendous computational cost, which is especially inhibiting for regional climate modeling as model uncertainties can originate from both RCMs and the driving GCMs. Here we propose the Ensemble-based Reconstructed Forcings (ERF) approach to regional climate modeling that achieves a similar level of bias reduction at a fraction of cost compared with the conventional MME approach. The new method constructs a single set of initial and boundary conditions (IBCs) by averaging the IBCs of multiple GCMs, and drives the RCM with this ensemble average of IBCs to conduct a single run. Using a regional climate model (RegCM4.3.4-CLM4.5), we tested the method over West Africa for multiple combination of (up to six) GCMs. Our results indicate that the performance of the ERF method is comparable to that of the MME average in simulating the mean climate. The bias reduction seen in ERF simulations is achieved by using more realistic IBCs in solving the system of equations underlying the RCM physics and dynamics. This endows the new method with a theoretical advantage in addition to reducing computational cost. The ERF output is an unaltered solution of the RCM as opposed to a climate state that might not be physically plausible due to the averaging of multiple solutions with the conventional MME approach. The ERF approach should be considered for use in major international efforts such as CORDEX. Key words: Multi-model ensemble, ensemble analysis, ERF, regional climate modeling

Stated Preferences for Forest Conservation in Southern Finland

DEFF Research Database (Denmark)

Lehtonen, Emmi; Kuuluvainen, Jari; Pouta, Eija

% support decreasing forest conservation. An average willingness-to-pay for increased biodiversity conservation was 60-212 € per household per year, depending on the described project and measurement method. In addition to costs per household, the number of conserved biotopes and endangered plant and animal......This study analyses Finnish citizens’ valuations and attitudes towards a forest conservation programme for southern Finland and the Pohjanmaa region. In particular, Finnish households’ willingness to accept expenses through increased taxation to guarantee a certain level of biodiversity...... conservation was investigated. Contingent valuation (CV) and choice experiment (CE) methods were applied. According to the CV results, 74% of respondents are prepared to pay for increased conservation and 16% support increased conservation but are not willing to pay for it. A further 5% are indifferent and 5...

Unsplit schemes for hyperbolic conservation laws with source terms in one space dimension

International Nuclear Information System (INIS)

Papalexandris, M.V.; Leonard, A.; Dimotakis, P.E.

1997-01-01

The present work is concerned with an application of the theory of characteristics to conservation laws with source terms in one space dimension, such as the Euler equations for reacting flows. Space-time paths are introduced on which the flow/chemistry equations decouple to a characteristic set of ODE's for the corresponding homogeneous laws, thus allowing the introduction of functions analogous to the Riemann invariants in classical theory. The geometry of these paths depends on the spatial gradients of the solution. This particular decomposition can be used in the design of efficient unsplit algorithms for the numerical integration of the equations. As a first step, these ideas are implemented for the case of a scalar conservation law with a nonlinear source term. The resulting algorithm belongs to the class of MUSCL-type, shock-capturing schemes. Its accuracy and robustness are checked through a series of tests. The stiffness of the source term is also studied. Then, the algorithm is generalized for a system of hyperbolic equations, namely the Euler equations for reacting flows. A numerical study of unstable detonations is performed. 57 refs

Duality properties of Gorringe Leach equations

Science.gov (United States)

Grandati, Yves; Bérard, Alain; Mohrbach, Hervé

2009-02-01

In the category of motions preserving the angular momentum direction, Gorringe and Leach exhibited two classes of differential equations having elliptical orbits. After enlarging slightly these classes, we show that they are related by a duality correspondence of the Arnold Vassiliev type. The specific associated conserved quantities (Laplace Runge Lenz vector and Fradkin Jauch Hill tensor) are then dual reflections of each other.

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.

Continuum analysis of biological systems conserved quantities, fluxes and forces

CERN Document Server

Suraishkumar, G K

2014-01-01

This book addresses the analysis, in the continuum regime, of biological systems at various scales, from the cellular level to the industrial one. It presents both fundamental conservation principles (mass, charge, momentum and energy) and relevant fluxes resulting from appropriate driving forces, which are important for the analysis, design and operation of biological systems. It includes the concept of charge conservation, an important principle for biological systems that is not explicitly covered in any other book of this kind. The book is organized in five parts: mass conservation; charge conservation; momentum conservation; energy conservation; and multiple conservations simultaneously applied. All mathematical aspects are presented step by step, allowing any reader with a basic mathematical background (calculus, differential equations, linear algebra, etc.) to follow the text with ease. The book promotes an intuitive understanding of all the relevant principles and in so doing facilitates their applica...

A spectral element-FCT method for the compressible Euler equations

International Nuclear Information System (INIS)

Giannakouros, J.; Karniadakis, G.E.

1994-01-01

A new algorithm based on spectral element discretizations and flux-corrected transport concepts is developed for the solution of the Euler equations of inviscid compressible fluid flow. A conservative formulation is proposed based on one- and two-dimensional cell-averaging and reconstruction procedures, which employ a staggered mesh of Gauss-Chebyshev and Gauss-Lobatto-Chebyshev collocation points. Particular emphasis is placed on the construction of robust boundary and interfacial conditions in one- and two-dimensions. It is demonstrated through shock-tube problems and two-dimensional simulations that the proposed algorithm leads to stable, non-oscillatory solutions of high accuracy. Of particular importance is the fact that dispersion errors are minimal, as show through experiments. From the operational point of view, casting the method in a spectral element formulation provides flexibility in the discretization, since a variable number of macro-elements or collocation points per element can be employed to accomodate both accuracy and geometric requirements

Discretely Conservative Finite-Difference Formulations for Nonlinear Conservation Laws in Split Form: Theory and Boundary Conditions

Science.gov (United States)

Fisher, Travis C.; Carpenter, Mark H.; Nordstroem, Jan; Yamaleev, Nail K.; Swanson, R. Charles

2011-01-01

Simulations of nonlinear conservation laws that admit discontinuous solutions are typically restricted to discretizations of equations that are explicitly written in divergence form. This restriction is, however, unnecessary. Herein, linear combinations of divergence and product rule forms that have been discretized using diagonal-norm skew-symmetric summation-by-parts (SBP) operators, are shown to satisfy the sufficient conditions of the Lax-Wendroff theorem and thus are appropriate for simulations of discontinuous physical phenomena. Furthermore, special treatments are not required at the points that are near physical boundaries (i.e., discrete conservation is achieved throughout the entire computational domain, including the boundaries). Examples are presented of a fourth-order, SBP finite-difference operator with second-order boundary closures. Sixth- and eighth-order constructions are derived, and included in E. Narrow-stencil difference operators for linear viscous terms are also derived; these guarantee the conservative form of the combined operator.

A high-order relaxation method with projective integration for solving nonlinear systems of hyperbolic conservation laws

Science.gov (United States)

Lafitte, Pauline; Melis, Ward; Samaey, Giovanni

2017-07-01

We present a general, high-order, fully explicit relaxation scheme which can be applied to any system of nonlinear hyperbolic conservation laws in multiple dimensions. The scheme consists of two steps. In a first (relaxation) step, the nonlinear hyperbolic conservation law is approximated by a kinetic equation with stiff BGK source term. Then, this kinetic equation is integrated in time using a projective integration method. After taking a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution, the time derivative is estimated and used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time step restriction on the outer time step is similar to the CFL condition for the hyperbolic conservation law. Moreover, the number of inner time steps is also independent of the stiffness of the BGK source term. We discuss stability and consistency, and illustrate with numerical results (linear advection, Burgers' equation and the shallow water and Euler equations) in one and two spatial dimensions.

Nuclear fuel management via fuel quality factor averaging

International Nuclear Information System (INIS)

Mingle, J.O.

1978-01-01

The numerical procedure of prime number averaging is applied to the fuel quality factor distribution of once and twice-burned fuel in order to evolve a fuel management scheme. The resulting fuel shuffling arrangement produces a near optimal flat power profile both under beginning-of-life and end-of-life conditions. The procedure is easily applied requiring only the solution of linear algebraic equations. (author)

Collisional drift fluid equations and implications for drift waves

International Nuclear Information System (INIS)

Pfirsch, Dieter; Correa-Restrepo, Dario

1996-01-01

The usual theoretical description of drift-wave turbulence (considered to be one possible cause of anomalous transport in a plasma), e.g. the Hasegawa-Wakatani theory, makes use of various approximations, the effects of which are extremely difficult to assess. This concerns in particular the conservation laws for energy and momentum. The latter law is important in relation to charge separation and the resulting electric fields, which are possibly related to the L-H transition. Energy conservation is crucial to the stability behaviour, it will be discussed by means of an example. New collisional multi-species drift-fluid equations were derived by a new method which yields, in a transparent way, conservation of energy and total angular momentum and the law for energy dissipation. Both electrostatic and electromagnetic field variations are considered. The only restriction involved is the validity of the drift approximation; in particular, there are no assumptions restricting the geometry of the system. The method is based primarily on a Lagrangian for dissipationless fluids in the drift approximation with isotropic pressures. The dissipative terms are introduced by adding corresponding terms to the ideal equations of motion and of the pressures. The equations of motion, of course, no longer result from a Lagrangian via Hamilton's principle. However, their relation to the ideal equations also implies a relation to the ideal Lagrangian, which can be used to advantage. Instead of introducing heat conduction one can also assume isothermal behaviour, e.g. T v (x) = constant. Assumptions of this kind are often made in the literature. The new method of introducing dissipation is not restricted to the present kind of theory; it can equally well be applied to theories such as multi-fluid theories without using the drift approximation of the present paper. (author)

A direct Arbitrary-Lagrangian-Eulerian ADER-WENO finite volume scheme on unstructured tetrahedral meshes for conservative and non-conservative hyperbolic systems in 3D

Science.gov (United States)

Boscheri, Walter; Dumbser, Michael

2014-10-01

In this paper we present a new family of high order accurate Arbitrary-Lagrangian-Eulerian (ALE) one-step ADER-WENO finite volume schemes for the solution of nonlinear systems of conservative and non-conservative hyperbolic partial differential equations with stiff source terms on moving tetrahedral meshes in three space dimensions. A WENO reconstruction technique is used to achieve high order of accuracy in space, while an element-local space-time Discontinuous Galerkin finite element predictor on moving curved meshes is used to obtain a high order accurate one-step time discretization. Within the space-time predictor the physical element is mapped onto a reference element using a high order isoparametric approach, where the space-time basis and test functions are given by the Lagrange interpolation polynomials passing through a predefined set of space-time nodes. Since our algorithm is cell-centered, the final mesh motion is computed by using a suitable node solver algorithm. A rezoning step as well as a flattener strategy are used in some of the test problems to avoid mesh tangling or excessive element deformations that may occur when the computation involves strong shocks or shear waves. The ALE algorithm presented in this article belongs to the so-called direct ALE methods because the final Lagrangian finite volume scheme is based directly on a space-time conservation formulation of the governing PDE system, with the rezoned geometry taken already into account during the computation of the fluxes. We apply our new high order unstructured ALE schemes to the 3D Euler equations of compressible gas dynamics, for which a set of classical numerical test problems has been solved and for which convergence rates up to sixth order of accuracy in space and time have been obtained. We furthermore consider the equations of classical ideal magnetohydrodynamics (MHD) as well as the non-conservative seven-equation Baer-Nunziato model of compressible multi-phase flows with

A Novel Partial Differential Algebraic Equation (PDAE) Solver

DEFF Research Database (Denmark)

Lim, Young-il; Chang, Sin-Chung; Jørgensen, Sten Bay

2004-01-01

For solving partial differential algebraic equations (PDAEs), the space-time conservation element/solution element (CE/SE) method is addressed in this study. The method of lines (MOL) using an implicit time integrator is compared with the CE/SE method in terms of computational efficiency, solution...... or nonlinear adsorption isotherm are solved by the two methods. The CE/SE method enforces both local and global flux conservation in space and time, and uses a simple stencil structure (two points at the previous time level and one point at the present time level). Thus, accurate and computationally...

Experimental Verification of Current Shear Design Equations for HSRC Beams

Directory of Open Access Journals (Sweden)

Attaullah Shah

2012-07-01

Full Text Available Experimental research on the shear capacity of HSRC (High Strength Reinforced Concrete beams is relatively very limited as compared to the NSRC (Normal Strength Reinforced Concrete beams. Most of the Building Codes determine the shear strength of HSRC with the help of empirical equations based on experimental work of NSRC beams and hence these equations are generally regarded as un-conservative for HSRC beams particularly at low level of longitudinal reinforcement. In this paper, 42 beams have been tested in two sets, such that in 21 beams no transverse reinforcement has been used, whereas in the remaining 21 beams, minimum transverse reinforcement has been used as per ACI-318 (American Concrete Institute provisions. Two values of compressive strength 52 and 61 MPa, three values of longitudinal steel ratio and seven values of shear span to depth ratio have been have been used. The beams were tested under concentrated load at the mid span. The results are compared with the equations proposed by different international building codes like ACI, AASHTO LRFD, EC (Euro Code, Canadian Code and Japanese Code for shear strength of HSRC beams.From comparison, it has been observed that some codes are less conservative for shear design of HSRC beams and further research is required to rationalize these equations.

Analysis of the average daily radon variations in the soil air

International Nuclear Information System (INIS)

Holy, K.; Matos, M.; Boehm, R.; Stanys, T.; Polaskova, A.; Hola, O.

1998-01-01

In this contribution the search of the relation between the daily variations of the radon concentration and the regular daily oscillations of the atmospheric pressure are presented. The deviation of the radon activity concentration in the soil air from the average daily value reaches only a few percent. For the dry summer months the average daily course of the radon activity concentration can be described by the obtained equation. The analysis of the average daily courses could give the information concerning the depth of the gas permeable soil layer. The soil parameter is determined by others method with difficulty

Fully Electromagnetic Nonlinear Gyrokinetic Equations for Tokamak Edge Turbulence

International Nuclear Information System (INIS)

Hahm, T.S.; Wang, Lu; Madsen, J.

2008-01-01

An energy conserving set of the fully electromagnetic nonlinear gyrokinetic Vlasov equation and Maxwell's equations, which is applicable to both L-mode turbulence with large amplitude and H-mode turbulence in the presence of high E x B shear has been derived. The phase-space action variational Lie perturbation method ensures the preservation of the conservation laws of the underlying Vlasov-Maxwell system. Our generalized ordering takes ρ i θi ∼ L E ∼ L p i is the thermal ion Larmor radius and ρ θi = B/B θ ρ i ), as typically observed in the tokamak H-mode edge, with L E and L p being the radial electric field and pressure gradient lengths. We take k # perpendicular# ρ i ∼ 1 for generality, and keep the relative fluctuation amplitudes e(delta)φ/T i ∼ (delta)B/B up to the second order. Extending the electrostatic theory in the presence of high E x B shear [Hahm, Phys. Plasmas 3, 4658 (1996)], contributions of electromagnetic fluctuations to the particle charge density and current are explicitly evaluated via pull-back transformation from the gyrocenter distribution function in the gyrokinetic Maxwell's equation

Estimating order-picking times for return heuristic - equations and simulations

Directory of Open Access Journals (Sweden)

Grzegorz Tarczyński

2015-09-01

Full Text Available Background: A key element of the evaluation of warehouse operation is the average order-picking time. In warehouses where the order-picking process is carried out according to the "picker-to-part" rule the order-picking time is usually proportional to the distance covered by the picker while picking items. This distance can by estimated by simulations or using mathematical equations. In the paper only the best described in the literature one-block rectangular warehouses are considered. Material and methods: For the one-block rectangular warehouses there are well known five routing heuristics. In the paper the author considers the return heuristic in two variants. The paper presents well known Hall's and De Koster's equations for the average distance traveled by the picker while completing items from one pick list. The author presents own proposals for calculating the expected distance. Results: the results calculated by the use of mathematical equations (the formulas of Hall, De Koster and own propositions were compared with the average values obtained using computer simulations. For the most cases the average error does not exceed 1% (except for Hall's equations. To carry out simulation the computer software Warehouse Real-Time Simulator was used. Conclusions: the order-picking time is a function of many variables and its optimization is not easy. It can be done in two stages: firstly using mathematical equations the set of the potentially best variants is established, next the results are verified using simulations. The results calculated by the use of equations are not precise, but possible to achieve immediately. The simulations are more time-consuming, but allow to analyze the order-picking process more accurately.

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

Why not energy conservation?

International Nuclear Information System (INIS)

Carlson, Shawn

2016-01-01

Energy conservation is a deep principle that is obeyed by all of the fundamental forces of nature. It puts stringent constraints on all systems, particularly systems that are ‘isolated,’ meaning that no energy can enter or escape. Notwithstanding the success of the principle of stationary action, it is fair to wonder to what extent physics can be formulated from the principle of stationary energy. We show that if one interprets mechanical energy as a state function, then its stationarity leads to a novel formulation of classical mechanics. However, unlike Lagrangian and Hamiltonian mechanics, which deliver their state functions via algebraic proscriptions (i.e., the Lagrangian is always the difference between a system’s kinetic and potential energies), this new formalism identifies its state functions as the solutions to a differential equation. This is an important difference because differential equations can generate more general solutions than algebraic recipes. When applied to Newtonian systems for which the energy function is separable, these state functions are always the mechanical energy. However, while the stationary state function for a charged particle moving in an electromagnetic field proves not to be energy, the function nevertheless correctly encodes the dynamics of the system. Moreover, the stationary state function for a free relativistic particle proves not to be the energy either. Rather, our differential equation yields the relativistic free-particle Lagrangian (plus a non-dynamical constant) in its correct dynamical context. To explain how this new formalism can consistently deliver stationary state functions that give the correct dynamics but that are not always the mechanical energy, we propose that energy conservation is a specific realization of a deeper principle of stationarity that governs both relativistic and non-relativistic mechanics. (paper)

Ordinary differential equations a graduate text

CERN Document Server

Bhamra, K S

2015-01-01

ORDINARY DIFFERENTIAL EQUATIONS: A Graduate Text presents a systematic and comprehensive introduction to ODEs for graduate and postgraduate students. The systematic organized text on differential inequalities, Gronwall's inequality, Nagumo's theorems, Osgood's criteria and applications of different equations of first order is dealt with in a greater depth. The book discusses qualitative and quantitative aspects of the Strum - Liouville problems, Green's function, integral equations, Laplace transform and is supported by a number of worked-out examples in each lesson to make the concepts clear. A lot of stress on stability theory is laid down, especially on Lyapunov and Poincare stability theory. A numerous figures in various lessons (in particular lessons dealing with stability theory) have been added to clarify the key concepts in DE theory. Nonlinear oscillation in conservative systems and Hamiltonian systems highlights basic nature of the systems considered. Perturbation techniques lesson deals in fairly d...

Land degradation and adoption of soil conservation technologies ...

African Journals Online (AJOL)

The study investigates the causes of land degradation, and adoption of soil conservation practices using a two-stage decision making process. The data for the study were collected with the aid of structured questionnaire and analyzed with descriptive analysis, difference regression equation and simultaneous probit model.

Waves, conservation laws and symmetries of a third-order nonlinear ...

African Journals Online (AJOL)

order is under consideration. Important properties concerning advanced character such like conservation laws and the equation of continuity are given. Characteristic wave properties such like dispersion relations and both the group and phase ...

An improved front tracking method for the Euler equations

NARCIS (Netherlands)

Witteveen, J.A.S.; Koren, B.; Bakker, P.G.

2007-01-01

An improved front tracking method for hyperbolic conservation laws is presented. The improved method accurately resolves discontinuities as well as continuous phenomena. The method is based on an improved front interaction model for a physically more accurate modeling of the Euler equations, as

Energy conservation in Newmark based time integration algorithms

DEFF Research Database (Denmark)

Krenk, Steen

2006-01-01

Energy balance equations are established for the Newmark time integration algorithm, and for the derived algorithms with algorithmic damping introduced via averaging, the so-called a-methods. The energy balance equations form a sequence applicable to: Newmark integration of the undamped equations...... of motion, an extended form including structural damping, and finally the generalized form including structural as well as algorithmic damping. In all three cases the expression for energy, appearing in the balance equation, is the mechanical energy plus some additional terms generated by the discretization...

Materials with memory initial-boundary value problems for constitutive equations with internal variables

CERN Document Server

Alber, Hans-Dieter

1998-01-01

This book contributes to the mathematical theory of systems of differential equations consisting of the partial differential equations resulting from conservation of mass and momentum, and of constitutive equations with internal variables. The investigations are guided by the objective of proving existence and uniqueness, and are based on the idea of transforming the internal variables and the constitutive equations. A larger number of constitutive equations from the engineering sciences are presented. The book is therefore suitable not only for specialists, but also for mathematicians seeking for an introduction in the field, and for engineers with a sound mathematical background.

Generalized reduced MHD equations

International Nuclear Information System (INIS)

Kruger, S.E.; Hegna, C.C.; Callen, J.D.

1998-07-01

A new derivation of reduced magnetohydrodynamic (MHD) equations is presented. A multiple-time-scale expansion is employed. It has the advantage of clearly separating the three time scales of the problem associated with (1) MHD equilibrium, (2) fluctuations whose wave vector is aligned perpendicular to the magnetic field, and (3) those aligned parallel to the magnetic field. The derivation is carried out without relying on a large aspect ratio assumption; therefore this model can be applied to any general toroidal configuration. By accounting for the MHD equilibrium and constraints to eliminate the fast perpendicular waves, equations are derived to evolve scalar potential quantities on a time scale associated with the parallel wave vector (shear-alfven wave time scale), which is the time scale of interest for MHD instability studies. Careful attention is given in the derivation to satisfy energy conservation and to have manifestly divergence-free magnetic fields to all orders in the expansion parameter. Additionally, neoclassical closures and equilibrium shear flow effects are easily accounted for in this model. Equations for the inner resistive layer are derived which reproduce the linear ideal and resistive stability criterion of Glasser, Greene, and Johnson

Hyperbolic Conservation Laws and Related Analysis with Applications

CERN Document Server

Holden, Helge; Karlsen, Kenneth

2014-01-01

This book presents thirteen papers, representing the most significant advances and current trends in nonlinear hyperbolic conservation laws and related analysis with applications. Topics covered include a survey on multidimensional systems of conservation laws as well as novel results  on liquid crystals, conservation laws with discontinuous flux functions, and applications to sedimentation.  Also included are articles on recent advances in the Euler equations and the Navier-Stokes-Fourier-Poisson system, in addition to new results on collective phenomena described by the Cucker-Smale model.    The Workshop on Hyperbolic Conservation Laws and Related Analysis with Applications at the International Centre for Mathematical Sciences (Edinburgh, UK) held in Edinburgh, September 2011, produced this fine collection of original research and survey articles. Many leading mathematicians attended the event and submitted their contributions for this volume. It is addressed to researchers and graduate students inter...

Hamilton's equations for a fluid membrane

International Nuclear Information System (INIS)

Capovilla, R; Guven, J; Rojas, E

2005-01-01

Consider a homogeneous fluid membrane described by the Helfrich-Canham energy, quadratic in the mean curvature of the membrane surface. The shape equation that determines equilibrium configurations is fourth order in derivatives and cubic in the mean curvature. We introduce a Hamiltonian formulation of this equation which dismantles it into a set of coupled first-order equations. This involves interpreting the Helfrich-Canham energy as an action; equilibrium surfaces are generated by the evolution of space curves. Two features complicate the implementation of a Hamiltonian framework. (i) The action involves second derivatives. This requires treating the velocity as a phase-space variable and the introduction of its conjugate momentum. The canonical Hamiltonian is constructed on this phase space. (ii) The action possesses a local symmetry-reparametrization invariance. The two labels we use to parametrize points on the surface are themselves physically irrelevant. This symmetry implies primary constraints, one for each label, that need to be implemented within the Hamiltonian. The two Lagrange multipliers associated with these constraints are identified as the components of the acceleration tangential to the surface. The conservation of the primary constraints implies two secondary constraints, fixing the tangential components of the momentum conjugate to the position. Hamilton's equations are derived and the appropriate initial conditions on the phase-space variables are identified. Finally, it is shown how the shape equation can be reconstructed from these equations

Taguchi method for partial differential equations with application in tumor growth.

Science.gov (United States)

Ilea, M; Turnea, M; Rotariu, M; Arotăriţei, D; Popescu, Marilena

2014-01-01

The growth of tumors is a highly complex process. To describe this process, mathematical models are needed. A variety of partial differential mathematical models for tumor growth have been developed and studied. Most of those models are based on the reaction-diffusion equations and mass conservation law. A variety of modeling strategies have been developed, each focusing on tumor growth. Systems of time-dependent partial differential equations occur in many branches of applied mathematics. The vast majority of mathematical models in tumor growth are formulated in terms of partial differential equations. We propose a mathematical model for the interactions between these three cancer cell populations. The Taguchi methods are widely used by quality engineering scientists to compare the effects of multiple variables, together with their interactions, with a simple and manageable experimental design. In Taguchi's design of experiments, variation is more interesting to study than the average. First, Taguchi methods are utilized to search for the significant factors and the optimal level combination of parameters. Except the three parameters levels, other factors levels other factors levels would not be considered. Second, cutting parameters namely, cutting speed, depth of cut, and feed rate are designed using the Taguchi method. Finally, the adequacy of the developed mathematical model is proved by ANOVA. According to the results of ANOVA, since the percentage contribution of the combined error is as small. Many mathematical models can be quantitatively characterized by partial differential equations. The use of MATLAB and Taguchi method in this article illustrates the important role of informatics in research in mathematical modeling. The study of tumor growth cells is an exciting and important topic in cancer research and will profit considerably from theoretical input. Interpret these results to be a permanent collaboration between math's and medical oncologists.

On Conservation Forms and Invariant Solutions for Classical Mechanics Problems of Liénard Type

Directory of Open Access Journals (Sweden)

Gülden Gün Polat

2014-01-01

Full Text Available In this study we apply partial Noether and λ-symmetry approaches to a second-order nonlinear autonomous equation of the form y′′+fyy′+g(y=0, called Liénard equation corresponding to some important problems in classical mechanics field with respect to f(y and g(y functions. As a first approach we utilize partial Lagrangians and partial Noether operators to obtain conserved forms of Liénard equation. Then, as a second approach, based on the λ-symmetry method, we analyze λ-symmetries for the case that λ-function is in the form of λ(x,y,y′=λ1(x,yy′+λ2(x,y. Finally, a classification problem for the conservation forms and invariant solutions are considered.

Construction of Difference Equations Using Lie Groups

International Nuclear Information System (INIS)

Axford, R.A.

1998-01-01

The theory of prolongations of the generators of groups of point transformations to the grid point values of dependent variables and grid spacings is developed and applied to the construction of group invariant numerical algorithms. The concepts of invariant difference operators and generalized discrete sources are introduced for the discretization of systems of inhomogeneous differential equations and shown to produce exact difference equations. Invariant numerical flux functions are constructed from the general solutions of first order partial differential equations that come out of the evaluation of the Lie derivatives of conservation forms of difference schemes. It is demonstrated that invariant numerical flux functions with invariant flux or slope limiters can be determined to yield high resolution difference schemes. The introduction of an invariant flux or slope limiter can be done so as not to break the symmetry properties of a numerical flux-function

Assessment of soil erosion and conservation on agricultural sloping lands using plot data in the semi-arid hilly loess region of China

Directory of Open Access Journals (Sweden)

T.X. Zhu

2014-11-01

New hydrological insights for the region: The results revealed that runoff per unit area slightly increased with slope angle on SSP, but reached a maximum at 15° and then decreased with slope angle on LSP. Soil loss per unit area increased with slope angle on both SSP and LSP. An average of 36.4% less runoff but only 3.6% less soil loss per unit area was produced on LSP than on SSP. The S factor calculated using the slope factor equations in USLE/RUSLE was significantly greater than that estimated from the measured soil loss on the plots. Rainstorms with recurrence intervals greater than 2 years were responsible for more than two thirds of the total soil and water loss. The effectiveness in reducing surface runoff by five types of conservation practices was mixed. However, all the conservation practices yielded much less soil loss than cropland.

Massively parallel computation of conservation laws

Energy Technology Data Exchange (ETDEWEB)

Garbey, M [Univ. Claude Bernard, Villeurbanne (France); Levine, D [Argonne National Lab., IL (United States)

1990-01-01

The authors present a new method for computing solutions of conservation laws based on the use of cellular automata with the method of characteristics. The method exploits the high degree of parallelism available with cellular automata and retains important features of the method of characteristics. It yields high numerical accuracy and extends naturally to adaptive meshes and domain decomposition methods for perturbed conservation laws. They describe the method and its implementation for a Dirichlet problem with a single conservation law for the one-dimensional case. Numerical results for the one-dimensional law with the classical Burgers nonlinearity or the Buckley-Leverett equation show good numerical accuracy outside the neighborhood of the shocks. The error in the area of the shocks is of the order of the mesh size. The algorithm is well suited for execution on both massively parallel computers and vector machines. They present timing results for an Alliant FX/8, Connection Machine Model 2, and CRAY X-MP.

Green`s function of Maxwell`s equations and corresponding implications for iterative methods

Energy Technology Data Exchange (ETDEWEB)

Singer, B.S. [Macquarie Univ., Sydney (Australia); Fainberg, E.B. [Inst. of Physics of the Earth, Moscow (Russian Federation)

1996-12-31

Energy conservation law imposes constraints on the norm and direction of the Hilbert space vector representing a solution of Maxwell`s equations. In this paper, we derive these constrains and discuss the corresponding implications for the Green`s function of Maxwell`s equations in a dissipative medium. It is shown that Maxwell`s equations can be reduced to an integral equation with a contracting kernel. The equation can be solved using simple iterations. Software based on this algorithm have successfully been applied to a wide range of problems dealing with high contrast models. The matrix corresponding to the integral equation has a well defined spectrum. The equation can be symmetrized and solved using different approaches, for instance one of the conjugate gradient methods.

Coarse-grained forms for equations describing the microscopic motion of particles in a fluid.

Science.gov (United States)

Das, Shankar P; Yoshimori, Akira

2013-10-01

Exact equations of motion for the microscopically defined collective density ρ(x,t) and the momentum density ĝ(x,t) of a fluid have been obtained in the past starting from the corresponding Langevin equations representing the dynamics of the fluid particles. In the present work we average these exact equations of microscopic dynamics over the local equilibrium distribution to obtain stochastic partial differential equations for the coarse-grained densities with smooth spatial and temporal dependence. In particular, we consider Dean's exact balance equation for the microscopic density of a system of interacting Brownian particles to obtain the basic equation of the dynamic density functional theory with noise. Our analysis demonstrates that on thermal averaging the dependence of the exact equations on the bare interaction potential is converted to dependence on the corresponding thermodynamic direct correlation functions in the coarse-grained equations.

Relative efforts of countries to conserve world’s megafauna

Directory of Open Access Journals (Sweden)

Peter A. Lindsey

2017-04-01

Full Text Available Surprisingly little attention has been paid to variation among countries in contributions to conservation. As a first step, we developed a Megafauna Conservation Index (MCI that assesses the spatial, ecological and financial contributions of 152 nations towards conservation of the world’s terrestrial megafauna. We chose megafauna because they are particularly valuable in economic, ecological and societal terms, and are challenging and expensive to conserve. We categorised these 152 countries as being above- or below-average performers based on whether their contribution to megafauna conservation was higher or lower than the global mean; ‘major’ performers or underperformers were those whose contribution exceeded 1 SD over or under the mean, respectively. Ninety percent of countries in North/Central America and 70% of countries in Africa were classified as major or above-average performers, while approximately one-quarter of countries in Asia (25% and Europe (21% were identified as major underperformers. We present our index to emphasise the need for measuring conservation performance, to help nations identify how best they could improve their efforts, and to present a starting point for the development of more robust and inclusive measures (noting how the IUCN Red List evolved over time. Our analysis points to three approaches that countries could adopt to improve their contribution to global megafauna conservation, depending on their circumstances: (1 upgrading or expanding their domestic protected area networks, with a particular emphasis on conserving large carnivore and herbivore habitat, (2 increase funding for conservation at home or abroad, or (3 ‘rewilding’ their landscapes. Once revised and perfected, we recommend publishing regular conservation rankings in the popular media to recognise major-performers, foster healthy pride and competition among nations, and identify ways for governments to improve their performance.

Computation techniques and computer programs to analyze Stirling cycle engines using characteristic dynamic energy equations

Science.gov (United States)

Larson, V. H.

1982-01-01

The basic equations that are used to describe the physical phenomena in a Stirling cycle engine are the general energy equations and equations for the conservation of mass and conversion of momentum. These equations, together with the equation of state, an analytical expression for the gas velocity, and an equation for mesh temperature are used in this computer study of Stirling cycle characteristics. The partial differential equations describing the physical phenomena that occurs in a Stirling cycle engine are of the hyperbolic type. The hyperbolic equations have real characteristic lines. By utilizing appropriate points along these curved lines the partial differential equations can be reduced to ordinary differential equations. These equations are solved numerically using a fourth-fifth order Runge-Kutta integration technique.

Variational principle for nonlinear gyrokinetic Vlasov--Maxwell equations

International Nuclear Information System (INIS)

Brizard, Alain J.

2000-01-01

A new variational principle for the nonlinear gyrokinetic Vlasov--Maxwell equations is presented. This Eulerian variational principle uses constrained variations for the gyrocenter Vlasov distribution in eight-dimensional extended phase space and turns out to be simpler than the Lagrangian variational principle recently presented by H. Sugama [Phys. Plasmas 7, 466 (2000)]. A local energy conservation law is then derived explicitly by the Noether method. In future work, this new variational principle will be used to derive self-consistent, nonlinear, low-frequency Vlasov--Maxwell bounce-gyrokinetic equations, in which the fast gyromotion and bounce-motion time scales have been eliminated

International Conference on Differential Equations and Mathematical Physics

CERN Document Server

Saitō, Yoshimi

1987-01-01

The meeting in Birmingham, Alabama, provided a forum for the discussion of recent developments in the theory of ordinary and partial differential equations, both linear and non-linear, with particular reference to work relating to the equations of mathematical physics. The meeting was attended by about 250 mathematicians from 22 countries. The papers in this volume all involve new research material, with at least outline proofs; some papers also contain survey material. Topics covered include: Schrödinger theory, scattering and inverse scattering, fluid mechanics (including conservative systems and inertial manifold theory attractors), elasticity, non-linear waves, and feedback control theory.

Application of autoregressive moving average model in reactor noise analysis

International Nuclear Information System (INIS)

Tran Dinh Tri

1993-01-01

The application of an autoregressive (AR) model to estimating noise measurements has achieved many successes in reactor noise analysis in the last ten years. The physical processes that take place in the nuclear reactor, however, are described by an autoregressive moving average (ARMA) model rather than by an AR model. Consequently more correct results could be obtained by applying the ARMA model instead of the AR model to reactor noise analysis. In this paper the system of the generalised Yule-Walker equations is derived from the equation of an ARMA model, then a method for its solution is given. Numerical results show the applications of the method proposed. (author)

A closed set of conservation laws and the evolution of the electron magnetic moment in the collisionless solar wind

International Nuclear Information System (INIS)

Alexander, P.

1993-01-01

A hydromagnetic equation system for the interplanetary collisionless solar wind is used to derive a set of conservation laws for that medium. It is found that every equation of the original system, including the closure relation, is related to one conservation law. The set that has been derived does not only include the traditional laws, but also a new one for the magnetic moment of the electrons. The conservation set is then used to obtain the space constants for the solar coronal expansion. The new law yields a constant that has not been predicted by other models

Conservative flight with a varying load factor and closed form ...

Indian Academy of Sciences (India)

Conservative flight performance of an aircraft with constant load factor was analysed by ... Within the frame work of flat earth hypotheses the equations of motion of an aircraft as obtained by ..... load factor function if this inequality holds good.

Angular momentum conservation for uniformly expanding flows

International Nuclear Information System (INIS)

Hayward, Sean A

2007-01-01

Angular momentum has recently been defined as a surface integral involving an axial vector and a twist 1-form, which measures the twisting around the spacetime due to a rotating mass. The axial vector is chosen to be a transverse, divergence-free, coordinate vector, which is compatible with any initial choice of axis and integral curves. Then a conservation equation expresses the rate of the change of angular momentum along a uniformly expanding flow as a surface integral of angular momentum densities, with the same form as the standard equation for an axial Killing vector, apart from the inclusion of an effective energy tensor for gravitational radiation

Higher conservation laws for ten-dimensional supersymmetric Yang-Mills theories

International Nuclear Information System (INIS)

Abdalla, E.; Forger, M.; Freiburg Univ.; Jacques, M.

1988-01-01

It is shown that ten-dimensional supersymmetric Yang-Mills theories are integrable systems, in the (weak) sense of admitting a (superspace) Lax representation for their equations of motion. This is achieved by means of an explicit proof that the equations of motion are not only a consequence of but in fact fully equivalent to the superspace constraint F αβ =0. Moreover, a procedure for deriving infinite series of non-local conservation laws is outlined. (orig.)

Ensemble averaged two-phase flow numerical simulation in vertical ducts for the void-studying behavior in BWRs

International Nuclear Information System (INIS)

Mohsen Sharifpur; Mahmoud Salehi; Ali Nouri Brojerdi; Ali Arefmanesh

2003-01-01

Investigation upon generation of vapor in the two-phase flow and predication of its behaviour is an important problem in nuclear industries. Here, the use of the ensemble averaging is to drive the governing equations for each phase in the bubbly two phase flow (two fluid model) and to simulate the water channel inside the four fuel rods along the vertical line. The governing equations will be simplified by having the experience on BWRs and data, which are obtained to find the distribution of void fraction, velocity and other parameters for each phase along the tube. Finally, we compare the results with the simulated results obtained from RELAP 5 Mode 2. The advantage of this work is to offer a new technique to solve the ensemble averaged two-phase flow by imposing the energy balance equation rather than to use the ordinary energy equations. (author)

Ex situ conservation of plant diversity in the world's botanic gardens.

Science.gov (United States)

Mounce, Ross; Smith, Paul; Brockington, Samuel

2017-10-01

Botanic gardens conserve plant diversity ex situ and can prevent extinction through integrated conservation action. Here we quantify how that diversity is conserved in ex situ collections across the world's botanic gardens. We reveal that botanic gardens manage at least 105,634 species, equating to 30% of all plant species diversity, and conserve over 41% of known threatened species. However, we also reveal that botanic gardens are disproportionately temperate, with 93% of species held in the Northern Hemisphere. Consequently, an estimated 76% of species absent from living collections are tropical in origin. Furthermore, phylogenetic bias ensures that over 50% of vascular genera, but barely 5% of non-vascular genera, are conserved ex situ. While botanic gardens are discernibly responding to the threat of species extinction, just 10% of network capacity is devoted to threatened species. We conclude that botanic gardens play a fundamental role in plant conservation, but identify actions to enhance future conservation of biodiversity.

Oscillation of a class of fractional differential equations with damping term.

Science.gov (United States)

Qin, Huizeng; Zheng, Bin

2013-01-01

We investigate the oscillation of a class of fractional differential equations with damping term. Based on a certain variable transformation, the fractional differential equations are converted into another differential equations of integer order with respect to the new variable. Then, using Riccati transformation, inequality, and integration average technique, some new oscillatory criteria for the equations are established. As for applications, oscillation for two certain fractional differential equations with damping term is investigated by the use of the presented results.

Optimal conservation outcomes require both restoration and protection.

Science.gov (United States)

Possingham, Hugh P; Bode, Michael; Klein, Carissa J

2015-01-01

Conservation outcomes are principally achieved through the protection of intact habitat or the restoration of degraded habitat. Restoration is generally considered a lower priority action than protection because protection is thought to provide superior outcomes, at lower costs, without the time delay required for restoration. Yet while it is broadly accepted that protected intact habitat safeguards more biodiversity and generates greater ecosystem services per unit area than restored habitat, conservation lacks a theory that can coherently compare the relative outcomes of the two actions. We use a dynamic landscape model to integrate these two actions into a unified conservation theory of protection and restoration. Using nonlinear benefit functions, we show that both actions are crucial components of a conservation strategy that seeks to optimise either biodiversity conservation or ecosystem services provision. In contrast to conservation orthodoxy, in some circumstances, restoration should be strongly preferred to protection. The relative priority of protection and restoration depends on their costs and also on the different time lags that are inherent to both protection and restoration. We derive a simple and easy-to-interpret heuristic that integrates these factors into a single equation that applies equally to biodiversity conservation and ecosystem service objectives. We use two examples to illustrate the theory: bird conservation in tropical rainforests and coastal defence provided by mangrove forests.

Optimal conservation outcomes require both restoration and protection.

Directory of Open Access Journals (Sweden)

Hugh P Possingham

2015-01-01

Full Text Available Conservation outcomes are principally achieved through the protection of intact habitat or the restoration of degraded habitat. Restoration is generally considered a lower priority action than protection because protection is thought to provide superior outcomes, at lower costs, without the time delay required for restoration. Yet while it is broadly accepted that protected intact habitat safeguards more biodiversity and generates greater ecosystem services per unit area than restored habitat, conservation lacks a theory that can coherently compare the relative outcomes of the two actions. We use a dynamic landscape model to integrate these two actions into a unified conservation theory of protection and restoration. Using nonlinear benefit functions, we show that both actions are crucial components of a conservation strategy that seeks to optimise either biodiversity conservation or ecosystem services provision. In contrast to conservation orthodoxy, in some circumstances, restoration should be strongly preferred to protection. The relative priority of protection and restoration depends on their costs and also on the different time lags that are inherent to both protection and restoration. We derive a simple and easy-to-interpret heuristic that integrates these factors into a single equation that applies equally to biodiversity conservation and ecosystem service objectives. We use two examples to illustrate the theory: bird conservation in tropical rainforests and coastal defence provided by mangrove forests.

A Hamiltonian functional for the linearized Einstein vacuum field equations

International Nuclear Information System (INIS)

Rosas-RodrIguez, R

2005-01-01

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

Peculiar symmetry structure of some known discrete nonautonomous equations

International Nuclear Information System (INIS)

Garifullin, R N; Habibullin, I T; Yamilov, R I

2015-01-01

We study the generalized symmetry structure of three known discrete nonautonomous equations. One of them is the semidiscrete dressing chain of Shabat. Two others are completely discrete equations defined on the square lattice. The first one is a discrete analogue of the dressing chain introduced by Levi and Yamilov. The second one is a nonautonomous generalization of the potential discrete KdV equation or, in other words, the H1 equation of the well-known Adler−Bobenko−Suris list. We demonstrate that these equations have generalized symmetries in both directions if and only if their coefficients, depending on the discrete variables, are periodic. The order of the simplest generalized symmetry in at least one direction depends on the period and may be arbitrarily high. We substantiate this picture by some theorems in the case of small periods. In case of an arbitrarily large period, we show that it is possible to construct two hierarchies of generalized symmetries and conservation laws. The same picture should take place in case of any nonautonomous equation of the Adler−Bobenko−Suris list. (paper)

The definition and computation of average neutron lifetimes

International Nuclear Information System (INIS)

Henry, A.F.

1983-01-01

A precise physical definition is offered for a class of average lifetimes for neutrons in an assembly of materials, either multiplying or not, or if the former, critical or not. A compact theoretical expression for the general member of this class is derived in terms of solutions to the transport equation. Three specific definitions are considered. Particular exact expressions for these are derived and reduced to simple algebraic formulas for one-group and two-group homogeneous bare-core models

Discrete exterior calculus discretization of incompressible Navier–Stokes equations over surface simplicial meshes

KAUST Repository

Mohamed, Mamdouh S.

2016-02-11

A conservative discretization of incompressible Navier–Stokes equations is developed based on discrete exterior calculus (DEC). A distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product. The governing equations are first rewritten using the exterior calculus notation, replacing vector calculus differential operators by the exterior derivative, Hodge star and wedge product operators. The discretization is then carried out by substituting with the corresponding discrete operators based on the DEC framework. Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured-triangular meshes, and first order accuracy for otherwise unstructured meshes. By construction, the method is conservative in that both mass and vorticity are conserved up to machine precision. The relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step.

Discrete exterior calculus discretization of incompressible Navier-Stokes equations over surface simplicial meshes

Science.gov (United States)

Mohamed, Mamdouh S.; Hirani, Anil N.; Samtaney, Ravi

2016-05-01

A conservative discretization of incompressible Navier-Stokes equations is developed based on discrete exterior calculus (DEC). A distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product. The governing equations are first rewritten using the exterior calculus notation, replacing vector calculus differential operators by the exterior derivative, Hodge star and wedge product operators. The discretization is then carried out by substituting with the corresponding discrete operators based on the DEC framework. Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured-triangular meshes, and first order accuracy for otherwise unstructured meshes. By construction, the method is conservative in that both mass and vorticity are conserved up to machine precision. The relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step.

Simulations of incompressible Navier Stokes equations on curved surfaces using discrete exterior calculus

Science.gov (United States)

Samtaney, Ravi; Mohamed, Mamdouh; Hirani, Anil

2015-11-01

We present examples of numerical solutions of incompressible flow on 2D curved domains. The Navier-Stokes equations are first rewritten using the exterior calculus notation, replacing vector calculus differential operators by the exterior derivative, Hodge star and wedge product operators. A conservative discretization of Navier-Stokes equations on simplicial meshes is developed based on discrete exterior calculus (DEC). The discretization is then carried out by substituting the corresponding discrete operators based on the DEC framework. By construction, the method is conservative in that both the discrete divergence and circulation are conserved up to machine precision. The relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step. Numerical examples include Taylor vortices on a sphere, Stuart vortices on a sphere, and flow past a cylinder on domains with varying curvature. Supported by the KAUST Office of Competitive Research Funds under Award No. URF/1/1401-01.

Frequency modulation at a moving material interface and a conservation law for wave number. [acoustic wave reflection and transmission

Science.gov (United States)

Kleinstein, G. G.; Gunzburger, M. D.

1976-01-01

An integral conservation law for wave numbers is considered. In order to test the validity of the proposed conservation law, a complete solution for the reflection and transmission of an acoustic wave impinging normally on a material interface moving at a constant speed is derived. The agreement between the frequency condition thus deduced from the dynamic equations of motion and the frequency condition derived from the jump condition associated with the integral equation supports the proposed law as a true conservation law. Additional comparisons such as amplitude discontinuities and Snells' law in a moving media further confirm the stated proposition. Results are stated concerning frequency and wave number relations across a shock front as predicted by the proposed conservation law.