Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.

Study of carotid arterial plaque stress for symptomatic and asymptomatic patients

Stroke is one of the leading causes of death in the world, resulting mostly from the sudden ruptures of atherosclerosis carotid plaques. Until now, the exact plaque rupture mechanism has not been fully understood, and also the plaque rupture risk stratification. The advanced multi-spectral magnetic resonance imaging (MRI) has allowed the plaque components to be visualized in-vivo and reconstructed by computational modeling. In the study, plaque stress analysis using fully coupled fluid structure interaction was applied to 20 patients (12 symptomatic and 8 asymptomatic) reconstructed from in-vivo MRI, followed by a detailed biomechanics analysis, and morphological feature study. The locally extreme stress conditions can be found in the fibrous cap region, 85% at the plaque shoulder based on the present study cases. Local maximum stress values predicted in the plaque region were found to be significantly higher in symptomatic patients than that in asymptomatic patients (200±43. kPa vs. 127±37. kPa, p=0.001). Plaque stress level, defined by excluding 5% highest stress nodes in the fibrous cap region based on the accumulative histogram of stress experienced on the computational nodes in the fibrous cap, was also significantly higher in symptomatic patients than that in asymptomatic patients (154±32. kPa vs. 111±23. kPa, p<0.05). Although there was no significant difference in lipid core size between the two patient groups, symptomatic group normally had a larger lipid core and a significantly thinner fibrous cap based on the reconstructed plaques using 3D interpolation from stacks of 2D contours. Plaques with a higher stenosis were more likely to have extreme stress conditions upstream of plaque throat. The combined analyses of plaque MR image and plaque stress will advance our understanding of plaque rupture, and provide a useful tool on assessing plaque rupture risk.

Solid volume fraction estimation of bone:marrow replica models using ultrasound transit time spectroscopy

The acceptance of broadband ultrasound attenuation (BUA) for the assessment of osteoporosis suffers from a limited understanding of both ultrasound wave propagation through cancellous bone and its exact dependence upon the material and structural properties. It has recently been proposed that ultrasound wave propagation in cancellous bone may be described by a concept of parallel sonic rays; the transit time of each ray defined by the proportion of bone and marrow propagated. A Transit Time Spectrum (TTS) describes the proportion of sonic rays having a particular transit time, effectively describing the lateral inhomogeneity of transit times over the surface aperture of the receive ultrasound transducer. The aim of this study was to test the hypothesis that the solid volume fraction (SVF) of simplified bone:marrow replica models may be reliably estimated from the corresponding ultrasound transit time spectrum. Transit time spectra were derived via digital deconvolution of the experimentally measured input and output ultrasonic signals, and compared to predicted TTS based on the parallel sonic ray concept, demonstrating agreement in both position and amplitude of spectral peaks. Solid volume fraction was calculated from the TTS; agreement between true (geometric calculation) with predicted (computer simulation) and experimentally-derived values were R2=99.9% and R2=97.3% respectively. It is therefore envisaged that ultrasound transit time spectroscopy (UTTS) offers the potential to reliably estimate bone mineral density and hence the established T-score parameter for clinical osteoporosis assessment.

Influence of time step in the simulation modelling of evapotranspiration

Three simulations of evapotranspiration were done with two values of time step,viz 10 min and one day. Inputs to the model were weather data, including directly measured upward and downward radiation, and soil characteristics. Three soils were used for each simulation. Analysis of the results shows that the time step has a direct influence on the prediction of potential evapotranspiration, but a complex interaction of this effect with the soil moisture characteristic, rate of increase of ground cover and bare soil evaporation determines the actual transpiration predicted. The results indicate that as small a time step as possible should be used in the simulation.

Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment This fraction is the result of undesirable genotype-by-environment interactions (G x E) and measured by the genetic correlation (r(g)) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of G x E over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels

Analytical studies of gain optimization in CO2-N2 gasdynamic lasers employing two-dimensional wedge nozzles

An analytical method has been proposed to optimise the small-signaloptical gain of CO2-N2 gasdynamic lasers (gdl) employing two-dimensional (2D) wedge nozzles. Following our earlier work the equations governing the steady, inviscid, quasi-one-dimensional flow in the wedge nozzle of thegdl are reduced to a universal form so that their solutions depend on a single unifying parameter. These equations are solved numerically to obtain similar solutions for the various flow quantities, which variables are subsequently used to optimize the small-signal-gain. The corresponding optimum values like reservoir pressure and temperature and 2D nozzle area ratio also have been predicted and graphed for a wide range of laser gas compositions, with either H2O or He as the catalyst. A large number of graphs are presented which may be used to obtain the optimum values of small signal gain for a wide range of laser compositions without further computations.

