The decomposition of some RDX and HMX based materials in the one dimensional time to explosion apparatus. Part 2. Estimating the violence of the cook-off event

Various methods of assessment have been applied to the One Dimensional Time to Explosion (ODTX) apparatus and experiments with the aim of allowing an estimate of the comparative violence of the explosion event to be made. Non-mechanical methods used were a simple visual inspection, measuring the increase in the void volume of the anvils following an explosion and measuring the velocity of the sound produced by the explosion over 1 metre. Mechanical methods used included monitoring piezo-electric devices inserted in the frame of the machine and measuring the rotational velocity of a rotating bar placed on the top of the anvils after it had been displaced by the shock wave. This last method, which resembles original Hopkinson Bar experiments, seemed the easiest to apply and analyse, giving relative rankings of violence and the possibility of the calculation of a “detonation” pressure.

Bistable systems with stochastic noise: virtues and limits of effective one-dimensional Langevin equations

The understanding of the statistical properties and of the dynamics of multistable systems is gaining more and more importance in a vast variety of scientific fields. This is especially relevant for the investigation of the tipping points of complex systems. Sometimes, in order to understand the time series of given observables exhibiting bimodal distributions, simple one-dimensional Langevin models are fitted to reproduce the observed statistical properties, and used to investing-ate the projected dynamics of the observable. This is of great relevance for studying potential catastrophic changes in the properties of the underlying system or resonant behaviours like those related to stochastic resonance-like mechanisms. In this paper, we propose a framework for encasing this kind of studies, using simple box models of the oceanic circulation and choosing as observable the strength of the thermohaline circulation. We study the statistical properties of the transitions between the two modes of operation of the thermohaline circulation under symmetric boundary forcings and test their agreement with simplified one-dimensional phenomenological theories. We extend our analysis to include stochastic resonance-like amplification processes. We conclude that fitted one-dimensional Langevin models, when closely scrutinised, may result to be more ad-hoc than they seem, lacking robustness and/or well-posedness. They should be treated with care, more as an empiric descriptive tool than as methodology with predictive power.

Quasilimiting behavior for one-dimensional diffusions with killing

This paper extends and clarifies results of Steinsaltz and Evans [Trans. Amer. Math. Soc. 359 (2007) 1285–1234], which found conditions for convergence of a killed one-dimensional diffusion conditioned on survival, to a quasistationary distribution whose density is given by the principal eigenfunction of the generator. Under the assumption that the limit of the killing at infinity differs from the principal eigenvalue we prove that convergence to quasistationarity occurs if and only if the principal eigenfunction is integrable. When the killing at ∞ is larger than the principal eigenvalue, then the eigenfunction is always integrable. When the killing at ∞ is smaller, the eigenfunction is integrable only when the unkilled process is recurrent; otherwise, the process conditioned on survival converges to 0 density on any bounded interval.

The unsteady flow of a weakly compressible fluid in a thin porous layer II: three-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a three-dimensional layer, composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Numerical solution of this three-dimensional evolution problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l, a situation which occurs frequently in the application to oil and gas reservoir recovery and which leads to significant stiffness in the numerical problem. Under the assumption that $\epsilon\propto h/l\ll 1$, we show that, to leading order in $\epsilon$, the pressure field varies only in the horizontal directions away from the wells (the outer region). We construct asymptotic expansions in $\epsilon$ in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive expressions for all significant process quantities. The only computations required are for the solution of non-stiff linear, elliptic, two-dimensional boundary-value, and eigenvalue problems. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the layer, $\epsilon$, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighbourhood of wells and away from wells.

Comparison of polar ozone loss rates simulated by one-dimensional and three-dimensional models with Match observations in recent Antarctic and Arctic winters

Simulations of ozone loss rates using a three-dimensional chemical transport model and a box model during recent Antarctic and Arctic winters are compared with experimental loss rates. The study focused on the Antarctic winter 2003, during which the first Antarctic Match campaign was organized, and on Arctic winters 1999/2000, 2002/2003. The maximum ozone loss rates retrieved by the Match technique for the winters and levels studied reached 6 ppbv/sunlit hour and both types of simulations could generally reproduce the observations at 2-sigma error bar level. In some cases, for example, for the Arctic winter 2002/2003 at 475 K level, an excellent agreement within 1-sigma standard deviation level was obtained. An overestimation was also found with the box model simulation at some isentropic levels for the Antarctic winter and the Arctic winter 1999/2000, indicating an overestimation of chlorine activation in the model. Loss rates in the Antarctic show signs of saturation in September, which have to be considered in the comparison. Sensitivity tests were performed with the box model in order to assess the impact of kinetic parameters of the ClO-Cl2O2 catalytic cycle and total bromine content on the ozone loss rate. These tests resulted in a maximum change in ozone loss rates of 1.2 ppbv/sunlit hour, generally in high solar zenith angle conditions. In some cases, a better agreement was achieved with fastest photolysis of Cl2O2 and additional source of total inorganic bromine but at the expense of overestimation of smaller ozone loss rates derived later in the winter.

