Unification of the Two-Parameter Equation of State and the Principle of Corresponding States

DEFF Research Database (Denmark)

Mollerup, Jørgen

1998-01-01

A two-parameter equation of state is a two-parameter corresponding states model. A two-parameter corresponding states model is composed of two scale factor correlations and a reference fluid equation of state. In a two-parameter equation of state the reference equation of state is the two-paramet...

A family of conjugate gradient methods for large-scale nonlinear equations

Directory of Open Access Journals (Sweden)

Dexiang Feng

2017-09-01

Full Text Available Abstract In this paper, we present a family of conjugate gradient projection methods for solving large-scale nonlinear equations. At each iteration, it needs low storage and the subproblem can be easily solved. Compared with the existing solution methods for solving the problem, its global convergence is established without the restriction of the Lipschitz continuity on the underlying mapping. Preliminary numerical results are reported to show the efficiency of the proposed method.

Integral equation methods for electromagnetics

CERN Document Server

Volakis, John

2012-01-01

This text/reference is a detailed look at the development and use of integral equation methods for electromagnetic analysis, specifically for antennas and radar scattering. Developers and practitioners will appreciate the broad-based approach to understanding and utilizing integral equation methods and the unique coverage of historical developments that led to the current state-of-the-art. In contrast to existing books, Integral Equation Methods for Electromagnetics lays the groundwork in the initial chapters so students and basic users can solve simple problems and work their way up to the mo

Analysis of spectral methods for the homogeneous Boltzmann equation

KAUST Repository

Filbet, Francis; Mouhot, Clé ment

2011-01-01

The development of accurate and fast algorithms for the Boltzmann collision integral and their analysis represent a challenging problem in scientific computing and numerical analysis. Recently, several works were devoted to the derivation of spectrally accurate schemes for the Boltzmann equation, but very few of them were concerned with the stability analysis of the method. In particular there was no result of stability except when the method was modified in order to enforce the positivity preservation, which destroys the spectral accuracy. In this paper we propose a new method to study the stability of homogeneous Boltzmann equations perturbed by smoothed balanced operators which do not preserve positivity of the distribution. This method takes advantage of the "spreading" property of the collision, together with estimates on regularity and entropy production. As an application we prove stability and convergence of spectral methods for the Boltzmann equation, when the discretization parameter is large enough (with explicit bound). © 2010 American Mathematical Society.

Analysis of spectral methods for the homogeneous Boltzmann equation

KAUST Repository

Filbet, Francis

2011-04-01

The development of accurate and fast algorithms for the Boltzmann collision integral and their analysis represent a challenging problem in scientific computing and numerical analysis. Recently, several works were devoted to the derivation of spectrally accurate schemes for the Boltzmann equation, but very few of them were concerned with the stability analysis of the method. In particular there was no result of stability except when the method was modified in order to enforce the positivity preservation, which destroys the spectral accuracy. In this paper we propose a new method to study the stability of homogeneous Boltzmann equations perturbed by smoothed balanced operators which do not preserve positivity of the distribution. This method takes advantage of the "spreading" property of the collision, together with estimates on regularity and entropy production. As an application we prove stability and convergence of spectral methods for the Boltzmann equation, when the discretization parameter is large enough (with explicit bound). © 2010 American Mathematical Society.

Neutron Star masses from the Field Correlator Method Equation of State

Directory of Open Access Journals (Sweden)

Zappalà D.

2014-04-01