Use of S1 nuclease and formamide in combination for the reassociation studies on GC-rich DNA

The optimal parameters in the use of nuclease S1 in DNA reassociation kinetics in the presence of formamide have been determined. The conditions are especially suitable for the study of DNA rich in mole percent GC. A 10-fold dilution of the reassociation samples leading to a decrease in both NaCl and formamide concentrations, consequently resulting in a lowering of Tm by only 1.5°C, and the S1 digestion at temperatures identical to the reassociation assay in order to retain the stability of the duplex, are two important aspects of this system. Under these conditions, the kinetics of reassociation followed the theoretically predicted pattern, while the earlier reported methods have shown lower values.

Stock-recruitment relationships of the tiger prawns (Penaeus esculentus and Penaeus semisulcatus) in the Australian Northern Prawn Fishery

This paper investigates the stock-recruitment and equilibrium yield dynamics for the two species of tiger prawns (Penaeus esculentus and Penaeus semisulcatus) in Australia's most productive prawn fishery: the Northern Prawn Fishery. Commercial trawl logbooks for 1970-93 and research surveys are used to develop population models for these prawns. A population model that incorporates continuous recruitment is developed. Annual spawning stock and recruitment indices are then estimated from the population model. Spawning stock indices represent the abundance of female prawns that are likely to spawn; recruitment indices represent the abundance of all prawns less than a certain size. The relationships between spawning stock and subsequent recruitment (SRR), between recruitment and subsequent spawning stock (RSR), and between recruitment and commercial catch were estimated through maximum-likelihood models that incorporated autoregressive terms. Yield as a function of fishing effort was estimated by constraining to equilibrium the SRR and RSR. The resulting production model was then used to determine maximum sustainable yield (MSY) and its corresponding fishing effort (f(MSY)). Long-term yield estimates for the two tiger prawn species range between 3700 and 5300 t. The fishing effort at present is close to the level that should produce MSY for both species of tiger prawns. However, current landings, recruitment and spawning stock are below the equilibrium values predicted by the models. This may be because of uncertainty in the spawning stock-recruitment relationships, a change in carrying capacity, biased estimates of fishing effort, unreliable catch statistics, or simplistic assumptions about stock structure. Although our predictions of tiger prawn yields are uncertain, management will soon have to consider new measures to counteract the effects of future increases in fishing effort.

Phenotypic and genotypic characterisation of Pasteurella multocida isolated from pigs at slaughter in New Zealand

AIMS: To examine pigs at slaughter in New Zealand for the presence of Pasteurella multocida, and to determine for isolates, their biochemical profi les, somatic and capsular types, and the presence or absence of the HSB and toxA genes, associated with haemorrhagic septicaemia (HS) and progressive atrophic rhinitis (PAR), respectively. METHODS: Swabs from 173 lungs, 158 palatine tonsils and 82 nasal passages of pigs at two abattoirs in New Zealand were cultured for P. multocida using conventional techniques, and isolated colonies were subjected to biochemical tests for identi- fi cation of biovars. Somatic serotyping was conducted using an agar gel immunodiffusion (AGID) test. Polymerase chain reaction (PCR) assays were used to confi rm phenotypic identifi cation of colonies using species-specifi c primers, capsule type using serogroup-specifi c primers and multiplex PCR, and to test for the presence of HSB and toxA genes. RESULTS: Pasteurella multocida was isolated from 11/173 (6.4%) lung, 32/158 (20.2%) palatine tonsil and 5/82 (6.1 %) nasal swab samples, a total of 48 isolates from 413 samples (11.6%). Isolation rates per farm ranged from 1–53% of tissue samples collected from pigs 5–6 months of age. On phenotypic characterisation, isolates were allocated to seven main biovars, viz 1, 2, 3, 5, 9, 12, and a dulcitol-negative variant of Biovar 8, the majority (30/48) being Biovar 3. Of the 42 isolates for which somatic serotyping was conducted, 10% were Serovar 1, 79% were Serovar 3, 2% were Serovar 6,1, 2% were Serovar 12, and 7% could not be typed. All 48 isolates were confi rmed as P. multocida using a species-specifi c PCR. In the capsular multiplex PCR, 92% of isolates were Capsular (Cap) type A, 2% were Cap D, and 6% could not be typed. None of the samples were positive for the HSB or toxA genes. CONCLUSION: Serovars or capsular types of P. multocida associated with HS or PAR in pigs were not detected. Establishment of species-specifi c, capsular and toxin PCR assays allowed the rapid screening of isolates of P. multocida, while serotyping provided an additional tool for epidemiological and tracing purposes.

