The Y-procedure: How to extract the chemical transformation rate from reaction-diffusion data with no assumption on the kinetic model

The paper presents a new method to extract the chemical transformation rate from reaction–diffusion data with no assumption on the kinetic model (“kinetic model-free procedure”). It is a new non-steady-state kinetic characterization procedure for heterogeneous catalysts. The mathematical foundation of the Y-procedure is a Laplace-domain analysis of the two inert zones in a TZTR followed by transposition to the Fourier domain. When combined with time discretization and filtering the Y-procedure leads to an efficient practical method for reconstructing the concentration and reaction rate in the active zone. Using the Y-procedure the concentration and reaction rate of a non-steady state catalytic process can be determined without any pre-assumption regarding the type of kinetic dependence. The Y-procedure is the basis for advanced software for non-steady state kinetic data interpretation. The Y-procedure can be used to relate changes in the catalytic reaction rate and kinetic parameters to changes in the surface composition (storage) of a catalyst.

Numerical simulations of fluid-structure interactions in single-reed mouthpieces

Most single-reed woodwind instrument models rely on a quasistationary approximation to describe the relationship between the volume flow and. the pressure difference across the reed channel. Semiempirical models based on the quasistationary approximation are very useful in explaining the fundamental characteristics of this family of instruments such as self-sustained oscillations and threshold of blowing pressure. However, they fail at explaining more complex phenomena associated with the fluid-structure interaction during dynamic flow regimes, such as the transient and steady-state behavior of the system as a function. of the mouthpiece geometry. Previous studies have discussed the accuracy of the quasistationary approximation but the amount of literature on the subject is sparse, mainly due to the difficulties involved in the measurement of dynamic flows in channels with an oscillating reed. In this paper, a numerical technique based on the lattice Boltzmann method and a finite difference scheme is proposed in order to investigate the characteristics of fully coupled fluid-structure interaction in single-reed mouthpieces with different channel configurations. Results obtained for a stationary simulation with a static reed agree very well with those predicted by the literature based on the quasistationary approximation. However, simulations carried out for a dynamic regime with dn oscillating reed show that the phenomenon associated with flow detachment and reattachment diverges considerably frorn the theoretical assumptions. Furthermore, in the case of long reed channels, the results obtained for the vena contracta factor are in significant disagreement with those predicted by theory. For short channels, the assumption of constant vena contracta was found to be valid for only 40% of the duty cycle. (c) 2007 Acoustical Society of America.

Dynamics of Gompertzian tumour growth under environmental fluctuations

We investigate the effect of correlated additive and multiplicative Gaussian white noise oil the Gompertzian growth of tumours. Our results are obtained by Solving numerically the time-dependent Fokker-Planck equation (FPE) associated with the stochastic dynamics. In Our numerical approach we have adopted B-spline functions as a truncated basis to expand the approximated eigenfunctions. The eigenfunctions and eigenvalues obtained using this method are used to derive approximate solutions of the dynamics under Study. We perform simulations to analyze various aspects, of the probability distribution. of the tumour cell populations in the transient- and steady-state regimes. More precisely, we are concerned mainly with the behaviour of the relaxation time (tau) to the steady-state distribution as a function of (i) of the correlation strength (lambda) between the additive noise and the multiplicative noise and (ii) as a function of the multiplicative noise intensity (D) and additive noise intensity (alpha). It is observed that both the correlation strength and the intensities of additive and multiplicative noise, affect the relaxation time.

Transonic Aeroelastic Stability Analysis Using a Kriging-Based Schur Complement Formulation

A method is described to allow searches for transonic aeroelastic instability of realistically sized aircraft models in multidimensional parameter spaces when computational fluid dynamics are used to model the aerodynamics. Aeroelastic instability is predicted from a small nonlinear eigenvalue problem. The approximation of the computationally expensive interaction term modeling the fluid response is formulated to allow the automated and blind search for aeroelastic instability. The approximation uses a kriging interpolation of exact numerical samples covering the parameter space. The approach, demonstrated for the Goland wing and the multidisciplinary optimization transport wing, results in stability analyses over whole flight envelopes at an equivalent cost of several steady-state simulations.

