A tool for predicting the dynamic response of biotrickling filters for VOC removal

This article presents the development of a MATLAB® computer program to simulate the performance of biotrickling filters. Since these filters behave differently during spraying and nonspraying cycles, the presented simulation tool is built on top of a mathematical description of each situation. The resulting variable-structure model is then used as the basis for simulation experiments. The model presented herein represents the first attempt to take into account the variable spraying pattern usually found in industrial installations. Overall, the software is flexible and easy to use, allowing the user to specify the emission concentration pattern, the gas concentration pattern, as well as the spraying cycle periods for up to two different emission patterns per day. The model is able to predict experimental data from a biotrickling filter treating isopropanol under intermittent conditions of loading and spraying. Simulation examples are then provided to study the effect of variable inlet concentrations and gas flow rates

Hartree-Fock-Bogoliubov Model and Simulation of Attractive and Repulsive Bose-Einstein Condensates

We describe a model of dynamic Bose-Einstein condensates near a Feshbach resonance that is computationally feasible under assumptions of spherical or cylindrical symmetry. Simulations in spherical symmetry approximate the experimentally measured time to collapse of an unstably attractive condensate only when the molecular binding energy in the model is correct, demonstrating that the quantum fluctuations and atom-molecule pairing included in the model are the dominant mechanisms during collapse. Simulations of condensates with repulsive interactions find some quantitative disagreement, suggesting that pairing and quantum fluctuations are not the only significant factors for condensate loss or burst formation. Inclusion of three-body recombination was found to be inconsequential in all of our simulations, though we do not consider recent experiments [1] conducted at higher densities

Method of lines solution to the transient SBS equations for nanosecond Stokes pulses

The spectral and temporal evolution of distributed sensing based on stimulated Brillouin scattering (SBS) in optical fibers for severalnanosecondStokes pulses is demonstrated by using the method of lines (MOL) solution of the transient SBS equations. A superbee fluxlimiter is utilized to avoid numerical damping and dispersion that would otherwise be brought on by the approximation of spatial derivativesassociated with the partial differential equations (PDEs). In order to increase computational efficiency, an approach is adopted wherebythe sparse PDE Jacobian matrix integrator option of the ODE solver(s) is employed. Simulation examples of SBS-based sensing for fiberscontaining sections with different Brillouin frequencies are presented. To the best of our knowledge, this MOL solution is proposed for thefirst time for modeling of the transient SBS equations for nanosecond Stokes pulses with different waveforms in a SBS based fiber opticsensor

Relativistic Wavepackets in Classically Chaotic Quantum Cosmological Billiards

Close to a spacelike singularity, pure gravity and supergravity in four to eleven spacetime dimensions admit a cosmological billiard description based on hyperbolic Kac-Moody groups. We investigate the quantum cosmological billiards of relativistic wavepackets towards the singularity, employing flat and hyperbolic space descriptions for the quantum billiards. We find that the strongly chaotic classical billiard motion of four-dimensional pure gravity corresponds to a spreading wavepacket subject to successive redshifts and tending to zero as the singularity is approached. We discuss the possible implications of these results in the context of singularity resolution and compare them with those of known semiclassical approaches. As an aside, we obtain exact solutions for the one-dimensional relativistic quantum billiards with moving walls.Comment: 18 pages, 10 figure

Radiolysis of NaCl at high and low temperatures: development of size distribution of bubbles and colloids

New experimental results are presented on low temperature irradiation (18 °C) of rock-salt samples which had been exposed to initial doses up to 320 GRad at 100 °C. Differential scanning calorimetry (DSC) shows that the latent heat of melting (LHM) of sodium colloids decreases during subsequent low-temperature irradiation, whereas the stored energy (SE) increases slowly, indicating that the process of radiolysis continues. The decrease of the LHM is due to dissolution of large colloids, because the intensities of the melting peaks decrease during the second stage irradiation at low temperature. The model is formulated to describe the nucleation kinetics and the evolution of the size distribution of chlorine precipitates and sodium colloids in NaCl under high dose irradiation. It is shown that the mechanism of dissolution of large Na colloids during low temperature irradiation can be related to melting of sodium colloids.

Role of the Kelvin-Helmholtz instability in the evolution of magnetized relativistic sheared plasma flows