Runoff, sediment loss and water quality from forest roads in a southeast Queensland coastal plain Pinus plantation

Runoff and sediment loss from forest roads were monitored for a two-year period in a Pinus plantation in southeast Queensland. Two classes of road were investigated: a gravelled road, which is used as a primary daily haulage route for the logging area, and an ungravelled road, which provides the main access route for individual logging compartments and is intensively used as a haulage route only during the harvest of these areas (approximately every 30 years). Both roads were subjected to routine traffic loads and maintenance during the study. Surface runoff in response to natural rainfall was measured and samples taken for the determination of sediment and nutrient (total nitrogen, total phosphorus, dissolved organic carbon and total iron) loads from each road. Results revealed that the mean runoff coefficient (runoff depth/rainfall depth) was consistently higher from the gravelled road plot with 0.57, as compared to the ungravelled road with 0.38. Total sediment loss over the two-year period was greatest from the gravelled road plot at 5.7 t km−1 compared to the ungravelled road plot with 3.9 t km−1. Suspended solids contributed 86% of the total sediment loss from the gravelled road, and 72% from the ungravelled road over the two years. Nitrogen loads from the two roads were both relatively constant throughout the study, and averaged 5.2 and 2.9 kg km−1 from the gravelled and ungravelled road, respectively. Mean annual phosphorus loads were 0.6 kg km−1 from the gravelled road and 0.2 kg km−1 from the ungravelled road. Organic carbon and total iron loads increased in the second year of the study, which was a much wetter year, and are thought to reflect the breakdown of organic matter in roadside drains and increased sediment generation, respectively. When road and drain maintenance (grading) was performed runoff and sediment loss were increased from both road types. Additionally, the breakdown of the gravel road base due to high traffic intensity during wet conditions resulted in the formation of deep (10 cm) ruts which increased erosion. The Water Erosion Prediction Project (WEPP):Road model was used to compare predicted to observed runoff and sediment loss from the two road classes investigated. For individual rainfall events, WEPP:Road predicted output showed strong agreement with observed values of runoff and sediment loss. WEPP:Road predictions for annual sediment loss from the entire forestry road network in the study area also showed reasonable agreement with the extrapolated observed values.

Simulating crop-parasitic weed interactions using APSIM: Model evaluation and application

The parasitic weed Orobanche crenata inflicts major damage on faba bean, lentil, pea and other crops in Mediterranean environments. The development of methods to control O. crenata is to a large extent hampered by the complexity of host-parasite systems. Using a model of host-parasite interactions can help to explain and understand this intricacy. This paper reports on the evaluation and application of a model simulating host-parasite competition as affected by environment and management that was implemented in the framework of the Agricultural Production Systems Simulator (APSIM). Model-predicted faba bean and O. crenata growth and development were evaluated against independent data. The APSIM-Fababean and -Parasite modules displayed a good capability to reproduce effects of pedoclimatic conditions, faba bean sowing date and O. crenata infestation on host-parasite competition. The r(2) values throughout exceeded 0.84 (RMSD: 5.36 days) for phenological, 0.85 (RMSD: 223.00 g m(-2)) for host growth and 0.78 (RMSD: 99.82 g m(-2)) for parasite growth parameters. Inaccuracies of simulated faba bean root growth that caused some bias of predicted parasite number and host yield loss may be dealt with by more flexibly simulating vertical root distribution. The model was applied in simulation experiments to determine optimum sowing windows for infected and non-infected faba bean in Mediterranean environments. Simulation results proved realistic and testified to the capability of APSIM to contribute to the development of tactical approaches in parasitic weed control.

Diet quality and liveweight gain of steers grazing Leucaena-grass pasture estimated with faecal near infrared reflectance spectroscopy (F.NIRS)