Predicting die life from die temperature for high pressure dies casting

The objective of this research was to determine the surface temperature of a high pressure die casting die during casting conditions. This was achieved by instrumentation of an insert which was placed in the shotplate region of the die. This research overcame the challenge of directly measuring the die surface temperature during a HPDC production casting cycle and shows that this is an effective method to determine the die surface temperature during the casting cycle. The instrumentation results gave a peak and minimum temperature of 500 C and 240 C respectively during steady state running conditions with a molten aluminium casting temperature of 660 C. Stress analysis from the steady state measured temperature of the die surface was calculated through a simple FEA model and the resulting stress uctuation was applied to a fatigue equation for the die material, the predicted number of cycles for cracking to start was found to correlate well with observed die damage.

Transient Heat Conduction of Variable Stiffness Composite Laminate

Recent research on Variable Stiffness (VS) laminates, which are constructed by steering the fiber orientation as a spatial function of location, have shown to improve laminate performance under mechanical loads. Two distinct cases of stiffness variation can be achieved either by variation of the fiber orientation in the direction of the global x-axis, or perpendicular to it. In the present paper, thermal analysis of a VS laminate is performed to study the effect of steering fibers on transient heat conduction under uniform heat flux using finite element method. The goal of the present paper is a parametric study of the effect of variable stiffness properties on transient response including time to reach steady state and temperature profile. Also, stress resultants and maximum stress location are investigated under different boundary conditions. A FEM algorithm is applied to exactly incorporate the boundary conditions for stress resultant analysis.

Transient Heat Conduction of Variable Stiffness Composite Laminate

Recent research on Variable Stiffness (VS) laminates, which are constructed by steering the fiber orientation as a spatial function of location, have shown to improve laminate performance under mechanical loads. Two distinct cases of stiffness variation can be achieved either by variation of the fiber orientation in the direction of the global x-axis, or perpendicular to it. In the present paper, thermal analysis of VS laminate is performed to study the effect of steering fibers on transient heat conduction under uniform heat flux using finite element method. The goal of the present paper is a parametric study of the
effect of variable stiffness properties on transient response including time to reach steady state and temperature profile. Also, stress resultants and maximum stress location are investigated under different boundary conditions. A FEM algorithm is applied to exactly incorporate the boundary conditions.

Functional characterization and membrane topology of Escherichia coli WecA, a sugar-phosphate transferase initiating the biosynthesis of enterobacterial common antigen and O-antigen lipopolysaccharide

WecA is an integral membrane protein that initiates the biosynthesis of enterobacterial common antigen and O-antigen lipopolysaccharide (LPS) by catalyzing the transfer of N-acetylglucosamine (GlcNAc)-1-phosphate onto undecaprenyl phosphate (Und-P) to form Und-P-P-GlcNAc. WecA belongs to a large family of eukaryotic and prokaryotic prenyl sugar transferases. Conserved aspartic acids in putative cytoplasmic loops 2 (Asp90 and Asp91) and 3 (Asp156 and Asp159) were targeted for replacement mutagenesis with either glutamic acid or asparagine. We examined the ability of each mutant protein to complement O-antigen LPS synthesis in a wecA-deficient strain and also determined the steady-state kinetic parameters of the mutant proteins in an in vitro transfer assay. Apparent K(m) and V(max) values for UDP-GlcNAc, Mg(2+), and Mn(2+) suggest that Asp156 is required for catalysis, while Asp91 appears to interact preferentially with Mg(2+), possibly playing a role in orienting the substrates. Topological analysis using the substituted cysteine accessibility method demonstrated the cytosolic location of Asp90, Asp91, and Asp156 and provided a more refined overall topological map of WecA. Also, we show that cells expressing a WecA derivative C terminally fused with the green fluorescent protein exhibited a punctate distribution of fluorescence on the bacterial surface, suggesting that WecA localizes to discrete regions in the bacterial plasma membrane.

