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.

Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

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.

Concordance between dispersal and mitochondrial gene flow: isolation by distance in a tropical teleost, Lates calcarifer (Australian barramundi)

Patterns of population subdivision and the relationship between gene flow and geographical distance in the tropical estuarine fish Lares calcarifer (Centropomidae) were investigated using mtDNA control region sequences. Sixty-three putative haplotypes were resolved from a total of 270 individuals from nine localities within three geographical regions spanning the north Australian coastline. Despite a continuous estuarine distribution throughout the sampled range, no haplotypes were shared among regions. However, within regions, common haplotypes were often shared among localities. Both sequence-based (average Phi(ST)=0.328) and haplotype-based (average Phi(ST)=0.182) population subdivision analyses indicated strong geographical structuring. Depending on the method of calculation, geographical distance explained either 79 per cent (sequence-based) or 23 per cent (haplotype-based) of the variation in mitochondrial gene flow. Such relationships suggest that genetic differentiation of L. calcarifer has been generated via isolation-by-distance, possibly in a stepping-stone fashion. This pattern of genetic structure is concordant with expectations based on the life history of L. calcarifer and direct studies of its dispersal patterns. Mitochondrial DNA variation, although generally in agreement with patterns of allozyme variation, detected population subdivision at smaller spatial scales. Our analysis of mtDNA variation in L. calcarifer confirms that population genetic models can detect population structure of not only evolutionary significance but also of demographic significance. Further, it demonstrates the power of inferring such structure from hypervariable markers, which correspond to small effective population sizes.

Two dimensional computational fluid dynamic models for waste stabilisation ponds

Traditional waste stabilisation pond (WSP) models encounter problems predicting pond performance because they cannot account for the influence of pond features, such as inlet structure or pond geometry, on fluid hydrodynamics. In this study, two dimensional (2-D) computational fluid dynamics (CFD) models were compared to experimental residence time distributions (RTD) from literature. In one of the-three geometries simulated, the 2-D CFD model successfully predicted the experimental RTD. However, flow patterns in the other two geometries were not well described due to the difficulty of representing the three dimensional (3-D) experimental inlet in the 2-D CFD model, and the sensitivity of the model results to the assumptions used to characterise the inlet. Neither a velocity similarity nor geometric similarity approach to inlet representation in 2-D gave results correlating with experimental data. However. it was shown that 2-D CFD models were not affected by changes in values of model parameters which are difficult to predict, particularly the turbulent inlet conditions. This work suggests that 2-D CFD models cannot be used a priori to give an adequate description of the hydrodynamic patterns in WSP. (C) 1998 Elsevier Science Ltd. All rights reserved.

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.

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.

Analysis of small signal stability margins using genetic optimization

Power system small signal stability analysis aims to explore different small signal stability conditions and controls, namely: (1) exploring the power system security domains and boundaries in the space of power system parameters of interest, including load flow feasibility, saddle node and Hopf bifurcation ones; (2) finding the maximum and minimum damping conditions; and (3) determining control actions to provide and increase small signal stability. These problems are presented in this paper as different modifications of a general optimization to a minimum/maximum, depending on the initial guesses of variables and numerical methods used. In the considered problems, all the extreme points are of interest. Additionally, there are difficulties with finding the derivatives of the objective functions with respect to parameters. Numerical computations of derivatives in traditional optimization procedures are time consuming. In this paper, we propose a new black-box genetic optimization technique for comprehensive small signal stability analysis, which can effectively cope with highly nonlinear objective functions with multiple minima and maxima, and derivatives that can not be expressed analytically. The optimization result can then be used to provide such important information such as system optimal control decision making, assessment of the maximum network's transmission capacity, etc. (C) 1998 Elsevier Science S.A. All rights reserved.

A state-contingent production approach to principal-agent problems with an application to point-source pollution control

