Efficient simulation of unsaturated flow in heterogeneous soils using an exponential integrator

- Provided a practical variable-stepsize implementation of the exponential Euler method (EEM). - Introduced a new second-order variant of the scheme that enables the local error to be estimated at the cost of a single additional function evaluation. - New EEM implementation outperformed sophisticated implementations of the backward differentiation formulae (BDF) of order 2 and was competitive with BDF of order 5 for moderate to high tolerances.

Topographical, structural and geophysical characterization of fracture zones: implications for groundwater flow and vulnerability

The main objective of this study is to evaluate selected geophysical, structural and topographic methods on regional, local, and tunnel and borehole scales, as indicators of the properties of fracture zones or fractures relevant to groundwater flow. Such information serves, for example, groundwater exploration and prediction of the risk of groundwater inflow in underground construction. This study aims to address how the features detected by these methods link to groundwater flow in qualitative and semi-quantitative terms and how well the methods reveal properties of fracturing affecting groundwater flow in the studied sites. The investigated areas are: (1) the Päijänne Tunnel for water-conveyance whose study serves as a verification of structures identified on regional and local scales; (2) the Oitti fuel spill site, to telescope across scales and compare geometries of structural assessment; and (3) Leppävirta, where fracturing and hydrogeological environment have been studied on the scale of a drilled well. The methods applied in this study include: the interpretation of lineaments from topographic data and their comparison with aeromagnetic data; the analysis of geological structures mapped in the Päijänne Tunnel; borehole video surveying; groundwater inflow measurements; groundwater level observations; and information on the tunnel s deterioration as demonstrated by block falls. The study combined geological and geotechnical information on relevant factors governing groundwater inflow into a tunnel and indicators of fracturing, as well as environmental datasets as overlays for spatial analysis using GIS. Geophysical borehole logging and fluid logging were used in Leppävirta to compare the responses of different methods to fracturing and other geological features on the scale of a drilled well. Results from some of the geophysical measurements of boreholes were affected by the large diameter (gamma radiation) or uneven surface (caliper) of these structures. However, different anomalies indicating more fractured upper part of the bedrock traversed by well HN4 in Leppävirta suggest that several methods can be used for detecting fracturing. Fracture trends appear to align similarly on different scales in the zone of the Päijänne Tunnel. For example, similarities of patterns were found between the regional magnetic trends, correlating with orientations of topographic lineaments interpreted as expressions of fracture zones. The same structural orientations as those of the larger structures on local or regional scales were observed in the tunnel, even though a match could not be made in every case. The size and orientation of the observation space (patch of terrain at the surface, tunnel section, or borehole), the characterization method, with its typical sensitivity, and the characteristics of the location, influence the identification of the fracture pattern. Through due consideration of the influence of the sampling geometry and by utilizing complementary fracture characterization methods in tandem, some of the complexities of the relationship between fracturing and groundwater flow can be addressed. The flow connections demonstrated by the response of the groundwater level in monitoring wells to pressure decrease in the tunnel and the transport of MTBE through fractures in bedrock in Oitti, highlight the importance of protecting the tunnel water from a risk of contamination. In general, the largest values of drawdown occurred in monitoring wells closest to the tunnel and/or close to the topographically interpreted fracture zones. It seems that, to some degree, the rate of inflow shows a positive correlation with the level of reinforcement, as both are connected with the fracturing in the bedrock. The following geological features increased the vulnerability of tunnel sections to pollution, especially when several factors affected the same locations: (1) fractured bedrock, particularly with associated groundwater inflow; (2) thin or permeable overburden above fractured rock; (3) a hydraulically conductive layer underneath the surface soil; and (4) a relatively thin bedrock roof above the tunnel. The observed anisotropy of the geological media should ideally be taken into account in the assessment of vulnerability of tunnel sections and eventually for directing protective measures.

Uniformity of film thickness on rotating planetary planar substrates