Full Text Available We analyse the hadron-quark phase transition in neutron stars by confronting the hadronic Equation of State (EoS obtained according to the microscopic Brueckner-Hartree-Fock many body theory, with the quark matter EoS derived within the Field Correlator Method. In particular, the latter EoS is only parametrized in terms of the gluon condensate and the large distance quark-antiquark potential, so that the comparison of the results of this analysis with the most recent measurements of heavy neutron star masses provides some physical constraints on these two parameters.

Analysis of a time fractional wave-like equation with the homotopy analysis method

International Nuclear Information System (INIS)

Xu Hang; Cang Jie

2008-01-01

The time fractional wave-like differential equation with a variable coefficient is studied analytically. By using a simple transformation, the governing equation is reduced to two fractional ordinary differential equations. Then the homotopy analysis method is employed to derive the solutions of these equations. The accurate series solutions are obtained. Especially, when h f =h g =-1, these solutions are exactly the same as those results given by the Adomian decomposition method. The present work shows the validity and great potential of the homotopy analysis method for solving nonlinear fractional differential equations. The basic idea described in this Letter is expected to be further employed to solve other similar nonlinear problems in fractional calculus

Some aspects of equations of state

International Nuclear Information System (INIS)

Frisch, H.L.

1979-02-01

Some elementary properties of the equation of state of molecules repulsing each other as point centers of force are developed briefly. An inequality for the Lennard--Jones gas is presented. The scaled particle theory equation of state of hard spheres is also reviewed briefly. Means of possibly applying these concepts to represent thermodynamic data on model detonating gases are suggested

Predictive equation of state method for heavy materials based on the Dirac equation and density functional theory

Science.gov (United States)

Wills, John M.; Mattsson, Ann E.

2012-02-01

Density functional theory (DFT) provides a formally predictive base for equation of state properties. Available approximations to the exchange/correlation functional provide accurate predictions for many materials in the periodic table. For heavy materials however, DFT calculations, using available functionals, fail to provide quantitative predictions, and often fail to be even qualitative. This deficiency is due both to the lack of the appropriate confinement physics in the exchange/correlation functional and to approximations used to evaluate the underlying equations. In order to assess and develop accurate functionals, it is essential to eliminate all other sources of error. In this talk we describe an efficient first-principles electronic structure method based on the Dirac equation and compare the results obtained with this method with other methods generally used. Implications for high-pressure equation of state of relativistic materials are demonstrated in application to Ce and the light actinides. Sandia National Laboratories is a multi-program laboratory managed andoperated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Analysis of a renormalization group method and normal form theory for perturbed ordinary differential equations

Science.gov (United States)

DeVille, R. E. Lee; Harkin, Anthony; Holzer, Matt; Josić, Krešimir; Kaper, Tasso J.

2008-06-01

For singular perturbation problems, the renormalization group (RG) method of Chen, Goldenfeld, and Oono [Phys. Rev. E. 49 (1994) 4502-4511] has been shown to be an effective general approach for deriving reduced or amplitude equations that govern the long time dynamics of the system. It has been applied to a variety of problems traditionally analyzed using disparate methods, including the method of multiple scales, boundary layer theory, the WKBJ method, the Poincaré-Lindstedt method, the method of averaging, and others. In this article, we show how the RG method may be used to generate normal forms for large classes of ordinary differential equations. First, we apply the RG method to systems with autonomous perturbations, and we show that the reduced or amplitude equations generated by the RG method are equivalent to the classical Poincaré-Birkhoff normal forms for these systems up to and including terms of O(ɛ2), where ɛ is the perturbation parameter. This analysis establishes our approach and generalizes to higher order. Second, we apply the RG method to systems with nonautonomous perturbations, and we show that the reduced or amplitude equations so generated constitute time-asymptotic normal forms, which are based on KBM averages. Moreover, for both classes of problems, we show that the main coordinate changes are equivalent, up to translations between the spaces in which they are defined. In this manner, our results show that the RG method offers a new approach for deriving normal forms for nonautonomous systems, and it offers advantages since one can typically more readily identify resonant terms from naive perturbation expansions than from the nonautonomous vector fields themselves. Finally, we establish how well the solution to the RG equations approximates the solution of the original equations on time scales of O(1/ɛ).

Solving the discrete KdV equation with homotopy analysis method

International Nuclear Information System (INIS)

Zou, L.; Zong, Z.; Wang, Z.; He, L.

2007-01-01

In this Letter, we apply the homotopy analysis method to differential-difference equations. We take the discrete KdV equation as an example, and successfully obtain double periodic wave solutions and solitary wave solutions. It illustrates the validity and the great potential of the homotopy analysis method in solving discrete KdV equation. Comparisons are made between the results of the proposed method and exact solutions. The results reveal that the proposed method is very effective and convenient

On the Effective Equation of State of Dark Energy

DEFF Research Database (Denmark)

Sloth, Martin Snoager

2010-01-01

In an effective field theory model with an ultraviolet momentum cutoff, there is a relation between the effective equation of state of dark energy and the ultraviolet cutoff scale. It implies that a measure of the equation of state of dark energy different from minus one, does not rule out vacuum...... energy as dark energy. It also indicates an interesting possibility that precise measurements of the infrared properties of dark energy can be used to probe the ultraviolet cutoff scale of effective quantum field theory coupled to gravity. In a toy model with a vacuum energy dominated universe...... with a Planck scale cutoff, the dark energy effective equation of state is -0.96....

Error Analysis of Galerkin's Method for Semilinear Equations

Directory of Open Access Journals (Sweden)

Tadashi Kawanago

2012-01-01

Full Text Available We establish a general existence result for Galerkin's approximate solutions of abstract semilinear equations and conduct an error analysis. Our results may be regarded as some extension of a precedent work (Schultz 1969. The derivation of our results is, however, different from the discussion in his paper and is essentially based on the convergence theorem of Newton’s method and some techniques for deriving it. Some of our results may be applicable for investigating the quality of numerical verification methods for solutions of ordinary and partial differential equations.

Comparative analysis of approximations used in the methods of Faddeev equations and hyperspherical harmonics

International Nuclear Information System (INIS)

Mukhtarova, M.I.

1988-01-01

Comparative analysis of approximations, used in the methods of Faddeev equations and hyperspherical harmonics (MHH) was conducted. The differences in solutions of these methods, related with introduction of approximation of sufficient partial states into the three-nucleon problem, is shown. MHH method is preferred. It is shown that MHH advantage can be manifested clearly when studying new classes of interactions: three-particle, Δ-isobar, nonlocal and other interactions

Multi-scale method for the resolution of the neutronic kinetics equations

International Nuclear Information System (INIS)

Chauvet, St.

2008-10-01

In this PhD thesis and in order to improve the time/precision ratio of the numerical simulation calculations, we investigate multi-scale techniques for the resolution of the reactor kinetics equations. We choose to focus on the mixed dual diffusion approximation and the quasi-static methods. We introduce a space dependency for the amplitude function which only depends on the time variable in the standard quasi-static context. With this new factorization, we develop two mixed dual problems which can be solved with Cea's solver MINOS. An algorithm is implemented, performing the resolution of these problems defined on different scales (for time and space). We name this approach: the Local Quasi-Static method. We present here this new multi-scale approach and its implementation. The inherent details of amplitude and shape treatments are discussed and justified. Results and performances, compared to MINOS, are studied. They illustrate the improvement on the time/precision ratio for kinetics calculations. Furthermore, we open some new possibilities to parallelize computations with MINOS. For the future, we also introduce some improvement tracks with adaptive scales. (author)

Stability theory for dynamic equations on time scales

CERN Document Server

Martynyuk, Anatoly A

2016-01-01

This monograph is a first in the world to present three approaches for stability analysis of solutions of dynamic equations. The first approach is based on the application of dynamic integral inequalities and the fundamental matrix of solutions of linear approximation of dynamic equations. The second is based on the generalization of the direct Lyapunovs method for equations on time scales, using scalar, vector and matrix-valued auxiliary functions. The third approach is the application of auxiliary functions (scalar, vector, or matrix-valued ones) in combination with differential dynamic inequalities. This is an alternative comparison method, developed for time continuous and time discrete systems. In recent decades, automatic control theory in the study of air- and spacecraft dynamics and in other areas of modern applied mathematics has encountered problems in the analysis of the behavior of solutions of time continuous-discrete linear and/or nonlinear equations of perturbed motion. In the book “Men of Ma...

Analysis of equations of state for ICF

International Nuclear Information System (INIS)

Ogando, F.; Velarde, P.

1998-01-01

In the field of inertial confinement fusion (ICF), numerical simulation plays a very important role reproducing real results with a good accuracy. Once considered the needed simplificative approximations, it is obtained and equation system that has to be solved in a more or less complex way. Anyway these simulation codes have to fit the behaviour of real materials with the help of some flexible parameters. In the case of fluid dynamics, the key to fit the equation system to the real material behaviour are the equations of state (EOS) involved in the problem. This is why there is an actual need of using accurate data in order to obtain results that match the real behaviour of materials. Nowadays there are several available sources of real EOS data. These sources are based either in theoretical models fixed with experimental results or in analytical approximations to those models. In this article a comparative analysis is made between two of the most used EOS data sources in simulations. The aim of this work is to check the quality of these data that is commonly used and on which the accuracy of the results depend. (Author) 3 refs

Generating a Multiphase Equation of State with Swarm Intelligence

Science.gov (United States)

Cox, Geoffrey

2017-06-01

Hydrocode calculations require knowledge of the variation of pressure of a material with density and temperature, which is given by the equation of state. An accurate model needs to account for discontinuities in energy, density and properties of a material across a phase boundary. When generating a multiphase equation of state the modeller attempts to balance the agreement between the available data for compression, expansion and phase boundary location. However, this can prove difficult because minor adjustments in the equation of state for a single phase can have a large impact on the overall phase diagram. Recently, Cox and Christie described a method for combining statistical-mechanics-based condensed matter physics models with a stochastic analysis technique called particle swarm optimisation. The models produced show good agreement with experiment over a wide range of pressure-temperature space. This talk details the general implementation of this technique, shows example results, and describes the types of analysis that can be performed with this method.

A relaxation-projection method for compressible flows. Part I: The numerical equation of state for the Euler equations

International Nuclear Information System (INIS)

Saurel, Richard; Franquet, Erwin; Daniel, Eric; Le Metayer, Olivier

2007-01-01

A new projection method is developed for the Euler equations to determine the thermodynamic state in computational cells. It consists in the resolution of a mechanical relaxation problem between the various sub-volumes present in a computational cell. These sub-volumes correspond to the ones traveled by the various waves that produce states with different pressures, velocities, densities and temperatures. Contrarily to Godunov type schemes the relaxed state corresponds to mechanical equilibrium only and remains out of thermal equilibrium. The pressure computation with this relaxation process replaces the use of the conventional equation of state (EOS). A simplified relaxation method is also derived and provides a specific EOS (named the Numerical EOS). The use of the Numerical EOS gives a cure to spurious pressure oscillations that appear at contact discontinuities for fluids governed by real gas EOS. It is then extended to the computation of interface problems separating fluids with different EOS (liquid-gas interface for example) with the Euler equations. The resulting method is very robust, accurate, oscillation free and conservative. For the sake of simplicity and efficiency the method is developed in a Lagrange-projection context and is validated over exact solutions. In a companion paper [F. Petitpas, E. Franquet, R. Saurel, A relaxation-projection method for compressible flows. Part II: computation of interfaces and multiphase mixtures with stiff mechanical relaxation. J. Comput. Phys. (submitted for publication)], the method is extended to the numerical approximation of a non-conservative hyperbolic multiphase flow model for interface computation and shock propagation into mixtures

Novel algorithm of large-scale simultaneous linear equations

International Nuclear Information System (INIS)

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

2010-01-01

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

Methods of mathematical modelling continuous systems and differential equations

CERN Document Server

Witelski, Thomas

2015-01-01

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

Steady-state transport equation resolution by particle methods, and numerical results

International Nuclear Information System (INIS)

Mercier, B.

1985-10-01

A method to solve steady-state transport equation has been given. Principles of the method are given. The method is studied in two different cases; estimations given by the theory are compared to numerical results. Results got in 1-D (spherical geometry) and in 2-D (axisymmetric geometry) are given [fr

A multi scale approximation solution for the time dependent Boltzmann-transport equation

International Nuclear Information System (INIS)

Merk, B.

2004-03-01

The basis of all transient simulations for nuclear reactor cores is the reliable calculation of the power production. The local power distribution is generally calculated by solving the space, time, energy and angle dependent neutron transport equation known as Boltzmann equation. The computation of exact solutions of the Boltzmann equation is very time consuming. For practical numerical simulations approximated solutions are usually unavoidable. The objective of this work is development of an effective multi scale approximation solution for the Boltzmann equation. Most of the existing methods are based on separation of space and time. The new suggested method is performed without space-time separation. This effective approximation solution is developed on the basis of an expansion for the time derivative of different approximations to the Boltzmann equation. The method of multiple scale expansion is used for the expansion of the time derivative, because the problem of the stiff time behaviour can't be expressed by standard expansion methods. This multiple scale expansion is used in this work to develop approximation solutions for different approximations of the Boltzmann equation, starting from the expansion of the point kinetics equations. The resulting analytic functions are used for testing the applicability and accuracy of the multiple scale expansion method for an approximation solution with 2 delayed neutron groups. The results are tested versus the exact analytical results for the point kinetics equations. Very good agreement between both solutions is obtained. The validity of the solution with 2 delayed neutron groups to approximate the behaviour of the system with 6 delayed neutron groups is demonstrated in an additional analysis. A strategy for a solution with 4 delayed neutron groups is described. A multiple scale expansion is performed for the space-time dependent diffusion equation for one homogenized cell with 2 delayed neutron groups. The result is

The quasi-steady state of all-vanadium redox flow batteries: A scale analysis

International Nuclear Information System (INIS)

Sharma, A.K.; Vynnycky, M.; Ling, C.Y.; Birgersson, E.; Han, M.

2014-01-01

Highlights: • We present a transient 2D model for a VRFB (conservation of species and charge); • Carry out scale analysis of the species conservation equation; • Derive the condition characterizing the quasi-steadiness of VRFB operation; • Verify it by comparing charge-discharge curve with transient simulations. - Abstract: In general, mathematical models for all-vanadium redox flow batteries (VRFB) that seek to capture the transport phenomena are transient in nature. In this paper, we carry out scale analysis of VRFB operation and derive the conditions when it can be assumed to be quasi-steady state in nature, i.e., time-dependence only through a boundary condition. We find that it is true for typical tank volume and flow rate employed for VRFBs. The proposed analysis is generic and can also be employed for other types of redox flow batteries

Spectral representations of neutron-star equations of state

International Nuclear Information System (INIS)

Lindblom, Lee

2010-01-01

Methods are developed for constructing spectral representations of cold (barotropic) neutron-star equations of state. These representations are faithful in the sense that every physical equation of state has a representation of this type and conversely every such representation satisfies the minimal thermodynamic stability criteria required of any physical equation of state. These spectral representations are also efficient, in the sense that only a few spectral coefficients are generally required to represent neutron-star equations of state quiet accurately. This accuracy and efficiency is illustrated by constructing spectral fits to a large collection of 'realistic' neutron-star equations of state.

Scattering theory methods for bound state problems

International Nuclear Information System (INIS)

Raphael, R.B.; Tobocman, W.

1978-01-01

For the analysis of the properties of a bound state system one may use in place of the Schroedinger equation the Lippmann-Schwinger (LS) equation for the wave function or the LS equation for the reactance operator. Use of the LS equation for the reactance operator constrains the solution to have correct asymptotic behaviour, so this approach would appear to be desirable when the bound state wave function is to be used to calculate particle transfer form factors. The Schroedinger equation based N-level analysis of the s-wave bound states of a square well is compared to the ones based on the LS equation. It is found that the LS equation methods work better than the Schroedinger equation method. The method that uses the LS equation for the wave function gives the best results for the wave functions while the method that uses the LS equation for the reactance operator gives the best results for the binding energies. The accuracy of the reactance operator based method is remarkably insensitive to changes in the oscillator constant used for the harmonic oscillator function basis set. It is also remarkably insensitive to the number of nodes in the bound state wave function. (Auth.)

Solution of the Schroedinger equation by a spectral method

International Nuclear Information System (INIS)

Feit, M.D.; Fleck, J.A. Jr.; Steiger, A.

1982-01-01

A new computational method for determining the eigenvalues and eigenfunctions of the Schroedinger equation is described. Conventional methods for solving this problem rely on diagonalization of a Hamiltonian matrix or iterative numerical solutions of a time independent wave equation. The new method, in contrast, is based on the spectral properties of solutions to the time-dependent Schroedinger equation. The method requires the computation of a correlation function from a numerical solution psi(r, t). Fourier analysis of this correlation function reveals a set of resonant peaks that correspond to the stationary states of the system. Analysis of the location of these peaks reveals the eigenvalues with high accuracy. Additional Fourier transforms of psi(r, t) with respect to time generate the eigenfunctions. The effectiveness of the method is demonstrated for a one-dimensional asymmetric double well potential and for the two-dimensional Henon--Heiles potential

A gradual update method for simulating the steady-state solution of stiff differential equations in metabolic circuits.

Science.gov (United States)

Shiraishi, Emi; Maeda, Kazuhiro; Kurata, Hiroyuki

2009-02-01

Numerical simulation of differential equation systems plays a major role in the understanding of how metabolic network models generate particular cellular functions. On the other hand, the classical and technical problems for stiff differential equations still remain to be solved, while many elegant algorithms have been presented. To relax the stiffness problem, we propose new practical methods: the gradual update of differential-algebraic equations based on gradual application of the steady-state approximation to stiff differential equations, and the gradual update of the initial values in differential-algebraic equations. These empirical methods show a high efficiency for simulating the steady-state solutions for the stiff differential equations that existing solvers alone cannot solve. They are effective in extending the applicability of dynamic simulation to biochemical network models.

Stochastic Galerkin methods for the steady-state Navier–Stokes equations

Energy Technology Data Exchange (ETDEWEB)

Sousedík, Bedřich, E-mail: sousedik@umbc.edu [Department of Mathematics and Statistics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250 (United States); Elman, Howard C., E-mail: elman@cs.umd.edu [Department of Computer Science and Institute for Advanced Computer Studies, University of Maryland, College Park, MD 20742 (United States)

2016-07-01

We study the steady-state Navier–Stokes equations in the context of stochastic finite element discretizations. Specifically, we assume that the viscosity is a random field given in the form of a generalized polynomial chaos expansion. For the resulting stochastic problem, we formulate the model and linearization schemes using Picard and Newton iterations in the framework of the stochastic Galerkin method, and we explore properties of the resulting stochastic solutions. We also propose a preconditioner for solving the linear systems of equations arising at each step of the stochastic (Galerkin) nonlinear iteration and demonstrate its effectiveness for solving a set of benchmark problems.

Mathematical and computational methods for semiclassical Schrödinger equations

KAUST Repository

Jin, Shi

2011-04-28

We consider time-dependent (linear and nonlinear) Schrödinger equations in a semiclassical scaling. These equations form a canonical class of (nonlinear) dispersive models whose solutions exhibit high-frequency oscillations. The design of efficient numerical methods which produce an accurate approximation of the solutions, or at least of the associated physical observables, is a formidable mathematical challenge. In this article we shall review the basic analytical methods for dealing with such equations, including WKB asymptotics, Wigner measure techniques and Gaussian beams. Moreover, we shall give an overview of the current state of the art of numerical methods (most of which are based on the described analytical techniques) for the Schrödinger equation in the semiclassical regime. © 2011 Cambridge University Press.

Entropy viscosity method applied to Euler equations

International Nuclear Information System (INIS)

Delchini, M. O.; Ragusa, J. C.; Berry, R. A.

2013-01-01

The entropy viscosity method [4] has been successfully applied to hyperbolic systems of equations such as Burgers equation and Euler equations. The method consists in adding dissipative terms to the governing equations, where a viscosity coefficient modulates the amount of dissipation. The entropy viscosity method has been applied to the 1-D Euler equations with variable area using a continuous finite element discretization in the MOOSE framework and our results show that it has the ability to efficiently smooth out oscillations and accurately resolve shocks. Two equations of state are considered: Ideal Gas and Stiffened Gas Equations Of State. Results are provided for a second-order time implicit schemes (BDF2). Some typical Riemann problems are run with the entropy viscosity method to demonstrate some of its features. Then, a 1-D convergent-divergent nozzle is considered with open boundary conditions. The correct steady-state is reached for the liquid and gas phases with a time implicit scheme. The entropy viscosity method correctly behaves in every problem run. For each test problem, results are shown for both equations of state considered here. (authors)

Advances in iterative methods for nonlinear equations

CERN Document Server

Busquier, Sonia

2016-01-01

This book focuses on the approximation of nonlinear equations using iterative methods. Nine contributions are presented on the construction and analysis of these methods, the coverage encompassing convergence, efficiency, robustness, dynamics, and applications. Many problems are stated in the form of nonlinear equations, using mathematical modeling. In particular, a wide range of problems in Applied Mathematics and in Engineering can be solved by finding the solutions to these equations. The book reveals the importance of studying convergence aspects in iterative methods and shows that selection of the most efficient and robust iterative method for a given problem is crucial to guaranteeing a good approximation. A number of sample criteria for selecting the optimal method are presented, including those regarding the order of convergence, the computational cost, and the stability, including the dynamics. This book will appeal to researchers whose field of interest is related to nonlinear problems and equations...

Stochastic fractional differential equations: Modeling, method and analysis

International Nuclear Information System (INIS)

Pedjeu, Jean-C.; Ladde, Gangaram S.

2012-01-01

By introducing a concept of dynamic process operating under multi-time scales in sciences and engineering, a mathematical model described by a system of multi-time scale stochastic differential equations is formulated. The classical Picard–Lindelöf successive approximations scheme is applied to the model validation problem, namely, existence and uniqueness of solution process. Naturally, this leads to the problem of finding closed form solutions of both linear and nonlinear multi-time scale stochastic differential equations of Itô–Doob type. Finally, to illustrate the scope of ideas and presented results, multi-time scale stochastic models for ecological and epidemiological processes in population dynamic are outlined.

Comparative analysis of solution methods of the punctual kinetic equations

International Nuclear Information System (INIS)

Hernandez S, A.

2003-01-01

The following one written it presents a comparative analysis among different analytical solutions for the punctual kinetics equation, which present two variables of interest: a) the temporary behavior of the neutronic population, and b) The temporary behavior of the different groups of precursors of delayed neutrons. The first solution is based on a method that solves the transfer function of the differential equation for the neutronic population, in which intends to obtain the different poles that give the stability of this transfer function. In this section it is demonstrated that the temporary variation of the reactivity of the system can be managed as it is required, since the integration time for this method doesn't affect the result. However, the second solution is based on an iterative method like that of Runge-Kutta or the Euler method where the algorithm was only used to solve first order differential equations giving this way solution to each differential equation that conforms the equations of punctual kinetics. In this section it is demonstrated that only it can obtain a correct temporary behavior of the neutronic population when it is integrated on an interval of very short time, forcing to the temporary variation of the reactivity to change very quick way without one has some control about the time. In both methods the same change is used so much in the reactivity of the system like in the integration times, giving validity to the results graph the one the temporary behavior of the neutronic population vs. time. (Author)

Excited-state lifetime measurements: Linearization of the Foerster equation by the phase-plane method

International Nuclear Information System (INIS)

Love, J.C.; Demas, J.N.

1983-01-01

The Foerster equation describes excited-state decay curves involving resonance intermolecular energy transfer. A linearized solution based on the phase-plane method has been developed. The new method is quick, insensitive to the fitting region, accurate, and precise

A methodology for uncertainty analysis of reference equations of state

DEFF Research Database (Denmark)

Cheung, Howard; Frutiger, Jerome; Bell, Ian H.

We present a detailed methodology for the uncertainty analysis of reference equations of state (EOS) based on Helmholtz energy. In recent years there has been an increased interest in uncertainties of property data and process models of thermal systems. In the literature there are various...... for uncertainty analysis is suggested as a tool for EOS. The uncertainties of the EOS properties are calculated from the experimental values and the EOS model structure through the parameter covariance matrix and subsequent linear error propagation. This allows reporting the uncertainty range (95% confidence...

Dynamic analysis of isotropic nanoplates subjected to moving load using state-space method based on nonlocal second order plate theory

Energy Technology Data Exchange (ETDEWEB)

Nami, Mohammad Rahim [Shiraz University, Shiraz, Iran (Iran, Islamic Republic of); Janghorban, Maziar [Islamic Azad University, Marvdash (Iran, Islamic Republic of)

2015-06-15

In this work, dynamic analysis of rectangular nanoplates subjected to moving load is presented. In order to derive the governing equations of motion, second order plate theory is used. To capture the small scale effects, the nonlocal elasticity theory is adopted. It is assumed that the nanoplate is subjected to a moving concentrated load with the constant velocity V in the x direction. To solve the governing equations, state-space method is used to find the deflections of rectangular nanoplate under moving load. The results obtained here reveal that the nonlocality has significant effect on the deflection of rectangular nanoplate subjected to moving load.

Sensitivity of rocky planet structures to the equation of state

International Nuclear Information System (INIS)

Swift, D.C.

2009-01-01

Structures were calculated for Mercury, Venus, Earth, the Moon, and Mars, using a core-mantle model and adjusting the core radius to reproduce the observed mass and diameter of each body. Structures were calculated using Fe and basalt equations of state of different degrees of sophistication for the core and mantle. The choice of equation of state had a significant effect on the inferred structure. For each structure, the moment of inertia ratio was calculated and compared with observed values. Linear Grueneisen equations of state fitted to limited portions of shock data reproduced the observed moments of inertia significantly better than did more detailed equations of state incorporating phase transitions, presumably reflecting the actual compositions of the bodies. The linear Grueneisen equations of state and corresponding structures seem however to be a reasonable starting point for comparative simulations of large-scale astrophysical impacts.

Approximate Solutions of Nonlinear Partial Differential Equations by Modified q-Homotopy Analysis Method

Directory of Open Access Journals (Sweden)

Shaheed N. Huseen

2013-01-01

Full Text Available A modified q-homotopy analysis method (mq-HAM was proposed for solving nth-order nonlinear differential equations. This method improves the convergence of the series solution in the nHAM which was proposed in (see Hassan and El-Tawil 2011, 2012. The proposed method provides an approximate solution by rewriting the nth-order nonlinear differential equation in the form of n first-order differential equations. The solution of these n differential equations is obtained as a power series solution. This scheme is tested on two nonlinear exactly solvable differential equations. The results demonstrate the reliability and efficiency of the algorithm developed.

BADGER v1.0: A Fortran equation of state library

Science.gov (United States)

Heltemes, T. A.; Moses, G. A.

2012-12-01

The BADGER equation of state library was developed to enable inertial confinement fusion plasma codes to more accurately model plasmas in the high-density, low-temperature regime. The code had the capability to calculate 1- and 2-T plasmas using the Thomas-Fermi model and an individual electron accounting model. Ion equation of state data can be calculated using an ideal gas model or via a quotidian equation of state with scaled binding energies. Electron equation of state data can be calculated via the ideal gas model or with an adaptation of the screened hydrogenic model with ℓ-splitting. The ionization and equation of state calculations can be done in local thermodynamic equilibrium or in a non-LTE mode using a variant of the Busquet equivalent temperature method. The code was written as a stand-alone Fortran library for ease of implementation by external codes. EOS results for aluminum are presented that show good agreement with the SESAME library and ionization calculations show good agreement with the FLYCHK code. Program summaryProgram title: BADGERLIB v1.0 Catalogue identifier: AEND_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEND_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 41 480 No. of bytes in distributed program, including test data, etc.: 2 904 451 Distribution format: tar.gz Programming language: Fortran 90. Computer: 32- or 64-bit PC, or Mac. Operating system: Windows, Linux, MacOS X. RAM: 249.496 kB plus 195.630 kB per isotope record in memory Classification: 19.1, 19.7. Nature of problem: Equation of State (EOS) calculations are necessary for the accurate simulation of high energy density plasmas. Historically, most EOS codes used in these simulations have relied on an ideal gas model. This model is inadequate for low

An analytic equation of state for Ising-like models

International Nuclear Information System (INIS)

O'Connor, Denjoe; Santiago, J A; Stephens, C R

2007-01-01

Using an environmentally friendly renormalization we derive, from an underlying field theory representation, a formal expression for the equation of state, y = f(x), that exhibits all desired asymptotic and analyticity properties in the three limits x → 0, x → ∞ and x → -1. The only necessary inputs are the Wilson functions γ λ , γ ψ and γ φ 2 , associated with a renormalization of the transverse vertex functions. These Wilson functions exhibit a crossover between the Wilson-Fisher fixed point and the fixed point that controls the coexistence curve. Restricting to the case N = 1, we derive a one-loop equation of state for 2 < d < 4 naturally parameterized by a ratio of nonlinear scaling fields. For d = 3 we show that a non-parameterized analytic form can be deduced. Various asymptotic amplitudes are calculated directly from the equation of state in all three asymptotic limits of interest and comparison made with known results. By positing a scaling form for the equation of state inspired by the one-loop result, but adjusted to fit the known values of the critical exponents, we obtain better agreement with known asymptotic amplitudes

Mathematical Analysis of Vehicle Delivery Scale of Bike-Sharing Rental Nodes

Science.gov (United States)

Zhai, Y.; Liu, J.; Liu, L.

2018-04-01

Aiming at the lack of scientific and reasonable judgment of vehicles delivery scale and insufficient optimization of scheduling decision, based on features of the bike-sharing usage, this paper analyses the applicability of the discrete time and state of the Markov chain, and proves its properties to be irreducible, aperiodic and positive recurrent. Based on above analysis, the paper has reached to the conclusion that limit state (steady state) probability of the bike-sharing Markov chain only exists and is independent of the initial probability distribution. Then this paper analyses the difficulty of the transition probability matrix parameter statistics and the linear equations group solution in the traditional solving algorithm of the bike-sharing Markov chain. In order to improve the feasibility, this paper proposes a "virtual two-node vehicle scale solution" algorithm which considered the all the nodes beside the node to be solved as a virtual node, offered the transition probability matrix, steady state linear equations group and the computational methods related to the steady state scale, steady state arrival time and scheduling decision of the node to be solved. Finally, the paper evaluates the rationality and accuracy of the steady state probability of the proposed algorithm by comparing with the traditional algorithm. By solving the steady state scale of the nodes one by one, the proposed algorithm is proved to have strong feasibility because it lowers the level of computational difficulty and reduces the number of statistic, which will help the bike-sharing companies to optimize the scale and scheduling of nodes.

MATHEMATICAL ANALYSIS OF VEHICLE DELIVERY SCALE OF BIKE-SHARING RENTAL NODES

Directory of Open Access Journals (Sweden)

Y. Zhai

2018-04-01

Full Text Available Aiming at the lack of scientific and reasonable judgment of vehicles delivery scale and insufficient optimization of scheduling decision, based on features of the bike-sharing usage, this paper analyses the applicability of the discrete time and state of the Markov chain, and proves its properties to be irreducible, aperiodic and positive recurrent. Based on above analysis, the paper has reached to the conclusion that limit state (steady state probability of the bike-sharing Markov chain only exists and is independent of the initial probability distribution. Then this paper analyses the difficulty of the transition probability matrix parameter statistics and the linear equations group solution in the traditional solving algorithm of the bike-sharing Markov chain. In order to improve the feasibility, this paper proposes a "virtual two-node vehicle scale solution" algorithm which considered the all the nodes beside the node to be solved as a virtual node, offered the transition probability matrix, steady state linear equations group and the computational methods related to the steady state scale, steady state arrival time and scheduling decision of the node to be solved. Finally, the paper evaluates the rationality and accuracy of the steady state probability of the proposed algorithm by comparing with the traditional algorithm. By solving the steady state scale of the nodes one by one, the proposed algorithm is proved to have strong feasibility because it lowers the level of computational difficulty and reduces the number of statistic, which will help the bike-sharing companies to optimize the scale and scheduling of nodes.

Assessment of UF6 Equation of State

Energy Technology Data Exchange (ETDEWEB)

Brady, P; Chand, K; Warren, D; Vandersall, J

2009-02-11

A common assumption in the mathematical analysis of flows of compressible fluids is to treat the fluid as a perfect gas. This is an approximation, as no real fluid obeys the perfect gas relationships over all temperature and pressure conditions. An assessment of the validity of treating the UF{sub 6} gas flow field within a gas centrifuge with perfect gas relationships has been conducted. The definition of a perfect gas is commonly stated in two parts: (1) the gas obeys the thermal equation of state, p = {rho}RT (thermally perfect), and, (2) the gas specific heats are constant (calorically perfect). Analysis indicates the thermally perfect assumption is valid for all flow conditions within the gas centrifuge, including shock fields. The low operating gas pressure is the primary factor in the suitability of the thermally perfect equation of state for gas centrifuge computations. UF{sub 6} is not calorically perfect, as the specific heats vary as a function of temperature. This effect is insignificant within the bulk of the centrifuge gas field, as gas temperatures vary over a narrow range. The exception is in the vicinity of shock fields, where temperature, pressure, and density gradients are large, and the variation of specific heats with temperature should be included in the technically detailed analyses. Results from a normal shock analysis incorporating variable specific heats is included herein, presented in the conventional form of shock parameters as a function of inlet Mach Number. The error introduced by assuming constant specific heats is small for a nominal UF{sub 6} shock field, such that calorically perfect shock relationships can be used for scaling and initial analyses. The more rigorous imperfect gas analysis should be used for detailed analyses.

A asymptotic numerical method for the steady-state convection diffusion equation

International Nuclear Information System (INIS)

Wu Qiguang

1988-01-01

In this paper, A asymptotic numerical method for the steady-state Convection diffusion equation is proposed, which need not take very fine mesh size in the neighbourhood of the boundary layer. Numerical computation for model problem show that we can obtain the numerical solution in the boundary layer with moderate step size

An original state equation for argillite

International Nuclear Information System (INIS)

Cariou, S.; Dormieux, L.; Skoczylas, F.

2010-01-01

Document available in extended abstract form only. Argillite is available in large quantities in France and exhibits a very small permeability: that is the reason why it is an excellent candidate for the storage of nuclear waste depositories. This study focuses particularly on argillite from Meuse-Haute-Marne (East of France) located 500 m deep. Following the observation that Biot theory does not permit us to explain some of our experiments, we propose in this contribution to use the tools of continuum micro-mechanics in order to develop a macroscopic state equation that would be consistent with our experimental data. In order to apply homogenization tools to argillite, the first step is to describe the microstructure of this material. Argillite can be modelled as a multi-scale material: - At the macro-scale, argillite looks like a homogenized material. - At the meso-scale, a shale matrix surrounding some inclusions can be distinguished. - At the micro-scale, the shale matrix reveals to be made of a solid phase and of a partially saturated porous network; capillary effects are taken into account at this scale. - At the nano-scale, the solid phase is not homogeneous but is a set of particles that are parallel platelets with an electrolyte in-between; the solid phase is not physically inert with regard to water in the porous network; first, liquid pressure in the electrolyte exists and is equal to liquid pressure in the porous network; then, the unbalance between the ion concentration in-between the platelets building up the solid phase and the ion concentration in the porous network implies an overpressure in the electrolyte. This complex model of argillite leads us though several homogenization steps to an original macroscopic state equation. We propose to compare this state equation with experimental data. Results of a poro-mechanical test run on a partially saturated sample of argillite are shown. One can read that the deformations of the sample are identical

Numerical methods and analysis of the nonlinear Vlasov equation on unstructured meshes of phase space

International Nuclear Information System (INIS)

Besse, Nicolas

2003-01-01

This work is dedicated to the mathematical and numerical studies of the Vlasov equation on phase-space unstructured meshes. In the first part, new semi-Lagrangian methods are developed to solve the Vlasov equation on unstructured meshes of phase space. As the Vlasov equation describes multi-scale phenomena, we also propose original methods based on a wavelet multi-resolution analysis. The resulting algorithm leads to an adaptive mesh-refinement strategy. The new massively-parallel computers allow to use these methods with several phase-space dimensions. Particularly, these numerical schemes are applied to plasma physics and charged particle beams in the case of two-, three-, and four-dimensional Vlasov-Poisson systems. In the second part we prove the convergence and give error estimates for several numerical schemes applied to the Vlasov-Poisson system when strong and classical solutions are considered. First we show the convergence of a semi-Lagrangian scheme on an unstructured mesh of phase space, when the regularity hypotheses for the initial data are minimal. Then we demonstrate the convergence of classes of high-order semi-Lagrangian schemes in the framework of the regular classical solution. In order to reconstruct the distribution function, we consider symmetrical Lagrange polynomials, B-Splines and wavelets bases. Finally we prove the convergence of a semi-Lagrangian scheme with propagation of gradients yielding a high-order and stable reconstruction of the solution. (author) [fr

Solution of two group neutron diffusion equation by using homotopy analysis method

International Nuclear Information System (INIS)

Cavdar, S.

2010-01-01

The Homotopy Analysis Method (HAM), proposed in 1992 by Shi Jun Liao and has been developed since then, is based on differential geometry as well as homotopy which is a fundamental concept in topology. It has proved to be useful for obtaining series solutions of many such problems involving algebraic, linear/non-linear, ordinary/partial differential equations, differential-integral equations, differential-difference equations, and coupled equations of them. Briefly, through HAM, it is possible to construct a continuous mapping of an initial guess approximation to the exact solution of the equation of concern. An auxiliary linear operator is chosen to construct such kind of a continuous mapping and an auxiliary parameter is used to ensure the convergence of series solution. We present the solutions of two-group neutron diffusion equation through HAM in this work. We also compare the results with that obtained by other well-known solution analytical and numeric methods.

Tailored parameter optimization methods for ordinary differential equation models with steady-state constraints.

Science.gov (United States)

Fiedler, Anna; Raeth, Sebastian; Theis, Fabian J; Hausser, Angelika; Hasenauer, Jan

2016-08-22

Ordinary differential equation (ODE) models are widely used to describe (bio-)chemical and biological processes. To enhance the predictive power of these models, their unknown parameters are estimated from experimental data. These experimental data are mostly collected in perturbation experiments, in which the processes are pushed out of steady state by applying a stimulus. The information that the initial condition is a steady state of the unperturbed process provides valuable information, as it restricts the dynamics of the process and thereby the parameters. However, implementing steady-state constraints in the optimization often results in convergence problems. In this manuscript, we propose two new methods for solving optimization problems with steady-state constraints. The first method exploits ideas from optimization algorithms on manifolds and introduces a retraction operator, essentially reducing the dimension of the optimization problem. The second method is based on the continuous analogue of the optimization problem. This continuous analogue is an ODE whose equilibrium points are the optima of the constrained optimization problem. This equivalence enables the use of adaptive numerical methods for solving optimization problems with steady-state constraints. Both methods are tailored to the problem structure and exploit the local geometry of the steady-state manifold and its stability properties. A parameterization of the steady-state manifold is not required. The efficiency and reliability of the proposed methods is evaluated using one toy example and two applications. The first application example uses published data while the second uses a novel dataset for Raf/MEK/ERK signaling. The proposed methods demonstrated better convergence properties than state-of-the-art methods employed in systems and computational biology. Furthermore, the average computation time per converged start is significantly lower. In addition to the theoretical results, the

Modeling Pore-Scale Oil-Gas Systems Using Gradient Theory with Peng-Robinson Equation of State

KAUST Repository

Fan, Xiaolin

2016-06-01

This research addresses a sequential convex splitting method for numerical simulation of multicomponent two-phase fluids mixture in a single-pore at constant temperature, which is modeled by the gradient theory with Peng-Robinson equation of state. The gradient theory of thermodynamics and variational calculus are utilized to obtain a system of chemical equilibrium equations which are transformed into a transient system as a numerical strategy on which the numerical scheme is based. The proposed numerical algorithm avoids computing Hessian matrix arising from the second-order derivative of homogeneous contribution of free energy; it is also quite robust. This scheme is proved to be unconditionally component-wise energy stable. The Raviart-Thomas mixed finite element method is applied to spatial discretization.

Modeling Pore-Scale Oil-Gas Systems Using Gradient Theory with Peng-Robinson Equation of State

KAUST Repository

Fan, Xiaolin; Kou, Jisheng; Qiao, Zhonghua; Sun, Shuyu

2016-01-01

This research addresses a sequential convex splitting method for numerical simulation of multicomponent two-phase fluids mixture in a single-pore at constant temperature, which is modeled by the gradient theory with Peng-Robinson equation of state. The gradient theory of thermodynamics and variational calculus are utilized to obtain a system of chemical equilibrium equations which are transformed into a transient system as a numerical strategy on which the numerical scheme is based. The proposed numerical algorithm avoids computing Hessian matrix arising from the second-order derivative of homogeneous contribution of free energy; it is also quite robust. This scheme is proved to be unconditionally component-wise energy stable. The Raviart-Thomas mixed finite element method is applied to spatial discretization.

Fractal-Based Methods and Inverse Problems for Differential Equations: Current State of the Art

Directory of Open Access Journals (Sweden)

Herb E. Kunze

2014-01-01

Full Text Available We illustrate, in this short survey, the current state of the art of fractal-based techniques and their application to the solution of inverse problems for ordinary and partial differential equations. We review several methods based on the Collage Theorem and its extensions. We also discuss two innovative applications: the first one is related to a vibrating string model while the second one considers a collage-based approach for solving inverse problems for partial differential equations on a perforated domain.

Multi-scale approximation of Vlasov equation

International Nuclear Information System (INIS)

Mouton, A.

2009-09-01

One of the most important difficulties of numerical simulation of magnetized plasmas is the existence of multiple time and space scales, which can be very different. In order to produce good simulations of these multi-scale phenomena, it is recommended to develop some models and numerical methods which are adapted to these problems. Nowadays, the two-scale convergence theory introduced by G. Nguetseng and G. Allaire is one of the tools which can be used to rigorously derive multi-scale limits and to obtain new limit models which can be discretized with a usual numerical method: this procedure is so-called a two-scale numerical method. The purpose of this thesis is to develop a two-scale semi-Lagrangian method and to apply it on a gyrokinetic Vlasov-like model in order to simulate a plasma submitted to a large external magnetic field. However, the physical phenomena we have to simulate are quite complex and there are many questions without answers about the behaviour of a two-scale numerical method, especially when such a method is applied on a nonlinear model. In a first part, we develop a two-scale finite volume method and we apply it on the weakly compressible 1D isentropic Euler equations. Even if this mathematical context is far from a Vlasov-like model, it is a relatively simple framework in order to study the behaviour of a two-scale numerical method in front of a nonlinear model. In a second part, we develop a two-scale semi-Lagrangian method for the two-scale model developed by E. Frenod, F. Salvarani et E. Sonnendrucker in order to simulate axisymmetric charged particle beams. Even if the studied physical phenomena are quite different from magnetic fusion experiments, the mathematical context of the one-dimensional paraxial Vlasov-Poisson model is very simple for establishing the basis of a two-scale semi-Lagrangian method. In a third part, we use the two-scale convergence theory in order to improve M. Bostan's weak-* convergence results about the finite

A method for the selection of a functional form for a thermodynamic equation of state using weighted linear least squares stepwise regression

Science.gov (United States)

Jacobsen, R. T.; Stewart, R. B.; Crain, R. W., Jr.; Rose, G. L.; Myers, A. F.

1976-01-01

A method was developed for establishing a rational choice of the terms to be included in an equation of state with a large number of adjustable coefficients. The methods presented were developed for use in the determination of an equation of state for oxygen and nitrogen. However, a general application of the methods is possible in studies involving the determination of an optimum polynomial equation for fitting a large number of data points. The data considered in the least squares problem are experimental thermodynamic pressure-density-temperature data. Attention is given to a description of stepwise multiple regression and the use of stepwise regression in the determination of an equation of state for oxygen and nitrogen.

Iterative Splitting Methods for Differential Equations

CERN Document Server

Geiser, Juergen

2011-01-01

Iterative Splitting Methods for Differential Equations explains how to solve evolution equations via novel iterative-based splitting methods that efficiently use computational and memory resources. It focuses on systems of parabolic and hyperbolic equations, including convection-diffusion-reaction equations, heat equations, and wave equations. In the theoretical part of the book, the author discusses the main theorems and results of the stability and consistency analysis for ordinary differential equations. He then presents extensions of the iterative splitting methods to partial differential

Transport equation solving methods

International Nuclear Information System (INIS)

Granjean, P.M.

1984-06-01

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

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

Implicit methods for equation-free analysis: convergence results and analysis of emergent waves in microscopic traffic models

DEFF Research Database (Denmark)

Marschler, Christian; Sieber, Jan; Berkemer, Rainer

2014-01-01

We introduce a general formulation for an implicit equation-free method in the setting of slow-fast systems. First, we give a rigorous convergence result for equation-free analysis showing that the implicitly defined coarse-level time stepper converges to the true dynamics on the slow manifold...... against the direction of traffic. Equation-free analysis enables us to investigate the behavior of the microscopic traffic model on a macroscopic level. The standard deviation of cars' headways is chosen as the macroscopic measure of the underlying dynamics such that traveling wave solutions correspond...... to equilibria on the macroscopic level in the equation-free setup. The collapse of the traffic jam to the free flow then corresponds to a saddle-node bifurcation of this macroscopic equilibrium. We continue this bifurcation in two parameters using equation-free analysis....

Numerical method for the nonlinear Fokker-Planck equation

International Nuclear Information System (INIS)

Zhang, D.S.; Wei, G.W.; Kouri, D.J.; Hoffman, D.K.

1997-01-01

A practical method based on distributed approximating functionals (DAFs) is proposed for numerically solving a general class of nonlinear time-dependent Fokker-Planck equations. The method relies on a numerical scheme that couples the usual path-integral concept to the DAF idea. The high accuracy and reliability of the method are illustrated by applying it to an exactly solvable nonlinear Fokker-Planck equation, and the method is compared with the accurate K-point Stirling interpolation formula finite-difference method. The approach is also used successfully to solve a nonlinear self-consistent dynamic mean-field problem for which both the cumulant expansion and scaling theory have been found by Drozdov and Morillo [Phys. Rev. E 54, 931 (1996)] to be inadequate to describe the occurrence of a long-lived transient bimodality. The standard interpretation of the transient bimodality in terms of the flat region in the kinetic potential fails for the present case. An alternative analysis based on the effective potential of the Schroedinger-like Fokker-Planck equation is suggested. Our analysis of the transient bimodality is strongly supported by two examples that are numerically much more challenging than other examples that have been previously reported for this problem. copyright 1997 The American Physical Society

Multiple time scale methods in tokamak magnetohydrodynamics

International Nuclear Information System (INIS)

Jardin, S.C.

1984-01-01

Several methods are discussed for integrating the magnetohydrodynamic (MHD) equations in tokamak systems on other than the fastest time scale. The dynamical grid method for simulating ideal MHD instabilities utilizes a natural nonorthogonal time-dependent coordinate transformation based on the magnetic field lines. The coordinate transformation is chosen to be free of the fast time scale motion itself, and to yield a relatively simple scalar equation for the total pressure, P = p + B 2 /2μ 0 , which can be integrated implicitly to average over the fast time scale oscillations. Two methods are described for the resistive time scale. The zero-mass method uses a reduced set of two-fluid transport equations obtained by expanding in the inverse magnetic Reynolds number, and in the small ratio of perpendicular to parallel mobilities and thermal conductivities. The momentum equation becomes a constraint equation that forces the pressure and magnetic fields and currents to remain in force balance equilibrium as they evolve. The large mass method artificially scales up the ion mass and viscosity, thereby reducing the severe time scale disparity between wavelike and diffusionlike phenomena, but not changing the resistive time scale behavior. Other methods addressing the intermediate time scales are discussed

Symmetries of the Euler compressible flow equations for general equation of state

Energy Technology Data Exchange (ETDEWEB)

Boyd, Zachary M. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Ramsey, Scott D. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Baty, Roy S. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2015-10-15

The Euler compressible flow equations exhibit different Lie symmetries depending on the equation of state (EOS) of the medium in which the flow occurs. This means that, in general, different types of similarity solution will be available in different flow media. We present a comprehensive classification of all EOS’s to which the Euler equations apply, based on the Lie symmetries admitted by the corresponding flow equations, restricting to the case of 1-D planar, cylindrical, or spherical geometry. The results are conveniently summarized in tables. This analysis also clarifies past work by Axford and Ovsiannikov on symmetry classification.

Solutions of hyperbolic equations with the CIP-BS method

International Nuclear Information System (INIS)

Utsumi, Takayuki; Koga, James; Yamagiwa, Mitsuru; Yabe, Takashi; Aoki, Takayuki

2004-01-01

In this paper, we show that a new numerical method, the Constrained Interpolation Profile - Basis Set (CIP-BS) method, can solve general hyperbolic equations efficiently. This method uses a simple polynomial basis set that is easily extendable to any desired higher-order accuracy. The interpolating profile is chosen so that the subgrid scale solution approaches the local real solution owing to the constraints from the spatial derivatives of the master equations. Then, introducing scalar products, the linear and nonlinear partial differential equations are uniquely reduced to the ordinary differential equations for values and spatial derivatives at the grid points. The method gives stable, less diffusive, and accurate results. It is successfully applied to the continuity equation, the Burgers equation, the Korteweg-de Vries equation, and one-dimensional shock tube problems. (author)

Equation of state with scale-invariant hidden local symmetry and gravitational waves

Directory of Open Access Journals (Sweden)

Lee Hyun Kyu

2018-01-01

Full Text Available The equation of state (EoS for the effective theory proposed recently in the frame work of the scale-invariant hidden local symmetry is discussed briefly. The EoS is found to be relatively stiffer at lower density and but relatively softer at higher density. The particular features of EoS on the gravitational waves are discussed. A relatively stiffer EoS for the neutron stars with the lower density induces a larger deviation of the gravitational wave form from the point-particle-approximation. On the other hand, a relatively softer EoS for the merger remnant of the higher density inside might invoke a possibility of the immediate formation of a black hole for short gamma ray bursts or the appearance of the higher peak frequency for gravitational waves from remnant oscillations. It is anticipated that this particular features could be probed in detail by the detections of gravitational waves from the binary neutron star mergers.

Partial differential equations with variable exponents variational methods and qualitative analysis

CERN Document Server

Radulescu, Vicentiu D

2015-01-01

Partial Differential Equations with Variable Exponents: Variational Methods and Qualitative Analysis provides researchers and graduate students with a thorough introduction to the theory of nonlinear partial differential equations (PDEs) with a variable exponent, particularly those of elliptic type. The book presents the most important variational methods for elliptic PDEs described by nonhomogeneous differential operators and containing one or more power-type nonlinearities with a variable exponent. The authors give a systematic treatment of the basic mathematical theory and constructive meth

Analysis and applications of a group contribution sPC-SAFT equation of state

DEFF Research Database (Denmark)

Tihic, Amra; von Solms, Nicolas; Michelsen, Michael Locht

2009-01-01

A group contribution (GC) method for estimating pure compound parameters for the molecular-based perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state (EoS) is proposed in a previous work [A. Tihic, G.M. Kontogeorgis, N. von Solms, M.L Michelsen, L Constantinou. Ind. En...

Analysis of equations of state and temperature dependence of thermal expansivity and bulk modulus for silicon

International Nuclear Information System (INIS)

Pandya, Tushar C; Bhatt, Apoorva D; Thakar, Nilesh A

2012-01-01

In the present paper an attempt has been made for the comparative study of different equations of state for silicon (Phase-1, cubic diamond structure) in the pressure range of 0-11 GPa. We compare the results of different equations of state (EOS) with available experimental data. The Kwon and Kim EOS is found to give far better agreement with the available experimental data. Results obtained by Poirier-Tarantola, Vinet, Tait and Suzuki's equations of state are not giving satisfactory agreement with the available experimental data. In the present study simple methods based on thermodynamic functions are presented to investigate the temperature dependence of thermal expansivity and bulk modulus for silicon. The results are reported for silicon. The calculated values of thermal expansivity are in good agreement with experimental data.

Equation of state, universal profiles, scaling and macroscopic quantum effects in warm dark matter galaxies

Energy Technology Data Exchange (ETDEWEB)

Vega, H.J. de [Sorbonne Universites, Universite Pierre et Marie Curie UPMC Paris VI, LPTHE CNRS UMR 7589, Paris Cedex 05 (France); Sanchez, N.G. [Observatoire de Paris PSL Research University, Sorbonne Universites UPMC Paris VI, Observatoire de Paris, LERMA CNRS UMR 8112, Paris (France)

2017-02-15

The Thomas-Fermi approach to galaxy structure determines self-consistently and non-linearly the gravitational potential of the fermionic warm dark matter (WDM) particles given their quantum distribution function f(E). This semiclassical framework accounts for the quantum nature and high number of DM particles, properly describing gravitational bounded and quantum macroscopic systems as neutron stars, white dwarfs and WDM galaxies. We express the main galaxy magnitudes as the halo radius r{sub h}, mass M{sub h}, velocity dispersion and phase space density in terms of the surface density which is important to confront to observations. From these expressions we derive the general equation of state for galaxies, i.e., the relation between pressure and density, and provide its analytic expression. Two regimes clearly show up: (1) Large diluted galaxies for M{sub h} >or similar 2.3 x 10{sup 6} M {sub CircleDot} and effective temperatures T{sub 0} > 0.017 K described by the classical self-gravitating WDM Boltzman gas with a space-dependent perfect gas equation of state, and (2) Compact dwarf galaxies for 1.6 x 10{sup 6} M {sub CircleDot} >or similar M{sub h} >or similar M{sub h,min} ≅ 3.10 x 10{sup 4} (2 keV/m){sup (16)/(5)} M {sub CircleDot}, T{sub 0} < 0.011 K described by the quantum fermionic WDM regime with a steeper equation of state close to the degenerate state. In particular, the T{sub 0} = 0 degenerate or extreme quantum limit yields the most compact and smallest galaxy. In the diluted regime, the halo radius r{sub h}, the squared velocity v{sup 2}(r{sub h}) and the temperature T{sub 0} turn to exhibit square-root of M{sub h} scaling laws. The normalized density profiles ρ(r)/ρ(0) and the normalized velocity profiles v{sup 2}(r)/v{sup 2}(0) are universal functions of r/r{sub h} reflecting the WDM perfect gas behavior in this regime. These theoretical results contrasted to robust and independent sets of galaxy data remarkably reproduce the observations. For

Efficient propagation of the hierarchical equations of motion using the matrix product state method

Science.gov (United States)

Shi, Qiang; Xu, Yang; Yan, Yaming; Xu, Meng

2018-05-01

We apply the matrix product state (MPS) method to propagate the hierarchical equations of motion (HEOM). It is shown that the MPS approximation works well in different type of problems, including boson and fermion baths. The MPS method based on the time-dependent variational principle is also found to be applicable to HEOM with over one thousand effective modes. Combining the flexibility of the HEOM in defining the effective modes and the efficiency of the MPS method thus may provide a promising tool in simulating quantum dynamics in condensed phases.

A note on the convergence of the Zakharov-Kuznetsov equation by homotopy analysis method

Directory of Open Access Journals (Sweden)

Amir Fallahzadeh

2014-07-01

Full Text Available In this paper, the convergence of Zakharov-Kuznetsov (ZK equation by homotopy analysis method (HAM is investigated. A theorem is proved to guarantee the convergence of HAMand to find the series solution of this equation via a reliable algorithm.

Performance and scaling of locally-structured grid methods forpartial differential equations

Energy Technology Data Exchange (ETDEWEB)

Colella, Phillip; Bell, John; Keen, Noel; Ligocki, Terry; Lijewski, Michael; Van Straalen, Brian

2007-07-19

In this paper, we discuss some of the issues in obtaining high performance for block-structured adaptive mesh refinement software for partial differential equations. We show examples in which AMR scales to thousands of processors. We also discuss a number of metrics for performance and scalability that can provide a basis for understanding the advantages and disadvantages of this approach.

Uncertainty assessment of equations of state with application to an organic Rankine cycle

DEFF Research Database (Denmark)

Frutiger, Jerome; Bell, Ian; O’Connell, John P.

2017-01-01

Evaluations of equations of state (EoS) should include uncertainty. This study presents a genericmethod to analyse EoS from a detailed uncertainty analysis of the mathematical form and the dataused to obtain EoS parameter values. The method is illustrated by comparison of Soave–Redlich–Kwong (SRK......) cubic EoS with perturbed-chain statistical associating fluid theory (PC-SAFT) EoS for anorganic Rankine cycle (ORC) for heat recovery to power fromthe exhaust gas of a marine diesel engineusing cyclopentane as working fluid. Uncertainties of the EoS input parameters including......Evaluations of equations of state (EoS) should include uncertainty. This study presents a genericmethod to analyse EoS from a detailed uncertainty analysis of the mathematical form and the dataused to obtain EoS parameter values. The method is illustrated by comparison of Soave–Redlich–Kwong (SRK...

Two-scale approach to oscillatory singularly perturbed transport equations

CERN Document Server

Frénod, Emmanuel

2017-01-01

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

Quadratic inner element subgrid scale discretisation of the Boltzmann transport equation

International Nuclear Information System (INIS)

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

2012-01-01

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

A Hamiltonian-based derivation of Scaled Boundary Finite Element Method for elasticity problems

International Nuclear Information System (INIS)

Hu Zhiqiang; Lin Gao; Wang Yi; Liu Jun

2010-01-01

The Scaled Boundary Finite Method (SBFEM) is a semi-analytical solution approach for solving partial differential equation. For problem in elasticity, the governing equations can be obtained by mechanically based formulation, Scaled-boundary-transformation-based formulation and principle of virtual work. The governing equations are described in the frame of Lagrange system and the unknowns are displacements. But in the solution procedure, the auxiliary variables are introduced and the equations are solved in the state space. Based on the observation that the duality system to solve elastic problem proposed by W.X. Zhong is similar to the above solution approach, the discretization of the SBFEM and the duality system are combined to derive the governing equations in the Hamilton system by introducing the dual variables in this paper. The Precise Integration Method (PIM) used in Duality system is also an efficient method for the solution of the governing equations of SBFEM in displacement and boundary stiffness matrix especially for the case which results some numerical difficulties in the usually uses the eigenvalue method. Numerical examples are used to demonstrate the validity and effectiveness of the PIM for solution of boundary static stiffness.

Bound states of the Dirac equation with some physical potentials by the Nikiforov-Uvarov method

Energy Technology Data Exchange (ETDEWEB)

Setare, Mohammad R; Haidari, S [Department of Physics, University of Kurdistan, Pasdaran Avenue, Sanandaj (Iran, Islamic Republic of)], E-mail: rezakord@ipm.ir, E-mail: heidary.somayeh@gmail.com

2010-01-15

Exact analytical solutions for the s-wave Dirac equation with the reflectionless-type, Rosen-Morse and Manning-Rosen potentials are obtained, under the condition of spin symmetry. We obtained bound state energy eigenvalues and corresponding spinor wave function in the framework of the Nikiforov-Uvarov (NU) method.

Comparative analysis among several methods used to solve the point kinetic equations

International Nuclear Information System (INIS)

Nunes, Anderson L.; Goncalves, Alessandro da C.; Martinez, Aquilino S.; Silva, Fernando Carvalho da

2007-01-01

The main objective of this work consists on the methodology development for comparison of several methods for the kinetics equations points solution. The evaluated methods are: the finite differences method, the stiffness confinement method, improved stiffness confinement method and the piecewise constant approximations method. These methods were implemented and compared through a systematic analysis that consists basically of confronting which one of the methods consume smaller computational time with higher precision. It was calculated the relative which function is to combine both criteria in order to reach the goal. Through the analyses of the performance factor it is possible to choose the best method for the solution of point kinetics equations. (author)

Comparative analysis among several methods used to solve the point kinetic equations

Energy Technology Data Exchange (ETDEWEB)

Nunes, Anderson L.; Goncalves, Alessandro da C.; Martinez, Aquilino S.; Silva, Fernando Carvalho da [Universidade Federal, Rio de Janeiro, RJ (Brazil). Coordenacao dos Programas de Pos-graduacao de Engenharia. Programa de Engenharia Nuclear; E-mails: alupo@if.ufrj.br; agoncalves@con.ufrj.br; aquilino@lmp.ufrj.br; fernando@con.ufrj.br

2007-07-01

The main objective of this work consists on the methodology development for comparison of several methods for the kinetics equations points solution. The evaluated methods are: the finite differences method, the stiffness confinement method, improved stiffness confinement method and the piecewise constant approximations method. These methods were implemented and compared through a systematic analysis that consists basically of confronting which one of the methods consume smaller computational time with higher precision. It was calculated the relative which function is to combine both criteria in order to reach the goal. Through the analyses of the performance factor it is possible to choose the best method for the solution of point kinetics equations. (author)

Reproducible and Verifiable Equations of State Using Microfabricated Materials

Science.gov (United States)

Martin, J. F.; Pigott, J. S.; Panero, W. R.

2017-12-01

Accurate interpretation of observable geophysical data, relevant to the structure, composition, and evolution of planetary interiors, requires precise determination of appropriate equations of state. We present the synthesis of controlled-geometry nanofabricated samples and insulation layers for the laser-heated diamond anvil cell. We present electron-gun evaporation, sputter deposition, and photolithography methods to mass-produce Pt/SiO2/Fe/SiO2 stacks and MgO insulating disks to be used in LHDAC experiments to reduce uncertainties in equation of state measurements due to large temperature gradients. We present a reanalysis of published iron PVT data to establish a statistically-valid extrapolation of the equation of state to inner core conditions with quantified uncertainties, addressing the complication of covariance in equation of state parameters. We use this reanalysis, together with the synthesized samples, to propose a scheme for measurement and validation of high-precision equations of state relevant to the Earth and super-Earth exoplanets.

A reduction method for phase equilibrium calculations with cubic equations of state

Directory of Open Access Journals (Sweden)

D. V. Nichita

2006-09-01

Full Text Available In this work we propose a new reduction method for phase equilibrium calculations using a general form of cubic equations of state (CEOS. The energy term in the CEOS is a quadratic form, which is diagonalized by applying a linear transformation. The number of the reduction parameters is related to the rank of the matrix C with elements (1-Cij, where Cij denotes the binary interaction parameters (BIPs. The dimensionality of the problem depends only on the number of reduction parameters, and is independent of the number of components in the mixture.

Solutions of the Noh Problem for Various Equations of State Using Lie Groups

International Nuclear Information System (INIS)

Axford, R.A.

1998-01-01

A method for developing invariant equations of state for which solutions of the Noh problem will exist is developed. The ideal gas equation of state is shown to be a special case of the general method. Explicit solutions of the Noh problem in planar, cylindrical and spherical geometry are determined for a Mie-Gruneisen and the stiff gas equation of state

Finite difference method for inner-layer equations in the resistive MagnetoHydroDynamic stability analysis

International Nuclear Information System (INIS)

Tokuda, Shinji; Watanabe, Tomoko.

1996-08-01

The matching problem in resistive MagnetoHydroDynamic stability analysis by the asymptotic matching method has been reformulated as an initial-boundary value problem for the inner-layer equations describing the plasma dynamics in the thin layer around a rational surface. The third boundary conditions at boundaries of a finite interval are imposed on the inner layer equations in the formulation instead of asymptotic conditions at infinities. The finite difference method for this problem has been applied to model equations whose solutions are known in a closed form. It has been shown that the initial value problem and the associated eigenvalue problem for the model equations can be solved by the finite difference method with numerical stability. The formulation presented here enables the asymptotic matching method to be a practical method for the resistive MHD stability analysis. (author)

Information theoretical methods as discerning quantifiers of the equations of state of neutron stars

Energy Technology Data Exchange (ETDEWEB)

Avellar, M.G.B. de, E-mail: mgb.avellar@iag.usp.br [Instituto de Astronomia, Geofísica e Ciências Atmosféricas – Universidade de São Paulo, Rua do Matão 1226, Cidade Universitária, 05508-090, São Paulo, SP (Brazil); Souza, R.A. de, E-mail: rodrigo.souza@usp.br [Instituto de Astronomia, Geofísica e Ciências Atmosféricas – Universidade de São Paulo, Rua do Matão 1226, Cidade Universitária, 05508-090, São Paulo, SP (Brazil); Horvath, J.E., E-mail: foton@iag.usp.br [Instituto de Astronomia, Geofísica e Ciências Atmosféricas – Universidade de São Paulo, Rua do Matão 1226, Cidade Universitária, 05508-090, São Paulo, SP (Brazil); Paret, D.M., E-mail: dmanreza@fisica.uh.cu [Facultad de Física, Universidad de la Habana, San Lázaro y L, Vedado La Habana, 10400 (Cuba)

2014-11-07

In this work we use the statistical measures of information entropy, disequilibrium and complexity to discriminate different approaches and parametrizations for different equations of state for quark stars. We confirm the usefulness of such quantities to quantify the role of interactions in such stars. We find that within this approach, a quark matter equation of state such as SU(2) NJL with vectorial coupling and phase transition is slightly favoured and deserves deeper studies. - Highlights: • We used information theory tools to discern different compositions for compact stars. • Hadronic and quark stars analogues behave differently when analyzed with these tools. • The effects of different equations of state are singled out in this work.

Solvation quantities from a COSMO-RS equation of state

International Nuclear Information System (INIS)

Panayiotou, C.; Tsivintzelis, I.; Aslanidou, D.; Hatzimanikatis, V.

2015-01-01

Highlights: • Extension of the successful COSMO-RS model to an equation-of-state model. • Two scaling constants, obtained from atom-specific contributions. • Overall estimation of the solvation quantities and contributions. - Abstract: This work focuses on the extension of the successful COSMO-RS model of mixtures into an equation-of-state model of fluids and its application for the estimation of solvation/hydration quantities of a variety of chemical substances. These quantities include free-energies, enthalpies and entropies of hydration as well as the separate contributions to each of them. Emphasis is given on the estimation of contributions from the conformational changes of solutes upon solvation and the associated restructuring of solvent in its immediate neighborhood. COSMO-RS is a quantum-mechanics based group/segment contribution model in which the Quasi-Chemical (QC) approach is used for the description of the non-random distribution of interacting segments in the system. Thus, the equation-of-state development is done through such a QC framework. The new model will not need any adjustable parameters for the strong specific interactions, such as hydrogen bonds, since they will be provided by the quantum-mechanics based cosmo-files – a key feature of COSMO-RS model. It will need, however, one volumetric and one energy parameter per fluid, which are scaling constants or molecular descriptors of the fluid and are obtained from rather easily available data such as densities, boiling points, vapor pressures, heats of vaporization or second virial coefficients. The performance and the potential of the new equation-of-state model to become a fully predictive model are critically discussed

FLRW Cosmology with Horava-Lifshitz Gravity: Impacts of Equations of State

Science.gov (United States)

Tawfik, A.; Abou El Dahab, E.

2017-07-01

Inspired by Lifshitz theory for quantum critical phenomena in condensed matter, Horava proposed a theory for quantum gravity with an anisotropic scaling in ultraviolet. In Horava-Lifshitz gravity (HLG), we have studied the impacts of six types of equations of state on the evolution of various cosmological parameters such as Hubble parameters and scale factor. From the comparison of the general relativity gravity with the HLG with detailed and without with non-detailed balance conditions, remarkable differences are found. Also, a noticeable dependence of singular and non-singular Big Bang on the equations of state is observed. We conclude that HLG explains various epochs in the early universe and might be able to reproduce the entire cosmic history with and without singular Big Bang.

The equation of state of liquid Flibe

International Nuclear Information System (INIS)

Chen, Xiang M.; Schrock, V.E.; Peterson, P.F.

1991-01-01

Flibe (Li 2 BeF 4 ) is a candidate material for the liquid blanket in the HYLIFE-2 fusion reactor. The thermodynamic properties of the material are important for the study of thermohydraulic behavior of the concept design, including the compressible analysis of the blanket isochoric heating problem and resulting jet breakup. The equation of state provides the relationship between all the thermodynamic properties. Previously, a soft sphere model of liquid equation of state was used for describing a number of liquid metals. In this paper we have fitted the available experimental data for liquid Flibe with a modified soft sphere model. 5 refs

Application of cellular neural network (CNN) method to the nuclear reactor dynamics equations

International Nuclear Information System (INIS)

Hadad, K.; Piroozmand, A.

2007-01-01

This paper describes the application of a multilayer cellular neural network (CNN) to model and solve the nuclear reactor dynamic equations. An equivalent electrical circuit is analyzed and the governing equations of a bare, homogeneous reactor core are modeled via CNN. The validity of the CNN result is compared with numerical solution of the system of nonlinear governing partial differential equations (PDE) using MATLAB. Steady state as well as transient simulations, show very good comparison between the two methods. We used our CNN model to simulate space-time response of different reactivity excursions in a typical nuclear reactor. On line solution of reactor dynamic equations is used as an aid to reactor operation decision making. The complete algorithm could also be implemented using very large scale integrated circuit (VLSI) circuitry. The efficiency of the calculation method makes it useful for small size nuclear reactors such as the ones used in space missions

Local density approximation for a perturbative equation of state

International Nuclear Information System (INIS)

Astrakharchik, G. E.

2005-01-01

Knowledge of a series expansion of the equation of state provides a deep insight into the physical nature of a quantum system. Starting from a generic 'perturbative' equation of state of a homogeneous ultracold gas we make predictions for the properties of the gas in the presence of harmonic confinement. The local density approximation is used to obtain the chemical potential, total and release energies, Thomas-Fermi size, and density profile of a trapped system in three-, two-, and one-dimensional geometries. The frequencies of the lowest breathing modes are calculated using scaling and sum-rule approaches and could be used in an experiment as a high-precision tool for obtaining the expansion terms of the equation of state. The derived formalism is applied to dilute Bose and Fermi gases in different dimensions and to integrable one-dimensional models. The physical meaning of the expansion terms in a number of systems is discussed

A Hartree-Fock-Slater-Boltzmann-Saha method for detailed atomic structure and equation of state of plasmas

International Nuclear Information System (INIS)

Jiang Minhao; Meng Xujun

2005-01-01

The effect of the free electron background in plasmas is introduced in Hartree-Fock-Slater self-consistent field atomic model to correct the single electron energies for each electron configuration, and to provide accurate atomic data for Boltzmann-Saha equation. In the iteration process chemical potential is adjusted to change the free electron background to satisfy simultaneously the conservation of the free electrons in Saha equation as well as in Hartree-Fock-Slater self-consistent field atomic model. As examples the equations of state of the carbon and aluminum plasmas are calculated to show the applicability of this method. (authors)

Numerical Analysis of Partial Differential Equations

CERN Document Server

Lui, S H

2011-01-01

A balanced guide to the essential techniques for solving elliptic partial differential equations Numerical Analysis of Partial Differential Equations provides a comprehensive, self-contained treatment of the quantitative methods used to solve elliptic partial differential equations (PDEs), with a focus on the efficiency as well as the error of the presented methods. The author utilizes coverage of theoretical PDEs, along with the nu merical solution of linear systems and various examples and exercises, to supply readers with an introduction to the essential concepts in the numerical analysis

A method for statistical steady state thermal analysis of reactor cores

International Nuclear Information System (INIS)

Whetton, P.A.

1980-01-01

This paper presents a method for performing a statistical steady state thermal analysis of a reactor core. The technique is only outlined here since detailed thermal equations are dependent on the core geometry. The method has been applied to a pressurised water reactor core and the results are presented for illustration purposes. Random hypothetical cores are generated using the Monte-Carlo method. The technique shows that by splitting the parameters into two types, denoted core-wise and in-core, the Monte Carlo method may be used inexpensively. The idea of using extremal statistics to characterise the low probability events (i.e. the tails of a distribution) is introduced together with a method of forming the final probability distribution. After establishing an acceptable probability of exceeding a thermal design criterion, the final probability distribution may be used to determine the corresponding thermal response value. If statistical and deterministic (i.e. conservative) thermal response values are compared, information on the degree of pessimism in the deterministic method of analysis may be inferred and the restrictive performance limitations imposed by this method relieved. (orig.)

Salzburger State Reactance Scale (SSR Scale): Validation of a Scale Measuring State Reactance.

Science.gov (United States)

Sittenthaler, Sandra; Traut-Mattausch, Eva; Steindl, Christina; Jonas, Eva

This paper describes the construction and empirical evaluation of an instrument for measuring state reactance, the Salzburger State Reactance (SSR) Scale. The results of a confirmatory factor analysis supported a hypothesized three-factor structure: experience of reactance, aggressive behavioral intentions, and negative attitudes. Correlations with divergent and convergent measures support the validity of this structure. The SSR Subscales were strongly related to the other state reactance measures. Moreover, the SSR Subscales showed modest positive correlations with trait measures of reactance. The SSR Subscales correlated only slightly or not at all with neighboring constructs (e.g., autonomy, experience of control). The only exception was fairness scales, which showed moderate correlations with the SSR Subscales. Furthermore, a retest analysis confirmed the temporal stability of the scale. Suggestions for further validation of this questionnaire are discussed.

STEADY STATE AND PSEUDO-TRANSIENT ELECTRIC POTENTIAL USING THE POISSONBOLTZMANN EQUATION

Directory of Open Access Journals (Sweden)

L. C. dos Santos

2015-03-01

Full Text Available A method for analysis of the electric potential profile in saline solutions was developed for systems with one or two infinite flat plates. A modified Poisson-Boltzmann equation, taking into account nonelectrostatic interactions between ions and surfaces, was used. To solve the stated problem in the steady-state approach the finite-difference method was used. For the formulated pseudo-transient problem, we solved the set of ordinary differential equations generated from the algebraic equations of the stationary case. A case study was also carried out in relation to temperature, solution concentration, surface charge and salt-type. The results were validated by the stationary problem solution, which had also been used to verify the ionic specificity for different salts. The pseudo-transient approach allowed a better understanding of the dynamic behavior of the ion-concentration profile and other properties due to the surface charge variation.

Time Scale in Least Square Method

Directory of Open Access Journals (Sweden)

Özgür Yeniay

2014-01-01

Full Text Available Study of dynamic equations in time scale is a new area in mathematics. Time scale tries to build a bridge between real numbers and integers. Two derivatives in time scale have been introduced and called as delta and nabla derivative. Delta derivative concept is defined as forward direction, and nabla derivative concept is defined as backward direction. Within the scope of this study, we consider the method of obtaining parameters of regression equation of integer values through time scale. Therefore, we implemented least squares method according to derivative definition of time scale and obtained coefficients related to the model. Here, there exist two coefficients originating from forward and backward jump operators relevant to the same model, which are different from each other. Occurrence of such a situation is equal to total number of values of vertical deviation between regression equations and observation values of forward and backward jump operators divided by two. We also estimated coefficients for the model using ordinary least squares method. As a result, we made an introduction to least squares method on time scale. We think that time scale theory would be a new vision in least square especially when assumptions of linear regression are violated.

Nuclear methods and the nuclear equation of state

CERN Document Server

1999-01-01

The theoretical study of the nuclear equation of state (EOS) is a field of research which deals with most of the fundamental problems of nuclear physics. This book gives an overview of the present status of the microscopic theory of the nuclear EOS. Its aim is essentially twofold: first, to serve as a textbook for students entering the field, by covering the different subjects as exhaustively and didactically as possible; second, to be a reference book for all researchers active in the theory of nuclear matter, by providing a report on the latest developments. Special emphasis is given to the

Thermodynamically constrained correction to ab initio equations of state

International Nuclear Information System (INIS)

French, Martin; Mattsson, Thomas R.

2014-01-01

We show how equations of state generated by density functional theory methods can be augmented to match experimental data without distorting the correct behavior in the high- and low-density limits. The technique is thermodynamically consistent and relies on knowledge of the density and bulk modulus at a reference state and an estimation of the critical density of the liquid phase. We apply the method to four materials representing different classes of solids: carbon, molybdenum, lithium, and lithium fluoride. It is demonstrated that the corrected equations of state for both the liquid and solid phases show a significantly reduced dependence of the exchange-correlation functional used.

Validation of the activity expansion method with ultrahigh pressure shock equations of state

Science.gov (United States)

Rogers, Forrest J.; Young, David A.

1997-11-01

Laser shock experiments have recently been used to measure the equation of state (EOS) of matter in the ultrahigh pressure region between condensed matter and a weakly coupled plasma. Some ultrahigh pressure data from nuclear-generated shocks are also available. Matter at these conditions has proven very difficult to treat theoretically. The many-body activity expansion method (ACTEX) has been used for some time to calculate EOS and opacity data in this region, for use in modeling inertial confinement fusion and stellar interior plasmas. In the present work, we carry out a detailed comparison with the available experimental data in order to validate the method. The agreement is good, showing that ACTEX adequately describes strongly shocked matter.

Scaling of differential equations

CERN Document Server

Langtangen, Hans Petter

2016-01-01

The book serves both as a reference for various scaled models with corresponding dimensionless numbers, and as a resource for learning the art of scaling. A special feature of the book is the emphasis on how to create software for scaled models, based on existing software for unscaled models. Scaling (or non-dimensionalization) is a mathematical technique that greatly simplifies the setting of input parameters in numerical simulations. Moreover, scaling enhances the understanding of how different physical processes interact in a differential equation model. Compared to the existing literature, where the topic of scaling is frequently encountered, but very often in only a brief and shallow setting, the present book gives much more thorough explanations of how to reason about finding the right scales. This process is highly problem dependent, and therefore the book features a lot of worked examples, from very simple ODEs to systems of PDEs, especially from fluid mechanics. The text is easily accessible and exam...

The Convergence Study of the Homotopy Analysis Method for Solving Nonlinear Volterra-Fredholm Integrodifferential Equations

Directory of Open Access Journals (Sweden)

Behzad Ghanbari

2014-01-01

Full Text Available We aim to study the convergence of the homotopy analysis method (HAM in short for solving special nonlinear Volterra-Fredholm integrodifferential equations. The sufficient condition for the convergence of the method is briefly addressed. Some illustrative examples are also presented to demonstrate the validity and applicability of the technique. Comparison of the obtained results HAM with exact solution shows that the method is reliable and capable of providing analytic treatment for solving such equations.

Toward enabling large-scale open-shell equation-of-motion coupled cluster calculations: triplet states of β-carotene

Energy Technology Data Exchange (ETDEWEB)

Hu, Hanshi; Bhaskaran-Nair, Kiran; Apra, Edoardo; Govind, Niranjan; Kowalski, Karol

2014-10-02

In this paper we discuss the application of novel parallel implementation of the coupled cluster (CC) and equation-of-motion coupled cluster methods (EOMCC) in calculations of excitation energies of triplet states in beta-carotene. Calculated excitation energies are compared with experimental data, where available. We also provide a detailed description of the new parallel algorithms for iterative CC and EOMCC models involving single and doubles excitations.

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

International Nuclear Information System (INIS)

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

1976-12-01

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

Numerical methods and analysis of multiscale problems

CERN Document Server

Madureira, Alexandre L

2017-01-01

This book is about numerical modeling of multiscale problems, and introduces several asymptotic analysis and numerical techniques which are necessary for a proper approximation of equations that depend on different physical scales. Aimed at advanced undergraduate and graduate students in mathematics, engineering and physics – or researchers seeking a no-nonsense approach –, it discusses examples in their simplest possible settings, removing mathematical hurdles that might hinder a clear understanding of the methods. The problems considered are given by singular perturbed reaction advection diffusion equations in one and two-dimensional domains, partial differential equations in domains with rough boundaries, and equations with oscillatory coefficients. This work shows how asymptotic analysis can be used to develop and analyze models and numerical methods that are robust and work well for a wide range of parameters.

on the properties of solutions and some applications on the TOV differential equation with a model of nuclear equation of state

International Nuclear Information System (INIS)

Esmail, S.F.H.

2006-01-01

the mathematical formulation of numerous physical problems results in differential equations actually non-linear differential equations . in our study we are interested in solutions of differential equations which describe the structure of neutron star in non-relativistic and relativistic cases. the aim of this work is to determine the mass and the radius of a neutron star, by solving the tolmann-oppenheimer-volkoff (TOV) differential equation using different models of the nuclear equation of state (EOS). analytically solutions are obtained for a simple form of the nuclear equation of state of Clayton model and poly trope model. for a more realistic equation of state the TOV differential equation is solved numerically using rung -Kutta method

Multiscale functions, scale dynamics, and applications to partial differential equations

Science.gov (United States)

Cresson, Jacky; Pierret, Frédéric

2016-05-01

Modeling phenomena from experimental data always begins with a choice of hypothesis on the observed dynamics such as determinism, randomness, and differentiability. Depending on these choices, different behaviors can be observed. The natural question associated to the modeling problem is the following: "With a finite set of data concerning a phenomenon, can we recover its underlying nature? From this problem, we introduce in this paper the definition of multi-scale functions, scale calculus, and scale dynamics based on the time scale calculus [see Bohner, M. and Peterson, A., Dynamic Equations on Time Scales: An Introduction with Applications (Springer Science & Business Media, 2001)] which is used to introduce the notion of scale equations. These definitions will be illustrated on the multi-scale Okamoto's functions. Scale equations are analysed using scale regimes and the notion of asymptotic model for a scale equation under a particular scale regime. The introduced formalism explains why a single scale equation can produce distinct continuous models even if the equation is scale invariant. Typical examples of such equations are given by the scale Euler-Lagrange equation. We illustrate our results using the scale Newton's equation which gives rise to a non-linear diffusion equation or a non-linear Schrödinger equation as asymptotic continuous models depending on the particular fractional scale regime which is considered.

Stochastic analysis of complex reaction networks using binomial moment equations.

Science.gov (United States)

Barzel, Baruch; Biham, Ofer

2012-09-01

The stochastic analysis of complex reaction networks is a difficult problem because the number of microscopic states in such systems increases exponentially with the number of reactive species. Direct integration of the master equation is thus infeasible and is most often replaced by Monte Carlo simulations. While Monte Carlo simulations are a highly effective tool, equation-based formulations are more amenable to analytical treatment and may provide deeper insight into the dynamics of the network. Here, we present a highly efficient equation-based method for the analysis of stochastic reaction networks. The method is based on the recently introduced binomial moment equations [Barzel and Biham, Phys. Rev. Lett. 106, 150602 (2011)]. The binomial moments are linear combinations of the ordinary moments of the probability distribution function of the population sizes of the interacting species. They capture the essential combinatorics of the reaction processes reflecting their stoichiometric structure. This leads to a simple and transparent form of the equations, and allows a highly efficient and surprisingly simple truncation scheme. Unlike ordinary moment equations, in which the inclusion of high order moments is prohibitively complicated, the binomial moment equations can be easily constructed up to any desired order. The result is a set of equations that enables the stochastic analysis of complex reaction networks under a broad range of conditions. The number of equations is dramatically reduced from the exponential proliferation of the master equation to a polynomial (and often quadratic) dependence on the number of reactive species in the binomial moment equations. The aim of this paper is twofold: to present a complete derivation of the binomial moment equations; to demonstrate the applicability of the moment equations for a representative set of example networks, in which stochastic effects play an important role.

On parametrised cold dense matter equation of state inference

Science.gov (United States)

Riley, Thomas E.; Raaijmakers, Geert; Watts, Anna L.

2018-04-01

Constraining the equation of state of cold dense matter in compact stars is a major science goal for observing programmes being conducted using X-ray, radio, and gravitational wave telescopes. We discuss Bayesian hierarchical inference of parametrised dense matter equations of state. In particular we generalise and examine two inference paradigms from the literature: (i) direct posterior equation of state parameter estimation, conditioned on observations of a set of rotating compact stars; and (ii) indirect parameter estimation, via transformation of an intermediary joint posterior distribution of exterior spacetime parameters (such as gravitational masses and coordinate equatorial radii). We conclude that the former paradigm is not only tractable for large-scale analyses, but is principled and flexible from a Bayesian perspective whilst the latter paradigm is not. The thematic problem of Bayesian prior definition emerges as the crux of the difference between these paradigms. The second paradigm should in general only be considered as an ill-defined approach to the problem of utilising archival posterior constraints on exterior spacetime parameters; we advocate for an alternative approach whereby such information is repurposed as an approximative likelihood function. We also discuss why conditioning on a piecewise-polytropic equation of state model - currently standard in the field of dense matter study - can easily violate conditions required for transformation of a probability density distribution between spaces of exterior (spacetime) and interior (source matter) parameters.

Equation of State Project Overview

Energy Technology Data Exchange (ETDEWEB)

Crockett, Scott [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2015-09-11

A general overview of the Equation of State (EOS) Project will be presented. The goal is to provide the audience with an introduction of what our more advanced methods entail (DFT, QMD, etc.. ) and how these models are being utilized to better constrain the thermodynamic models. These models substantially reduce our regions of interpolation between the various thermodynamic limits. I will also present a variety example of recent EOS work.

Validation of the activity expansion method with ultrahigh pressure shock equations of state

International Nuclear Information System (INIS)

Rogers, F.J.; Young, D.A.

1997-01-01

Laser shock experiments have recently been used to measure the equation of state (EOS) of matter in the ultrahigh pressure region between condensed matter and a weakly coupled plasma. Some ultrahigh pressure data from nuclear-generated shocks are also available. Matter at these conditions has proven very difficult to treat theoretically. The many-body activity expansion method (ACTEX) has been used for some time to calculate EOS and opacity data in this region, for use in modeling inertial confinement fusion and stellar interior plasmas. In the present work, we carry out a detailed comparison with the available experimental data in order to validate the method. The agreement is good, showing that ACTEX adequately describes strongly shocked matter. copyright 1997 The American Physical Society

Validation of the activity expansion method with ultrahigh pressure shock equations of state

Energy Technology Data Exchange (ETDEWEB)

Rogers, F.J.; Young, D.A. [Physics Department, Lawrence Livermore National Laboratory, P.O. Box 808, Livermore, California 94550 (United States)

1997-11-01

Laser shock experiments have recently been used to measure the equation of state (EOS) of matter in the ultrahigh pressure region between condensed matter and a weakly coupled plasma. Some ultrahigh pressure data from nuclear-generated shocks are also available. Matter at these conditions has proven very difficult to treat theoretically. The many-body activity expansion method (ACTEX) has been used for some time to calculate EOS and opacity data in this region, for use in modeling inertial confinement fusion and stellar interior plasmas. In the present work, we carry out a detailed comparison with the available experimental data in order to validate the method. The agreement is good, showing that ACTEX adequately describes strongly shocked matter. {copyright} {ital 1997} {ital The American Physical Society}

SOLVING NONLINEAR KLEIN-GORDON EQUATION WITH A QUADRATIC NONLINEAR TERM USING HOMOTOPY ANALYSIS METHOD

Directory of Open Access Journals (Sweden)

H. Jafari

2010-07-01

Full Text Available In this paper, nonlinear Klein-Gordon equation with quadratic term is solved by means of an analytic technique, namely the Homotopy analysis method (HAM.Comparisons are made between the Adomian decomposition method (ADM, the exact solution and homotopy analysis method. The results reveal that the proposed method is very effective and simple.

Equation of state in the presence of gravity

Science.gov (United States)

Kim, Hyeong-Chan; Kang, Gungwon

2016-11-01

We investigate how an equation of state for matter is affected when a gravity is present. For this purpose, we consider a box of ideal gas in the presence of Newtonian gravity. In addition to the ordinary thermodynamic quantities, a characteristic variable that represents a weight per unit area relative to the average pressure is required in order to describe a macroscopic state of the gas. Although the density and the pressure are not uniform due to the presence of gravity, the ideal gas law itself is satisfied for the thermodynamic quantities when averaged over the system. Assuming that the system follows an adiabatic process further, we obtain a new relation between the averaged pressure and density, which differs from the conventional equation of state for the ideal gas in the absence of gravity. Applying our results to a small volume in a Newtonian star, however, we find that the conventional one is reliable for most astrophysical situations when the characteristic scale is small. On the other hand, gravity effects become significant near the surface of a Newtonian star.

Mathematical Methods for Engineers and Scientists 3 Fourier Analysis, Partial Differential Equations and Variational Methods

CERN Document Server

Tang, Kwong-Tin

2007-01-01

Pedagogical insights gained through 30 years of teaching applied mathematics led the author to write this set of student oriented books. Topics such as complex analysis, matrix theory, vector and tensor analysis, Fourier analysis, integral transforms, ordinary and partial differential equations are presented in a discursive style that is readable and easy to follow. Numerous clearly stated, completely worked out examples together with carefully selected problem sets with answers are used to enhance students' understanding and manipulative skill. The goal is to make students comfortable and confident in using advanced mathematical tools in junior, senior, and beginning graduate courses.

Benchmarking and application of the state-of-the-art uncertainty analysis methods XSUSA and SHARK-X

International Nuclear Information System (INIS)

Aures, A.; Bostelmann, F.; Hursin, M.; Leray, O.

2017-01-01

Highlights: • Application of the uncertainty analysis methods XSUSA and SHARK-X. • Propagation of nuclear data uncertainty through PWR pin cell depletion calculation. • Uncertainty quantification of eigenvalue, nuclide densities and Doppler coefficient. • Top contributor to overall output uncertainty by sensitivity analysis. • Comparison with SAMPLER and TSUNAMI of the SCALE code package. - Abstract: This study presents collaborative work performed between GRS and PSI on benchmarking and application of the state-of-the-art uncertainty analysis methods XSUSA and SHARK-X. Applied to a PWR pin cell depletion calculation, both methods propagate input uncertainty from nuclear data to output uncertainty. The uncertainty of the multiplication factors, nuclide densities, and fuel temperature coefficients derived by both methods are compared at various burnup steps. Comparisons of these quantities are furthermore performed with the SAMPLER module of SCALE 6.2. The perturbation theory based TSUNAMI module of both SCALE 6.1 and SCALE 6.2 is additionally applied for comparisons of the reactivity coefficient.

Approximate method for stochastic chemical kinetics with two-time scales by chemical Langevin equations

International Nuclear Information System (INIS)

Wu, Fuke; Tian, Tianhai; Rawlings, James B.; Yin, George

2016-01-01

The frequently used reduction technique is based on the chemical master equation for stochastic chemical kinetics with two-time scales, which yields the modified stochastic simulation algorithm (SSA). For the chemical reaction processes involving a large number of molecular species and reactions, the collection of slow reactions may still include a large number of molecular species and reactions. Consequently, the SSA is still computationally expensive. Because the chemical Langevin equations (CLEs) can effectively work for a large number of molecular species and reactions, this paper develops a reduction method based on the CLE by the stochastic averaging principle developed in the work of Khasminskii and Yin [SIAM J. Appl. Math. 56, 1766–1793 (1996); ibid. 56, 1794–1819 (1996)] to average out the fast-reacting variables. This reduction method leads to a limit averaging system, which is an approximation of the slow reactions. Because in the stochastic chemical kinetics, the CLE is seen as the approximation of the SSA, the limit averaging system can be treated as the approximation of the slow reactions. As an application, we examine the reduction of computation complexity for the gene regulatory networks with two-time scales driven by intrinsic noise. For linear and nonlinear protein production functions, the simulations show that the sample average (expectation) of the limit averaging system is close to that of the slow-reaction process based on the SSA. It demonstrates that the limit averaging system is an efficient approximation of the slow-reaction process in the sense of the weak convergence.

An Evaluation of Kernel Equating: Parallel Equating with Classical Methods in the SAT Subject Tests[TM] Program. Research Report. ETS RR-09-06

Science.gov (United States)

Grant, Mary C.; Zhang, Lilly; Damiano, Michele

2009-01-01

This study investigated kernel equating methods by comparing these methods to operational equatings for two tests in the SAT Subject Tests[TM] program. GENASYS (ETS, 2007) was used for all equating methods and scaled score kernel equating results were compared to Tucker, Levine observed score, chained linear, and chained equipercentile equating…

A robust and accurate numerical method for transcritical turbulent flows at supercritical pressure with an arbitrary equation of state

International Nuclear Information System (INIS)

Kawai, Soshi; Terashima, Hiroshi; Negishi, Hideyo

2015-01-01

This paper addresses issues in high-fidelity numerical simulations of transcritical turbulent flows at supercritical pressure. The proposed strategy builds on a tabulated look-up table method based on REFPROP database for an accurate estimation of non-linear behaviors of thermodynamic and fluid transport properties at the transcritical conditions. Based on the look-up table method we propose a numerical method that satisfies high-order spatial accuracy, spurious-oscillation-free property, and capability of capturing the abrupt variation in thermodynamic properties across the transcritical contact surface. The method introduces artificial mass diffusivity to the continuity and momentum equations in a physically-consistent manner in order to capture the steep transcritical thermodynamic variations robustly while maintaining spurious-oscillation-free property in the velocity field. The pressure evolution equation is derived from the full compressible Navier–Stokes equations and solved instead of solving the total energy equation to achieve the spurious pressure oscillation free property with an arbitrary equation of state including the present look-up table method. Flow problems with and without physical diffusion are employed for the numerical tests to validate the robustness, accuracy, and consistency of the proposed approach.

Optical Bloch equations with multiply connected states

International Nuclear Information System (INIS)

Stacey, D N; Lucas, D M; Allcock, D T C; Szwer, D J; Webster, S C

2008-01-01

The optical Bloch equations, which give the time evolution of the elements of the density matrix of an atomic system subject to radiation, are generalized so that they can be applied when transitions between pairs of states can proceed by more than one stimulated route. The case considered is that for which the time scale of interest in the problem is long compared with that set by the differences in detuning of the radiation fields stimulating via the different routes. It is shown that the Bloch equations then reduce to the standard form of linear differential equations with constant coefficients. The theory is applied to a two-state system driven by two lasers with different intensities and frequencies and to a three-state Λ-system with one laser driving one transition and two driving the second. It is also shown that the theory reproduces well the observed response of a cold 40 Ca + ion when subject to a single laser frequency driving the 4S 1/2 -4P 1/2 transition and a laser with two strong sidebands driving 3D 3/2 -4P 1/2

xRage Equation of State

Energy Technology Data Exchange (ETDEWEB)

Grove, John W. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2016-08-16

The xRage code supports a variety of hydrodynamic equation of state (EOS) models. In practice these are generally accessed in the executing code via a pressure-temperature based table look up. This document will describe the various models supported by these codes and provide details on the algorithms used to evaluate the equation of state.

Scale calculus and the Schroedinger equation

International Nuclear Information System (INIS)

Cresson, Jacky

2003-01-01

This paper is twofold. In a first part, we extend the classical differential calculus to continuous nondifferentiable functions by developing the notion of scale calculus. The scale calculus is based on a new approach of continuous nondifferentiable functions by constructing a one parameter family of differentiable functions f(t,ε) such that f(t,ε)→f(t) when ε goes to zero. This led to several new notions as representations: fractal functions and ε-differentiability. The basic objects of the scale calculus are left and right quantum operators and the scale operator which generalizes the classical derivative. We then discuss some algebraic properties of these operators. We define a natural bialgebra, called quantum bialgebra, associated with them. Finally, we discuss a convenient geometric object associated with our study. In a second part, we define a first quantization procedure of classical mechanics following the scale relativity theory developed by Nottale. We obtain a nonlinear Schroedinger equation via the classical Newton's equation of dynamics using the scale operator. Under special assumptions we recover the classical Schroedinger equation and we discuss the relevance of these assumptions

Investigation of two and three parameter equations of state for cryogenic fluids

International Nuclear Information System (INIS)

Jenkins, S.L.; Majumdar, A.K.; Hendricks, R.C.

1990-01-01

Two-phase flows are a common occurrence in cryogenic engines and an accurate evaluation of the heat-transfer coefficient in two-phase flow is of significant importance in their analysis and design. The thermodynamic equation of state plays a key role in calculating the heat transfer coefficient which is a function of thermodynamic and thermophysical properties. An investigation has been performed to study the performance of two- and three-parameter equations of state to calculate the compressibility factor of cryogenic fluids along the saturation loci. The two-parameter equations considered here are van der Waals and Redlich-Kwong equations of state. The three-parameter equation represented here is the generalized Benedict-Webb-Rubin (BWR) equation of Lee and Kesler. Results have been compared with the modified BWR equation of Bender and the extended BWR equations of Stewart. Seven cryogenic fluids have been tested; oxygen, hydrogen, helium, nitrogen, argon, neon, and air. The performance of the generalized BWR equation is poor for hydrogen and helium. The van der Waals equation is found to be inaccurate for air near the critical point. For helium, all three equations of state become inaccurate near the critical point. 13 refs

Nonequilibrium Equation of State in Suspensions of Active Colloids

Directory of Open Access Journals (Sweden)

Félix Ginot

2015-01-01

Full Text Available Active colloids constitute a novel class of materials composed of colloidal-scale particles locally converting chemical energy into motility, mimicking micro-organisms. Evolving far from equilibrium, these systems display structural organizations and dynamical properties distinct from thermalized colloidal assemblies. Harvesting the potential of this new class of systems requires the development of a conceptual framework to describe these intrinsically nonequilibrium systems. We use sedimentation experiments to probe the nonequilibrium equation of state of a bidimensional assembly of active Janus microspheres and conduct computer simulations of a model of self-propelled hard disks. Self-propulsion profoundly affects the equation of state, but these changes can be rationalized using equilibrium concepts. We show that active colloids behave, in the dilute limit, as an ideal gas with an activity-dependent effective temperature. At finite density, increasing the activity is similar to increasing adhesion between equilibrium particles. We quantify this effective adhesion and obtain a unique scaling law relating activity and effective adhesion in both experiments and simulations. Our results provide a new and efficient way to understand the emergence of novel phases of matter in active colloidal suspensions.

Solutions of Heat-Like and Wave-Like Equations with Variable Coefficients by Means of the Homotopy Analysis Method

International Nuclear Information System (INIS)

Alomari, A. K.; Noorani, M. S. M.; Nazar, R.

2008-01-01

We employ the homotopy analysis method (HAM) to obtain approximate analytical solutions to the heat-like and wave-like equations. The HAM contains the auxiliary parameter ħ, which provides a convenient way of controlling the convergence region of series solutions. The analysis is accompanied by several linear and nonlinear heat-like and wave-like equations with initial boundary value problems. The results obtained prove that HAM is very effective and simple with less error than the Adomian decomposition method and the variational iteration method

Analysis and development of adjoint-based h-adaptive direct discontinuous Galerkin method for the compressible Navier-Stokes equations

Science.gov (United States)

Cheng, Jian; Yue, Huiqiang; Yu, Shengjiao; Liu, Tiegang

2018-06-01

In this paper, an adjoint-based high-order h-adaptive direct discontinuous Galerkin method is developed and analyzed for the two dimensional steady state compressible Navier-Stokes equations. Particular emphasis is devoted to the analysis of the adjoint consistency for three different direct discontinuous Galerkin discretizations: including the original direct discontinuous Galerkin method (DDG), the direct discontinuous Galerkin method with interface correction (DDG(IC)) and the symmetric direct discontinuous Galerkin method (SDDG). Theoretical analysis shows the extra interface correction term adopted in the DDG(IC) method and the SDDG method plays a key role in preserving the adjoint consistency. To be specific, for the model problem considered in this work, we prove that the original DDG method is not adjoint consistent, while the DDG(IC) method and the SDDG method can be adjoint consistent with appropriate treatment of boundary conditions and correct modifications towards the underlying output functionals. The performance of those three DDG methods is carefully investigated and evaluated through typical test cases. Based on the theoretical analysis, an adjoint-based h-adaptive DDG(IC) method is further developed and evaluated, numerical experiment shows its potential in the applications of adjoint-based adaptation for simulating compressible flows.

Dual-scale Galerkin methods for Darcy flow

Science.gov (United States)

Wang, Guoyin; Scovazzi, Guglielmo; Nouveau, Léo; Kees, Christopher E.; Rossi, Simone; Colomés, Oriol; Main, Alex

2018-02-01

The discontinuous Galerkin (DG) method has found widespread application in elliptic problems with rough coefficients, of which the Darcy flow equations are a prototypical example. One of the long-standing issues of DG approximations is the overall computational cost, and many different strategies have been proposed, such as the variational multiscale DG method, the hybridizable DG method, the multiscale DG method, the embedded DG method, and the Enriched Galerkin method. In this work, we propose a mixed dual-scale Galerkin method, in which the degrees-of-freedom of a less computationally expensive coarse-scale approximation are linked to the degrees-of-freedom of a base DG approximation. We show that the proposed approach has always similar or improved accuracy with respect to the base DG method, with a considerable reduction in computational cost. For the specific definition of the coarse-scale space, we consider Raviart-Thomas finite elements for the mass flux and piecewise-linear continuous finite elements for the pressure. We provide a complete analysis of stability and convergence of the proposed method, in addition to a study on its conservation and consistency properties. We also present a battery of numerical tests to verify the results of the analysis, and evaluate a number of possible variations, such as using piecewise-linear continuous finite elements for the coarse-scale mass fluxes.

Information Equation of State

Directory of Open Access Journals (Sweden)

M. Paul Gough

2008-07-01

Full Text Available Landauer’s principle is applied to information in the universe. Once stars began forming there was a constant information energy density as the increasing proportion of matter at high stellar temperatures exactly compensated for the expanding universe. The information equation of state was close to the dark energy value, w = -1, for a wide range of redshifts, 10 > z > 0.8, over one half of cosmic time. A reasonable universe information bit content of only 1087 bits is sufficient for information energy to account for all dark energy. A time varying equation of state with a direct link between dark energy and matter, and linked to star formation in particular, is clearly relevant to the cosmic coincidence problem. In answering the ‘Why now?’ question we wonder ‘What next?’ as we expect the information equation of state to tend towards w = 0 in the future.c

Particle number fluctuations for the van der Waals equation of state

International Nuclear Information System (INIS)

Vovchenko, V; Anchishkin, D V; Gorenstein, M I

2015-01-01

The van der Waals (VDW) equation of state describes a thermal equilibrium in system of particles, where both repulsive and attractive interactions between them are included. This equation predicts the existence of the first order liquid–gas phase transition and the critical point. The standard form of the VDW equation is given by the pressure function in a canonical ensemble (CE) with a fixed number of particles. In this paper the VDW equation is derived within the grand canonical ensemble (GCE) formulation. We argue that this procedure can be useful for new physical applications, in particular, the fluctuations of the number of particles, which are absent in the CE, can be studied in the GCE. For the VDW equation of state in the GCE the particle number fluctuations are calculated for the whole phase diagram, both outside and inside the liquid–gas mixed phase region. It is shown that the scaled variance of these fluctuations remains finite within the mixed phase and goes to infinity at the critical point. The GCE formulation of the VDW equation of state can also be an important step for its application in the statistical description of hadronic systems, where numbers of different particle species are usually not conserved. (paper)

Extreme compression behaviour of equations of state

International Nuclear Information System (INIS)

Shanker, J.; Dulari, P.; Singh, P.K.

2009-01-01

The extreme compression (P→∞) behaviour of various equations of state with K' ∞ >0 yields (P/K) ∞ =1/K' ∞ , an algebraic identity found by Stacey. Here P is the pressure, K the bulk modulus, K ' =dK/dP, and K' ∞ , the value of K ' at P→∞. We use this result to demonstrate further that there exists an algebraic identity also between the higher pressure derivatives of bulk modulus which is satisfied at extreme compression by different types of equations of state such as the Birch-Murnaghan equation, Poirier-Tarantola logarithmic equation, generalized Rydberg equation, Keane's equation and the Stacey reciprocal K-primed equation. The identity has been used to find a relationship between λ ∞ , the third-order Grueneisen parameter at P→∞, and pressure derivatives of bulk modulus with the help of the free-volume formulation without assuming any specific form of equation of state.

Equation of limiting plasticity of the metal upon complex stress state

International Nuclear Information System (INIS)

Tin'gaev, A.K.

2002-01-01

A method for evaluation of the limiting plasticity of the metal in the zones of complex 3D stress state is presented. An analytic equation is derived for limiting plasticity. Parameters of the equation are expresses through the standard characteristics of the mechanical properties determined upon static tension of the smooth sample. Introduced into the obtained analytical equation is a universal fracture constant which indirectly characterizes the state of the material from the point of view of its capacity for elastic overstrain relaxation in the form of plastic flow or fracture. The new equation makes it possible to estimate the limiting plasticity of the metal in a state of complex stress on the basis of traditional characteristics of mechanical properties, which are not difficult to determine [ru

The limiting temperature of hot nuclei from microscopic equation of state

International Nuclear Information System (INIS)

Baldo, M.; Ferreira, L.S.; Nicotra, O.E.

2004-01-01

The limiting temperature T lim of a series of nuclei is calculated employing a set of microscopic nuclear equations of state (EOS's). It is shown that the value of T lim is sensitive to the nuclear matter equation of state used. Comparison with the values extracted in recent phenomenological analysis appears to favor a definite selection of EOS's. On the basis of this phenomenological analysis, it therefore seems possible to check the microscopic calculations of the nuclear EOS at finite temperature, which is hardly accessible through other experimental information

An equation of state for purely kinetic k-essence inspired by cosmic topological defects

Energy Technology Data Exchange (ETDEWEB)

Cordero, Ruben; Gonzalez, Eduardo L.; Queijeiro, Alfonso [Instituto Politecnico Nacional, Departamento de Fisica, Escuela Superior de Fisica y Matematicas, Ciudad de Mexico (Mexico)

2017-06-15

We investigate the physical properties of a purely kinetic k-essence model with an equation of state motivated in superconducting membranes. We compute the equation of state parameter w and discuss its physical evolution via a nonlinear equation of state. Using the adiabatic speed of sound and energy density, we restrict the range of parameters of the model in order to have an acceptable physical behavior. We study the evolution of the scale factor and address the question of the possible existence of finite-time future singularities. Furthermore, we analyze the evolution of the luminosity distance d{sub L} with redshift z by comparing (normalizing) it with the ΛCDM model. Since the equation of state parameter is z-dependent the evolution of the luminosity distance is also analyzed using the Alcock-Paczynski test. (orig.)

State Equation Determination of Cow Dung Biogas

Science.gov (United States)

Marzuki, A.; Wicaksono, L. B.

2017-08-01

A state function is a thermodynamic function which relates various macroscopically measurable properties of a system (state variable) describing the state of matter under a given set of physical conditions. A good understanding of a biogas state function plays a very important role in an effort to maximize biogas processes and to help predicting combation performance. This paper presents a step by step process of an experimental study aimed at determining the equation of state of cow dung biogas. The equation was derived from the data obtained from the experimental results of compressibility (κ) and expansivity (β) following the general form of gas state equation dV = βdT + κdP. In this equation, dV is gas volume variation, dT is temperature variation, and dP is pressure variation. From these results, we formulated a unique state equation from which the biogas critical temperature (Tc) and critical pressure were then determined (Tc = 266.7 K, Pc = 5096647.5 Pa).

Workshop on Numerical Methods for Ordinary Differential Equations

CERN Document Server

Gear, Charles; Russo, Elvira

1989-01-01

Developments in numerical initial value ode methods were the focal topic of the meeting at L'Aquila which explord the connections between the classical background and new research areas such as differental-algebraic equations, delay integral and integro-differential equations, stability properties, continuous extensions (interpolants for Runge-Kutta methods and their applications, effective stepsize control, parallel algorithms for small- and large-scale parallel architectures). The resulting proceedings address many of these topics in both research and survey papers.

On petroleum fluid characterization with the PC-SAFT equation of state

DEFF Research Database (Denmark)

Liang, Xiaodong; Yan, Wei; Thomsen, Kaj

2014-01-01

The perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state has shown promising results for describing complex phase behaviors and high pressure properties of various systems. It has been proposed as an alternative to the classical cubic equations of state in the petroleum...... industry. It is, however, far from a simple task to develop a sophisticated oil characterization method for the PC-SAFT EOS. In this work, in order to answer some fundamental questions of developing new characterization methods for PC-SAFT, six methods are proposed to estimate the model parameters...

Study on TVD parameters sensitivity of a crankshaft using multiple scale and state space method considering quadratic and cubic non-linearities

Directory of Open Access Journals (Sweden)

R. Talebitooti

Full Text Available In this paper the effect of quadratic and cubic non-linearities of the system consisting of the crankshaft and torsional vibration damper (TVD is taken into account. TVD consists of non-linear elastomer material used for controlling the torsional vibration of crankshaft. The method of multiple scales is used to solve the governing equations of the system. Meanwhile, the frequency response of the system for both harmonic and sub-harmonic resonances is extracted. In addition, the effects of detuning parameters and other dimensionless parameters for a case of harmonic resonance are investigated. Moreover, the external forces including both inertia and gas forces are simultaneously applied into the model. Finally, in order to study the effectiveness of the parameters, the dimensionless governing equations of the system are solved, considering the state space method. Then, the effects of the torsional damper as well as all corresponding parameters of the system are discussed.

Large-scale Granger causality analysis on resting-state functional MRI

Science.gov (United States)

D'Souza, Adora M.; Abidin, Anas Zainul; Leistritz, Lutz; Wismüller, Axel

2016-03-01

We demonstrate an approach to measure the information flow between each pair of time series in resting-state functional MRI (fMRI) data of the human brain and subsequently recover its underlying network structure. By integrating dimensionality reduction into predictive time series modeling, large-scale Granger Causality (lsGC) analysis method can reveal directed information flow suggestive of causal influence at an individual voxel level, unlike other multivariate approaches. This method quantifies the influence each voxel time series has on every other voxel time series in a multivariate sense and hence contains information about the underlying dynamics of the whole system, which can be used to reveal functionally connected networks within the brain. To identify such networks, we perform non-metric network clustering, such as accomplished by the Louvain method. We demonstrate the effectiveness of our approach to recover the motor and visual cortex from resting state human brain fMRI data and compare it with the network recovered from a visuomotor stimulation experiment, where the similarity is measured by the Dice Coefficient (DC). The best DC obtained was 0.59 implying a strong agreement between the two networks. In addition, we thoroughly study the effect of dimensionality reduction in lsGC analysis on network recovery. We conclude that our approach is capable of detecting causal influence between time series in a multivariate sense, which can be used to segment functionally connected networks in the resting-state fMRI.

Equation-of-State Modeling of Phase Equilibria in Petroleum Fluids

DEFF Research Database (Denmark)

Jørgensen, Marianne

1996-01-01

The Soave-Redlich-Kwong (SRK) equation of state was used to investigate and develop several aspects of the modeling of natural petroleum fluids.A new method was presented for numerical evaluation of PVT experiments. This method was used in the estimation of binary interaction parameters. A comphr......The Soave-Redlich-Kwong (SRK) equation of state was used to investigate and develop several aspects of the modeling of natural petroleum fluids.A new method was presented for numerical evaluation of PVT experiments. This method was used in the estimation of binary interaction parameters....... A comphrensive study of pseudoization procedures is presented. It is concluded that the compared methods exhibit results of comparable accuracy, and that six to eight pseudocomponents are needed for optimal representation of petroleum fluids.Finally, it is investigated how well the EOS can represent the VLLE...

Numerical analysis for multi-group neutron-diffusion equation using Radial Point Interpolation Method (RPIM)

International Nuclear Information System (INIS)

Kim, Kyung-O; Jeong, Hae Sun; Jo, Daeseong

2017-01-01

Highlights: • Employing the Radial Point Interpolation Method (RPIM) in numerical analysis of multi-group neutron-diffusion equation. • Establishing mathematical formation of modified multi-group neutron-diffusion equation by RPIM. • Performing the numerical analysis for 2D critical problem. - Abstract: A mesh-free method is introduced to overcome the drawbacks (e.g., mesh generation and connectivity definition between the meshes) of mesh-based (nodal) methods such as the finite-element method and finite-difference method. In particular, the Point Interpolation Method (PIM) using a radial basis function is employed in the numerical analysis for the multi-group neutron-diffusion equation. The benchmark calculations are performed for the 2D homogeneous and heterogeneous problems, and the Multiquadrics (MQ) and Gaussian (EXP) functions are employed to analyze the effect of the radial basis function on the numerical solution. Additionally, the effect of the dimensionless shape parameter in those functions on the calculation accuracy is evaluated. According to the results, the radial PIM (RPIM) can provide a highly accurate solution for the multiplication eigenvalue and the neutron flux distribution, and the numerical solution with the MQ radial basis function exhibits the stable accuracy with respect to the reference solutions compared with the other solution. The dimensionless shape parameter directly affects the calculation accuracy and computing time. Values between 1.87 and 3.0 for the benchmark problems considered in this study lead to the most accurate solution. The difference between the analytical and numerical results for the neutron flux is significantly increased in the edge of the problem geometry, even though the maximum difference is lower than 4%. This phenomenon seems to arise from the derivative boundary condition at (x,0) and (0,y) positions, and it may be necessary to introduce additional strategy (e.g., the method using fictitious points and

Nodal spectrum method for solving neutron diffusion equation

International Nuclear Information System (INIS)

Sanchez, D.; Garcia, C. R.; Barros, R. C. de; Milian, D.E.

1999-01-01

Presented here is a new numerical nodal method for solving static multidimensional neutron diffusion equation in rectangular geometry. Our method is based on a spectral analysis of the nodal diffusion equations. These equations are obtained by integrating the diffusion equation in X, Y directions and then considering flat approximations for the current. These flat approximations are the only approximations that are considered in this method, as a result the numerical solutions are completely free from truncation errors. We show numerical results to illustrate the methods accuracy for coarse mesh calculations

Navier-Stokes equations by the finite element method

International Nuclear Information System (INIS)

Portella, P.E.

1984-01-01

A computer program to solve the Navier-Stokes equations by using the Finite Element Method is implemented. The solutions variables investigated are stream-function/vorticity in the steady case and velocity/pressure in the steady state and transient cases. For steady state flow the equations are solved simultaneously by the Newton-Raphson method. For the time dependent formulation, a fractional step method is employed to discretize in time and artificial viscosity is used to preclude spurious oscilations in the solution. The element used is the three node triangle. Some numerical examples are presented and comparisons are made with applications already existent. (Author) [pt

Comparison of Kernel Equating and Item Response Theory Equating Methods

Science.gov (United States)

Meng, Yu

2012-01-01

The kernel method of test equating is a unified approach to test equating with some advantages over traditional equating methods. Therefore, it is important to evaluate in a comprehensive way the usefulness and appropriateness of the Kernel equating (KE) method, as well as its advantages and disadvantages compared with several popular item…

Multicomponent equations of state for electrolytes

DEFF Research Database (Denmark)

Lin, Yi; Thomsen, Kaj; Hemptinne, Jean-Charles de

2007-01-01

. The parameters in the equations of state were fitted to experimental data consisting of apparent molar volumes, osmotic coefficients, mean ionic activity coefficients, and solid-liquid equilibrium data. The results of the parameter fitting are presented. The ability of the equations of state to reproduce...

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

International Nuclear Information System (INIS)

Clouet, J.F.; Samba, G.

2005-01-01

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

Functional analytic methods in complex analysis and applications to partial differential equations

International Nuclear Information System (INIS)

Mshimba, A.S.A.; Tutschke, W.

1990-01-01

The volume contains 24 lectures given at the Workshop on Functional Analytic Methods in Complex Analysis and Applications to Partial Differential Equations held in Trieste, Italy, between 8-19 February 1988, at the ICTP. A separate abstract was prepared for each of these lectures. Refs and figs

Finite element analysis of multi-material models using a balancing domain decomposition method combined with the diagonal scaling preconditioner

International Nuclear Information System (INIS)

Ogino, Masao

2016-01-01

Actual problems in science and industrial applications are modeled by multi-materials and large-scale unstructured mesh, and the finite element analysis has been widely used to solve such problems on the parallel computer. However, for large-scale problems, the iterative methods for linear finite element equations suffer from slow or no convergence. Therefore, numerical methods having both robust convergence and scalable parallel efficiency are in great demand. The domain decomposition method is well known as an iterative substructuring method, and is an efficient approach for parallel finite element methods. Moreover, the balancing preconditioner achieves robust convergence. However, in case of problems consisting of very different materials, the convergence becomes bad. There are some research to solve this issue, however not suitable for cases of complex shape and composite materials. In this study, to improve convergence of the balancing preconditioner for multi-materials, a balancing preconditioner combined with the diagonal scaling preconditioner, called Scaled-BDD method, is proposed. Some numerical results are included which indicate that the proposed method has robust convergence for the number of subdomains and shows high performances compared with the original balancing preconditioner. (author)

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

Energy Technology Data Exchange (ETDEWEB)

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

2015-04-16

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

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

An Equation of State for the Thermodynamic Properties of Cyclohexane

Energy Technology Data Exchange (ETDEWEB)

Zhou, Yong, E-mail: Yong.Zhou2@honeywell.com; Liu, Jun [Honeywell Integrated Technology China Co. Ltd., 430 Li Bing Road, Zhangjiang Hi-Tech Park, Shanghai 201203 (China); Penoncello, Steven G. [Center for Applied Thermodynamic Studies, College of Engineering, University of Idaho, Moscow, Idaho 83844 (United States); Lemmon, Eric W. [Applied Chemicals and Materials Division, National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80305 (United States)

2014-12-15

An equation of state for cyclohexane has been developed using the Helmholtz energy as the fundamental property with independent variables of density and temperature. Multi-property fitting technology was used to fit the equation of state to data for pρT, heat capacities, sound speeds, virial coefficients, vapor pressures, and saturated densities. The equation of state was developed to conform to the Maxwell criteria for two-phase vapor-liquid equilibrium states, and is valid from the triple-point temperature to 700 K, with pressures up to 250 MPa and densities up to 10.3 mol dm{sup −3}. In general, the uncertainties (k = 2, indicating a level of confidence of 95%) in density for the equation of state are 0.1% (liquid and vapor) up to 500 K, and 0.2% above 500 K, with higher uncertainties within the critical region. Between 283 and 473 K with pressures lower than 30 MPa, the uncertainty is as low as 0.03% in density in the liquid phase. The uncertainties in the speed of sound are 0.2% between 283 and 323 K in the liquid, and 1% elsewhere. Other uncertainties are 0.05% in vapor pressure and 2% in heat capacities. The behavior of the equation of state is reasonable within the region of validity and at higher and lower temperatures and pressures. A detailed analysis has been performed in this article.

Nested variant of the method of moments of coupled cluster equations for vertical excitation energies and excited-state potential energy surfaces.

Science.gov (United States)

Kowalski, Karol

2009-05-21

In this article we discuss the problem of proper balancing of the noniterative corrections to the ground- and excited-state energies obtained with approximate coupled cluster (CC) and equation-of-motion CC (EOMCC) approaches. It is demonstrated that for a class of excited states dominated by single excitations and for states with medium doubly excited component, the newly introduced nested variant of the method of moments of CC equations provides mathematically rigorous way of balancing the ground- and excited-state correlation effects. The resulting noniterative methodology accounting for the effect of triples is tested using its parallel implementation on the systems, for which iterative CC/EOMCC calculations with full inclusion of triply excited configurations or their most important subset are numerically feasible.

On Solution of a Fractional Diffusion Equation by Homotopy Transform Method

International Nuclear Information System (INIS)

Salah, A.; Hassan, S.S.A.

2012-01-01

The homotopy analysis transform method (HATM) is applied in this work in order to find the analytical solution of fractional diffusion equations (FDE). These equations are obtained from standard diffusion equations by replacing a second-order space derivative by a fractional derivative of order α and a first order time derivative by a fractional derivative. Furthermore, some examples are given. Numerical results show that the homotopy analysis transform method is easy to implement and accurate when applied to a fractional diffusion equations.

Fundamentals of equations of state

CERN Document Server

Eliezer, Shalom; Hora, Heinrich

2002-01-01

The equation of state was originally developed for ideal gases, and proved central to the development of early molecular and atomic physics. Increasingly sophisticated equations of state have been developed to take into account molecular interactions, quantization, relativistic effects, etc. Extreme conditions of matter are encountered both in nature and in the laboratory, for example in the centres of stars, in relativistic collisions of heavy nuclei, in inertial confinement fusion (where a temperature of 10 9 K and a pressure exceeding a billion atmospheres can be achieved). A sound knowledg

Numerical Analysis of an H1-Galerkin Mixed Finite Element Method for Time Fractional Telegraph Equation

Directory of Open Access Journals (Sweden)

Jinfeng Wang

2014-01-01

Full Text Available We discuss and analyze an H1-Galerkin mixed finite element (H1-GMFE method to look for the numerical solution of time fractional telegraph equation. We introduce an auxiliary variable to reduce the original equation into lower-order coupled equations and then formulate an H1-GMFE scheme with two important variables. We discretize the Caputo time fractional derivatives using the finite difference methods and approximate the spatial direction by applying the H1-GMFE method. Based on the discussion on the theoretical error analysis in L2-norm for the scalar unknown and its gradient in one dimensional case, we obtain the optimal order of convergence in space-time direction. Further, we also derive the optimal error results for the scalar unknown in H1-norm. Moreover, we derive and analyze the stability of H1-GMFE scheme and give the results of a priori error estimates in two- or three-dimensional cases. In order to verify our theoretical analysis, we give some results of numerical calculation by using the Matlab procedure.

Direct Yaw Control of Vehicle using State Dependent Riccati Equation with Integral Terms

Directory of Open Access Journals (Sweden)

SANDHU, F.

2016-05-01

Full Text Available Direct yaw control of four-wheel vehicles using optimal controllers such as the linear quadratic regulator (LQR and the sliding mode controller (SMC either considers only certain parameters constant in the nonlinear equations of vehicle model or totally neglect their effects to obtain simplified models, resulting in loss of states for the system. In this paper, a modified state-dependent Ricatti equation method obtained by the simplification of the vehicle model is proposed. This method overcomes the problem of the lost states by including state integrals. The results of the proposed system are compared with the sliding mode slip controller and state-dependent Ricatti equation method using high fidelity vehicle model in the vehicle simulation software package, Carsim. Results show 38% reduction in the lateral velocity, 34% reduction in roll and 16% reduction in excessive yaw by only increasing the fuel consumption by 6.07%.

On numerical solution of Burgers' equation by homotopy analysis method

International Nuclear Information System (INIS)

Inc, Mustafa

2008-01-01

In this Letter, we present the Homotopy Analysis Method (shortly HAM) for obtaining the numerical solution of the one-dimensional nonlinear Burgers' equation. The initial approximation can be freely chosen with possible unknown constants which can be determined by imposing the boundary and initial conditions. Convergence of the solution and effects for the method is discussed. The comparison of the HAM results with the Homotopy Perturbation Method (HPM) and the results of [E.N. Aksan, Appl. Math. Comput. 174 (2006) 884; S. Kutluay, A. Esen, Int. J. Comput. Math. 81 (2004) 1433; S. Abbasbandy, M.T. Darvishi, Appl. Math. Comput. 163 (2005) 1265] are made. The results reveal that HAM is very simple and effective. The HAM contains the auxiliary parameter h, which provides us with a simple way to adjust and control the convergence region of solution series. The numerical solutions are compared with the known analytical and some numerical solutions

Dark energy cosmology with generalized linear equation of state

International Nuclear Information System (INIS)

Babichev, E; Dokuchaev, V; Eroshenko, Yu

2005-01-01

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

National scale biomass estimators for United States tree species

Science.gov (United States)

Jennifer C. Jenkins; David C. Chojnacky; Linda S. Heath; Richard A. Birdsey

2003-01-01

Estimates of national-scale forest carbon (C) stocks and fluxes are typically based on allometric regression equations developed using dimensional analysis techniques. However, the literature is inconsistent and incomplete with respect to large-scale forest C estimation. We compiled all available diameter-based allometric regression equations for estimating total...

A new method for the solution of the Schroedinger equation

International Nuclear Information System (INIS)

Amore, Paolo; Aranda, Alfredo; De Pace, Arturo

2004-01-01

We present a new method for the solution of the Schroedinger equation applicable to problems of a non-perturbative nature. The method works by identifying three different scales in the problem, which then are treated independently: an asymptotic scale, which depends uniquely on the form of the potential at large distances; an intermediate scale, still characterized by an exponential decay of the wavefunction; and, finally, a short distance scale, in which the wavefunction is sizable. The notion of optimized perturbation is then used in the last two regimes. We apply the method to the quantum anharmonic oscillator and find it suitable to treat both energy eigenvalues and wavefunctions, even for strong couplings

Hybrid massively parallel fast sweeping method for static Hamilton–Jacobi equations

Energy Technology Data Exchange (ETDEWEB)

Detrixhe, Miles, E-mail: mdetrixhe@engineering.ucsb.edu [Department of Mechanical Engineering (United States); University of California Santa Barbara, Santa Barbara, CA, 93106 (United States); Gibou, Frédéric, E-mail: fgibou@engineering.ucsb.edu [Department of Mechanical Engineering (United States); University of California Santa Barbara, Santa Barbara, CA, 93106 (United States); Department of Computer Science (United States); Department of Mathematics (United States)

2016-10-01

The fast sweeping method is a popular algorithm for solving a variety of static Hamilton–Jacobi equations. Fast sweeping algorithms for parallel computing have been developed, but are severely limited. In this work, we present a multilevel, hybrid parallel algorithm that combines the desirable traits of two distinct parallel methods. The fine and coarse grained components of the algorithm take advantage of heterogeneous computer architecture common in high performance computing facilities. We present the algorithm and demonstrate its effectiveness on a set of example problems including optimal control, dynamic games, and seismic wave propagation. We give results for convergence, parallel scaling, and show state-of-the-art speedup values for the fast sweeping method.

Hybrid massively parallel fast sweeping method for static Hamilton–Jacobi equations

International Nuclear Information System (INIS)

Detrixhe, Miles; Gibou, Frédéric

2016-01-01

The fast sweeping method is a popular algorithm for solving a variety of static Hamilton–Jacobi equations. Fast sweeping algorithms for parallel computing have been developed, but are severely limited. In this work, we present a multilevel, hybrid parallel algorithm that combines the desirable traits of two distinct parallel methods. The fine and coarse grained components of the algorithm take advantage of heterogeneous computer architecture common in high performance computing facilities. We present the algorithm and demonstrate its effectiveness on a set of example problems including optimal control, dynamic games, and seismic wave propagation. We give results for convergence, parallel scaling, and show state-of-the-art speedup values for the fast sweeping method.

Homogeneous axisymmetric model with a limitting stiff equation of state

International Nuclear Information System (INIS)

Korkina, M.P.; Martynenko, V.G.

1976-01-01

A solution is obtained for Einstein's equations in which all metric coefficients are time functions for a limiting stiff equation of the substance state. Thr solution describes a homogeneous cosmological model with cylindrical symmetry. It is shown that the same metrics can be induced by a massless scalar only time-dependent field. Analysis of this solution is presented

On the application of finite element method in the solution of steady state diffusion equation

International Nuclear Information System (INIS)

Ono, S.

1982-01-01

The solution of the steady state neutron diffusion equation is obtained by using the finite element method. Specifically the variational approach is used for one dimensional problems and the weighted residual method (Galerkin) for one and two dimensional problems. The spatial domain is divided into retangular elements and the neutron flux is approximated by linear (one dimensional case), and bilinear (two-dimensional case) functions. Numerical results are obtained with a FORTRAN IV computer program and compared with those obtained by the finite difference CITATION code. The results show that linear or bilinear functions, do not satisfactorily describe the differential parameters in highly heterogeneous reactor cases, but provide good results for integral parameters such as multiplication factor. (Author) [pt

Equations-of-motion treatment of pairing correlations: Seniority-one states

International Nuclear Information System (INIS)

Andreozzi, F.; Covello, A.; Gargano, A.; Porrino, A.

1988-01-01

In prior work we have developed an equations-of-motion method for treating seniority-one states in pairing-force theory. Here we present a new and simpler version of that method. Some numerical applications to Sn isotopes show its considerable practical value

Equation of state study of Laser Megajoule capsules ablator materials

International Nuclear Information System (INIS)

Colin-Lalu, Pierre

2016-01-01

This PhD thesis enters the field of inertial confinement fusion studies. In particular, it focuses on the equation of state tables of ablator materials synthesized on LMJ capsules. This work is indeed aims at improving the theoretical models introduced into the equation of state tables. We focused in the Mbar-eV pressure-temperature range because it can be access on kJ-scale laser facilities.In order to achieve this, we used the QEOS model, which is simple to use, configurable, and easily modifiable.First, quantum molecular dynamics (QMD) simulations were performed to generate cold compression curve as well as shock compression curves along the principal Hugoniot. Simulations were compared to QEOS model and showed that atomic bond dissociation has an effect on the compressibility. Results from these simulations are then used to parametrize the Grueneisen parameter in order to generate a tabulated equation of state that includes dissociation. It allowed us to show its influence on shock timing in a hydrodynamic simulation.Second, thermodynamic states along the Hugoniot were measured during three experimental campaigns upon the LULI2000 and GEKKO XII laser facilities. Experimental data confirm QMD simulations.This study was performed on two ablator materials which are an undoped polymer CHO, and a silicon-doped polymer CHOSi. Results showed universal shock compression properties. (author) [fr

Comparison of methods for calculating thermodynamic properties of binary mixtures in the sub and super critical state: Lee-Kesler and cubic equations of state for binary mixtures containing either CO2 or H2S

International Nuclear Information System (INIS)

Yang, Jyisy; Griffiths, Peter R.; Goodwin, Anthony R.H.

2003-01-01

The (ρ,T,p) and (vapor + liquid) equilibria for fluid mixtures containing either CO 2 or H 2 S have been determined from 13 equations of state. The estimated values have been compared with published experimental results. CO 2 and H 2 S were used to represent non-polar and polar fluids, respectively. The equations of state investigated were as follows: (a) the Lee-Kesler equation; (b) two equations that included new reference fluids for the Lee-Kesler method; (c) three so-called extended equations of state; and (d) seven cubic equations of state. After adjustment of the binary interaction parameters the predicted values differed from the experimental data by about 0.8% for CO 2 mixtures while for H 2 S mixtures the uncertainty was about ±2.8%. Somewhat larger errors, although still lower than ±5%, were obtained for co-existing phase densities; the Lee-Kesler method provided results of the highest accuracy. The cubic equations proposed by Schmidt and Wenzel and Valderrama provide the most reliable predictions of both single and co-existing phase densities. Comparison of the predicted (vapor + liquid) equilibrium with experiment shows that each of the seven cubic equations provides results of similar accuracy and all within ±6%

Discontinuous Galerkin Subgrid Finite Element Method for Heterogeneous Brinkman’s Equations

KAUST Repository

Iliev, Oleg P.; Lazarov, Raytcho D.; Willems, Joerg

2010-01-01

We present a two-scale finite element method for solving Brinkman's equations with piece-wise constant coefficients. This system of equations model fluid flows in highly porous, heterogeneous media with complex topology of the heterogeneities. We

Simple functional-differential equations for the bound-state wave-function components

International Nuclear Information System (INIS)

Kamuntavicius, G.P.

1986-01-01

The author presents a new method of a direct derivation of differential equations for the wave-function components of identical-particles systems. The method generates in a simple manner all the possible variants of these equations. In some cases they are the differential equations of Faddeev or Yakubovskii. It is shown that the case of the bound states allows to formulate very simple equations for the components which are equivalent to the Schroedinger equation for the complete wave function. The components with a minimal antisymmetry are defined and the corresponding equations are derived. (Auth.)

Psychometric evaluation of the altered states of consciousness rating scale (OAV.

Directory of Open Access Journals (Sweden)

Erich Studerus

2010-08-01

Full Text Available The OAV questionnaire has been developed to integrate research on altered states of consciousness (ASC. It measures three primary and one secondary dimensions of ASC that are hypothesized to be invariant across ASC induction methods. The OAV rating scale has been in use for more than 20 years and applied internationally in a broad range of research fields, yet its factorial structure has never been tested by structural equation modeling techniques and its psychometric properties have never been examined in large samples of experimentally induced ASC.The present study conducted a psychometric evaluation of the OAV in a sample of psilocybin (n = 327, ketamine (n = 162, and MDMA (n = 102 induced ASC that was obtained by pooling data from 43 experimental studies. The factorial structure was examined by confirmatory factor analysis, exploratory structural equation modeling, hierarchical item clustering (ICLUST, and multiple indicators multiple causes (MIMIC modeling. The originally proposed model did not fit the data well even if zero-constraints on non-target factor loadings and residual correlations were relaxed. Furthermore, ICLUST suggested that the "oceanic boundlessness" and "visionary restructuralization" factors could be combined on a high level of the construct hierarchy. However, because these factors were multidimensional, we extracted and examined 11 new lower order factors. MIMIC modeling indicated that these factors were highly measurement invariant across drugs, settings, questionnaire versions, and sexes. The new factors were also demonstrated to have improved homogeneities, satisfactory reliabilities, discriminant and convergent validities, and to differentiate well among the three drug groups.The original scales of the OAV were shown to be multidimensional constructs. Eleven new lower order scales were constructed and demonstrated to have desirable psychometric properties. The new lower order scales are most likely better suited to

Thermodynamics and equations of state of matter from ideal gas to quark-gluon plasma

CERN Document Server

Fortov, Vladimir E

2016-01-01

The monograph presents a comparative analysis of different thermodynamic models of the equations of state. The basic ideological premises of the theoretical methods and the experiment are considered. The principal attention is on the description of states that are of greatest interest for the physics of high energy concentrations which are either already attained or can be reached in the near future in controlled terrestrial conditions, or are realized in astrophysical objects at different stages of their evolution. Ultra-extreme astrophysical and nuclear-physical applications are also analyzed where the thermodynamics of matter is affected substantially by relativism, high-power gravitational and magnetic fields, thermal radiation, transformation of nuclear particles, nucleon neutronization, and quark deconfinement. The book is intended for a wide range of specialists engaged in the study of the equations of state of matter and high energy density physics, as well as for senior students and postgraduates.

Differential and difference equations a comparison of methods of solution

CERN Document Server

Maximon, Leonard C

2016-01-01

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

KRYSI, Ordinary Differential Equations Solver with Sdirk Krylov Method

International Nuclear Information System (INIS)

Hindmarsh, A.C.; Norsett, S.P.

2001-01-01

1 - Description of program or function: KRYSI is a set of FORTRAN subroutines for solving ordinary differential equations initial value problems. It is suitable for both stiff and non-stiff systems. When solving the implicit stage equations in the stiff case, KRYSI uses a Krylov subspace iteration method called the SPIGMR (Scaled Preconditioned Incomplete Generalized Minimum Residual) method. No explicit Jacobian storage is required, except where used in pre- conditioning. A demonstration problem is included with a description of two pre-conditioners that are natural for its solution by KRYSI. 2 - Method of solution: KRYSI uses a three-stage, third-order singly diagonally implicit Runge-Kutta (SDIRK) method. In the stiff case, a preconditioned Krylov subspace iteration within a (so-called) inexact Newton iteration is used to solve the system of nonlinear algebraic equations

Scaling of the steady state and stability behaviour of single and two-phase natural circulation systems

International Nuclear Information System (INIS)

Vijayan, P.K.; Nayak, A.K.; Bade, M.H.; Kumar, N.; Saha, D.; Sinha, R.K.

2002-01-01

Scaling methods for both single-phase and two-phase natural circulation systems have been presented. For single-phase systems, simulation of the steady state flow can be achieved by preserving just one nondimensional parameter. For uniform diameter two-phase systems also, it is possible to simulate the steady state behaviour with just one non-dimensional parameter. Simulation of the stability behaviour requires geometric similarity in addition to the similarity of the physical parameters appearing in the governing equations. The scaling laws proposed have been tested with experimental data in case of single-phase natural circulation. (author)

The nuclear equation of state

International Nuclear Information System (INIS)

Kahana, S.

1986-01-01

The role of the nuclear equation of state in determining the fate of the collapsing cores of massive stars is examined in light of both recent theoretical advances in this subject and recent experimental measurements with relativistic heavy ions. The difficulties existing in attempts to bring the softer nuclear matter apparently required by the theory of Type II supernovae into consonance with the heavy ion data are discussed. Relativistic mean field theory is introduced as a candidate for derivation of the equation of state, and a simple form for the saturation compressibility is obtained. 28 refs., 4 figs., 1 tab

The nuclear equation of state

Energy Technology Data Exchange (ETDEWEB)

Kahana, S.

1986-01-01

The role of the nuclear equation of state in determining the fate of the collapsing cores of massive stars is examined in light of both recent theoretical advances in this subject and recent experimental measurements with relativistic heavy ions. The difficulties existing in attempts to bring the softer nuclear matter apparently required by the theory of Type II supernovae into consonance with the heavy ion data are discussed. Relativistic mean field theory is introduced as a candidate for derivation of the equation of state, and a simple form for the saturation compressibility is obtained. 28 refs., 4 figs., 1 tab.

Neutron stars, fast pulsars, supernovae and the equation of state of dense matter

International Nuclear Information System (INIS)

Glendening, N.K.

1989-01-01

We discuss the prospects for obtaining constraints on the equation of state from astrophysical sources. Neutron star masses although few are known at present, provide a very direct constraint in as much as the connection to the equation of state involves only the assumption that Einstein's general theory of relativity is correct at the macroscopic scale. If the millisecond pulses briefly observed in the remnant of SN1987A can be attributed to uniform rotation of a pulsar, then a very severe constraint is placed on the equation of state. The theory again is very secure. The precise nature of the constraint is not yet understood, but it appears that the equation of state must be neither too soft nor stiff, and it may be that there is information not only on the stiffness of the equation of state but on its shape. Supernovae simulations involve such a plethora of physical processes including those involved in the evolution of the precollapse configuration, not all of them known or understood, that they provide no constraint at the present time. Not even the broad category of mechanism for the explosion is agreed upon (prompt shock, delayed shock, or nuclear explosion). In connection with very fast pulsars, we include some speculations on pure quark matter stars, and on possible scenarios for understanding the disappearance of the fast pulsar in SN1987A. 47 refs., 16 figs., 1 tab

Neutron stars, fast pulsars, supernovae and the equation of state of dense matter

Energy Technology Data Exchange (ETDEWEB)

Glendening, N.K.

1989-06-01

We discuss the prospects for obtaining constraints on the equation of state from astrophysical sources. Neutron star masses although few are known at present, provide a very direct constraint in as much as the connection to the equation of state involves only the assumption that Einstein's general theory of relativity is correct at the macroscopic scale. If the millisecond pulses briefly observed in the remnant of SN1987A can be attributed to uniform rotation of a pulsar, then a very severe constraint is placed on the equation of state. The theory again is very secure. The precise nature of the constraint is not yet understood, but it appears that the equation of state must be neither too soft nor stiff, and it may be that there is information not only on the stiffness of the equation of state but on its shape. Supernovae simulations involve such a plethora of physical processes including those involved in the evolution of the precollapse configuration, not all of them known or understood, that they provide no constraint at the present time. Not even the broad category of mechanism for the explosion is agreed upon (prompt shock, delayed shock, or nuclear explosion). In connection with very fast pulsars, we include some speculations on pure quark matter stars, and on possible scenarios for understanding the disappearance of the fast pulsar in SN1987A. 47 refs., 16 figs., 1 tab.

A phenomenological one-parameter equation of state for osmotic pressures of PEG and other neutral flexible polymers in good solvents

DEFF Research Database (Denmark)

Cohen, J.A.; Podgornik, R; Hansen, Per Lyngs

2009-01-01

We present a phenomenological one-parameter scaling equation of state that accurately represents osmotic pressures of neutral flexible polymers in good solvents from the dilute through the semidilute regime. The equation comprises a sum of scaled van't Hoff and des Cloizeaux terms including a fit...

Fully Digital Chaotic Differential Equation-based Systems And Methods

KAUST Repository

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

2012-01-01

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

Fully Digital Chaotic Differential Equation-based Systems And Methods

KAUST Repository

Radwan, Ahmed Gomaa Ahmed

2012-09-06

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

Methods for Equating Mental Tests.

Science.gov (United States)

1984-11-01

1983) compared conventional and IRT methods for equating the Test of English as a Foreign Language ( TOEFL ) after chaining. Three conventional and...three IRT equating methods were examined in this study; two sections of TOEFL were each (separately) equated. The IRT methods included the following: (a...group. A separate base form was established for each of the six equating methods. Instead of equating the base-form TOEFL to itself, the last (eighth

Nonparametric reconstruction of the dark energy equation of state

Energy Technology Data Exchange (ETDEWEB)

Heitmann, Katrin [Los Alamos National Laboratory; Holsclaw, Tracy [Los Alamos National Laboratory; Alam, Ujjaini [Los Alamos National Laboratory; Habib, Salman [Los Alamos National Laboratory; Higdon, David [Los Alamos National Laboratory; Sanso, Bruno [UC SANTA CRUZ; Lee, Herbie [UC SANTA CRUZ

2009-01-01

The major aim of ongoing and upcoming cosmological surveys is to unravel the nature of dark energy. In the absence of a compelling theory to test, a natural approach is to first attempt to characterize the nature of dark energy in detail, the hope being that this will lead to clues about the underlying fundamental theory. A major target in this characterization is the determination of the dynamical properties of the dark energy equation of state w. The discovery of a time variation in w(z) could then lead to insights about the dynamical origin of dark energy. This approach requires a robust and bias-free method for reconstructing w(z) from data, which does not rely on restrictive expansion schemes or assumed functional forms for w(z). We present a new non parametric reconstruction method for the dark energy equation of state based on Gaussian Process models. This method reliably captures nontrivial behavior of w(z) and provides controlled error bounds. We demollstrate the power of the method on different sets of simulated supernova data. The GP model approach is very easily extended to include diverse cosmological probes.

Approximate analytical solution of diffusion equation with fractional time derivative using optimal homotopy analysis method

Directory of Open Access Journals (Sweden)

S. Das

2013-12-01

Full Text Available In this article, optimal homotopy-analysis method is used to obtain approximate analytic solution of the time-fractional diffusion equation with a given initial condition. The fractional derivatives are considered in the Caputo sense. Unlike usual Homotopy analysis method, this method contains at the most three convergence control parameters which describe the faster convergence of the solution. Effects of parameters on the convergence of the approximate series solution by minimizing the averaged residual error with the proper choices of parameters are calculated numerically and presented through graphs and tables for different particular cases.

Numerical bifurcation analysis of delay differential equations arising from physiological modeling.

Science.gov (United States)

Engelborghs, K; Lemaire, V; Bélair, J; Roose, D

2001-04-01

This paper has a dual purpose. First, we describe numerical methods for continuation and bifurcation analysis of steady state solutions and periodic solutions of systems of delay differential equations with an arbitrary number of fixed, discrete delays. Second, we demonstrate how these methods can be used to obtain insight into complex biological regulatory systems in which interactions occur with time delays: for this, we consider a system of two equations for the plasma glucose and insulin concentrations in a diabetic patient subject to a system of external assistance. The model has two delays: the technological delay of the external system, and the physiological delay of the patient's liver. We compute stability of the steady state solution as a function of two parameters, compare with analytical results and compute several branches of periodic solutions and their stability. These numerical results allow to infer two categories of diabetic patients for which the external system has different efficiency.

Equations of state for light water

International Nuclear Information System (INIS)

Rubin, G.A.; Granziera, M.R.

1983-01-01

The equations of state for light water were developed, based on the tables of Keenan and Keyes. Equations are presented, describing the specific volume, internal energy, enthalpy and entropy of saturated steam, superheated vapor and subcooled liquid as a function of pressure and temperature. For each property, several equations are shown, with different precisions and different degress of complexity. (Author) [pt

Comparison of scale analysis and numerical simulation for saturated zone convective mixing processes

International Nuclear Information System (INIS)

Oldenburg, C.M.

1998-01-01

Scale analysis can be used to predict a variety of quantities arising from natural systems where processes are described by partial differential equations. For example, scale analysis can be applied to estimate the effectiveness of convective missing on the dilution of contaminants in groundwater. Scale analysis involves substituting simple quotients for partial derivatives and identifying and equating the dominant terms in an order-of-magnitude sense. For free convection due to sidewall heating of saturated porous media, scale analysis shows that vertical convective velocity in the thermal boundary layer region is proportional to the Rayleigh number, horizontal convective velocity is proportional to the square root of the Rayleigh number, and thermal boundary layer thickness is proportional to the inverse square root of the Rayleigh number. These scale analysis estimates are corroborated by numerical simulations of an idealized system. A scale analysis estimate of mixing time for a tracer mixing by hydrodynamic dispersion in a convection cell also agrees well with numerical simulation for two different Rayleigh numbers. Scale analysis for the heating-from-below scenario produces estimates of maximum velocity one-half as large as the sidewall case. At small values of the Rayleigh number, this estimate is confirmed by numerical simulation. For larger Rayleigh numbers, simulation results suggest maximum velocities are similar to the sidewall heating scenario. In general, agreement between scale analysis estimates and numerical simulation results serves to validate the method of scale analysis. Application is to radioactive repositories

Robust periodic steady state analysis of autonomous oscillators based on generalized eigenvalues

NARCIS (Netherlands)

Mirzavand, R.; Maten, ter E.J.W.; Beelen, T.G.J.; Schilders, W.H.A.; Abdipour, A.

2011-01-01

In this paper, we present a new gauge technique for the Newton Raphson method to solve the periodic steady state (PSS) analysis of free-running oscillators in the time domain. To find the frequency a new equation is added to the system of equations. Our equation combines a generalized eigenvector

Robust periodic steady state analysis of autonomous oscillators based on generalized eigenvalues

NARCIS (Netherlands)

Mirzavand, R.; Maten, ter E.J.W.; Beelen, T.G.J.; Schilders, W.H.A.; Abdipour, A.; Michielsen, B.; Poirier, J.R.

2012-01-01

In this paper, we present a new gauge technique for the Newton Raphson method to solve the periodic steady state (PSS) analysis of free-running oscillators in the time domain. To find the frequency a new equation is added to the system of equations. Our equation combines a generalized eigenvector

A simple scaling law for the equation of state and the radial distribution functions calculated by density-functional theory molecular dynamics

Science.gov (United States)

Danel, J.-F.; Kazandjian, L.

2018-06-01

It is shown that the equation of state (EOS) and the radial distribution functions obtained by density-functional theory molecular dynamics (DFT-MD) obey a simple scaling law. At given temperature, the thermodynamic properties and the radial distribution functions given by a DFT-MD simulation remain unchanged if the mole fractions of nuclei of given charge and the average volume per atom remain unchanged. A practical interest of this scaling law is to obtain an EOS table for a fluid from that already obtained for another fluid if it has the right characteristics. Another practical interest of this result is that an asymmetric mixture made up of light and heavy atoms requiring very different time steps can be replaced by a mixture of atoms of equal mass, which facilitates the exploration of the configuration space in a DFT-MD simulation. The scaling law is illustrated by numerical results.

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

International Nuclear Information System (INIS)

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

1986-01-01

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

New method for solving three-dimensional Schroedinger equation

International Nuclear Information System (INIS)

Melezhik, V.S.

1992-01-01

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

Equation of state of warm condensed matter

Energy Technology Data Exchange (ETDEWEB)

Barbee, T.W., III; Young, D.A.; Rogers, F.J.

1998-03-01

Recent advances in computational condensed matter theory have yielded accurate calculations of properties of materials. These calculations have, for the most part, focused on the low temperature (T=0) limit. An accurate determination of the equation of state (EOS) at finite temperature also requires knowledge of the behavior of the electron and ion thermal pressure as a function of T. Current approaches often interpolate between calculated T=0 results and approximations valid in the high T limit. Plasma physics-based approaches are accurate in the high temperature limit, but lose accuracy below T{approximately}T{sub Fermi}. We seek to ``connect up`` these two regimes by using ab initio finite temperature methods (including linear-response[1] based phonon calculations) to derive an equation of state of condensed matter for T{<=}T{sub Fermi}. We will present theoretical results for the principal Hugoniot of shocked materials, including carbon and aluminum, up to pressures P>100 GPa and temperatures T>10{sup 4}K, and compare our results with available experimental data.

Stoichiometric network analysis and associated dimensionless kinetic equations. Application to a model of the Bray-Liebhafsky reaction.

Science.gov (United States)

Schmitz, Guy; Kolar-Anić, Ljiljana Z; Anić, Slobodan R; Cupić, Zeljko D

2008-12-25

The stoichiometric network analysis (SNA) introduced by B. L. Clarke is applied to a simplified model of the complex oscillating Bray-Liebhafsky reaction under batch conditions, which was not examined by this method earlier. This powerful method for the analysis of steady-states stability is also used to transform the classical differential equations into dimensionless equations. This transformation is easy and leads to a form of the equations combining the advantages of classical dimensionless equations with the advantages of the SNA. The used dimensionless parameters have orders of magnitude given by the experimental information about concentrations and currents. This simplifies greatly the study of the slow manifold and shows which parameters are essential for controlling its shape and consequently have an important influence on the trajectories. The effectiveness of these equations is illustrated on two examples: the study of the bifurcations points and a simple sensitivity analysis, different from the classical one, more based on the chemistry of the studied system.

Sufficient Descent Conjugate Gradient Methods for Solving Convex Constrained Nonlinear Monotone Equations

Directory of Open Access Journals (Sweden)

San-Yang Liu

2014-01-01

Full Text Available Two unified frameworks of some sufficient descent conjugate gradient methods are considered. Combined with the hyperplane projection method of Solodov and Svaiter, they are extended to solve convex constrained nonlinear monotone equations. Their global convergence is proven under some mild conditions. Numerical results illustrate that these methods are efficient and can be applied to solve large-scale nonsmooth equations.

Quintom models with an equation of state crossing -1

International Nuclear Information System (INIS)

Zhao Wen; Zhang Yang

2006-01-01

In this paper, we investigate a kind of special quintom model, which is made of a quintessence field φ 1 and a phantom field φ 2 , and the potential function has the form of V(φ 1 2 -φ 2 2 ). This kind of quintom field can be separated into two kinds: the hessence model, which has the state of φ 1 2 >φ 2 2 , and the hantom model with the state φ 1 2 2 2 . We discuss the evolution of these models in the ω-ω ' plane (ω is the state equation of the dark energy, and ω ' is its time derivative in units of Hubble time), and find that according to ω>-1 or ' plane can be divided into four parts. The late time attractor solution, if existing, is always quintessencelike or Λ-like for hessence field, so the big rip does not exist. But for hantom field, its late time attractor solution can be phantomlike or Λ-like, and sometimes, the big rip is unavoidable. Then we consider two special cases: one is the hessence field with an exponential potential, and the other is with a power law potential. We investigate their evolution in the ω-ω ' plane. We also develop a theoretical method of constructing the hessence potential function directly from the effective equation-of-state function ω(z). We apply our method to five kinds of parametrizations of equation-of-state parameter, where ω crossing -1 can exist, and find they all can be realized. At last, we discuss the evolution of the perturbations of the quintom field, and find the perturbations of the quintom δ Q and the metric Φ are all finite even at the state of ω=-1 and ω ' ≠0

Application of the N-quantum approximation method to bound state problems

International Nuclear Information System (INIS)

Raychaudhuri, A.

1977-01-01

The N-quantum approximation (NQA) method is examined in the light of its application to bound state problems. Bound state wave functions are obtained as expansion coefficients in a truncated Haag expansion. From the equations of motion for the Heisenberg field and the NQA expansion, an equation satisfied by the wave function is derived. Two different bound state systems are considered. In one case, the bound state problem of two identical scalars by scalar exchange is analyzed using the NQA. An integral equation satisfied by the wave function is derived. In the nonrelativistic limit, the equation is shown to reduce to the Schroedinger equation. The equation is solved numerically, and the results compared with those obtained for this system by other methods. The NQA method is also applied to the bound state of two spin 1/2 particles with electromagnetic interaction. The integral equation for the wave function is shown to agree with the corresponding Bethe Salpeter equation in the nonrelativistic limit. Using the Dirac (4 x 4) matrices the wave function is expanded in terms of structure functions and the equation for the wave function is reduced to two disjoint sets of coupled equation for the structure functions

Evaluation of Maryland abutment scour equation through selected threshold velocity methods

Science.gov (United States)

Benedict, S.T.

2010-01-01

The U.S. Geological Survey, in cooperation with the Maryland State Highway Administration, used field measurements of scour to evaluate the sensitivity of the Maryland abutment scour equation to the critical (or threshold) velocity variable. Four selected methods for estimating threshold velocity were applied to the Maryland abutment scour equation, and the predicted scour to the field measurements were compared. Results indicated that performance of the Maryland abutment scour equation was sensitive to the threshold velocity with some threshold velocity methods producing better estimates of predicted scour than did others. In addition, results indicated that regional stream characteristics can affect the performance of the Maryland abutment scour equation with moderate-gradient streams performing differently from low-gradient streams. On the basis of the findings of the investigation, guidance for selecting threshold velocity methods for application to the Maryland abutment scour equation are provided, and limitations are noted.

Hot QCD equations of state and relativistic heavy ion collisions

Science.gov (United States)

Chandra, Vinod; Kumar, Ravindra; Ravishankar, V.

2007-11-01

We study two recently proposed equations of state obtained from high-temperature QCD and show how they can be adapted to use them for making predictions for relativistic heavy ion collisions. The method involves extracting equilibrium distribution functions for quarks and gluons from the equation of state (EOS), which in turn will allow a determination of the transport and other bulk properties of the quark gluon-plasma. Simultaneously, the method also yields a quasiparticle description of interacting quarks and gluons. The first EOS is perturbative in the QCD coupling constant and has contributions of O(g5). The second EOS is an improvement over the first, with contributions up to O[g6ln(1/g)]; it incorporates the nonperturbative hard thermal contributions. The interaction effects are shown to be captured entirely by the effective chemical potentials for the gluons and the quarks, in both cases. The chemical potential is seen to be highly sensitive to the EOS. As an application, we determine the screening lengths, which are, indeed, the most important diagnostics for QGP. The screening lengths are seen to behave drastically differently depending on the EOS considered and therefore yield a way to distinguish the two equations of state in heavy ion collisions.

Green's function-stochastic methods framework for probing nonlinear evolution problems: Burger's equation, the nonlinear Schroedinger's equation, and hydrodynamic organization of near-molecular-scale vorticity

International Nuclear Information System (INIS)

Keanini, R.G.

2011-01-01

Research highlights: → Systematic approach for physically probing nonlinear and random evolution problems. → Evolution of vortex sheets corresponds to evolution of an Ornstein-Uhlenbeck process. → Organization of near-molecular scale vorticity mediated by hydrodynamic modes. → Framework allows calculation of vorticity evolution within random strain fields. - Abstract: A framework which combines Green's function (GF) methods and techniques from the theory of stochastic processes is proposed for tackling nonlinear evolution problems. The framework, established by a series of easy-to-derive equivalences between Green's function and stochastic representative solutions of linear drift-diffusion problems, provides a flexible structure within which nonlinear evolution problems can be analyzed and physically probed. As a preliminary test bed, two canonical, nonlinear evolution problems - Burgers' equation and the nonlinear Schroedinger's equation - are first treated. In the first case, the framework provides a rigorous, probabilistic derivation of the well known Cole-Hopf ansatz. Likewise, in the second, the machinery allows systematic recovery of a known soliton solution. The framework is then applied to a fairly extensive exploration of physical features underlying evolution of randomly stretched and advected Burger's vortex sheets. Here, the governing vorticity equation corresponds to the Fokker-Planck equation of an Ornstein-Uhlenbeck process, a correspondence that motivates an investigation of sub-sheet vorticity evolution and organization. Under the assumption that weak hydrodynamic fluctuations organize disordered, near-molecular-scale, sub-sheet vorticity, it is shown that these modes consist of two weakly damped counter-propagating cross-sheet acoustic modes, a diffusive cross-sheet shear mode, and a diffusive cross-sheet entropy mode. Once a consistent picture of in-sheet vorticity evolution is established, a number of analytical results, describing the

Isothermal Multiphase Flash Calculations with the PC-SAFT Equation of State

International Nuclear Information System (INIS)

Justo-Garcia, Daimler N.; Garcia-Sanchez, Fernando; Romero-Martinez, Ascencion

2008-01-01

A computational approach for isothermal multiphase flash calculations with the PC-SAFT (Perturbed-Chain Statistical Associating Fluid Theory) equation of state is presented. In the framework of the study of fluid phase equilibria of multicomponent systems, the general multiphase problem is the single most important calculation which consists of finding the correct number and types of phases and their corresponding equilibrium compositions such that the Gibbs energy of the system is a minimum. For solving this problem, the system Gibbs energy was minimized using a rigorous method for thermodynamic stability analysis to find the most stable state of the system. The efficiency and reliability of the approach to predict and calculate complex phase equilibria are illustrated by solving three typical problems encountered in the petroleum industry

Meta-analysis a structural equation modeling approach

CERN Document Server

Cheung, Mike W-L

2015-01-01

Presents a novel approach to conducting meta-analysis using structural equation modeling. Structural equation modeling (SEM) and meta-analysis are two powerful statistical methods in the educational, social, behavioral, and medical sciences. They are often treated as two unrelated topics in the literature. This book presents a unified framework on analyzing meta-analytic data within the SEM framework, and illustrates how to conduct meta-analysis using the metaSEM package in the R statistical environment. Meta-Analysis: A Structural Equation Modeling Approach begins by introducing the impo

Relativistic equations of state at finite temperature

International Nuclear Information System (INIS)

Santos, A.M.S.; Menezes, D.P.

2004-01-01

In this work we study the effects of temperature on the equations of state obtained within a relativistic model with and without β equilibrium, over a wide range of densities. We integrate the TOV equations. We also compare the results of the equation of state, effective mass and strangeness fraction from the TM1, NL3 and GL sets of parameters, as well as investigating the importance of antiparticles in the treatment. The have checked that TM1 and NL3 are not appropriate for the description of neutron and protoneutron stars. (author)

Solitary wave solutions to the modified form of Camassa-Holm equation by means of the homotopy analysis method

International Nuclear Information System (INIS)

Abbasbandy, S.

2009-01-01

Solitary wave solutions to the modified form of Camassa-Holm (CH) equation are sought. In this work, the homotopy analysis method (HAM), one of the most effective method, is applied to obtain the soliton wave solutions with and without continuity of first derivatives at crest

Converging shock flows for a Mie-Grüneisen equation of state

Science.gov (United States)

Ramsey, Scott D.; Schmidt, Emma M.; Boyd, Zachary M.; Lilieholm, Jennifer F.; Baty, Roy S.

2018-04-01

Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equations admit a wide variety of scale-invariant solutions (including the famous Noh, Sedov, and Guderley shock solutions) when the included equation of state (EOS) closure model assumes a certain scale-invariant form. However, this scale-invariant EOS class does not include even simple models used for shock compression of crystalline solids, including many broadly applicable representations of Mie-Grüneisen EOS. Intuitively, this incompatibility naturally arises from the presence of multiple dimensional scales in the Mie-Grüneisen EOS, which are otherwise absent from scale-invariant models that feature only dimensionless parameters (such as the adiabatic index in the ideal gas EOS). The current work extends previous efforts intended to rectify this inconsistency, by using a scale-invariant EOS model to approximate a Mie-Grüneisen EOS form. To this end, the adiabatic bulk modulus for the Mie-Grüneisen EOS is constructed, and its key features are used to motivate the selection of a scale-invariant approximation form. The remaining surrogate model parameters are selected through enforcement of the Rankine-Hugoniot jump conditions for an infinitely strong shock in a Mie-Grüneisen material. Finally, the approximate EOS is used in conjunction with the 1D inviscid Euler equations to calculate a semi-analytical Guderley-like imploding shock solution in a metal sphere and to determine if and when the solution may be valid for the underlying Mie-Grüneisen EOS.

Application of conjugate gradient method to Commix-1B three-dimensional momentum equation

International Nuclear Information System (INIS)

King, J.B.; Domanus, H.

1987-01-01

Conjugate gradient method which is a special case of the variational method was implemented in the momentum section of the COMMIX-1B thermal hydraulics code. The comparisons between this method and the conventional iterative method of Successive Over Relation (S.O.R.) were made. Using COMMIX-1B, three steady state problems were analyzed. These problems were flow distribution in a scaled model of the Clinch River Fast Breeder Reactor outlet plenum, flow of coolant in the cold leg and downcomer of a PWR and isothermal air flow through a partially blocked pipe. It was found that if the conjugate gradient method is used, the execution time required to solve the resulting COMMIX-1B system of equations can be reduced by a factor of about 2 for the first two problems. For the isothermal air flow problem, the conjugate gradient method did not improve the execution time

Simple equation method for nonlinear partial differential equations and its applications

Directory of Open Access Journals (Sweden)

Taher A. Nofal

2016-04-01

Full Text Available In this article, we focus on the exact solution of the some nonlinear partial differential equations (NLPDEs such as, Kodomtsev–Petviashvili (KP equation, the (2 + 1-dimensional breaking soliton equation and the modified generalized Vakhnenko equation by using the simple equation method. In the simple equation method the trial condition is the Bernoulli equation or the Riccati equation. It has been shown that the method provides a powerful mathematical tool for solving nonlinear wave equations in mathematical physics and engineering problems.

Partial differential equations with numerical methods

CERN Document Server

Larsson, Stig

2003-01-01

The book is suitable for advanced undergraduate and beginning graduate students of applied mathematics and engineering. The main theme is the integration of the theory of linear PDEs and the numerical solution of such equations. For each type of PDE, elliptic, parabolic, and hyperbolic, the text contains one chapter on the mathematical theory of the differential equation, followed by one chapter on finite difference methods and one on finite element methods. As preparation, the two-point boundary value problem and the initial-value problem for ODEs are discussed in separate chapters. There is also one chapter on the elliptic eigenvalue problem and eigenfunction expansion. The presentation does not presume a deep knowledge of mathematical and functional analysis. Some background on linear functional analysis and Sobolev spaces, and also on numerical linear algebra, is reviewed in two appendices.

Fractional dynamic calculus and fractional dynamic equations on time scales

CERN Document Server

Georgiev, Svetlin G

2018-01-01

Pedagogically organized, this monograph introduces fractional calculus and fractional dynamic equations on time scales in relation to mathematical physics applications and problems. Beginning with the definitions of forward and backward jump operators, the book builds from Stefan Hilger’s basic theories on time scales and examines recent developments within the field of fractional calculus and fractional equations. Useful tools are provided for solving differential and integral equations as well as various problems involving special functions of mathematical physics and their extensions and generalizations in one and more variables. Much discussion is devoted to Riemann-Liouville fractional dynamic equations and Caputo fractional dynamic equations.  Intended for use in the field and designed for students without an extensive mathematical background, this book is suitable for graduate courses and researchers looking for an introduction to fractional dynamic calculus and equations on time scales. .

Reduction of the allotropic transition temperature in nanocrystalline zirconium: Predicted by modified equation of state (MEOS) method and molecular dynamics simulation

Science.gov (United States)

Salati, Amin; Mokhtari, Esmail; Panjepour, Masoud; Aryanpour, Gholamreza

2013-04-01

The temperature at which polymorphic phase transformation occurs in nanocrystalline (NC) materials is different from that of coarse-grained specimens. This anomaly has been related to the role of grain boundary component in these materials and can be predicted by a dilated crystal model. In this study, based on this model, a modified equation of state (MEOS) method (instead of equation of state, EOS, method) is used to calculate the total Gibbs free energy of each phase (β-Zr or α-Zr) in NC Zr. Thereupon, the change in the total Gibbs free energy for β-Zr to α-Zr phase transformation (ΔGβ→α) via the grain size is calculated by this method. Similar to polymorphic transformation in other NC materials (Fe, Nb, Co, TiO2, Al2O3 and ZnS), it is found that the estimated transformation temperature in NC Zr (β→α) is reduced with decreasing grain size. Finally, a molecular dynamics (MD) simulation is employed to confirm the theoretical results.

Bound states of quarks calculated with stochastic integration of the Bethe-Salpeter equation

International Nuclear Information System (INIS)

Salomon, M.

1992-07-01

We have computed the masses, wave functions and sea quark content of mesons in their ground state by integrating the Bethe-Salpeter equation with a stochastic algorithm. This method allows the inclusion of a large set of diagrams. Inspection of the kernel of the equation shows that q-q-bar pairs with similar constituent masses in a singlet spin state exhibit a high bound state which is not present in other pairs. The pion, kaon and eta belongs to this category. 19 refs., 2 figs., 2 tabs

An exact and consistent adjoint method for high-fidelity discretization of the compressible flow equations

Science.gov (United States)

Subramanian, Ramanathan Vishnampet Ganapathi

Methods and computing hardware advances have enabled accurate predictions of complex compressible turbulence phenomena, such as the generation of jet noise that motivates the present effort. However, limited understanding of underlying physical mechanisms restricts the utility of such predictions since they do not, by themselves, indicate a route to design improvement. Gradient-based optimization using adjoints can circumvent the flow complexity to guide designs. Such methods have enabled sensitivity analysis and active control of turbulence at engineering flow conditions by providing gradient information at computational cost comparable to that of simulating the flow. They accelerate convergence of numerical design optimization algorithms, though this is predicated on the availability of an accurate gradient of the discretized flow equations. This is challenging to obtain, since both the chaotic character of the turbulence and the typical use of discretizations near their resolution limits in order to efficiently represent its smaller scales will amplify any approximation errors made in the adjoint formulation. Formulating a practical exact adjoint that avoids such errors is especially challenging if it is to be compatible with state-of-the-art simulation methods used for the turbulent flow itself. Automatic differentiation (AD) can provide code to calculate a nominally exact adjoint, but existing general-purpose AD codes are inefficient to the point of being prohibitive for large-scale turbulence simulations. We analyze the compressible flow equations as discretized using the same high-order workhorse methods used for many high-fidelity compressible turbulence simulations, and formulate a practical space--time discrete-adjoint method without changing the basic discretization. A key step is the definition of a particular discrete analog of the continuous norm that defines our cost functional; our selection leads directly to an efficient Runge--Kutta-like scheme

A Comparison of Kernel Equating and Traditional Equipercentile Equating Methods and the Parametric Bootstrap Methods for Estimating Standard Errors in Equipercentile Equating

Science.gov (United States)

Choi, Sae Il

2009-01-01

This study used simulation (a) to compare the kernel equating method to traditional equipercentile equating methods under the equivalent-groups (EG) design and the nonequivalent-groups with anchor test (NEAT) design and (b) to apply the parametric bootstrap method for estimating standard errors of equating. A two-parameter logistic item response…

Two-state approximation of the Fadeev-Hahn equations

International Nuclear Information System (INIS)

Brener, S.E.

1993-01-01

The equations have been chosen which allow both to solve the scattering problems and to calculate the parameters of bound states of three particles with Coulomb interaction when the system energy is below the decay to three separate particles. The method of constructing of equations which are most suitable for concrete problems is considered. Different numerical schemes to calculate the low energy scattering cross sections with two-particle clusterization in 'in' and 'out' collision's channels have been developed. The bounds of applied approaches were determined and the peculiarities connected with differently defined reaction amplitudes under these approaches have been considered. The interpretation of obtained results at different definitions of reaction amplitudes was demonstrated, and the elastic, inelastic cross sections and muon transfer rates in hydrogen isotope mesic atom collisions have been calculated using Fadeev-Hahn equations. (author)

Common intersection points in dense fluids via equations of state

International Nuclear Information System (INIS)

Parsafar, G. A.; Noorian, R.

2001-01-01

Some new of state which are derived for dense fluids in recent years, namely the linear isotherm regularity, the dense system equation of state, Ihm-Song-Mason equation of state, and a newly derived semi-empirical equation of state have used to investigate the common intersection point of isobaric expansivity (α p ) in dense fluids. We have shown that the accuracy of these equations of state in predicting such a common intersection point is reduced from the new semi-imperial equation of state, dense system equation of state, linear isotherm regularity, to Ihm-Song-Mason equation of state. respectively. Form physical point of view, the van der Waals equation of state is used to investigate such an intersection point. It is shown that the van der Waals repulsion forces and temperature dependency of the effective molecular diameter are important for existence of this common point. Finally, we have shown that the common intersection points of the isotherms of thermal pressure coefficient, the isotherms of heat capacity at constant volume, and the iso chores of internal pressure for a fluid are related to each other. Also, the common intersection points of the reduced bulk modulus and 1/(Tα p ) for isotherms of a fluid both appear at the same density

A multilevel correction adaptive finite element method for Kohn-Sham equation

Science.gov (United States)

Hu, Guanghui; Xie, Hehu; Xu, Fei

2018-02-01

In this paper, an adaptive finite element method is proposed for solving Kohn-Sham equation with the multilevel correction technique. In the method, the Kohn-Sham equation is solved on a fixed and appropriately coarse mesh with the finite element method in which the finite element space is kept improving by solving the derived boundary value problems on a series of adaptively and successively refined meshes. A main feature of the method is that solving large scale Kohn-Sham system is avoided effectively, and solving the derived boundary value problems can be handled efficiently by classical methods such as the multigrid method. Hence, the significant acceleration can be obtained on solving Kohn-Sham equation with the proposed multilevel correction technique. The performance of the method is examined by a variety of numerical experiments.

Born approximation to a perturbative numerical method for the solution of the Schrodinger equation

International Nuclear Information System (INIS)

Adam, Gh.

1978-05-01

A perturbative numerical (PN) method is given for the solution of a regular one-dimensional Cauchy problem arising from the Schroedinger equation. The present method uses a step function approximation for the potential. Global, free of scaling difficulty, forward and backward PN algorithms are derived within first order perturbation theory (Born approximation). A rigorous analysis of the local truncation errors is performed. This shows that the order of accuracy of the method is equal to four. In between the mesh points, the global formula for the wavefunction is accurate within O(h 4 ), while that for the first order derivative is accurate within O(h 3 ). (author)

Structural equation and log-linear modeling: a comparison of methods in the analysis of a study on caregivers' health

Directory of Open Access Journals (Sweden)

Rosenbaum Peter L

2006-10-01

Full Text Available Abstract Background In this paper we compare the results in an analysis of determinants of caregivers' health derived from two approaches, a structural equation model and a log-linear model, using the same data set. Methods The data were collected from a cross-sectional population-based sample of 468 families in Ontario, Canada who had a child with cerebral palsy (CP. The self-completed questionnaires and the home-based interviews used in this study included scales reflecting socio-economic status, child and caregiver characteristics, and the physical and psychological well-being of the caregivers. Both analytic models were used to evaluate the relationships between child behaviour, caregiving demands, coping factors, and the well-being of primary caregivers of children with CP. Results The results were compared, together with an assessment of the positive and negative aspects of each approach, including their practical and conceptual implications. Conclusion No important differences were found in the substantive conclusions of the two analyses. The broad confirmation of the Structural Equation Modeling (SEM results by the Log-linear Modeling (LLM provided some reassurance that the SEM had been adequately specified, and that it broadly fitted the data.

Robust Scale Transformation Methods in IRT True Score Equating under Common-Item Nonequivalent Groups Design

Science.gov (United States)

He, Yong

2013-01-01

Common test items play an important role in equating multiple test forms under the common-item nonequivalent groups design. Inconsistent item parameter estimates among common items can lead to large bias in equated scores for IRT true score equating. Current methods extensively focus on detection and elimination of outlying common items, which…

A novel numerical flux for the 3D Euler equations with general equation of state

KAUST Repository

Toro, Eleuterio F.

2015-09-30

Here we extend the flux vector splitting approach recently proposed in (E F Toro and M E Vázquez-Cendón. Flux splitting schemes for the Euler equations. Computers and Fluids. Vol. 70, Pages 1-12, 2012). The scheme was originally presented for the 1D Euler equations for ideal gases and its extension presented in this paper is threefold: (i) we solve the three-dimensional Euler equations on general meshes; (ii) we use a general equation of state; and (iii) we achieve high order of accuracy in both space and time through application of the semi-discrete ADER methodology on general meshes. The resulting methods are systematically assessed for accuracy, robustness and efficiency on a carefully selected suite of test problems. Formal high accuracy is assessed through convergence rates studies for schemes of up to 4th order of accuracy in both space and time on unstructured meshes.

Dynamic Modeling and Analysis of the Large-Scale Rotary Machine with Multi-Supporting

Directory of Open Access Journals (Sweden)

Xuejun Li

2011-01-01

Full Text Available The large-scale rotary machine with multi-supporting, such as rotary kiln and rope laying machine, is the key equipment in the architectural, chemistry, and agriculture industries. The body, rollers, wheels, and bearings constitute a chain multibody system. Axis line deflection is a vital parameter to determine mechanics state of rotary machine, thus body axial vibration needs to be studied for dynamic monitoring and adjusting of rotary machine. By using the Riccati transfer matrix method, the body system of rotary machine is divided into many subsystems composed of three elements, namely, rigid disk, elastic shaft, and linear spring. Multiple wheel-bearing structures are simplified as springs. The transfer matrices of the body system and overall transfer equation are developed, as well as the response overall motion equation. Taken a rotary kiln as an instance, natural frequencies, modal shape, and response vibration with certain exciting axis line deflection are obtained by numerical computing. The body vibration modal curves illustrate the cause of dynamical errors in the common axis line measurement methods. The displacement response can be used for further measurement dynamical error analysis and compensation. The response overall motion equation could be applied to predict the body motion under abnormal mechanics condition, and provide theory guidance for machine failure diagnosis.

Critical comparison between equation of motion-Green's function methods and configuration interaction methods: analysis of methods and applications

International Nuclear Information System (INIS)

Freed, K.F.; Herman, M.F.; Yeager, D.L.

1980-01-01

A description is provided of the common conceptual origins of many-body equations of motion and Green's function methods in Liouville operator formulations of the quantum mechanics of atomic and molecular electronic structure. Numerical evidence is provided to show the inadequacies of the traditional strictly perturbative approaches to these methods. Nonperturbative methods are introduced by analogy with techniques developed for handling large configuration interaction calculations and by evaluating individual matrix elements to higher accuracy. The important role of higher excitations is exhibited by the numerical calculations, and explicit comparisons are made between converged equations of motion and configuration interaction calculations for systems where a fundamental theorem requires the equality of the energy differences produced by these different approaches. (Auth.)

State-dependent differential Riccati equation to track control of time-varying systems with state and control nonlinearities.

Science.gov (United States)

Korayem, M H; Nekoo, S R

2015-07-01

This work studies an optimal control problem using the state-dependent Riccati equation (SDRE) in differential form to track for time-varying systems with state and control nonlinearities. The trajectory tracking structure provides two nonlinear differential equations: the state-dependent differential Riccati equation (SDDRE) and the feed-forward differential equation. The independence of the governing equations and stability of the controller are proven along the trajectory using the Lyapunov approach. Backward integration (BI) is capable of solving the equations as a numerical solution; however, the forward solution methods require the closed-form solution to fulfill the task. A closed-form solution is introduced for SDDRE, but the feed-forward differential equation has not yet been obtained. Different ways of solving the problem are expressed and analyzed. These include BI, closed-form solution with corrective assumption, approximate solution, and forward integration. Application of the tracking problem is investigated to control robotic manipulators possessing rigid or flexible joints. The intention is to release a general program for automatic implementation of an SDDRE controller for any manipulator that obeys the Denavit-Hartenberg (D-H) principle when only D-H parameters are received as input data. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Initial state dependence of nonlinear kinetic equations: The classical electron gas

International Nuclear Information System (INIS)

Marchetti, M.C.; Cohen, E.G.D.; Dorfman, J.R.; Kirkpatrick, T.R.

1985-01-01

The method of nonequilibrium cluster expansion is used to study the decay to equilibrium of a weakly coupled inhomogeneous electron gas prepared in a local equilibrium state at the initial time, t=0. A nonlinear kinetic equation describing the long time behavior of the one-particle distribution function is obtained. For consistency, initial correlations have to be taken into account. The resulting kinetic equation-differs from that obtained when the initial state of the system is assumed to be factorized in a product of one-particle functions. The question of to what extent correlations in the initial state play an essential role in determining the form of the kinetic equation at long times is discussed. To that end, the present calculations are compared wih results obtained before for hard sphere gases and in general with strong short-range forces. A partial answer is proposed and some open questions are indicated

Stochastic wave-function unravelling of the generalized Lindblad equation using correlated states

International Nuclear Information System (INIS)

Moodley, Mervlyn; Nsio Nzundu, T; Paul, S

2012-01-01

We perform a stochastic wave-function unravelling of the generalized Lindblad master equation using correlated states, a combination of the system state vectors and the environment population. The time-convolutionless projection operator method using correlated projection superoperators is applied to a two-state system, a qubit, that is coupled to an environment consisting of two energy bands which are both populated. These results are compared to the data obtained from Monte Carlo wave-function simulations based on the unravelling of the master equation. We also show a typical quantum trajectory and the average time evolution of the state vector on the Bloch sphere. (paper)

Auxiliary equation method for solving nonlinear partial differential equations

International Nuclear Information System (INIS)

Sirendaoreji,; Jiong, Sun

2003-01-01

By using the solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct several kinds of exact travelling wave solutions for some nonlinear partial differential equations. By this method some physically important nonlinear equations are investigated and new exact travelling wave solutions are explicitly obtained with the aid of symbolic computation

Least squares shadowing sensitivity analysis of a modified Kuramoto–Sivashinsky equation

International Nuclear Information System (INIS)

Blonigan, Patrick J.; Wang, Qiqi

2014-01-01

Highlights: •Modifying the Kuramoto–Sivashinsky equation and changing its boundary conditions make it an ergodic dynamical system. •The modified Kuramoto–Sivashinsky equation exhibits distinct dynamics for three different ranges of system parameters. •Least squares shadowing sensitivity analysis computes accurate gradients for a wide range of system parameters. - Abstract: Computational methods for sensitivity analysis are invaluable tools for scientists and engineers investigating a wide range of physical phenomena. However, many of these methods fail when applied to chaotic systems, such as the Kuramoto–Sivashinsky (K–S) equation, which models a number of different chaotic systems found in nature. The following paper discusses the application of a new sensitivity analysis method developed by the authors to a modified K–S equation. We find that least squares shadowing sensitivity analysis computes accurate gradients for solutions corresponding to a wide range of system parameters

Similarity, Clustering, and Scaling Analyses for the Foreign Exchange Market ---Comprehensive Analysis on States of Market Participants with High-Frequency Financial Data---

Science.gov (United States)

Sato, A.; Sakai, H.; Nishimura, M.; Holyst, J. A.

This article proposes mathematical methods to quantify states of marketparticipants in the foreign exchange market (FX market) and conduct comprehensive analysis on behavior of market participants by means of high-frequency financial data. Based on econophysics tools and perspectives we study similarity measures for both rate movements and quotation activities among various currency pairs. We perform also clustering analysis on market states for observation days, and find scaling relationship between mean values of quotation activities and their standard deviations. Using these mathematical methods we can visualize states of the FX market comprehensively. Finally we conclude that states of market participants temporally vary due to both external and internal factors.

Development of a large-scale general purpose two-phase flow analysis code

International Nuclear Information System (INIS)

Terasaka, Haruo; Shimizu, Sensuke

2001-01-01

A general purpose three-dimensional two-phase flow analysis code has been developed for solving large-scale problems in industrial fields. The code uses a two-fluid model to describe the conservation equations for two-phase flow in order to be applicable to various phenomena. Complicated geometrical conditions are modeled by FAVOR method in structured grid systems, and the discretization equations are solved by a modified SIMPLEST scheme. To reduce computing time a matrix solver for the pressure correction equation is parallelized with OpenMP. Results of numerical examples show that the accurate solutions can be obtained efficiently and stably. (author)

A Solution Space for a System of Null-State Partial Differential Equations: Part 2

Science.gov (United States)

Flores, Steven M.; Kleban, Peter

2015-01-01

This article is the second of four that completely and rigorously characterize a solution space for a homogeneous system of 2 N + 3 linear partial differential equations in 2 N variables that arises in conformal field theory (CFT) and multiple Schramm-Löwner evolution (SLE). The system comprises 2 N null-state equations and three conformal Ward identities which govern CFT correlation functions of 2 N one-leg boundary operators. In the first article (Flores and Kleban, Commun Math Phys, arXiv:1212.2301, 2012), we use methods of analysis and linear algebra to prove that dim , with C N the Nth Catalan number. The analysis of that article is complete except for the proof of a lemma that it invokes. The purpose of this article is to provide that proof. The lemma states that if every interval among ( x 2, x 3), ( x 3, x 4),…,( x 2 N-1, x 2 N ) is a two-leg interval of (defined in Flores and Kleban, Commun Math Phys, arXiv:1212.2301, 2012), then F vanishes. Proving this lemma by contradiction, we show that the existence of such a nonzero function implies the existence of a non-vanishing CFT two-point function involving primary operators with different conformal weights, an impossibility. This proof (which is rigorous in spite of our occasional reference to CFT) involves two different types of estimates, those that give the asymptotic behavior of F as the length of one interval vanishes, and those that give this behavior as the lengths of two intervals vanish simultaneously. We derive these estimates by using Green functions to rewrite certain null-state PDEs as integral equations, combining other null-state PDEs to obtain Schauder interior estimates, and then repeatedly integrating the integral equations with these estimates until we obtain optimal bounds. Estimates in which two interval lengths vanish simultaneously divide into two cases: two adjacent intervals and two non-adjacent intervals. The analysis of the latter case is similar to that for one vanishing

Hydrostatic Equilibria of Rotating Stars with Realistic Equation of State

Science.gov (United States)

Yasutake, Nobutoshi; Fujisawa, Kotaro; Okawa, Hirotada; Yamada, Shoichi

Stars rotate generally, but it is a non-trivial issue to obtain hydrostatic equilibria for rapidly rotating stars theoretically, especially for baroclinic cases, in which the pressure depends not only on the density, but also on the temperature and compositions. It is clear that the stellar structures with realistic equation of state are the baroclinic cases, but there are not so many studies for such equilibria. In this study, we propose two methods to obtain hydrostatic equilibria considering rotation and baroclinicity, namely the weak-solution method and the strong-solution method. The former method is based on the variational principle, which is also applied to the calculation of the inhomogeneous phases, known as the pasta structures, in crust of neutron stars. We found this method might break the balance equation locally, then introduce the strong-solution method. Note that our method is formulated in the mass coordinate, and it is hence appropriated for the stellar evolution calculations.

Pion production and the nuclear equation of state

International Nuclear Information System (INIS)

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

1984-10-01

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

A novel scale for measuring mixed states in bipolar disorder.

Science.gov (United States)

Cavanagh, Jonathan; Schwannauer, Matthias; Power, Mick; Goodwin, Guy M

2009-01-01

Conventional descriptions of bipolar disorder tend to treat the mixed state as something of an afterthought. There is no scale that specifically measures the phenomena of the mixed state. This study aimed to test a novel scale for mixed state in a clinical and community population of bipolar patients. The scale included clinically relevant symptoms of both mania and depression in a bivariate scale. Recovered respondents were asked to recall their last manic episode. The scale allowed endorsement of one or more of the manic and depressive symptoms. Internal consistency analyses were carried out using Cronbach alpha. Factor analysis was carried out using a standard Principal Components Analysis followed by Varimax Rotation. A confirmatory factor analytic method was used to validate the scale structure in a representative clinical sample. The reliability analysis gave a Cronbach alpha value of 0.950, with a range of corrected-item-total-scale correlations from 0.546 (weight change) to 0.830 (mood). The factor analysis revealed a two-factor solution for the manic and depressed items which accounted for 61.2% of the variance in the data. Factor 1 represented physical activity, verbal activity, thought processes and mood. Factor 2 represented eating habits, weight change, passage of time and pain sensitivity. This novel scale appears to capture the key features of mixed states. The two-factor solution fits well with previous models of bipolar disorder and concurs with the view that mixed states may be more than the sum of their parts.

Solution of generalized control system equations at steady state

International Nuclear Information System (INIS)

Vilim, R.B.

1987-01-01

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

Iteratively-coupled propagating exterior complex scaling method for electron-hydrogen collisions

International Nuclear Information System (INIS)

Bartlett, Philip L; Stelbovics, Andris T; Bray, Igor

2004-01-01

A newly-derived iterative coupling procedure for the propagating exterior complex scaling (PECS) method is used to efficiently calculate the electron-impact wavefunctions for atomic hydrogen. An overview of this method is given along with methods for extracting scattering cross sections. Differential scattering cross sections at 30 eV are presented for the electron-impact excitation to the n = 1, 2, 3 and 4 final states, for both PECS and convergent close coupling (CCC), which are in excellent agreement with each other and with experiment. PECS results are presented at 27.2 eV and 30 eV for symmetric and asymmetric energy-sharing triple differential cross sections, which are in excellent agreement with CCC and exterior complex scaling calculations, and with experimental data. At these intermediate energies, the efficiency of the PECS method with iterative coupling has allowed highly accurate partial-wave solutions of the full Schroedinger equation, for L ≤ 50 and a large number of coupled angular momentum states, to be obtained with minimal computing resources. (letter to the editor)

Equation-of-motion coupled cluster method for high spin double electron attachment calculations

Energy Technology Data Exchange (ETDEWEB)

Musiał, Monika, E-mail: musial@ich.us.edu.pl; Lupa, Łukasz; Kucharski, Stanisław A. [Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice (Poland)

2014-03-21

The new formulation of the equation-of-motion (EOM) coupled cluster (CC) approach applicable to the calculations of the double electron attachment (DEA) states for the high spin components is proposed. The new EOM equations are derived for the high spin triplet and quintet states. In both cases the new equations are easier to solve but the substantial simplification is observed in the case of quintets. Out of 21 diagrammatic terms contributing to the standard DEA-EOM-CCSDT equations for the R{sub 2} and R{sub 3} amplitudes only four terms survive contributing to the R{sub 3} part. The implemented method has been applied to the calculations of the excited states (singlets, triplets, and quintets) energies of the carbon and silicon atoms and potential energy curves for selected states of the Na{sub 2} (triplets) and B{sub 2} (quintets) molecules.

Systematic analysis of scaling properties in deep inelastic scattering

International Nuclear Information System (INIS)

Beuf, Guillaume; Peschanski, Robi; Royon, Christophe; Salek, David

2008-01-01

Using the 'quality factor' method, we analyze the scaling properties of deep inelastic processes at the accelerator HERA and fixed target experiments for x≤10 -2 . We look for scaling formulas of the form σ γ * p (τ), where τ(L=logQ 2 ,Y) is a scaling variable suggested by the asymptotic properties of QCD evolution equations with rapidity Y. We consider four cases: 'fixed coupling', corresponding to the original geometric scaling proposal and motivated by the asymptotic properties of the Balitsky-Kovchegov equation with fixed QCD coupling constant; two versions, 'running coupling I, II,' of the scaling suggested by the Balitsky-Kovchegov equation with running coupling; and 'diffusive scaling' suggested by the QCD evolution equation with Pomeron loops. The quality factors, quantifying the phenomenological validity of the candidate scaling variables, are fitted on the total and deeply virtual Compton scattering cross-section data from HERA and predictions are made for the elastic vector meson and for the diffractive cross sections at fixed small x P or β. The first three scaling formulas have comparably good quality factors while the fourth one is disfavored. Adjusting initial conditions gives a significant improvement of the running coupling II scaling.

Proposed equation of state experimental program at NOVA/NIF

International Nuclear Information System (INIS)

Holmes, N.C.; Cauble, R.; Celliers, P.; Collins, G.; DaSilva, L.; Hammel, B.; Stewart, R.; Strand, O.; Sullivan, A.

1996-03-01

The high pressure equations of state of materials must be accurately known for confident design, simulation, and interpretation of ICF target experiments and for weapons. Here, the authors sketch out a program to perform precise and accurate equation of state (EOS) experiments using large, high-power lasers. The program is divided into three phases: (1) a driver qualification program which will determine the necessary drive and target design parameters; (2) characterization of low and high-Z standard materials; (3) accurate impedance-match experiments to determine equations of state of materials of interest relative to the qualified standards. Novel methods are proposed for the first phase which are an order of magnitude more sensitive than those applied previously. They identify and describe the choices for standards and the samples slated for initial studies. In all cases, their goal is to obtain data that is of sufficiently high quality to serve for design purposes and to validate theories of material response at high pressures of the order of 1--10 TPa. This is particularly important right now, for materials such as CH plastics and hydrogen (DT). In addition, they will suggest a program for EOS measurements for P > 10 TPa, the region of maximum theoretical EOS uncertainty for many materials

Instantaneous equations for multiphase flow in porous media without length-scale restrictions using a non-local averaging volume

International Nuclear Information System (INIS)

Espinosa-Paredes, Gilberto

2010-01-01

The aim of this paper is to propose a framework to obtain a new formulation for multiphase flow conservation equations without length-scale restrictions, based on the non-local form of the averaged volume conservation equations. The simplification of the local averaging volume of the conservation equations to obtain practical equations is subject to the following length-scale restrictions: d << l << L, where d is the characteristic length of the dispersed phases, l is the characteristic length of the averaging volume, and L is the characteristic length of the physical system. If the foregoing inequality does not hold, or if the scale of the problem of interest is of the order of l, the averaging technique and therefore, the macroscopic theories of multiphase flow should be modified in order to include appropriate considerations and terms in the corresponding equations. In these cases the local form of the averaged volume conservation equations are not appropriate to describe the multiphase system. As an example of the conservation equations without length-scale restrictions, the natural circulation boiling water reactor was consider to study the non-local effects on the thermal-hydraulic core performance during steady-state and transient behaviors, and the results were compared with the classic local averaging volume conservation equations.

The Noble-Abel Stiffened-Gas equation of state

Science.gov (United States)

Le Métayer, Olivier; Saurel, Richard

2016-04-01

Hyperbolic two-phase flow models have shown excellent ability for the resolution of a wide range of applications ranging from interfacial flows to fluid mixtures with several velocities. These models account for waves propagation (acoustic and convective) and consist in hyperbolic systems of partial differential equations. In this context, each phase is compressible and needs an appropriate convex equation of state (EOS). The EOS must be simple enough for intensive computations as well as boundary conditions treatment. It must also be accurate, this being challenging with respect to simplicity. In the present approach, each fluid is governed by a novel EOS named "Noble Abel stiffened gas," this formulation being a significant improvement of the popular "Stiffened Gas (SG)" EOS. It is a combination of the so-called "Noble-Abel" and "stiffened gas" equations of state that adds repulsive effects to the SG formulation. The determination of the various thermodynamic functions and associated coefficients is the aim of this article. We first use thermodynamic considerations to determine the different state functions such as the specific internal energy, enthalpy, and entropy. Then we propose to determine the associated coefficients for a liquid in the presence of its vapor. The EOS parameters are determined from experimental saturation curves. Some examples of liquid-vapor fluids are examined and associated parameters are computed with the help of the present method. Comparisons between analytical and experimental saturation curves show very good agreement for wide ranges of temperature for both liquid and vapor.

Homotopy analysis solutions of point kinetics equations with one delayed precursor group

International Nuclear Information System (INIS)

Zhu Qian; Luo Lei; Chen Zhiyun; Li Haofeng

2010-01-01

Homotopy analysis method is proposed to obtain series solutions of nonlinear differential equations. Homotopy analysis method was applied for the point kinetics equations with one delayed precursor group. Analytic solutions were obtained using homotopy analysis method, and the algorithm was analysed. The results show that the algorithm computation time and precision agree with the engineering requirements. (authors)

Discontinuous Galerkin Subgrid Finite Element Method for Heterogeneous Brinkman’s Equations

KAUST Repository

Iliev, Oleg P.

2010-01-01

We present a two-scale finite element method for solving Brinkman\\'s equations with piece-wise constant coefficients. This system of equations model fluid flows in highly porous, heterogeneous media with complex topology of the heterogeneities. We make use of the recently proposed discontinuous Galerkin FEM for Stokes equations by Wang and Ye in [12] and the concept of subgrid approximation developed for Darcy\\'s equations by Arbogast in [4]. In order to reduce the error along the coarse-grid interfaces we have added a alternating Schwarz iteration using patches around the coarse-grid boundaries. We have implemented the subgrid method using Deal.II FEM library, [7], and we present the computational results for a number of model problems. © 2010 Springer-Verlag Berlin Heidelberg.

Structural equation modeling methods and applications

CERN Document Server

Wang, Jichuan

2012-01-01

A reference guide for applications of SEM using Mplus Structural Equation Modeling: Applications Using Mplus is intended as both a teaching resource and a reference guide. Written in non-mathematical terms, this book focuses on the conceptual and practical aspects of Structural Equation Modeling (SEM). Basic concepts and examples of various SEM models are demonstrated along with recently developed advanced methods, such as mixture modeling and model-based power analysis and sample size estimate for SEM. The statistical modeling program, Mplus, is also featured and provides researchers with a

Mathematical Methods for Engineers and Scientists 2 Vector Analysis, Ordinary Differential Equations and Laplace Transforms

CERN Document Server

Tang, Kwong-Tin

2007-01-01

Pedagogical insights gained through 30 years of teaching applied mathematics led the author to write this set of student-oriented books. Topics such as complex analysis, matrix theory, vector and tensor analysis, Fourier analysis, integral transforms, ordinary and partial differential equations are presented in a discursive style that is readable and easy to follow. Numerous clearly stated, completely worked out examples together with carefully selected problem sets with answers are used to enhance students' understanding and manipulative skill. The goal is to make students comfortable and confident in using advanced mathematical tools in junior, senior, and beginning graduate courses.

A new equation of state for better liquid density prediction of natural gas systems

Science.gov (United States)

Nwankwo, Princess C.

Equations of state formulations, modifications and applications have remained active research areas since the success of van der Waal's equation in 1873. The need for better reservoir fluid modeling and characterization is of great importance to petroleum engineers who deal with thermodynamic related properties of petroleum fluids at every stage of the petroleum "life span" from its drilling, to production through the wellbore, to transportation, metering and storage. Equations of state methods are far less expensive (in terms of material cost and time) than laboratory or experimental forages and the results are interestingly not too far removed from the limits of acceptable accuracy. In most cases, the degree of accuracy obtained, by using various EOS's, though not appreciable, have been acceptable when considering the gain in time. The possibility of obtaining an equation of state which though simple in form and in use, could have the potential of further narrowing the present existing bias between experimentally determined and popular EOS estimated results spurred the interest that resulted in this study. This research study had as its chief objective, to develop a new equation of state that would more efficiently capture the thermodynamic properties of gas condensate fluids, especially the liquid phase density, which is the major weakness of other established and popular cubic equations of state. The set objective was satisfied by a new semi analytical cubic three parameter equation of state, derived by the modification of the attraction term contribution to pressure of the van der Waal EOS without compromising either structural simplicity or accuracy of estimating other vapor liquid equilibria properties. The application of new EOS to single and multi-component light hydrocarbon fluids recorded far lower error values than does the popular two parameter, Peng-Robinson's (PR) and three parameter Patel-Teja's (PT) equations of state. Furthermore, this research

Observational constraints on cosmological models with Chaplygin gas and quadratic equation of state

International Nuclear Information System (INIS)

Sharov, G.S.

2016-01-01

Observational manifestations of accelerated expansion of the universe, in particular, recent data for Type Ia supernovae, baryon acoustic oscillations, for the Hubble parameter H ( z ) and cosmic microwave background constraints are described with different cosmological models. We compare the ΛCDM, the models with generalized and modified Chaplygin gas and the model with quadratic equation of state. For these models we estimate optimal model parameters and their permissible errors with different approaches to calculation of sound horizon scale r s ( z d ). Among the considered models the best value of χ 2 is achieved for the model with quadratic equation of state, but it has 2 additional parameters in comparison with the ΛCDM and therefore is not favored by the Akaike information criterion.

A method of orbital analysis for large-scale first-principles simulations

Energy Technology Data Exchange (ETDEWEB)

Ohwaki, Tsukuru [Advanced Materials Laboratory, Nissan Research Center, Nissan Motor Co., Ltd., 1 Natsushima-cho, Yokosuka, Kanagawa 237-8523 (Japan); Otani, Minoru [Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8568 (Japan); Ozaki, Taisuke [Research Center for Simulation Science (RCSS), Japan Advanced Institute of Science and Technology (JAIST), 1-1 Asahidai, Nomi, Ishikawa 923-1292 (Japan)

2014-06-28

An efficient method of calculating the natural bond orbitals (NBOs) based on a truncation of the entire density matrix of a whole system is presented for large-scale density functional theory calculations. The method recovers an orbital picture for O(N) electronic structure methods which directly evaluate the density matrix without using Kohn-Sham orbitals, thus enabling quantitative analysis of chemical reactions in large-scale systems in the language of localized Lewis-type chemical bonds. With the density matrix calculated by either an exact diagonalization or O(N) method, the computational cost is O(1) for the calculation of NBOs associated with a local region where a chemical reaction takes place. As an illustration of the method, we demonstrate how an electronic structure in a local region of interest can be analyzed by NBOs in a large-scale first-principles molecular dynamics simulation for a liquid electrolyte bulk model (propylene carbonate + LiBF{sub 4})

A method of orbital analysis for large-scale first-principles simulations

International Nuclear Information System (INIS)

Ohwaki, Tsukuru; Otani, Minoru; Ozaki, Taisuke

2014-01-01

An efficient method of calculating the natural bond orbitals (NBOs) based on a truncation of the entire density matrix of a whole system is presented for large-scale density functional theory calculations. The method recovers an orbital picture for O(N) electronic structure methods which directly evaluate the density matrix without using Kohn-Sham orbitals, thus enabling quantitative analysis of chemical reactions in large-scale systems in the language of localized Lewis-type chemical bonds. With the density matrix calculated by either an exact diagonalization or O(N) method, the computational cost is O(1) for the calculation of NBOs associated with a local region where a chemical reaction takes place. As an illustration of the method, we demonstrate how an electronic structure in a local region of interest can be analyzed by NBOs in a large-scale first-principles molecular dynamics simulation for a liquid electrolyte bulk model (propylene carbonate + LiBF 4 )