Scattering by infinite one-dimensional rough surfaces

We consider the Dirichlet boundary-value problem for the Helmholtz equation in a non-locally perturbed half-plane. This problem models time-harmonic electromagnetic scattering by a one-dimensional, infinite, rough, perfectly conducting surface; the same problem arises in acoustic scattering by a sound-soft surface. ChandlerWilde & Zhang have suggested a radiation condition for this problem, a generalization of the Rayleigh expansion condition for diffraction gratings, and uniqueness of solution has been established. Recently, an integral equation formulation of the problem has also been proposed and, in the special case when the whole boundary is both Lyapunov and a small perturbation of a flat boundary, the unique solvability of this integral equation has been shown by Chandler-Wilde & Ross by operator perturbation arguments. In this paper we study the general case, with no limit on surface amplitudes or slopes, and show that the same integral equation has exactly one solution in the space of bounded and continuous functions for all wavenumbers. As an important corollary we prove that, for a variety of incident fields including the incident plane wave, the Dirichlet boundary-value problem for the scattered field has a unique solution.

Eigenvalues of a one-dimensional Dirac operator pencil

We study the spectrum of a one-dimensional Dirac operator pencil, with a coupling constant in front of the potential considered as the spectral parameter. Motivated by recent investigations of graphene waveguides, we focus on the values of the coupling constant for which the kernel of the Dirac operator contains a square integrable function. In physics literature such a function is called a confined zero mode. Several results on the asymptotic distribution of coupling constants giving rise to zero modes are obtained. In particular, we show that this distribution depends in a subtle way on the sign variation and the presence of gaps in the potential. Surprisingly, it also depends on the arithmetic properties of certain quantities determined by the potential. We further observe that variable sign potentials may produce complex eigenvalues of the operator pencil. Some examples and numerical calculations illustrating these phenomena are presented.

Solvothermal synthesis of one-dimensional chalcogenides containing group 13 elements

The solvothermal synthesis and characterisation of [C6H16N2][GaS2]2 (1), [C6H16N2][Ga2Se3(Se2)] (2), and mixed-metal phases with composition [C6H16N2][Ga2–xInxSe3(Se2)] (0 < x < 2)(3–5), is described. These materials have been characterised by single-crystal and powder X-ray diffraction, thermogravimetric analysis and UV/Vis diffuse reflectance spectroscopy. The materials contain one-dimensional anionic chains. In 1, these chains consist of edge-linked GaS4 tetrahedra, whilst in 2–5, the chains contain perselenide (Se2)2– units and comprise alternating four-membered [M2Se2] and five-membered [M2Se3] rings (where M = Ga, In). Compounds 3–5 represent the first examples of ternary mixed-metal [M2Se3(Se2)]2– chains.

Hydrothermal synthesis of [C6H16N2][In2Se3(Se2)]: a new one-dimensional indium selenide

A new organically templated indium selenide, [C6H16N2][In2Se3(Se2)], has been prepared hydrothermally from the reaction of indium, selenium and trans-1,4-diaminocyclohexane in water at 170 °C. This material was characterised by single-crystal and powder X-ray diffraction, thermogravimetric analysis, UV–vis diffuse reflectance spectroscopy, FT-IR and elemental analysis. The compound crystallises in the monoclinic space group C2/c (a=12.0221(16) Å, b=11.2498(15) Å, c=12.8470(17) Å, β=110.514(6)°). The crystal structure of [C6H16N2][In2Se3(Se2)] contains anionic chains of stoichiometry [In2Se3(Se2)]2−, which are aligned parallel to the [1 0 1] direction, and separated by diprotonated trans-1,4-diaminocyclohexane cations. The [In2Se3(Se2)]2− chains, which consist of alternating four-membered [In2Se2] and five-membered [In2Se3] rings, contain perselenide (Se2)2− units. UV–vis diffuse reflectance spectroscopy indicates that [C6H16N2][In2Se3(Se2)] has a band gap of 2.23(1) eV