Field test of multiple 1/10 scale tidal turbines in steady flows

Queen's University Belfast and Wave Barrier Ltd have developed a tidal testing platform to test hydrokinetic turbines at medium scale. Multiple turbines can be pushed through still water conditions, in steady-state pushing tests. Experiments were conducted to evaluate the interactions between two identical, mono-strut, horizontal axis tidal turbines (HATTs) of 1.5 m diameter (D) rotor. Their relative performance when located individually, in-plane and in-line are investigated. The data shows a high consistency in the power curves at different flow speeds, which indicates high repeatability in this Reynolds range. For an individual turbine, there is no performance difference when the rotor is mounted either upstream or downstream of the supporting structure. When placed in-plane, the turbines have no adverse effect on one another. When spaced in-line with 2D separation, there is a 63% reduction in the performance of the downstream turbine. At 6D downstream this performance reduction is still 59%, indicating some wake recovery between 2D and 6D, though the influence from the upstream rotor persists to at least 6D downstream of the first device. In contrast the performance of the downstream turbine when placed at 1.5D offset of the upstream device at 6D downstream is approximately recovered to the individual turbine performance.

Modelling of Organic Rankine Cycle for Waste Heat Recovery Process in Supercritical Condition

Organic Rankine Cycle (ORC) is the most commonly used method for recovering energy from small sources of heat. The investigation of the ORC in supercritical condition is a new research area as it has a potential to generate high power and thermal efficiency in a waste heat recovery system. This paper presents a steady state ORC model in supercritical condition and its simulations with a real engine’s exhaust data. The key component of ORC, evaporator, is modelled using finite volume method, modelling of all other components of the waste heat recovery system such as pump, expander and condenser are also presented. The aim of this paper is to investigate the effects of mass flow rate and evaporator outlet temperature on the efficiency of the waste heat recovery process. Additionally, the necessity of maintaining an optimum evaporator outlet temperature is also investigated. Simulation results show that modification of mass flow rate is the key to changing the operating temperature at the evaporator outlet.

Study on Hybrid Heat-Gas-Power Flow Algorithm for Integrated Community Energy System: 区域综合能源系统电/气/热混合潮流算法研究

With the integration of combined heat and power (CHP) units, air-conditioners and gas boilers, power, gas, and heat systems are becoming tightly linked to each other in the integrated community energy system (ICES). Interactions among the three systems are not well captured by traditional methods. To address this issue, a hybrid power-gas-heat flow calculation method was developed in this paper. In the proposed method, an energy hub model was presented to describe interactions among the three systems incorporating various CHP operating modes. In addition, three operating modes were proposed for the ICES including fully decoupled, partially coupled, and fully coupled. Numerical results indicated that the proposed algorithm can be used in the steady-state analysis of the ICES and reflect interactions among various energy systems.

The Type IIb SN 2011dh: Two years of observations and modelling of the lightcurves

