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.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Effect of load distance on self-selected manual lifting technique

The aim of this research is to determine the effects of constraining the horizontal distance of the feet from the load on the posture adopted at the start of the lift. Kinematic data were collected while each of 24 subjects lifted 3, 6, and 9 kg loads from a starting height 18 cm above the ground. The position of the feet was controlled relative to the load such that the horizontal distance from the hand to the ankle at the start of extension was either 20, 40, or 60 cm. Subjects performed 20 trials in each of six combinations of load and ankle-load distance chosen to provide three sets of equivilent load moment pairs. The initial horizontal distance from the load to the ankle had a large influence on the posture adopted to lift the load. Ankle and knee flexion, in particular, were reduced when the ankle-load distance was smaller, and particularly so when the distance was reduced to 20 cm. Hip flexion was reduced to a smaller extent, while lumbar vertebral flexion remained relatively unchanged. The inclination of the trunk at the start of the lift was unchanged when the ankle-load distance was 60 or 40 cm, but was 10 degrees greater when the load was 20 cm from the ankles, indicating that subjects adopted a posture closer to a stoop when the ankle-load distance was small. Comparison of conditions of equal load moment (but different load mass and ankle-load distance) revealed differences which mirrored the effects of ankle-load distance alone, suggesting that the effects of ankle-load distance on the posture adopted at the start of extension were largely independent of the load moment. While the forces and torques required to lift a load must be to some extent dependent on the load moment, rather than load or ankle-load distance per se, the posture adopted to lift the load is not.

Modeling ex vivo hematopoiesis using chemical engineering metaphors

Ex vivo hematopoiesis is increasingly used for clinical applications. Models of ex vivo hematopoiesis are required to better understand the complex dynamics and to optimize hematopoietic culture processes. A general mathematical modeling framework is developed which uses traditional chemical engineering metaphors to describe the complex hematopoietic dynamics. Tanks and tubular reactors are used to describe the (pseudo-) stochastic and deterministic elements of hematopoiesis, respectively. Cells at any point in the differentiation process can belong to either an immobilized, inert phase (quiescent cells) or a mobile, active phase (cycling cells). The model describes five processes: (1) flow (differentiation), (2) autocatalytic formation (growth),(3) degradation (death), (4) phase transition from immobilized to mobile phase (quiescent to cycling transition), and (5) phase transition from mobile to immobilized phase (cycling to quiescent transition). The modeling framework is illustrated with an example concerning the effect of TGF-beta 1 on erythropoiesis. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

Flush air data system calibration using numerical simulation

The use of computational fluid dynamics simulations for calibrating a flush air data system is described, In particular, the flush air data system of the HYFLEX hypersonic vehicle is used as a case study. The HYFLEX air data system consists of nine pressure ports located flush with the vehicle nose surface, connected to onboard pressure transducers, After appropriate processing, surface pressure measurements can he converted into useful air data parameters. The processing algorithm requires an accurate pressure model, which relates air data parameters to the measured pressures. In the past, such pressure models have been calibrated using combinations of flight data, ground-based experimental results, and numerical simulation. We perform a calibration of the HYFLEX flush air data system using computational fluid dynamics simulations exclusively, The simulations are used to build an empirical pressure model that accurately describes the HYFLEX nose pressure distribution ol cr a range of flight conditions. We believe that computational fluid dynamics provides a quick and inexpensive way to calibrate the air data system and is applicable to a broad range of flight conditions, When tested with HYFLEX flight data, the calibrated system is found to work well. It predicts vehicle angle of attack and angle of sideslip to accuracy levels that generally satisfy flight control requirements. Dynamic pressure is predicted to within the resolution of the onboard inertial measurement unit. We find that wind-tunnel experiments and flight data are not necessary to accurately calibrate the HYFLEX flush air data system for hypersonic flight.

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].

Frequency stabilised grating feedback laser diode for atom cooling applications

