Theory manual to OctVCE – a cartesian cell CFD code with special application to blast wave problems, Department of Mechanical Engineering Report 2007/12

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries, in particular, from explosions. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel for its simplicity and sufficiency for practical engineering design problems. The code uses a finite-volume formulation of the unsteady Euler equations with a second order explicit Runge-Kutta Godonov (MUSCL) scheme. Gradients are calculated using a least-squares method with a minmod limiter. Flux solvers used are AUSM, AUSMDV and EFM. No fluid-structure coupling or chemical reactions are allowed, but gas models can be perfect gas and JWL or JWLB for the explosive products. This report also describes the code’s ‘octree’ mesh adaptive capability and point-inclusion query procedures for the VCE geometry engine. Finally, some space will also be devoted to describing code parallelization using the shared-memory OpenMP paradigm. The user manual to the code is to be found in the companion report 2007/13.

Effects of memory load and distraction on performance and event-related slow potentials in a visuospatial working memory task

Brain electrical activity related to working memory was recorded at 15 scalp electrodes during a visuospatial delayed response task. Participants (N = 18) touched the remembered position of a target on a computer screen after either a 1 or 8 sec delay. These memory trials were compared to sensory trials in which the target remained present throughout the delay and response periods. Distracter stimuli identical to the target were briefly presented during the delay on 30% of trials. Responses were less accurate in memory than sensory trials, especially after the long delay. During the delay slow potentials developed that were significantly more negative in memory than sensory trials. The difference between memory and sensory trials was greater at anterior than posterior electrodes. On trials with distracters, the slow potentials generated by memory trials showed further enhancement of negativity whereas there were minimal effects on accuracy of performance. The results provide evidence that engagement of visuospatial working memory generates slow wave negativity with a timing and distribution consistent with frontal activation. Enhanced brain activity associated with working memory is required to maintain performance in the presence of distraction. © 1997 by the Massachusetts Institute of Technology

Characteristics of undular hydraulic jumps: Experiments and analysis

The writers measured velocity, pressure and energy distributions, wavelengths, and wave amplitudes along undular jumps in a smooth rectangular channel 0.25 m wide. In each case the upstream flow was a fully developed shear flow. Analysis of the data shows that the jump has strong three-dimensional features and that the aspect ratio of the channel is an important parameter. Energy dissipation on the centerline is far from negligible and is largely constrained to the reach between the start of the lateral shock waves and the first wave crest of the jump, in which the boundary layer develops under a strong adverse pressure gradient. A Boussinesq-type solution of the free-surface profile, velocity, and energy and pressure distributions is developed and compared with the data. Limitations of the two-dimensional analysis are discussed.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

A computational study of shock speeds in high-performance shock tubes

This paper describes U2DE, a finite-volume code that numerically solves the Euler equations. The code was used to perform multi-dimensional simulations of the gradual opening of a primary diaphragm in a shock tube. From the simulations, the speed of the developing shock wave was recorded and compared with other estimates. The ability of U2DE to compute shock speed was confirmed by comparing numerical results with the analytic solution for an ideal shock tube. For high initial pressure ratios across the diaphragm, previous experiments have shown that the measured shock speed can exceed the shock speed predicted by one-dimensional models. The shock speeds computed with the present multi-dimensional simulation were higher than those estimated by previous one-dimensional models and, thus, were closer to the experimental measurements. This indicates that multi-dimensional flow effects were partly responsible for the relatively high shock speeds measured in the experiments.

Velocimetry and thermometry of supersonic flow around a cylindrical body

The supersonic flow around a cylindrical body has been studied using two optical techniques. For both sets of measurements, the cylinder was mounted from the side of the tunnel, allowing investigation of the bow shock region as well as in the wake. A new technique, laser-enhanced ionization flow tagging, was used for streamwise velocity determinations behind the body. From these measurements, it was found that the downstream velocity outside the wake was (1.90 +/- 0.06) km/s, whereas inside the wake the velocity was about 0-500 m/s in the upstream direction. Planar laser induced fluorescence of nitric oxide was employed for temperature determinations. It was established that the freestream temperature was (2120 +/- 100) K, decreasing to around (1550 +/- 400) K in the wake.

