Analysis and Insights from a Dynamical Model of Nuclear Plant Safety Risk

In this paper, we expand upon previously reported results of a dynamical systems model for the impact of plant processes and programmatic performance on nuclear plant safety risk. We utilize both analytical techniques and numerical simulations typical of the analysis of nonlinear dynamical systems to obtain insights important for effective risk management. This includes use of bifurcation diagrams to show that period doubling bifurcations and regions of chaotic dynamics can occur. We also investigate the impact of risk mitigating functions (equipment reliability and loss prevention) on plant safety risk and demonstrate that these functions are capable of improving risk to levels that are better than those that are represented in a traditional risk assessment. Next, we analyze the system response to the presence of external noise and obtain some conclusions with respect to the allocation of resources to ensure that safety is maintained at optimal levels. In particular, we demonstrate that the model supports the importance of management and regulator attention to plants that have demonstrated poor performance by providing an external stimulus to obtain desired improvements. Equally important, the model suggests that excessive intervention, by either plant management or regulatory authorities, can have a deleterious impact on safety for plants that are operating with very effective programs and processes. Finally, we propose a modification to the model that accounts for the impact of plant risk culture on process performance and plant safety risk. We then use numerical simulations to demonstrate the important safety benefits of a strong risk culture.Nonlinear Dynamical Systems, Process Model, Risk Management

An interpolation problem in the Denjoy–Carleman classes

Inspired by some iterative algorithms useful for proving the real analyticity (or the Gevrey regularity) of a solution of a linear partial differential equation with real-analytic coefficients, we consider the following question. Given a smooth function defined on [a,b]subset of R$[a,b]\subset {\mathbb {R}}$ and given an increasing divergent sequence dn$d_n$ of positive integers such that the derivative of order dn$d_n$ of f has a growth of the type Mdn$M_{d_n}$, when can we deduce that f is a function in the Denjoy-Carleman class CM([a,b])C<^>M([a,b])? We provide a positive result and show that a suitable condition on the gaps between the terms of the sequence dn$d_n$ is needed

Study of the one-dimensional off-lattice hot-monomer reaction model

Hot monomers are particles having a transient mobility (a ballistic flight) prior to being definitely absorbed on a surface. After arriving at a surface, the excess energy coming from the kinetic energy in the gas phase is dissipated through degrees of freedom parallel to the surface plane. In this paper we study the hot monomer-monomer adsorption-reaction process on a continuum (off-lattice) one-dimensional space by means of Monte Carlo simulations. The system exhibits second-order irreversible phase transition between a reactive and saturated (absorbing) phases which belong to the directed percolation (DP) universality class. This result is interpreted by means of a coarse-grained Langevin description which allows as to extend the DP conjecture to transitions occurring in continuous media.Comment: 13 pages, 5 figures, final version to appear in J. Phys.

Numerical study of a first-order irreversible phase transition in a CO+NO catalyzed reaction model

The first-order irreversible phase transitions (IPT) of the Yaldran-Khan model (Yaldran-Khan, J. Catal. 131, 369, 1991) for the CO+NO reaction is studied using the constant coverage (CC) ensemble and performing epidemic simulations. The CC method allows the study of hysteretic effects close to coexistence as well as the location of both the upper spinodal point and the coexistence point. Epidemic studies show that at coexistence the number of active sites decreases according to a (short-time) power law followed by a (long-time) exponential decay. It is concluded that first-order IPT's share many characteristic of their reversible counterparts, such as the development of short ranged correlations, hysteretic effects, metastabilities, etc.Comment: 17 pages, 10 figure

Analytic and Gevrey Hypoellipticity for Perturbed Sums of Squares Operators

We prove a couple of results concerning pseudodifferential perturbations of differential operators being sums of squares of vector fields and satisfying H\"ormander's condition. The first is on the minimal Gevrey regularity: if a sum of squares with analytic coefficients is perturbed with a pseudodifferential operator of order strictly less than its subelliptic index it still has the Gevrey minimal regularity. We also prove a statement concerning real analytic hypoellipticity for the same type of pseudodifferential perturbations, provided the operator satisfies to some extra conditions (see Theorem 1.2 below) that ensure the analytic hypoellipticity

Measurements of the Yield Stress in Frictionless Granular Systems

