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.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

A consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements

To translate and transfer solution data between two totally different meshes (i.e. mesh 1 and mesh 2), a consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements is presented in this paper. The proposed algorithm has the following significant advantages: (1) The use of a point-searching strategy allows a point in one mesh to be accurately related to an element (containing this point) in another mesh. Thus, to translate/transfer the solution of any particular point from mesh 2 td mesh 1, only one element in mesh 2 needs to be inversely mapped. This certainly minimizes the number of elements, to which the inverse mapping is applied. In this regard, the present algorithm is very effective and efficient. (2) Analytical solutions to the local co ordinates of any point in a four-node quadrilateral element, which are derived in a rigorous mathematical manner in the context of this paper, make it possible to carry out an inverse mapping process very effectively and efficiently. (3) The use of consistent interpolation enables the interpolated solution to be compatible with an original solution and, therefore guarantees the interpolated solution of extremely high accuracy. After the mathematical formulations of the algorithm are presented, the algorithm is tested and validated through a challenging problem. The related results from the test problem have demonstrated the generality, accuracy, effectiveness, efficiency and robustness of the proposed consistent point-searching algorithm. Copyright (C) 1999 John Wiley & Sons, Ltd.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Effect of removing user fees on attendance for curative and preventive primary health care services in rural South Africa

User fees are used to recover costs and discourage unnecessary attendance at primary care clinics in many developing countries. In South Africa, user fees for children aged under 6 yea rs and pregnant women were removed in 1994, and in 1997 all user fees at all primary health care clinics were abolished. The intention of these policy changes was to improve access to health services for previously disadvantaged communities. We investigated the impact of these changes on clinic attendance patterns in Hlabisa health district. Average quarterly new registrations and total attendances for preventive services (antenatal care, immunization, growth monitoring) and curative services (treatment of ailments) at a mobile primary health care unit were studied from 1992 to 1998. Regression analysis was undertaken to assess whether trends were statistically significant. There was a sustained increase in new registrations (P = 0.0001) and total attendances (P = 0.0001)for curative services, and a fall in new registrations (P = 0.01) and total attendances for immunization and growth monitoring (P = 0.0002) over the study period. The upturn in demand for curative services started at the time of the first policy change. The decreases in antenatal registrations (P = 0.07) and attendances (P = 0.09) were not statistically significant The number of new registrations for immunization and growth monitoring increased following the first policy change but declined thereafter. We found no evidence that the second policy change influenced underlying trends. The removal of user fees improved access to curative services but this may have happened at the expense of some preventive services. Governments should remain vigilant about the effects of new health policies in order to ensure that objectives are being met.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Pt-RuO(2) electrodes prepared by thermal decomposition of polymeric precursors as catalysts for direct methanol fuel cell applications

The electrocatalytic activity of Pt and RuO(2) mixed electrodes of different compositions towards methanol oxidation was investigated. The catalysts were prepared by thermal decomposition of polymeric precursors and characterized by energy dispersive X-ray, scanning electronic microscopy, X-ray diffraction and cyclic voltammetry. This preparation method allowed obtaining uniform films with controlled stoichiometry and high surface area. Cyclic voltammetry experiments in the presence of methanol showed that mixed electrodes decreased the potential peak of methanol oxidation by approximately 100 mV (RHE) when compared to the electrode containing only Pt. In addition, voltammetric experiments indicated that the Pt(0.6)Ru(0.4)O(y) electrode led to higher oxidation current densities at lower potentials. Chronoamperometry experiments confirmed the contribution of RuO(2) to the catalytic activity as well as the better performance of the Pt(0.6)Ru(0.4)O(y) electrode composition. Formic acid and CO(2) were identified as being the reaction products formed in the electrolysis performed at 400 and 600 mV. The relative formation of CO(2) was favored in the electrolysis performed at 400 mV (RHE) with the Pt(0.6)Ru(0.4)O(y) electrode. The presence of RuO(2) in Pt-Ru-based electrodes is important for improving the catalytic activity towards methanol electrooxidation. Moreover, the thermal decomposition of polymeric precursors seems to be a promising route for the production of catalysts applicable to DMFC. (C) 2009 International Association for Hydrogen Energy. Published by Elsevier Ltd. All rights reserved.

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.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Splanchnic non-hepatic hemodynamics and metabolism during liver transplantation

Background/Aims: The aim of this study is to compare the splanchnic non-hepatic hemodynamics and the metabolic changes during orthotopic liver transplantation between the conventional with bypass and the piggyback methods. Methodology: A prospective, consecutive series of 59 primary transplants were analyzed. Oxygen consumption, glucose, potassium, and lactate metabolism were quantitatively estimated from blood samples from the radial artery and portal vein, collected up to 120 minutes after graft reperfusion. Mean arterial pressure, portal venous pressure, portal venous blood flow, and splanchnic vascular resistance were also measured or calculated at postreperfusion collection times. Results: There was a greater increase in portal venous blood flow (p=0.05) and lower splanchnic vascular resistance (p=0.04) in the piggyback group. Mean arterial pressure and portal venous pressure were similar for both groups. Oxygen, glucose and potassium consumption were higher in the piggyback group, but none of the metabolic parameters differed significantly between groups. Conclusions: In conclusion, the study detected a higher portal venous blood flow and a lower and splanchnic vascular resistance associated with the piggyback technique. After graft reperfusion, no difference in the splanchnic non-hepatic metabolic parameters was observed between the conventional with bypass and the piggyback methods of orthotopic liver transplantation.