We explore, via analytical and numerical methods, the Kelvin-Helmholtz (KH) instability in relativistic magnetized plasmas, with applications to astrophysical jets. We solve the single-fluid relativistic magnetohydrodynamic (RMHD) equations in conservative form using a scheme which is fourth order in space and time. To recover the primitive RMHD variables, we use a highly accurate, rapidly convergent algorithm which improves upon such schemes as the Newton-Raphson method. Although the exact RMHD equations are marginally stable, numerical discretization renders them unstable. We include numerical viscosity to restore numerical stability. In relativistic flows, diffusion can lead to a mathematical anomaly associated with frame transformations. However, in our KH studies, we remain in the rest frame of the system, and therefore do not encounter this anomaly. We use a two-dimensional slab geometry with periodic boundary conditions in both directions. The initial unperturbed velocity peaks along the central axis and vanishes asymptotically at the transverse boundaries. Remaining unperturbed quantities are uniform, with a flow-aligned unperturbed magnetic field. The early evolution in the nonlinear regime corresponds to the formation of counter-rotating vortices, connected by filaments, which persist in the absence of a magnetic field. A magnetic field inhibits the vortices through a series of stages, namely, field amplification, vortex disruption, turbulent breakdown, and an approach to a flow-aligned equilibrium configuration. Similar stages have been discussed in MHD literature. We examine how and to what extent these stages manifest in RMHD for a set of representative field strengths. To characterize field strength, we define a relativistic extension of the Alfvénic Mach number M(A). We observe close complementarity between flow and magnetic field behavior. Weaker fields exhibit more vortex rotation, magnetic reconnection, jet broadening, and intermediate turbulence. Sufficiently strong fields (M(A)<6) completely suppress vortex formation. Maximum jet deceleration, and viscous dissipation, occur for intermediate vortex-disruptive fields, while electromagnetic energy is maximized for the strongest fields which allow vortex formation. Highly relativistic flows destabilize the system, supporting modes with near-maximum growth at smaller wavelengths than the shear width of the velocity. This helps to explain early numerical breakdown of highly relativistic simulations using numerical viscosity, a long-standing problem. While magnetic fields generally stabilize the system, we have identified many features of the complex and turbulent reorganization that occur for sufficiently weak fields in RMHD flows, and have described the transition from disruptive to stabilizing fields at M(A)≈6. Our results are qualitatively similar to observations of numerous jets, including M87, whose knots may exhibit vortex-like behavior. Furthermore, in both the linear and nonlinear analyses, we have successfully unified the HD, MHD, RHD, and RMHD regimes

INSGFP/w human embryonic stem cells facilitate isolation of in vitro derived insulin-producing cells

AIMS/HYPOTHESIS: We aimed to generate human embryonic stem cell (hESC) reporter lines that would facilitate the characterisation of insulin-producing (INS⁺) cells derived in vitro. METHODS: Homologous recombination was used to insert sequences encoding green fluorescent protein (GFP) into the INS locus, to create reporter cell lines enabling the prospective isolation of viable INS⁺ cells. RESULTS: Differentiation of INS(GFP/w) hESCs using published protocols demonstrated that all GFP⁺ cells co-produced insulin, confirming the fidelity of the reporter gene. INS-GFP⁺ cells often co-produced glucagon and somatostatin, confirming conclusions from previous studies that early hESC-derived insulin-producing cells were polyhormonal. INS(GFP/w) hESCs were used to develop a 96-well format spin embryoid body (EB) differentiation protocol that used the recombinant protein-based, fully defined medium, APEL. Like INS-GFP⁺ cells generated with other methods, those derived using the spin EB protocol expressed a suite of pancreatic-related transcription factor genes including ISL1, PAX6 and NKX2.2. However, in contrast with previous methods, the spin EB protocol yielded INS-GFP⁺ cells that also co-expressed the beta cell transcription factor gene, NKX6.1, and comprised a substantial proportion of monohormonal INS⁺ cells. CONCLUSIONS/INTERPRETATION: INS(GFP/w) hESCs are a valuable tool for investigating the nature of early INS⁺ progenitors in beta cell ontogeny and will facilitate the development of novel protocols for generating INS⁺ cells from differentiating hESCs

Geometric numerical schemes for the KdV equation

Geometric discretizations that preserve certain Hamiltonian structures at the discrete level has been proven to enhance the accuracy of numerical schemes. In particular, numerous symplectic and multi-symplectic schemes have been proposed to solve numerically the celebrated Korteweg-de Vries (KdV) equation. In this work, we show that geometrical schemes are as much robust and accurate as Fourier-type pseudo-spectral methods for computing the long-time KdV dynamics, and thus more suitable to model complex nonlinear wave phenomena.Comment: 22 pages, 14 figures, 74 references. Other author's papers can be downloaded at http://www.lama.univ-savoie.fr/~dutykh