The free running linewidth of an external cavity grating feedback diode laser is on the order of a few megahertz and is limited by the mechanical and acoustic vibrations of the external cavity. Such frequency fluctuations can be removed by electronic feedback. We present a hybrid stabilisation technique that uses both a Fabry-Perot confocal cavity and an atomic resonance to achieve excellent short and long term frequency stability. The system has been shown to reduce the laser linewidth of an external cavity diode laser by an order of magnitude to 140 kHz, while limiting frequency excursions to 60 kHz relative to an absolute reference over periods of several hours. The scheme also presents a simple way to frequency offset two lasers many gigahertz apart which should find a use in atom cooling experiments, where hyperfine ground-state frequency separations are often required.

Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining was carried out to investigate microfracture, residual damage, and surface roughness associated with material removal and surface generation. Brittle versus plastic deformation was studied using Vickers indention and nano-indentation. To characterize the abrasive machining response, the 6H-SiC (0001) substrates were ground using diamond wheels with grit sizes of 25, 15 and 7 mum, and then polished with diamond suspensions of 3 and 0.05 mum. It is found that in indentation, there was a scale effect for brittle versus plastic deformation in 6H-SiC substrates. Also, in grinding, the scales of fracture and surface roughness of the substrates decreased with a decrease in diamond grit size. However, in polishing, a reduction in grit size of diamond suspensions gave no significant improvement in surface roughness. Furthermore, the results showed that fracture-free 6H-SiC (0001) surfaces were generated in polishing with the existence of the residual crystal defects, which were associated with the origin of defects in single crystal growth. (C) 2003 Elsevier Ltd. All rights reserved.

High speed versus conventional grinding in high removal rate machining of alumina and alumina-titania

High removal rate (up to 16.6 mm(3)/s per mm) grinding of alumina and alumina-titania was investigated with respect to material removal and basic grinding parameters using a resin-bond 160 mu m grit diamond wheel at the speeds of 40 and 160 m/s, respectively. The results show that the material removal for the single-phase polycrystalline alumina and the two-phase alumina-titania composite revealed identical mechanisms of microfracture and grain dislodgement under the grinding conditioned selected. There were no distinct differences in surface roughness and morphology for both materials ground at either conventional or high speed. An increase in material removal rate did not necessarily worsen the surface toughness for the two materials at both speeds. Also the grinding forces for the two ceramics demonstrated similar characteristics at any grinding speeds and specific removal rates. Both normal and tangential grinding forces and their force ratios at the high speed were lower than those at the conventional speed, regardless of removal rates. An increase in specific removal rate caused more rapid increases in normal and tangential forces obtained at the conventional grinding speed than those at the high speed. Furthermore, it is found that the high speed grinding at all the removal rates exerted a great amount of coolant-induced normal forces in grinding zone, which were 4-6 times higher than the pure normal grinding forces. (c) 2004 Elsevier Ltd. All rights reserved.

Ground-State Entanglement of the BCS Model

The concept of local concurrence is used to quantify the entanglement between a single qubit and the remainder of a multiqubit system. For the ground state of the BCS model in the thermodynamic limit the set of local concurrences completely describes the entanglement. As a measure for the entanglement of the full system we investigate the average local concurrence (ALC). We find that the ALC satisfies a simple relation with the order parameter. We then show that for finite systems with a fixed particle number, a relation between the ALC and the condensation energy exposes a threshold coupling. Below the threshold, entanglement measures besides the ALC are significant.

Machining of micro aspherical mould inserts

This paper reports a parametric investigation and development of grinding technologies for micro aspherical mould inserts using parallel grinding method. The parametric investigation revealed that at nanometric scale the undeformed chip thickness has little influence on the surface finish of ground inserts. The grinding trace spacing has a slightly larger influence on the surface finish. A new technique was developed to true and dress the resin bonded micro wheels with mesh size of #3000, which produced a satisfactory wheel form accuracy and relatively high grain packing density. A form error compensation technique was also developed, with which mould inserts of submicron form accuracy were consistently produced. Using the developed technologies, micro aspherical inserts of diameters ranging from 200 mu m to 1000 mu m with surface finish of around 10 nm and form error of similar to 0.2-0.4 mu m were successfully fabricated. (c) 2005 Elsevier Inc. All rights reserved.