For the purposes of obtaining a number of components with nearly identical thickness distributions over the substrate area and of minimizing the inhomogeneities of the film, it is logical to presume that a substrate rotating on its own axis and revolving around another axis will give more uniformity in film thickness than a substrate only revolving around one axis. In relation to the practical applications, an investigation has been undertaken to study the refinement that can be achieved by using a planar planetary substrate holder. It is shown theoretically that the use of the planetary substrate holder under ideal conditions of source and geometry does not offer any further improvement in uniformity of thickness over the conventional rotary work-holder. It is also shown that the geometrical parameters alone have little influence over the uniformity achieved on a planetary substrate, because of the complex cyclidal motion of any point on it. However, for any given geometry, a non-integral speed ratio of the planetary substrate and the work-holder shows considerably less variation in thickness over the substrate area.

Flow-graph techniques for epicyclic gear train analysis

Flow-graph techniques are applied in this article for the analysis of an epicyclic gear train. A gear system based on this is designed and constructed for use in Numerical Control Systems.

A parametric differentiation version with finite-difference scheme applicable to a class of problems in boundary layer flow with massive blowing

A numerical procedure, based on the parametric differentiation and implicit finite difference scheme, has been developed for a class of problems in the boundary-layer theory for saddle-point regions. Here, the results are presented for the case of a three-dimensional stagnation-point flow with massive blowing. The method compares very well with other methods for particular cases (zero or small mass blowing). Results emphasize that the present numerical procedure is well suited for the solution of saddle-point flows with massive blowing, which could not be solved by other methods.

Gene flow from transgenic plant populations: models and applications for risk assessment

The future use of genetically modified (GM) plants in food, feed and biomass production requires a careful consideration of possible risks related to the unintended spread of trangenes into new habitats. This may occur via introgression of the transgene to conventional genotypes, due to cross-pollination, and via the invasion of GM plants to new habitats. Assessment of possible environmental impacts of GM plants requires estimation of the level of gene flow from a GM population. Furthermore, management measures for reducing gene flow from GM populations are needed in order to prevent possible unwanted effects of transgenes on ecosystems. This work develops modeling tools for estimating gene flow from GM plant populations in boreal environments and for investigating the mechanisms of the gene flow process. To describe spatial dimensions of the gene flow, dispersal models are developed for the local and regional scale spread of pollen grains and seeds, with special emphasis on wind dispersal. This study provides tools for describing cross-pollination between GM and conventional populations and for estimating the levels of transgenic contamination of the conventional crops. For perennial populations, a modeling framework describing the dynamics of plants and genotypes is developed, in order to estimate the gene flow process over a sequence of years. The dispersal of airborne pollen and seeds cannot be easily controlled, and small amounts of these particles are likely to disperse over long distances. Wind dispersal processes are highly stochastic due to variation in atmospheric conditions, so that there may be considerable variation between individual dispersal patterns. This, in turn, is reflected to the large amount of variation in annual levels of cross-pollination between GM and conventional populations. Even though land-use practices have effects on the average levels of cross-pollination between GM and conventional fields, the level of transgenic contamination of a conventional crop remains highly stochastic. The demographic effects of a transgene have impacts on the establishment of trangenic plants amongst conventional genotypes of the same species. If the transgene gives a plant a considerable fitness advantage in comparison to conventional genotypes, the spread of transgenes to conventional population can be strongly increased. In such cases, dominance of the transgene considerably increases gene flow from GM to conventional populations, due to the enhanced fitness of heterozygous hybrids. The fitness of GM plants in conventional populations can be reduced by linking the selectively favoured primary transgene to a disfavoured mitigation transgene. Recombination between these transgenes is a major risk related to this technique, especially because it tends to take place amongst the conventional genotypes and thus promotes the establishment of invasive transgenic plants in conventional populations.

Experimental and numerical investigation of jet injection in a wall bounded supersonic flow

Numerical and experimental studies of a supersonic jet (Helium) inclined at 45 degrees to a oncoming Mach 2 flow have been carried out. The numerical study has been used to arrive at a geometry that could reduce an oncoming Mach 5.75 flow to Mach 2 flow and in determining the jet parameters. Experiments are carried out in the IISc. hypersonic shock tunnel HST2 at similar conditions obtained from numerical studies. Flow visualization studies carried out using Schlieren technique clearly show the presence of the bow shock in front of the jet exposed to supersonic cross flow. The jet Mach number is experimentally found to be approximate to 3. Visual observations show that the jet has penetrated up to 60% of the total height of the chamber.