We perform extensive molecular dynamics simulations of 2D frictionless granular materials to determine whether these systems can be characterized by a single static yield shear stress. We consider boundary-driven planar shear at constant volume and either constant shear force or constant shear velocity. Under steady flow conditions, these two ensembles give similar results for the average shear stress versus shear velocity. However, near jamming it is possible that the shear stress required to initiate shear flow can differ substantially from the shear stress required to maintain flow. We perform several measurements of the shear stress near the initiation and cessation of flow. At fixed shear velocity, we measure the average shear stress $\Sigma_{yv}$ in the limit of zero shear velocity. At fixed shear force, we measure the minimum shear stress $\Sigma_{yf}$ required to maintain steady flow at long times. We find that in finite-size systems $\Sigma_{yf} > \Sigma_{yv}$, which implies that there is a jump discontinuity in the shear velocity from zero to a finite value when these systems begin flowing at constant shear force. However, our simulations show that the difference $\Sigma_{yf} - \Sigma_{yv}$, and thus the discontinuity in the shear velocity, tend to zero in the infinite system size limit. Thus, our results indicate that in the large system limit, frictionless granular systems are characterized by a single static yield shear stress. We also monitor the short-time response of these systems to applied shear and show that the packing fraction of the system and shape of the velocity profile can strongly influence whether or not the shear stress at short times overshoots the long-time average value.Comment: 7 pages and 6 figure

Kolmogorov-Smirnov Test Distinguishes Attractors with Similar Dimensions

Recent advances in nonlinear dynamics have led to more informative characterizations of complex signals making it possible to probe correlations in data to which traditional linear statistical and spectral analyses were not sensitive. Many of these new tools require detailed knowledge of small scale structures of the attractor; knowledge that can be acquired only from relatively large amounts of precise data that are not contaminated by noise-not the kind of data one usually obtains from experiments. There is a need for tools that can take advantage of \u27\u27coarse-grained\u27\u27 information, but which nevertheless remain sensitive to higher-order correlations in the data. We propose that the correlation integral, now much used as an intermediate step in the calculation of dimensions and entropies, can be used as such a tool and that the Kolmogorov-Smirnov test is a convenient and reliable way of comparing correlation integrals quantitatively. This procedure makes it possible to distinguish between attractors with similar dimensions. For example, it can unambiguously distinguish (p \u3c 10(-8)) the Lorenz, Rossler, and Mackey-Glass (delay = 17) attractors whose correlation dimensions are within 1% of each other. We also show that the Kolmogorov-Smirnov test is a convenient way of comparing a data set with its surrogates

Optimization of the transmission of observable expectation values and observable statistics in Continuous Variable Teleportation

We analyze the statistics of observables in continuous variable quantum teleportation in the formalism of the characteristic function. We derive expressions for average values of output state observables in particular cumulants which are additive in terms of the input state and the resource of teleportation. Working with Squeezed Bell-like states, which may be optimized in a free parameter for better teleportation performance we discuss the relation between resources optimal for fidelity and for different observable averages. We obtain the values of the free parameter which optimize the central momenta and cumulants up to fourth order. For the cumulants the distortion between in and out states due to teleportation depends only on the resource. We obtain optimal parameters for the second and fourth order cumulants which do not depend on the squeezing of the resource. The second order central momenta which is equal to the second order cumulants and the photon number average are optimized by the same resource. We show that the optimal fidelity resource, found in reference (Phys. Rev. A {\bf 76}, 022301 (2007)) to depend also on the characteristics of input, tends for high squeezing to the resource which optimizes the second order momenta. A similar behavior is obtained for the resource which optimizes the photon statistics which is treated here using the sum of the squared differences in photon probabilities of input and output states as the distortion measure. This is interpreted to mean that the distortions associated to second order momenta dominates the behavior of the output state for large squeezing of the resource. Optimal fidelity and optimal photon statistics resources are compared and is shown that for mixtures of Fock states they are equivalent.Comment: 25 pages, 11 figure

Epidemic analysis of the second-order transition in the Ziff-Gulari-Barshad surface-reaction model

We study the dynamic behavior of the Ziff-Gulari-Barshad (ZGB) irreversible surface-reaction model around its kinetic second-order phase transition, using both epidemic and poisoning-time analyses. We find that the critical point is given by p_1 = 0.3873682 \pm 0.0000015, which is lower than the previous value. We also obtain precise values of the dynamical critical exponents z, \delta, and \eta which provide further numerical evidence that this transition is in the same universality class as directed percolation.Comment: REVTEX, 4 pages, 5 figures, Submitted to Physical Review

Damage Spreading in a Driven Lattice Gas Model

We studied damage spreading in a Driven Lattice Gas (DLG) model as a function of the temperature $T$, the magnitude of the external driving field $E$, and the lattice size. The DLG model undergoes an order-disorder second-order phase transition at the critical temperature $T_c(E)$, such that the ordered phase is characterized by high-density strips running along the direction of the applied field; while in the disordered phase one has a lattice-gas-like behaviour. It is found that the damage always spreads for all the investigated temperatures and reaches a saturation value $D_{sat}$ that depends only on $T$. $D_{sat}$ increases for $TT_c(E=\infty)$ and is free of finite-size effects. This behaviour can be explained as due to the existence of interfaces between the high-density strips and the lattice-gas-like phase whose roughness depends on $T$. Also, we investigated damage spreading for a range of finite fields as a function of $T$, finding a behaviour similar to that of the case with $E=\infty$.Comment: 13 pages, 7 figures. Submitted to "Journal of Statistical Mechanics: Theory and Experiment