Three drafts of Bos indicus cross steers (initially 178-216 kg) grazed Leucaena-grass pasture [Leucaena leucocephala subspecies glabrata cv. Cunningham with green panic (Panicum maximum cv. trichoglume)] from late winter through to autumn during three consecutive years in the Burnett region of south-east Queensland. Measured daily weight gain (DWGActual) of the steers was generally 0.7-1.1 kg/day during the summer months. Estimated intakes of metabolisable energy and dry matter (DM) were calculated from feeding standards as the intakes required by the steers to grow at the DWGActual. Diet attributes were predicted from near infrared reflectance spectroscopy spectra of faeces (F.NIRS) using established calibration equations appropriate for northern Australian forages. Inclusion of some additional reference samples from cattle consuming Leucaena diets into F.NIRS calibrations based on grass and herbaceous legume-grass pastures improved prediction of the proportion of Leucaena in the diet. Mahalanobis distance values supported the hypothesis that the F.NIRS predictions of diet crude protein concentration and DM digestibility (DMD) were acceptable. F.NIRS indicated that the percentage of Leucaena in the diet varied widely (10-99%). Diet crude protein concentration and DMD were usually high, averaging 12.4 and 62%, respectively, and were related asymptotically to the percentage of Leucaena in the diet (R2 = 0.48 and 0.33, respectively). F.NIRS calibrations for DWG were not satisfactory to predict this variable from an individual faecal sample since the s.e. of prediction were 0.33-0.40 kg/day. Cumulative steer liveweight (LW) predicted from F.NIRS DWG calibrations, which had been previously developed with tropical grass and grass-herbaceous legume pastures, greatly overestimated the measured steer LW; therefore, these calibrations were not useful. Cumulative steer LW predicted from a modified F.NIRS DWG calibration, which included data from the present study, was strongly correlated (R2 = 0.95) with steer LW but overestimated LW by 19-31 kg after 8 months. Additional reference data are needed to develop robust F.NIRS calibrations to encompass the diversity of Leucaena pastures of northern Australia. In conclusion, the experiment demonstrated that F.NIRS could improve understanding of diet quality and nutrient intake of cattle grazing Leucaena-grass pasture, and the relationships between nutrient supply and cattle growth.

Architectural modelling of maize under water stress

Two field experiments using maize (Pioneer 31H50) and three watering regimes [(i) irrigated for the whole crop cycle, until anthesis, (ii) not at all (experiment 1) and (iii) fully irrigated and rain grown for the whole crop cycle (experiment 2)] were conducted at Gatton, Australia, during the 2003-04 season. Data on crop ontogeny, leaf, sheath and internode lengths and leaf width, and senescence were collected at 1- to 3-day intervals. A glasshouse experiment during 2003 quantified the responses of leaf shape and leaf presentation to various levels of water stress. Data from experiment 1 were used to modify and parameterise an architectural model of maize (ADEL-Maize) to incorporate the impact of water stress on maize canopy characteristics. The modified model produced accurate fitted values for experiment 1 for final leaf area and plant height, but values during development for leaf area were lower than observed data. Crop duration was reasonably well fitted and differences between the fully irrigated and rain-grown crops were accurately predicted. Final representations of maize crop canopies were realistic. Possible explanations for low values of leaf area are provided. The model requires further development using data from the glasshouse study and before being validated using data from experiment 2 and other independent data. It will then be used to extend functionality in architectural models of maize. With further research and development, the model should be particularly useful in examining the response of maize production to water stress including improved prediction of total biomass and grain yield. This will facilitate improved simulation of plant growth and development processes allowing investigation of genotype by environment interactions under conditions of suboptimal water supply.

Seed persistence in the field may be predicted by laboratory-controlled aging

Weed management is complicated by the presence of soil seed banks. The complexity of soil-seed interactions means that seed persistence in the field is often difficult to measure, let alone predict. Field trials, although accurate in their context, are time-consuming and expensive to conduct for individual species. Some ex situ techniques for estimating seed life expectancy have been proposed, but these fail to simulate the environmental complexity of the field. Also, it has been questioned whether techniques such as the controlled aging test (CAT) are useful indicators of field persistence. This study aimed to test the validity of the standard CAT (seed aging at 45 C and 60% relative humidity) in use at the Royal Botanic Gardens, Kew, U.K., for predicting field seed-persistence. Comparison of seed persistence and CAT data for 27 northwest European species suggested a significant positive correlation of 0.31. Subsequently, 13 species of emerging and common weeds of Queensland were assessed for their seed longevity using the CAT. The seed longevity data of these species in the CAT were linked with field seed-persistence data according to three broad seed-persistence categories: <1 yr, 1 to 3 yr, and >3 yr. We discuss the scope for using the CAT as a tool for rapid assignment of species to these categories. There is a need for further studies that compare predictions of seed persistence based on the CAT with seed persistence in the field for a larger range of species and environments.

Interfacial area measurement in a gas - liquid ejector for a sodium chloride - air system

Interfacial area measurement has been carried out experimentally by measuring the bubble size and holdup for air-sodium chloride solution system. The size of the bubble is predominantly established by the air hold up. High speed photography technique for bubble size measurement and gamma ray attenuation method for holdup measurements are followed. The measured values are compared with the theoretically predicted values. Interracial area as a function of the liquid flow rate and also its distance from the nozzle of the ejector has been reported in this paper. The results obtained for this non-reactive system are also compared with those of air-water system.