We present optical and near-infrared (NIR) photometry and spectroscopy as well as modelling of the lightcurves of the Type IIb supernova (SN) 2011dh. Our extensive dataset, for which we present the observations obtained after day 100, spans two years, and complemented with Spitzer mid-infrared (MIR) data, we use it to build an optical-to-MIR bolometric lightcurve between days 3 and 732. To model the bolometric lightcurve before day 400 we use a grid of hydrodynamical SN models, which allows us to determine the errors in the derived quantities, and a bolometric correction determined with steady-state non-local thermodynamic equilibrium (NLTE) modelling. Using this method we find a helium core mass of 3.1^+0.7<inf>-0.4</inf> M<inf>⊙</inf> for SN 2011dh, consistent within error bars with previous results obtained using the bolometric lightcurve before day 80. We compute bolometric and broad-band lightcurves between days 100 and 500 from spectral steady-state NLTE models, presented and discussed in a companion paper. The preferred 12 M<inf>⊙</inf> (initial mass) model, previously found to agree well with the observed spectra, shows a good overall agreement with the observed lightcurves, although some discrepancies exist. Time-dependent NLTE modelling shows that after day ∼600 a steady-state assumption is no longer valid. The radioactive energy deposition in this phase is likely dominated by the positrons emitted in the decay of ⁵⁶Co, but seems insufficient to reproduce the lightcurves, and what energy source is dominating the emitted flux is unclear. We find an excess in the K and the MIR bands developing between days 100 and 250, during which an increase in the optical decline rate is also observed. A local origin of the excess is suggested by the depth of the He I 20 581 Å absorption. Steady-state NLTE models with a modest dust opacity in the core (τ = 0.44), turned on during this period, reproduce the observed behaviour, but an additional excess in the Spitzer 4.5 μm band remains. Carbon-monoxide (CO) first-overtone band emission is detected at day 206, and possibly at day 89, and assuming the additional excess to bedominated by CO fundamental band emission, we find fundamental to first-overtone band ratios considerably higher than observed in SN 1987A. The profiles of the [O i] 6300 Å and Mg i] 4571 Å lines show a remarkable similarity, suggesting that these lines originate from a common nuclear burning zone (O/Ne/Mg), and using small scale fluctuations in the line profiles we estimate a filling factor of ≲ 0.07 for the emitting material. This paper concludes our extensive observational and modelling work on SN 2011dh. The results from hydrodynamical modelling, steady-state NLTE modelling, and stellar evolutionary progenitor analysis are all consistent, and suggest an initial mass of ∼12 M<inf>⊙</inf> for the progenitor.

Helium line intensity ratio in microwave-generated plasmas

The line intensity ratio method provides a nonintrusive diagnostic for the measurement of electron temperature in microwave-generated plasmas. For optically thin plasmas of low density, a line intensity method using He I lines can often be used, and is based on the fact that the electron impact excitation rate coefficients for helium singlet and triplet states are insensitive to electron density but differ as a function of the electron temperature. Line intensity measurements are presented from microwave-generated helium plasmas. Both steady-state corona and collision-radiative theoretical models are used to evaluate the ground and excited state populations. The line ratio versus electron temperature obtained from both of these methods are compared with the results from measurements. However, it is not possible to diagnose the electron temperature from the line ratios alone due to the presence of significant opacity and nonnegligible 1s2s S-3 metastable fraction in the plasma. (C) 2004 American Institute of Physics.

Nonequilibrium processes from generalized Langevin equations: Realistic nanoscale systems connected to two thermal baths

We extend the generalized Langevin equation (GLE) method [L. Stella, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nosé-Hoover thermostats). We concentrate on the steady-state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e., ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.

Efficient simulations with electronic open boundaries

We present a reformulation of the hairy-probe method for introducing electronic open boundaries that is appropriate for steady-state calculations involving nonorthogonal atomic basis sets. As a check on the correctness of the method we investigate a perfect atomic wire of Cu atoms and a perfect nonorthogonal chain of H atoms. For both atom chains we find that the conductance has a value of exactly one quantum unit and that this is rather insensitive to the strength of coupling of the probes to the system, provided values of the coupling are of the same order as the mean interlevel spacing of the system without probes. For the Cu atom chain we find in addition that away from the regions with probes attached, the potential in the wire is uniform, while within them it follows a predicted exponential variation with position. We then apply the method to an initial investigation of the suitability of graphene as a contact material for molecular electronics. We perform calculations on a carbon nanoribbon to determine the correct coupling strength of the probes to the graphene and obtain a conductance of about two quantum units corresponding to two bands crossing the Fermi surface. We then compute the current through a benzene molecule attached to two graphene contacts and find only a very weak current because of the disruption of the π conjugation by the covalent bond between the benzene and the graphene. In all cases we find that very strong or weak probe couplings suppress the current.