Thermally developing Brinkman–Brinkman forced convection in rectangular ducts with isothermal walls

The Extended Weighted Residuals Method (EWRM) is applied to investigate the effects of viscous dissipation on the thermal development of forced convection in a porous-saturated duct of rectangular cross-section with isothermal boundary condition. The Brinkman flow model is employed for determination of the velocity field. The temperature in the flow field was computed by utilizing the Green’s function solution based on the EWRM. Following the computation of the temperature field, expressions are presented for the local Nusselt number and the bulk temperature as a function of the dimensionless longitudinal coordinate. In addition to the aspect ratio, the other parameters included in this computation are the Darcy number, viscosity ratio, and the Brinkman number.

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.

High-resolution bulk density images, using calibrated X-ray radiography of impregnated soil slices

Bulk density of undisturbed soil samples can be measured using computed tomography (CT) techniques with a spatial resolution of about 1 mm. However, this technique may not be readily accessible. On the other hand, x-ray radiographs have only been considered as qualitative images to describe morphological features. A calibration procedure was set up to generate two-dimensional, high-resolution bulk density images from x-ray radiographs made with a conventional x-ray diffraction apparatus. Test bricks were made to assess the accuracy of the method. Slices of impregnated soil samples were made using hardsetting seedbeds that had been gamma scanned at 5-mm depth increments in a previous study. The calibration procedure involved three stages: (i) calibration of the image grey levels in terms of glass thickness using a staircase made from glass cover slips, (ii) measurement of ratio between the soil and resin mass attenuation coefficients and the glass mass attenuation coefficient, using compacted bricks of known thickness and bulk density, and (iii) image correction accounting for the heterogeneity of the irradiation field. The procedure was simple, rapid, and the equipment was easily accessible. The accuracy of the bulk density determination was good (mean relative error 0.015), The bulk density images showed a good spatial resolution, so that many structural details could be observed. The depth functions were consistent with both the global shrinkage and the gamma probe data previously obtained. The suggested method would be easily applied to the new fuzzy set approach of soil structure, which requires generation of bulk density images. Also, it would be an invaluable tool for studies requiring high-resolution bulk density measurement, such as studies on soil surface crusts.

Calculating current densities and fields due to shielded bi-planar radio-frequency coils

A method for the accurate computation of the current densities produced in a wide-runged bi-planar radio-frequency coil is presented. The device has applications in magnetic resonance imaging. There is a set of opposing primary rungs, symmetrically placed on parallel planes and a similar arrangement of rungs on two parallel planes surrounding the primary serves as a shield. Current densities induced in these primary and shielding rungs are calculated to a high degree of accuracy using an integral-equation approach, combined with the inverse finite Hilbert transform. Once these densities are known, accurate electrical and magnetic fields are then computed without difficulty. Some test results are shown. The method is so rapid that it can be incorporated into optimization software. Some preliminary fields produced from optimized coils are presented.

Expokit: A software package for computing matrix exponentials

Expokit provides a set of routines aimed at computing matrix exponentials. More precisely, it computes either a small matrix exponential in full, the action of a large sparse matrix exponential on an operand vector, or the solution of a system of linear ODEs with constant inhomogeneity. The backbone of the sparse routines consists of matrix-free Krylov subspace projection methods (Arnoldi and Lanczos processes), and that is why the toolkit is capable of coping with sparse matrices of large dimension. The software handles real and complex matrices and provides specific routines for symmetric and Hermitian matrices. The computation of matrix exponentials is a numerical issue of critical importance in the area of Markov chains and furthermore, the computed solution is subject to probabilistic constraints. In addition to addressing general matrix exponentials, a distinct attention is assigned to the computation of transient states of Markov chains.

The Heisenberg antiferromagnet on an anisotropic triangular lattice: linear spin-wave theory

We consider the effect of quantum spin fluctuations on the ground-state properties of the Heisenberg antiferromagnet on an anisotropic triangular lattice using linear spin-wave (LSW) theory. This model should describe the magnetic properties of the insulating phase of the kappa-(BEDT-TTF)(2)X family of superconducting molecular crystals. The ground-state energy, the staggered magnetization, magnon excitation spectra, and spin-wave velocities are computed as functions of the ratio of the antiferromagnetic exchange between the second and first neighbours, J(2)/J(1). We find that near J(2)/J(1) = 0.5, i.e., in the region where the classical spin configuration changes from a Neel-ordered phase to a spiral phase, the staggered magnetization vanishes, suggesting the possibility of a quantum disordered state. in this region, the quantum correction to the magnetization is large but finite. This is in contrast to the case for the frustrated Heisenberg model on a square lattice, for which the quantum correction diverges logarithmically at the transition from the Neel to the collinear phase. For large J(2)/J(1), the model becomes a set of chains with frustrated interchain coupling. For J(2) > 4J(1), the quantum correction to the magnetization, within LSW theory, becomes comparable to the classical magnetization, suggesting the possibility of a quantum disordered state. We show that, in this regime, the quantum fluctuations are much larger than for a set of weakly coupled chains with non-frustrated interchain coupling.