Two basic representations of principal-agent relationships, the 'state-space' and 'parameterized distribution' formulations, have emerged. Although the state-space formulation appears more natural, analytical studies using this formulation have had limited success. This paper develops a state-space formulation of the moral-hazard problem using a general representation of production under uncertainty. A closed-form solution for the agency-cost problem is derived. Comparative-static results are deduced. Next we solve the principal's problem of selecting the optimal output given the agency-cost function. The analysis is applied to the problem of point-source pollution control. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Modeling the fluid-flow-induced stress and collapse in a dendritic network

An analytical approach to the stress development in the coherent dendritic network during solidification is proposed. Under the assumption that stresses are developed in the network as a result of the friction resisting shrinkage-induced interdendritic fluid flow, the model predicts the stresses in the solid. The calculations reflect the expected effects of postponed dendrite coherency, slower solidification conditions, and variations of eutectic volume fraction and shrinkage. Comparing the calculated stresses to the measured shear strength of equiaxed mushy zones shows that it is possible for the stresses to exceed the strength, thereby resulting in reorientation or collapse of the dendritic network.

Effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in deformable fluid-saturated porous media heated from below

In this paper, a solution method is presented to deal with fully coupled problems between medium deformation, pore-fluid flow and heat transfer in fluid-saturated porous media having supercritical Rayleigh numbers. To validate the present solution method, analytical solutions to a benchmark problem are derived for some special cases. After the solution method is validated, a numerical study is carried out to investigate the effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in fluid-saturated media when they are heated from below. The related numerical results have demonstrated that: (1) medium thermoelasticity has a little influence on the overall pattern of convective pore-fluid flow, but it has a considerable effect on the localization of medium deformation, pore-fluid flow, heat transfer and mineralization in a porous medium, especially when the porous medium is comprised of soft rock masses; (2) convective pore-fluid flow plays a very important role in the localization of medium deformation, heat transfer and mineralization in a porous medium. (C) 1999 Elsevier Science S.A. All rights reserved.

Hypoxia stimulates cerebral blood flow in the estuarine crocodile (Crocodylus porosus)

The effect of N-2 respiration on cerebral blood flow (CBF) velocity on the dorsal surface of cerebellum was examined in the estuarine crocodile, Crocodylus porosus, using epi-illumination microscopy. Twelve minutes of N-2 respiration resulted in a 126% increase in CBF velocity. N-2 respiration had no effect on blood pressure, indicating an underlying cerebral vasodilation. In addition, heart rate increased significantly. Systemic injections of aminophylline and the NO synthase (NOS) inhibitor nitro-L-arginine (L-NA) did not affect the hypoxia induced increase in CBF. We conclude that C. porosus responds to hypoxia with adenosine and nitric oxide (NO) independent cerebral vasodilation, and that this is likely to be a mechanism protecting the brain from energy deficiency during prolonged dives. (C) 1999 Elsevier Science Ireland Ltd. All rights reserved.

Targeting local tissues by transdermal application: Understanding drug physicochemical properties that best exploit protein binding and blood flow effects

The targeting of topically applied drug molecules into tissues below a site of application requires an understanding of the complex interrelationships between the drug, its formulation, the barrier properties of the skin, and the physiological processes occurring below the skin that are responsible for drug clearance from the site, tissue, and/or systemic distribution and eventual elimination. There is still a certain amount of controversy over the ability of topically applied drugs to penetrate into deeper tissues by diffusion or whether this occurs by redistribution in the systemic circulation. The major focus of our work in this area has been in determining how changes in drug structure and physicochemical properties, such as protein binding and lipophilicity, affect drug clearance into the local dermal microcirculation and lymphatics, as well as subsequent distribution into deeper tissues below an application site. The present study outlines our recent thinking on the drug molecule optimal physical attributes, in terms of plasma and tissue partitioning behaviour, that offer the greatest potential for deep tissue targeting. Drug Dev. Res. 46:309-315, 1999. (C) 1999 Wiley-Liss, Inc.