Synthesis, characterisation and magnetic properties of a one-dimensional Iron(II) coordination polymer

A new iron(II) coordination polymer, [FeCl2(NC7H9)2(N2C12H12)], has been synthesized under solvothermal conditions and structurally characterized by single-crystal X-ray diffraction. This material crystallizes in the monoclinic space group C2/c, with a = 11.2850(6), b = 13.8925(7), c = 17.0988(9) Å and β = 94.300(3)º (Z = 4). The crystal structure consists of neutral zig-zag chains, in which the iron(II) ions are octahedrally coordinated. The infinite polymer chains are packed into a three-dimensional structure through C–H···Cl interactions. Magnetic susceptibility measurements reveal the existence of weak antiferromagnetic interactions between the iron(II) ions. The effective magnetic moment, μ eff = 5.33 μ B , is consistent with a high-spin iron(II) configuration.

Coupling an aerosol box model with one-dimensional flow: a tool for understanding observations of new particle formation events

Field observations of new particle formation and the subsequent particle growth are typically only possible at a fixed measurement location, and hence do not follow the temporal evolution of an air parcel in a Lagrangian sense. Standard analysis for determining formation and growth rates requires that the time-dependent formation rate and growth rate of the particles are spatially invariant; air parcel advection means that the observed temporal evolution of the particle size distribution at a fixed measurement location may not represent the true evolution if there are spatial variations in the formation and growth rates. Here we present a zero-dimensional aerosol box model coupled with one-dimensional atmospheric flow to describe the impact of advection on the evolution of simulated new particle formation events. Wind speed, particle formation rates and growth rates are input parameters that can vary as a function of time and location, using wind speed to connect location to time. The output simulates measurements at a fixed location; formation and growth rates of the particle mode can then be calculated from the simulated observations at a stationary point for different scenarios and be compared with the ‘true’ input parameters. Hence, we can investigate how spatial variations in the formation and growth rates of new particles would appear in observations of particle number size distributions at a fixed measurement site. We show that the particle size distribution and growth rate at a fixed location is dependent on the formation and growth parameters upwind, even if local conditions do not vary. We also show that different input parameters used may result in very similar simulated measurements. Erroneous interpretation of observations in terms of particle formation and growth rates, and the time span and areal extent of new particle formation, is possible if the spatial effects are not accounted for.

STUDY OF THE EQUATORIAL ATLANTIC OCEAN MIXING LAYER USING A ONE-DIMENSIONAL TURBULENCE MODEL

The General Ocean Turbulence Model (GOTM) is applied to the diagnostic turbulence field of the mixing layer (ML) over the equatorial region of the Atlantic Ocean. Two situations were investigated: rainy and dry seasons, defined, respectively, by the presence of the intertropical convergence zone and by its northward displacement. Simulations were carried out using data from a PIRATA buoy located on the equator at 23 degrees W to compute surface turbulent fluxes and from the NASA/GEWEX Surface Radiation Budget Project to close the surface radiation balance. A data assimilation scheme was used as a surrogate for the physical effects not present in the one-dimensional model. In the rainy season, results show that the ML is shallower due to the weaker surface stress and stronger stable stratification; the maximum ML depth reached during this season is around 15 m, with an averaged diurnal variation of 7 m depth. In the dry season, the stronger surface stress and the enhanced surface heat balance components enable higher mechanical production of turbulent kinetic energy and, at night, the buoyancy acts also enhancing turbulence in the first meters of depth, characterizing a deeper ML, reaching around 60 m and presenting an average diurnal variation of 30 m.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Heuristics for the one-dimensional cutting stock problem with limited multiple stock lengths

This paper deals with the classical one-dimensional integer cutting stock problem, which consists of cutting a set of available stock lengths in order to produce smaller ordered items. This process is carried out in order to optimize a given objective function (e.g., minimizing waste). Our study deals with a case in which there are several stock lengths available in limited quantities. Moreover, we have focused on problems of low demand. Some heuristic methods are proposed in order to obtain an integer solution and compared with others. The heuristic methods are empirically analyzed by solving a set of randomly generated instances and a set of instances from the literature. Concerning the latter. most of the optimal solutions of these instances are known, therefore it was possible to compare the solutions. The proposed methods presented very small objective function value gaps. (C) 2008 Elsevier Ltd. All rights reserved.