Nonlinear cellular dynamics of the idealized detonation model: Regular cells

High-resolution numerical simulations of cellular detonations are performed using a parallelized adaptive grid solver, in the case where the channel width is very wide. In particular, the nonlinear response of a weakly unstable ZND detonation to two-dimensional perturbations is studied in the context of the idealized one-step chemistry model. For random perturbations, cells appear with a characteristic size in good agreement with that corresponding to the maximum growth rate from a linear stability analysis. However, the cells then grow and equilibrate at a larger size. It is also shown that the linear analysis predicts well the ratio of cell lengths to cell widths of the fully developed cells. The evolutionary dynamics of the growth are nonetheless quite slow, in that the detonation needs to run of the order of 1000 reaction lengths before the final size and equilibrium state is reached. For sinusoidal perturbations, it is found that there is a large band of wavelengths/cell sizes which can propagate over very long distances (~1000 reaction lengths). By perturbing the fully developed cells of each wavelength, it is found that smaller cells in this range are unstable to symmetry breaking, which again results in cellular growth to a larger final size. However, a range of larger cell sizes appear to be nonlinearly stable. As a result it is found that the final cell size of the model is non-unique, even for such a weakly unstable, regular cell case. Indeed, in the case studied, the equilibrium cell size varies by 100% with different initial conditions. Numerical dependencies of the cellular dynamics are also examined

Validating the FLASH Code: Vortex-Dominated Flows

As a component of the Flash Center's validation program, we compare FLASH simulation results with experimental results from Los Alamos National Laboratory. The flow of interest involves the lateral interaction between a planar Ma=1.2 shock wave with a cylinder of gaseous sulfur hexafluoride (SF_6) in air, and in particular the development of primary and secondary instabilities after the passage of the shock. While the overall evolution of the flow is comparable in the simulations and experiments, small-scale features are difficult to match. We focus on the sensitivity of numerical results to simulation parameters.Comment: 10 pages, 5 figures, presented at the 5th International Conference on High Energy Laboratory Astrophysics, Tucson, AZ, March 10-13, 200

Low Mach number effect in simulation of high Mach number flow

In this note, we relate the two well-known difficulties of Godunov schemes: the carbuncle phenomena in simulating high Mach number flow, and the inaccurate pressure profile in simulating low Mach number flow. We introduced two simple low-Mach-number modifications for the classical Roe flux to decrease the difference between the acoustic and advection contributions of the numerical dissipation. While the first modification increases the local numerical dissipation, the second decreases it. The numerical tests on the double-Mach reflection problem show that both modifications eliminate the kinked Mach stem suffered by the original flux. These results suggest that, other than insufficient numerical dissipation near the shock front, the carbuncle phenomena is strongly relevant to the non-comparable acoustic and advection contributions of the numerical dissipation produced by Godunov schemes due to the low Mach number effect.Comment: 9 pages, 1 figur

Repeat testing of low-level HIV-1 RNA: Assay performance and implementation in clinical trials

Objective: Assess the performance of HIV-1 RNA repeat testing of stored samples in cases of low-level viremia during clinical trials. Design: Prospective and retrospective analysis of randomized clinical trial samples and reference standards. Methods: To evaluate assay variability of the Cobas AmpliPrep/Cobas TaqMan HIV-1 Test, v2.0, three separate sources of samples were utilized: the World Health Organization (WHO) HIV reference standard (assayed using 50 independent measurements at six viral loads <200 copies/ml), retrospective analysis of four to six aliquots of plasma samples from four clinical trial participants, and prospective repeat testing of 120 samples from participants in randomized trials with low-level viremia. Results: The TaqMan assay on the WHO HIV-1 RNA standards at viral loads <200 copies/ml performed within the expected variability according to assay specifications. However, standards with low viral loads of 36 and 18 copies/ml reported values of ≥ 50 copies/ml in 66 and 18% of tests, respectively. In participants treated with antiretrovirals who had unexpected viremia of 50-200 copies/ml after achieving <50 copies/ml, retesting of multiple aliquots of stored plasma found <50 copies/ml in nearly all cases upon retesting (14/15; 93%). Repeat testing was prospectively implemented in four clinical trials for all samples with virologic rebound of 50-200 copies/ml (n = 120 samples from 92 participants) from which 42% (50/120) had a retest result of less than 50 copies/ml and 58% (70/120) retested ≥ 50 copies/ml. Conclusion: The TaqMan HIV-1 RNA assay shows variability around 50 copies/ml that affects clinical trial results and may impact clinical practice. In participants with a history of viral load suppression, unexpected low-level viremia may be because of assay variability rather than low drug adherence or true virologic failure. Retesting a stored aliquot of the same sample may differentiate between assay variability and virologic failure as the source of viremia. This retesting strategy could save time, money, and anxiety for patients and their providers, as well as decrease follow-up clinic visits without increasing the risk of virologic failure and resistance development

Spatial Representation and Navigation in a Bio-inspired Robot

A biologically inspired computational model of rodent repre-sentation?based (locale) navigation is presented. The model combines visual input in the form of realistic two dimensional grey-scale images and odometer signals to drive the firing of simulated place and head direction cells via Hebbian synapses. The space representation is built incrementally and on-line without any prior information about the environment and consists of a large population of location-sensitive units (place cells) with overlapping receptive fields. Goal navigation is performed using reinforcement learning in continuous state and action spaces, where the state space is represented by population activity of the place cells. The model is able to reproduce a number of behavioral and neuro-physiological data on rodents. Performance of the model was tested on both simulated and real mobile Khepera robots in a set of behavioral tasks and is comparable to the performance of animals in similar tasks

AMRITA -- A computational facility

Amrita is a software system for automating numerical investigations. The system is driven using its own powerful scripting language, Amrita, which facilitates both the composition and archiving of complete numerical investigations, as distinct from isolated computations. Once archived, an Amrita investigation can later be reproduced by any interested party, and not just the original investigator, for no cost other than the raw CPU time needed to parse the archived script. In fact, this entire lecture can be reconstructed in such a fashion. To do this, the script: constructs a number of shock-capturing schemes; runs a series of test problems, generates the plots shown; outputs the LATEX to typeset the notes; performs a myriad of behind-the-scenes tasks to glue everything together. Thus Amrita has all the characteristics of an operating system and should not be mistaken for a common-or-garden code