Peptidergic control of gastrointestinal blood flow in the estuarine crocodile, Crocodylus porosus

Peptidergic mechanisms influencing the resistance of the gastrointestinal vascular bed of the estuarine crocodile, Crocodylus porosus, were investigated. The gut was perfused in situ via the mesenteric and the celiac arteries, and the effects of different neuropeptides were tested using bolus injections. Effects on vascular resistance were recorded as changes in inflow pressures. Peptides found in sensory neurons [substance P, neurokinin A, and calcitonin gene-related peptide (CGRP)] all caused significant relaxation of the celiac vascular bed, as did vasoactive intestinal polypeptide (VIP), another well-known vasodilator. Except for VIP, the peptides also induced transitory gut contractions. Somatostatin and neuropeptide Y (NPY), which coexist in adrenergic neurons of the C. porosus, induced vasoconstriction in the celiac vascular bed without affecting the gut motility. Galanin caused vasoconstriction and occasionally activated the gut wall. To elucidate direct effects on individual vessels, the different peptides were tested on isolated ring preparations of the mesenteric and celiac arteries. Only CGRP and VIP relaxed the epinephrine-precontracted celiac artery, whereas the effects on the mesenteric artery were variable. Somatostatin and NPY did not affect the resting tonus of these vessels, but somatostatin potentiated the epinephrine-induced contraction of the celiac artery. Immunohistochemistry revealed the existence and localization of the above-mentioned peptides in nerve fibers innervating vessels of different sizes in the gut region. These data support the hypothesis of an important role for neuropeptides in the control of the vascular bed of the gastrointestinal tract in C. porosus.

Hepatic structure-pharmacokinetic relationships: The hepatic disposition and metabolite kinetics of a homologous series of O-acyl derivatives of salicylic acid

