Reconstruction of arbitrary biochemical reaction networks: A compressive sensing approach

Reconstruction of biochemical reaction networks (BRN) and genetic regulatory networks (GRN) in particular is a central topic in systems biology which raises crucial theoretical challenges in system identification. Nonlinear Ordinary Differential Equations (ODEs) that involve polynomial and rational functions are typically used to model biochemical reaction networks. Such nonlinear models make the problem of determining the connectivity of biochemical networks from time-series experimental data quite difficult. In this paper, we present a network reconstruction algorithm that can deal with ODE model descriptions containing polynomial and rational functions. Rather than identifying the parameters of linear or nonlinear ODEs characterised by pre-defined equation structures, our methodology allows us to determine the nonlinear ODEs structure together with their associated parameters. To solve the network reconstruction problem, we cast it as a compressive sensing (CS) problem and use sparse Bayesian learning (SBL) algorithms as a computationally efficient and robust way to obtain its solution. © 2012 IEEE.

An organizing center in a planar model of neuronal excitability

This paper studies the excitability properties of a generalized FitzHugh-Nagumo model. The model differs from the classical FitzHugh-Nagumo model in that it accounts for the effect of cooperative gating variables such as activation of calcium currents. Excitability is explored by unfolding a pitchfork bifurcation that is shown to organize five different types of excitability. In addition to the three classical types of neuronal excitability, two novel types are described and distinctly associated to the presence of cooperative variables. © 2012 Society for Industrial and Applied Mathematics.

Effect of95 241Am(n, γ) reaction branching ratio on fuel cycle and reactor design characteristics

The growing interest in innovative reactors and advanced fuel cycle designs requires more accurate prediction of various transuranic actinide concentrations during irradiation or following discharge because of their effect on reactivity or spent-fuel emissions, such as gamma and neutron activity and decay heat. In this respect, many of the important actinides originate from the 241Am(n,γ) reaction, which leads to either the ground or the metastable state of 242Am. The branching ratio for this reaction depends on the incident neutron energy and has very large uncertainty in the current evaluated nuclear data files. This study examines the effect of accounting for the energy dependence of the 241Am(n,γ) reaction branching ratio calculated from different evaluated data files for different reactor and fuel types on the reactivity and concentrations of some important actinides. The results of the study confirm that the uncertainty in knowing the 241Am(n,γ) reaction branching ratio has a negligible effect on the characteristics of conventional light water reactor fuel. However, in advanced reactors with large loadings of actinides in general, and 241Am in particular, the branching ratio data calculated from the different data files may lead to significant differences in the prediction of the fuel criticality and isotopic composition. Moreover, it was found that neutron energy spectrum weighting of the branching ratio in each analyzed case is particularly important and may result in up to a factor of 2 difference in the branching ratio value. Currently, most of the neutronic codes have a single branching ratio value in their data libraries, which is sometimes difficult or impossible to update in accordance with the neutron spectrum shape for the analyzed system.

Statistics of reaction progress variable and mixture fraction gradients from DNS of turbulent partially premixed flames

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios = 0.7 and 1.0 using three-dimensional, simplified chemistry-based direct numerical simulations (DNS). The simulation parameters are chosen such that the flames represent the thin reaction zones regime combustion. A random bimodal distribution of equivalence ratio is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity. The results suggest that the probability density functions (PDFs) of the mixture fraction gradient magnitude |Δξ| (i.e., P(|Δξ|)) can be reasonably approximated using a log-normal distribution. However, this presumed PDF distribution captures only the qualitative nature of the PDF of the reaction progress variable gradient magnitude |Δc| (i.e., P(|Δc|)). It has been found that a bivariate log-normal distribution does not sufficiently capture the quantitative behavior of the joint PDF of |Δξ| and |Δc| (i.e., P(|Δξ|, |Δc|)), and the agreement with the DNS data has been found to be poor in certain regions of the flame brush, particularly toward the burned gas side of the flame brush. Moreover, the variables |Δξ| and |Δc| show appreciable correlation toward the burned gas side of the flame brush. These findings are corroborated further using a DNS data of a lifted jet flame to study the flame geometry dependence of these statistics. © 2013 Copyright Taylor and Francis Group, LLC.

The effect of clearance on shrouded and unshrouded turbines at two levels of reaction

In this paper, the effect of seal clearance on the efficiency of a turbine with a shrouded rotor is compared with the effect of the tip clearance when the same turbine has an unshrouded rotor. The shrouded versus unshrouded comparison was undertaken for two turbine stage designs one having 50% reaction, the other having 24% reaction. Measurements for a range of clearances, including very small clearances, showed three important phenomena. Firstly, as the clearance is reduced, there is a "break-even clearance" at which both the shrouded turbine and the unshrouded turbine have the same efficiency. If the clearance is reduced further, the unshrouded turbine performs better than the shrouded turbine, with the difference at zero clearance termed the "offset loss." This is contrary to the traditional assumption that both shrouded and unshrouded turbines have the same efficiency at zero clearance. The physics of the break-even clearance and the offset loss are discussed. Secondly, the use of a lower reaction had the effect of reducing the tip leakage efficiency penalty for both the shrouded and the unshrouded turbines. In order to understand the effect of reaction on the tip leakage, an analytical model was used and it was found that the tip leakage efficiency penalty should be understood as the dissipated kinetic energy rather than either the tip leakage mass flow rate or the tip leakage loss coefficient. Thirdly, it was also observed that, at a fixed flow coefficient, the fractional change in the output power with clearance was approximately twice the fractional change in efficiency with clearance. This was explained by using an analytical model. © 2014 by ASME.

Direct epitaxial growth of θ-Ni2Si by reaction of a thin Ni(10 at.% Pt) film with Si(1 0 0) substrate

The reaction between an 11 nm Ni(10 at.% Pt) film on a Si substrate has been examined by in situ X-ray diffraction (XRD), atom probe tomography (APT) and transmission electron microscopy (TEM). In situ XRD experiments show the unusual formation of a phase without an XRD peak through consumption of the metal. According to APT, this phase has an Si concentration gradient in accordance with the θ-Ni2Si metastable phase. TEM analysis confirms the direct formation of θ-Ni2Si in epitaxy on Si(1 0 0) with two variants of the epitaxial relationship. © 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.