Investigations of separated flow over backward facing steps in IISc hypersonic shock tunnel

Two backward facing step (2 mm and 3 mm step height) models are selected for surface heat transfer measurements. The platinum thin film gauges are deposited on the Macor inserts using both hand paint and vacuum sputtering technique. Using the Eckert reference temperature method the heating rates has been theoretically calculated along the flat plate portion of the model and the theoretical estimates are compared with experimentally determined surface heat transfer rate. Theoretical analysis of heat flux distribution down stream of the backward facing step model has been carried out using Gai’s non-dimensional analysis. Based on the measured surface heating rates on the backward facing step, the reattachment distance is estimated for 2 and 3 mm step height at nominal Mach number of 7.6. It has been found from the present study that for 2 and 3 mm step height, it approximately takes about 10 and 8 step heights downstream of the model respectively for the flow to re-attach.

Magnetohydrodynamic flow between torsionally oscillating eccentric disks

The unsteady MHD flow of an incompressible, viscous electrically conducting fluid contained between two torsionally oscillating eccentric disks has been investigated. The state of uniform rotation of the central region visualised in the steady flow is seen to be absent in the case of oscillatory flow.

Gene flow from Corymbia hybrids in northern New South Wales

Some of the most productive taxa for forestry are interspecific F1 hybrids grown as exotics in the tropics and subtropics. Attributes of resilience, adaptability and vigour which engender the hybrids for wood production, may also exacerbate the risk they present from gene flow to native species gene pools or to local ecologies as weeds. To determine the biological and genetic factors that influence the extent of hybridisation, we examine the distribution and genealogy of wildlings surrounding plantings of locally-exotic Corymbia torelliana (Section Cadageria) near native C. henryi (Section Maculatae) in northern New South Wales. Our study showed pre-mating and pre- and post-zygotic barriers were incomplete, with in situ generation and natural establishment of both F1 hybrids (n = 3) and advanced generation hybrids under the disturbed conditions bordering native forest. As hybrids were located on alluvial flats exposed to frost, they also likely have an extended ecological range relative to native C. henryi. Despite the likely generation of large viable seed crops on F1 trees at the site over many years, establishment success and survival of advanced generation hybrids may be low, as only 5 immature and no mature advanced generation hybrids were identified. Propagation and genetic analysis of a seed crop from one F1 wildling showed early survival and vigour of seedlings in cultivation was high, and that at least for some F1 in some seasons, backcrossing to the recurrent native C. henryi parent is favoured (60%), whereas selfing (10%) and crossing with other F1 (30%) was less frequent. Transport of seed by stingless bees probably accounted for long distance dispersal from C. torelliana, but this mechanism does not appear to supplement gravity-dispersal of seed from the F1. Coupled with other evidence from studies of bee behaviour, controlled pollination in Corymbia sp., and long-term fitness in second generation eucalypt hybrids, we anticipate gene flow via pollen rather than seed will be the greater challenge for managing the risk of introgression of C. torelliana ancestry into native species from the planted F1 hybrid. If large sources of F1 pollen become available to compete with native pollen, gene flow will probably be frequent and hybrids may establish in disturbed conditions and in habitats beyond the ecological range of their native parent. Further study is needed to determine the degree to which outbreeding depression and poor survival inhibits on-going gene flow.

Unsteady laminar boundary layers in a compressible stagnation flow

The unsteady laminar compressible boundary-layer flow in the immediate vicinity of a two-dimensional stagnation point due to an incident stream whose velocity varies arbitrarily with time is considered. The governing partial differential equations, involving both time and the independent similarity variable, are transformed into new co-ordinates with finite ranges by means of a transformation which maps an infinite interval into a finite one. The resulting equations are solved by converting them into a matrix equation through the application of implicit finite-difference formulae. Computations have been carried out for two particular unsteady free-stream velocity distributions: (1) a constantly accelerating stream and (2) a fluctuating stream. The results show that in the former case both the skin-friction and the heat-transfer parameter increase steadily with time after a certain instant, while in the latter they oscillate thus responding to the fluctuations in the free-stream velocity.