A computerized tomography scan method for calculating the hernia sac and abdominal cavity volume in complex large incisional hernia with loss of domain

Preoperative progressive pneumoperitoneum (PPP) is a safe and effective procedure in the treatment of large incisional hernia (size > 10 cm in width or length) with loss of domain (LIHLD). There is no consensus in the literature on the amount of gas that must be insufflated in a PPP program or even how long it should be maintained. We describe a technique for calculating the hernia sac volume (HSV) and abdominal cavity volume (ACV) based on abdominal computerized tomography (ACT) scanning that eliminates the need for subjective criteria for inclusion in a PPP program and shows the amount of gas that must be insufflated into the abdominal cavity in the PPP program. Our technique is indicated for all patients with large or recurrent incisional hernias evaluated by a senior surgeon with suspected LIHLD. We reviewed our experience from 2001 to 2008 of 23 consecutive hernia surgical procedures of LIHLD undergoing preoperative evaluation with CT scanning and PPP. An ACT was required in all patients with suspected LIHLD in order to determine HSV and ACV. The PPP was performed only if the volume ratio HSV/ACV (VR = HSV/ACV) was a parts per thousand yen25% (VR a parts per thousand yen 25%). We have performed this procedure on 23 patients, with a mean age of 55.6 years (range 31-83). There were 16 women and 7 men with an average age of 55.6 years (range 31-83), and a mean BMI of 38.5 kg/m(2) (range 23-55.2). Almost all patients (21 of 23 patients-91.30%) were overweight; 43.5% (10 patients) were severely obese (obese class III). The mean calculated volumes for ACV and HSV were 9,410 ml (range 6,060-19,230 ml) and 4,500 ml (range 1,850-6,600 ml), respectively. The PPP is performed by permanent catheter placed in a minor surgical procedure. The total amount of CO(2) insufflated ranged from 2,000 to 7,000 ml (mean 4,000 ml). Patients required a mean of 10 PPP sessions (range 4-18) to achieve the desired volume of gas (that is the same volume that was calculated for the hernia sac). Since PPP sessions were performed once a day, 4-18 days were needed for preoperative preparation with PPP. The mean VR was 36% (ranged from 26 to 73%). We conclude that ACT provides objective data for volume calculation of both hernia sac and abdominal cavity and also for estimation of the volume of gas that should be insufflated into the abdominal cavity in PPP.

Vitamin D intoxication: a cause of hypocalcaemia and acute renal failure in a HIV patient

Hypercalcaemia in patients with HIV infection is usually associated with specific conditions such as lymphoma and granulomatous diseases. We described a case of severe hypercalcaemia consequent to vitamin D intoxication and secondary renal failure in a HIV patient under tenofovir using. Serum creatinine and calcium returned to near normal levels after vitamin D discontinuation, saline and furosemide administration. Some aspects of the drug-induced nephropathy are discussed.

One-Stage Revision of Infected Total Hip Arthroplasty with Bone Graft

There are many different opinions in the literature regarding the best procedure for revision of infected hip arthroplasty and hence in achieving long-term stabilization of a new implant. Thirty-two patients with 32 loose and infected total hip arthroplasties underwent revision with a bone graft in a 1-stage procedure. The bone graft was used in the acetabulum and femure in 25 patients, in the acetabulum alone in 4 patients and in the femur alone in 3 patients. A metal mesh was necessary in 15 patients to contain the morselized bone graft. At the time of surgical revision, 9 patients had a draining sinus, 6 had a closed sinus, and 17 had never had sinuses in the surgical wound. Antibiotic therapy was administered intravenously and orally for 6 months. Mean follow-up was 103 months (range, 63-183 months), and infection recurred in 2 (6.2%) cases. Further studies are necessary, and continuation of this method is justified.

Three-dimensional mapping of hippocampal anatomy in unmedicated and lithium-treated patients with bipolar disorder

Declarative memory impairments are common in patients with bipolar illness, suggesting underlying hippocampal pathology. However, hippocampal volume deficits are rarely observed in bipolar disorder. Here we used surface-based anatomic mapping to examine hippocampal anatomy in bipolar patients treated with lithium relative to matched control subjects and unmedicated patients with bipolar disorder. High-resolution brain magnetic resonance images were acquired from 33 patients with bipolar disorder ( 21 treated with lithium and 12 unmedicated), and 62 demographically matched healthy control subjects. Three-dimensional parametric mesh models were created from manual tracings of the hippocampal formation. Total hippocampal volume was significantly larger in lithium-treated bipolar patients compared with healthy controls (by 10.3%; p=0.001) and unmedicated bipolar patients ( by 13.9%; p=0.003). Statistical mapping results, confirmed by permutation testing, revealed localized deficits in the right hippocampus, in regions corresponding primarily to cornu ammonis vertical bar subfields, in unmedicated bipolar patients, as compared to both normal controls (p=0.01), and in lithium-treated bipolar patients (p=0.03). These findings demonstrate the sensitivity of these anatomic mapping methods for detecting subtle alterations in hippocampal structure in bipolar disorder. The observed reduction in subregions of the hippocampus in unmedicated bipolar patients suggests a possible neural correlate for memory deficits frequently reported in this illness. Moreover, increased hippocampal volume in lithium-treated bipolar patients may reflect postulated neurotrophic effects of this agent, a possibility warranting further study in longitudinal investigations.