Biomodel-guided stereotaxy

OBJECTIVES: To simplify the practice of stereotactic surgery by using an original method, apparatus, and solid anatomic replica for trajectory planning and to validate the method and apparatus in a laboratory and clinical trial. METHODS: The patient is marked with fiducials and scanned by using computed tomography or magnetic resonance imaging. The three-dimensional data are converted to a format acceptable to stereolithography. Stereolithography uses a laser to polymerize photosensitive resin into a solid plastic model (biomodel). Stereolithography can replicate brood vessels, soft tissue, tumor, and bone accurately (

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.

X-ray absorption and phase contrast imaging to study the interplay between plant roots and soil structure

Plant performance is, at least partly, linked to the location of roots with respect to soil structure features and the micro-environment surrounding roots. Measurements of root distributions from intact samples, using optical microscopy and field tracings have been partially successful but are imprecise and labour-intensive. Theoretically, X-ray computed micro-tomography represents an ideal solution for non-invasive imaging of plant roots and soil structure. However, before it becomes fast enough and affordable or easily accessible, there is still a need for a diagnostic tool to investigate root/soil interplay. Here, a method for detection of undisturbed plant roots and their immediate physical environment is presented. X-ray absorption and phase contrast imaging are combined to produce projection images of soil sections from which root distributions and soil structure can be analyzed. The clarity of roots on the X-ray film is sufficient to allow manual tracing on an acetate sheet fixed over the film. In its current version, the method suffers limitations mainly related to (i) the degree of subjectivity associated with manual tracing and (ii) the difficulty of separating live and dead roots. The method represents a simple and relatively inexpensive way to detect and quantify roots from intact samples and has scope for further improvements. In this paper, the main steps of the method, sampling, image acquisition and image processing are documented. The potential use of the method in an agronomic perspective is illustrated using surface and sub-surface soil samples from a controlled wheat trial. Quantitative characterization of root attributes, e.g. radius, length density, branching intensity and the complex interplay between roots and soil structure, is presented and discussed.

Magnetic resonance imaging (MRI) and diseases of the liver and biliary tract. Part 1. Basic principles, MRI in the assessment of diffuse and focal hepatic disease

Magnetic resonance imaging (MRI) relies on the physical properties of unpaired protons in tissues to generate images. Unpaired protons behave like tiny bar magnets and will align themselves in a magnetic field. Radiofrequency pulses will excite these aligned protons to higher energy states. As they return to their original state, they will release this energy as radio waves. The frequency of the radio waves depends on the local magnetic field and by varying this over a subject, it is possible to build the images we are familiar with. In general, MRI has not been sufficiently sensitive or specific in the assessment of diffuse liver disease for clinical use. However, because of the specific characteristics of fat and iron, it may be useful in the assessment of hepatic steatosis and iron overload. Magnetic resonance imaging is useful in the assessment of focal liver disease, particularly in conjunction with contrast agents. Haemangiomas have a characteristic bright appearance on T-2 weighted images because of the slow flowing blood in dilated sinusoids. Focal nodular hyperplasia (FNH) has a homogenous appearance, and enhances early in the arterial phase after gadolinium injection, while the central scar typically enhances late. Hepatic adenomas have a more heterogenous appearance and also enhance in the arterial phase, but less briskly than FNH. Hepatocellular carcinoma is similar to an adenoma, but typically occurs in a cirrhotic liver and has earlier washout of contrast. The appearance of metastases depends on the underlying primary malignancy. Overall, MRI appears more sensitive and specific than computed tomography with contrast for the detection and evaluation of malignant lesions. (C) 2000 Blackwell Science Asia Pty Ltd.

A computational analysis of substrate binding strength by phosphorylase kinase and protein kinase A

Protein kinases exhibit various degrees of substrate specificity. The large number of different protein kinases in the eukaryotic proteomes makes it impractical to determine the specificity of each enzyme experimentally. To test if it were possible to discriminate potential substrates from non-substrates by simple computational techniques, we analysed the binding enthalpies of modelled enzyme-substrate complexes and attempted to correlate it with experimental enzyme kinetics measurements. The crystal structures of phosphorylase kinase and cAMP-dependent protein kinase were used to generate models of the enzyme with a series of known peptide substrates and non-substrates, and the approximate enthalpy of binding assessed following energy minimization. We show that the computed enthalpies do not correlate closely with kinetic measurements, but the method can distinguish good substrates from weak substrates and non-substrates. Copyright (C) 2002 John Wiley Sons, Ltd.

Relationship of interleukin-6 and tumor necrosis factor-alpha with muscle mass and muscle strength in elderly men and women: The health ABC study

Background. A decline in muscle mass and muscle strength characterizes normal aging. As clinical and animal studies show it relationship between higher cytokine levels and low muscle mass, the aim of this study was to investigate whether markers, of inflammation are associated with muscle mass and strength in well-functioning elderly persons. Methods. We Used baseline data (1997-1998) of the Health, Aging, and Body Composition (Health ABC) Study on 3075 black and white men and women aged 70-79 years. Midthigh muscle cross-sectional area (computed tomography), appendicular muscle mass (dual-energy x-ray ab absorptiometry). isokinetic knee extensor strength (KinCom). and isometric inip strength were measured. plasma levels of interleukin-6 (IL-6) and tumor necrosis factor-alpha (TNF-alpha) were assessed by enzyme-linked immunosorbent assay (ELISA). Results. Higher cytokine levels were generally associated with lower muscle mass and lower muscle strength. The most consistent relationship across the gender and race groups was observed for IL-6 and grip strength: per SD increase in IL-6, grip strength was 1.1 to 2.4 kg lower (p < .05) after adjustment for age, clinic Site. health status, medications, physical activity. smoking. height. and body fat. Ail overall measure of elevated cytokine level was created by combining the levels of IL-6 and TNF-alpha. With the exception of white men, elderly persons having high levels of IL-6 (> 1.80 pg/ml) as well as high levels of TNF-alpha (> 3.20 pg/ml) had a smaller muscle area, less appendicular mass. a lower knee extensor strength. and a lower grip strength compared to those with low levels of both cytokines. Conclusions. Higher plasma concentrations of IL-6 and TNF-alpha are associated with lower muscle mass and lower muscle strength in well-functioning older men and women. Higher cytokine levels. as often observed in healthy older persons. may contribute to the loss Of muscle mass and strength that accompanies aging.

Ultrasound evaluation of polymer gel dosimeters

A new method for the evaluation of radiotherapy 3D polymer gel dosimeters has been developed using ultrasound to assess the significant structural changes that occur following irradiation of the dosimeters. The ultrasonic parameters of acoustic speed of propagation, attenuation and transmitted signal intensity were measured as a function of absorbed radiation dose. The dose sensitivities for each parameter were determined as 1.8 x 10(-4) s m(-1) Gy(-1), 3.9 dB m(-1) Gy(-1) and 3.2 V-1 Gy(-1) respectively. All parameters displayed a strong variation with absorbed dose that continued beyond absorbed doses of 15 Gy. The ultrasonic measurements demonstrated a significantly larger dynamic range in dose response curves than that achieved with previously published magnetic resonance imaging (MRI) dose response data. It is concluded that ultrasound shows great potential as a technique for the evaluation of polymer gel dosimeters.

Magnetization transfer imaging for polymer gel dosimetry

Off-resonance RF pre-saturation was used to obtain contrast in MRI images of polymer gel dosimeters irradiated to doses up to 50 Gy. Two different polymer gel dosimeters composed of 2-hydroxyethyl-acryl ate or methacrylic acid monomers mixed with N, N'-methylene-bisacrylamide (BIS), dispersed in an aqueous gelatin matrix were evaluated. Radiation-induced polymerization of the co-monomers generates a fast-relaxing insoluble polymer. Saturation of the polymer using off-resonance Gaussian RF pulses prior to a spin-echo read-out with a short echo time leads to contrast that is dependent on the absorbed dose. This contrast is attributed to magnetization transfer (MT) between free water and the polymer, and direct saturation of water was found to be negligible under the prevailing experimental conditions. The usefulness of MT imaging was assessed by computing the dose resolution obtained with this technique. We found a low value of dose resolution over a wide range of doses could be obtained with a single experiment. This is an advantage over multiple spin echo (MSE) experiments using a single echo spacing where an optimal dose resolution is achieved over only very limited ranges of doses. The results suggest MT imaging protocols may be developed into a useful tool for polymer gel dosimetry.

Assessment of image quality of a standard and two dose-reducing protocols in paediatric pelvic CT

Background: Concerns exist regarding the effect of radiation dose from paediatric pelvic CT scans and the potential later risk of radiation-induced neoplasm and teratogenic outcomes in these patients. Objective: To assess the diagnostic quality of CT images of the paediatric pelvis using either reduced mAs or increased pitch compared with standard settings. Materials and methods: A prospective study of pelvic CT scans of 105 paediatric patients was performed using one of three protocols: (1) 31 at a standard protocol of 200 mA with rotation time of 0.75 s at 120 kVp and a pitch factor approximating 1.4; (2) 31 at increased pitch factor approaching 2 and 200 mA; and (3) 43 at a reduced setting of 100 mA and a pitch factor of 1.4. All other settings remained the same in all three groups. Image quality was assessed by radiologists blinded to the protocol used in each scan. Results: No significant difference was found between the quality of images acquired at standard settings and those acquired at half the standard mAs. The use of increased pitch factor resulted in a higher proportion of poor images. Conclusions: Images acquired at 120 kVp using 75 mAs are equivalent in diagnostic quality to those acquired at 150 mAs. Reduced settings can provide useful imaging of the paediatric pelvis and should be considered as a standard protocol in these situations.