1 The hepatic disposition and metabolite kinetics of a homologous series of O-acyl (acetyl, propionyl, butanoyl, pentanoyl, hexanoyl and octanoyl) esters of salicylic acid (C2SA, C3SA, C4SA, C5SA, C6SA and C8SA, respectively) was determined using a single-pass, in-sills rat liver preparation. 2 The hepatic venous outflow profiles for the parent esters and the generated metabolite, salicylic acid (SA) were analysed by HPLC. Non-parametric moments analysis was used to determine the area under the curve (AUC'), mean transit time (MTT) and normalized variance (CV2) for the parent esters and generated SA. 3 Pregenerated SA ([C-14]-salicylic acid) was injected into each liver with the parent ester to determine its distribution characteristics. 4 The overall recovery of ester plus metabolite was 89% of the ester dose injected and independent of the ester carbon number, suggesting that ester extraction was due to hepatic metabolism to salicylic acid. 5 The metabolite AUC' value increased directly with the lipophilicity of the parent ester (from 0.12 for C2SA to 0.95 for C8SA). By contrast, the parent AUC' decreased with the lipophilicity (from 0.85 for C2SA to zero for C8SA). The metabolite MTT value also showed a trend to increase with the lipophilicity of the parent ester (from 15.72 s for C3SA to 61.97 s for C8SA). However, the parent MTT value shows no significant change across the series. 6 The two-compartment dispersion model was used to derive the kinetic parameters for parent ester, pregenerated SA and generated SA. Consequently, these parameters were used to estimate the values of AUG', MITT and CV2 for the parent ester and metabolite. The moments values obtained using the two-compartment dispersion model show similar trends to the corresponding moments values obtained from the outflow profiles using a non-parametric approach. 7 The more lipophilic aspirin analogues are more confined to the portal circulation after oral administration than aspirin due to their more extensive hepatic elimination avoiding systemic prostacyclin inhibition. Given that aspirin's selectivity as an anti-thrombotic agent has been postulated to be due to selective anti-platelet effects in the portal circulation, the more lipophilic and highly extracted analogues are potentially more selective anti-thrombotic agents than aspirin.

Proliferation in monocyte-derived dendritic cell cultures is caused by progenitor cells capable of myeloid differentiation

Dendritic cells (DC) can be generated by culture of adherent peripheral blood (PB) cells in the presence of granulocyte-macrophage colony-stimulating factor (GM-CSF) and interleukin-4 (IL-4). There is controversy as to whether these DC arise from proliferating precursors or simply from differentiation of monocytes. DC were generated from myeloid-enriched PB non-T cells or sorted monocytes. DC generated from either population functioned as potent antigen-presenting cells. Uptake of [H-3]-thymidine was observed in DC cultured from myeloid-enriched non-T cells. Addition of lipopolysaccharide or tumor necrosis factor-alpha led to maturation of the DC, but did not inhibit proliferation. Ki67(+) cells were observed in cytospins of these DC, and by double staining were CD3(-)CD19(-)CD11c(-)CD40(-) and myeloperoxidase(+), suggesting that they were myeloid progenitor cells. Analysis of the starting population by flow cytometry demonstrated small numbers of CD34(+)CD33(-)CD14(-) progenitor cells, and numerous granulocyte-macrophage colony-forming units were generated in standard assays. Thus, production of DC in vitro from adherent PB cells also enriches for progenitor cells that are capable of proliferation after exposure to GM-CSF. Of clinical importance, the yield of DC derived in the presence of GM-CSF and IL-4 cannot be expanded beyond the number of starting monocytes. (C) 1998 by The American Society of Hematology.

Proliferation in monocyte-derived dendritic cell cultures is due to cells capable of myeloid differentiation

Dendritic cells (DC) can be generated by culture of adherent peripheral blood (PB) cells in the presence of granulocyte-macrophage colony-stimulating factor (GM-CSF) and interleukin-4 (IL-4). There is controversy as to whether these DC arise from proliferating precursors or simply from differentiation of monocytes. DC were generated from myeloid-enriched PB non-T cells or sorted monocytes. DC generated from either population functioned as potent antigen-presenting cells. Uptake of [H-3]-thymidine was observed in DC cultured from myeloid-enriched non-T cells. Addition of lipopolysaccharide or tumor necrosis factor-alpha led to maturation of the DC, but did not inhibit proliferation. Ki67(+) cells were observed in cytospins of these DC, and by double staining were CD3(-)CD19(-)CD11c(-)CD40(-) and myeloperoxidase(+), suggesting that they were myeloid progenitor cells. Analysis of the starting population by flow cytometry demonstrated small numbers of CD34(+)CD33(-)CD14(-) progenitor cells, and numerous granulocyte-macrophage colony-forming units were generated in standard assays. Thus, production of DC in vitro from adherent PB cells also enriches for progenitor cells that are capable of proliferation after exposure to GM-CSF. Of clinical importance, the yield of DC derived in the presence of GM-CSF and IL-4 cannot be expanded beyond the number of starting monocytes. (C) 1998 by The American Society of Hematology.

Quantum oscillations in quasi-one-dimensional metals with spin-density-wave ground states

We consider the magnetoresistance oscillation phenomena in the Bechgaard salts (TMTSF)(2)X, where X = ClO4, PF6, and AsF6 in pulsed magnetic fields to 51 T. Of particular importance is the observation of a new magnetoresistance oscillation for X = ClO4 in its quenched state. In the absence of any Fermi-surface reconstruction due to anion order at low temperatures, all three materials exhibit nonmonotonic temperature dependence of the oscillation amplitude in the spin-density-wave (SDW) state. We discuss a model where, below a characteristic temperature T* within the SDW state, a magnetic breakdown gap opens. [S0163-1829(99)00904-2].