New graph-based algorithms to efficiently solve large scale open pit mining optimisation problems

In the mining optimisation literature, most researchers focused on two strategic-level and tactical-level open-pit mine optimisation problems, which are respectively termed ultimate pit limit (UPIT) or constrained pit limit (CPIT). However, many researchers indicate that the substantial numbers of variables and constraints in real-world instances (e.g., with 50-1000 thousand blocks) make the CPIT’s mixed integer programming (MIP) model intractable for use. Thus, it becomes a considerable challenge to solve the large scale CPIT instances without relying on exact MIP optimiser as well as the complicated MIP relaxation/decomposition methods. To take this challenge, two new graph-based algorithms based on network flow graph and conjunctive graph theory are developed by taking advantage of problem properties. The performance of our proposed algorithms is validated by testing recent large scale benchmark UPIT and CPIT instances’ datasets of MineLib in 2013. In comparison to best known results from MineLib, it is shown that the proposed algorithms outperform other CPIT solution approaches existing in the literature. The proposed graph-based algorithms leads to a more competent mine scheduling optimisation expert system because the third-party MIP optimiser is no longer indispensable and random neighbourhood search is not necessary.

Stability of pargasite during ultrahigh-temperature metamorphism:A consequence of titanium and REE partitioning

Orthopyroxene-clinopyroxene-plagioclase needles and symplectite along the cleavage planes and grain boundaries of fluorine-bearing titanian-ferroan pargasite from the Highland Complex, Sri Lanka, are interpreted as evidence for dehydration melting at ultrahigh-temperature conditions. High Ti (up to 0.4 pfu) and F (XF up to 0.56) content in pargasite extends its stability to higher temperatures, and the composition indicates the dehydration melting reaction may take place at ultrahigh-temperatures (~950 °C) at a pressure around 10 kbar, close to peak metamorphic conditions. The increase of Ti content close to the grain boundaries and cleavage planes in pargasite indicates titanium partitioning from the melt during dehydration melting enhanced the stability of the mineral toward ultrahigh-temperature conditions. The REE content in the pargasite shows a similar behavior to that of titanium. The cores with no breakdown assemblage consist of low and flat REE concentrations with respect to the high and Eu-depleted rim. Clinopyroxene in symplectite and needle-shaped lamellae within the pargasite porphyroblasts have similar REE patterns with slightly low-concentrations relative to that of pargasite. In the breakdown assemblage, LREEs are partitioned mainly into plagioclase while the HREEs are partitioned into orthopyroxene. The REE enrichment in the pargasite rims signals their relative partitioning between pargasite rims and melt. Modeling of the partitioning of Ti and REEs associated with pargasite breakdown demonstrates that its stability is greatly enhanced at UHT conditions. This investigation implies that the stability of hydrous minerals such as amphibole can be extended to UHT conditions, and expands our knowledge of metamorphism in the lower crust.

Automated topometric graph generation from floor plan analysis

The world is rich with information such as signage and maps to assist humans to navigate. We present a method to extract topological spatial information from a generic bitmap floor plan and build a topometric graph that can be used by a mobile robot for tasks such as path planning and guided exploration. The algorithm first detects and extracts text in an image of the floor plan. Using the locations of the extracted text, flood fill is used to find the rooms and hallways. Doors are found by matching SURF features and these form the connections between rooms, which are the edges of the topological graph. Our system is able to automatically detect doors and differentiate between hallways and rooms, which is important for effective navigation. We show that our method can extract a topometric graph from a floor plan and is robust against ambiguous cases most commonly seen in floor plans including elevators and stairwells.

Modeling of a pressure modulated desalination system using bond graph methodology

A desalination system is a complex multi energy domain system comprising power/energy flow across several domains such as electrical, thermal, and hydraulic. The dynamic modeling of a desalination system that comprehensively addresses all these multi energy domains is not adequately addressed in the literature. This paper proposes to address the issue of modeling the various energy domains for the case of a single stage flash evaporation desalination system. This paper presents a detailed bond graph modeling of a desalination unit with seamless integration of the power flow across electrical, thermal, and hydraulic domains. The paper further proposes a performance index function that leads to the tracking of the optimal chamber pressure giving the optimal flow rate for a given unit of energy expended. The model has been validated in steady state conditions by simulation and experimentation.

ESR measurements of the partitioning of some new spin probes in n-octanol-water

The stable free radical 1,1,3,3-tetramethylisoindolin-2-yloxyl (TMIO) has proved to be very suitable for use as a spin probe for a number of applications. Because it is soluble mainly in non-polar liquids, there is a need for new derivatives that can be used in a variety of environments. This has been done by introducing substituents in the 5-position of the aromatic ring, namely carboxyl (CTMIO), trimethylamino (TMTMIOI) and sodium sulphonate (NaTMIOS). An accurate ESR method was developed for the measurement of partition coefficients in n-octanol–water. For comparison purposes the method was also applied to some Tempo derivatives. The effect of temperature on the rotational correlation times and the nitrogen-14 hyperfine coupling constant of some of the spin probes was investigated. There is evidence for dimerization of CTMIO to form a biradical

Bond graph toolbox for handling complex variable

Bond graph is an apt modelling tool for any system working across multiple energy domains. Power electronics system modelling is usually the study of the interplay of energy in the domains of electrical, mechanical, magnetic and thermal. The usefulness of bond graph modelling in power electronic field has been realised by researchers. Consequently in the last couple of decades, there has been a steadily increasing effort in developing simulation tools for bond graph modelling that are specially suited for power electronic study. For modelling rotating magnetic fields in electromagnetic machine models, a support for vector variables is essential. Unfortunately, all bond graph simulation tools presently provide support only for scalar variables. We propose an approach to provide complex variable and vector support to bond graph such that it will enable modelling of polyphase electromagnetic and spatial vector systems. We also introduced a rotary gyrator element and use it along with the switched junction for developing the complex/vector variable's toolbox. This approach is implemented by developing a complex S-function tool box in Simulink inside a MATLAB environment This choice has been made so as to synthesise the speed of S-function, the user friendliness of Simulink and the popularity of MATLAB.

Ontogenetic depth partitioning by juvenile freshwater sawfish (Pristis microdon: Pristidae) in a riverine environment

The freshwater sawfish (Pristis microdon) is a critically endangered elasmobranch. Ontogenetic changes in the habitat use of juvenile P. microdon were studied using acoustic tracking in the Fitzroy River, Western Australia. Habitat partitioning was significant between 0+ (2007 year class) and larger 1+ (2006 year class) P. microdon. Smaller 0+ fish generally occupied shallower water (<0.6 m) compared with 1+ individuals, which mainly occurred in depths >0.6 m. Significant differences in hourly depth use were also revealed. The depth that 1+ P. microdon occupied was significantly influenced by lunar phase with these animals utilising a shallower and narrower depth range during the full moon compared with the new moon. This was not observed in 0+ individuals. Habitat partitioning was likely to be related to predator avoidance, foraging behaviours, and temperature and/or light regimes. The occurrence of 1+ P. microdon in deeper water may also result from a need for greater depths in which to manoeuvre. The present study demonstrates the utility of acoustic telemetry in monitoring P. microdon in a riverine environment. These results demonstrate the need to consider the habitat requirements of different P. microdon cohorts in the strategic planning of natural resources and will aid in the development of management strategies for this species.

Estimation and partitioning of polygenic variation captured by common SNPs for Alzheimer's disease, multiple sclerosis and endometriosis

Common diseases such as endometriosis (ED), Alzheimer's disease (AD) and multiple sclerosis (MS) account for a significant proportion of the health care burden in many countries. Genome-wide association studies (GWASs) for these diseases have identified a number of individual genetic variants contributing to the risk of those diseases. However, the effect size for most variants is small and collectively the known variants explain only a small proportion of the estimated heritability. We used a linear mixed model to fit all single nucleotide polymorphisms (SNPs) simultaneously, and estimated genetic variances on the liability scale using SNPs from GWASs in unrelated individuals for these three diseases. For each of the three diseases, case and control samples were not all genotyped in the same laboratory. We demonstrate that a careful analysis can obtain robust estimates, but also that insufficient quality control (QC) of SNPs can lead to spurious results and that too stringent QC is likely to remove real genetic signals. Our estimates show that common SNPs on commercially available genotyping chips capture significant variation contributing to liability for all three diseases. The estimated proportion of total variation tagged by all SNPs was 0.26 (SE 0.04) for ED, 0.24 (SE 0.03) for AD and 0.30 (SE 0.03) for MS. Further, we partitioned the genetic variance explained into five categories by a minor allele frequency (MAF), by chromosomes and gene annotation. We provide strong evidence that a substantial proportion of variation in liability is explained by common SNPs, and thereby give insights into the genetic architecture of the diseases.

DFL: A Data Flow Language

Many novel computer architectures like array and multiprocessors which achieve high performance through the use of concurrency exploit variations of the von Neumann model of computation. The effective utilization of the machines makes special demands on programmers and their programming languages, such as the structuring of data into vectors or the partitioning of programs into concurrent processes. In comparison, the data flow model of computation demands only that the principle of structured programming be followed. A data flow program, often represented as a data flow graph, is a program that expresses a computation by indicating the data dependencies among operators. A data flow computer is a machine designed to take advantage of concurrency in data flow graphs by executing data independent operations in parallel. In this paper, we discuss the design of a high level language (DFL: Data Flow Language) suitable for data flow computers. Some sample procedures in DFL are presented. The implementation aspects have not been discussed in detail since there are no new problems encountered. The language DFL embodies the concepts of functional programming, but in appearance closely resembles Pascal. The language is a better vehicle than the data flow graph for expressing a parallel algorithm. The compiler has been implemented on a DEC 1090 system in Pascal.

Improving decisions for invasive species management: Reformulation and extensions of the Panetta-Lawes eradication graph

Aim: Effective decisions for managing invasive species depend on feedback about the progress of eradication efforts. Panetta & Lawes. developed the eradograph, an intuitive graphical tool that summarizes the temporal trajectories of delimitation and extirpation to support decision-making. We correct and extend the tool, which was affected by incompatibilities in the units used to measure these features that made the axes impossible to interpret biologically. Location: Victoria, New South Wales and Queensland, Australia. Methods: Panetta and Lawes' approach represented delimitation with estimates of the changes in the area known to be infested and extirpation with changes in the mean time since the last detection. We retain the original structure but propose different metrics that improve biological interpretability. We illustrate the methods with a hypothetical example and real examples of invasion and treatment of branched broomrape (Orobanche ramosa L.) and the guava rust complex (Puccinia psidii (Winter 1884)) in Australia. Results: These examples illustrate the potential of the tool to guide decisions about the effectiveness of search and control activities. Main conclusions: The eradograph is a graphical data summary tool that provides insight into the progress of eradication. Our correction and extension of the tool make it easier to interpret and provide managers with better decision support. © 2013 John Wiley & Sons Ltd.

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.

Growth analysis and biomass partitioning of Cyperus iria in response to rice planting density and nitrogen rate

Cyperus iria is a weed of rice with widespread occurrence throughout the world. Because of concerns about excessive and injudicious use of herbicides, cultural weed management approaches that are safe and economical are needed. Developing such approaches will require a better understanding of weed biology and ecology, as well as of weed response to increases in crop density and nutrition. Knowledge of the effects of nitrogen (N) fertilizer on crop-weed competitive interactions could also help in the development of integrated weed management strategies. The present study was conducted in a screenhouse to determine the effects of rice planting density (0, 5, 10, and 20 plants pot−1) and N rate (0, 50, 100, and 150 kg ha−1) on the growth of C. iria. Tiller number per plant decreased by 73–88%, leaf number by 85–94%, leaf area by 85–98%, leaf biomass by 92–99%, and inflorescence biomass by 96–99% when weed plants were grown at 20 rice plants pot−1 (i.e., 400 plants m−2) compared with weed plants grown alone. All of these parameters increased when N rates were increased. On average, weed biomass increased by 118–389% and rice biomass by 121–275% with application of 50–150 kg N ha−1, compared to control. Addition of N favored weed biomass production relative to rice biomass. Increased N rates reduced the root-to-shoot weight ratio of C. iria. Rice interference reduced weed growth and biomass and completely suppressed C. iria when no N was applied at high planting densities (i.e., 20 plants pot−1). The weed showed phenotypic plasticity in response to N application, and the addition of N increased the competitive ability of the weed over rice at densities of 5 and 10 rice plants pot−1 compared with 20 plants pot−1. The results of the present study suggest that high rice density (i.e., 400 plants m−2) can help suppress C. iria growth even at high N rates (150 kg ha−1).

Denitrification losses from an intensively managed sub-tropical pasture: Impact of soil moisture on the partitioning of N2 and N2O emissions

Intensively managed pastures in subtropical Australia under dairy production are nitrogen (N) loaded agro-ecosystems, with an increased pool of N available for denitrification. The magnitude of denitrification losses and N2:N2O partitioning in these agro-ecosystems is largely unknown, representing a major uncertainty when estimating total N loss and replacement. This study investigated the influence of different soil moisture contents on N2 and N2O emissions from a subtropical dairy pasture in Queensland, Australia. Intact soil cores were incubated over 15 days at 80% and 100% water-filled pore space (WFPS), after the application of 15N labelled nitrate, equivalent to 50 kg N ha−1. This setup enabled the direct quantification of N2 and N2O emissions following fertilisation using the 15N gas flux method. The main product of denitrification in both treatments was N2. N2 emissions exceeded N2O emissions by a factor of 8 ± 1 at 80% WFPS and a factor of 17 ± 2 at 100% WFPS. The total amount of N-N2 lost over the incubation period was 21.27 kg ± 2.10 N2-N ha−1 at 80% WFPS and 25.26 kg ± 2.79 kg ha−1 at 100% WFPS respectively. N2 emissions remained high at 100% WFPS, while related N2O emissions decreased. At 80% WFPS, N2 emissions increased constantly over time while N2O fluxes declined. Consequently, N2/(N2 + N2O) product ratios increased over the incubation period in both treatments. N2/(N2 + N2O) product ratios responded significantly to soil moisture, confirming WFPS as a key driver of denitrification. The substantial amount of fertiliser lost as N2 reveals the agronomic significance of denitrification as a major pathway of N loss for sub-tropical pastures at high WFPS and may explain the low fertiliser N use efficiency observed for these agro-ecosystems.

Acoustic emission source location and damage detection in a metallic structure using a graph-theory-based geodesic pproach

A geodesic-based approach using Lamb waves is proposed to locate the acoustic emission (AE) source and damage in an isotropic metallic structure. In the case of the AE (passive) technique, the elastic waves take the shortest path from the source to the sensor array distributed in the structure. The geodesics are computed on the meshed surface of the structure using graph theory based on Dijkstra's algorithm. By propagating the waves in reverse virtually from these sensors along the geodesic path and by locating the first intersection point of these waves, one can get the AE source location. The same approach is extended for detection of damage in a structure. The wave response matrix of the given sensor configuration for the healthy and the damaged structure is obtained experimentally. The healthy and damage response matrix is compared and their difference gives the information about the reflection of waves from the damage. These waves are backpropagated from the sensors and the above method is used to locate the damage by finding the point where intersection of geodesics occurs. In this work, the geodesic approach is shown to be suitable to obtain a practicable source location solution in a more general set-up on any arbitrary surface containing finite discontinuities. Experiments were conducted on aluminum specimens of simple and complex geometry to validate this new method.

Using a glow graph to represent data flow and dependency in event logs

The idea of extracting knowledge in process mining is a descendant of data mining. Both mining disciplines emphasise data flow and relations among elements in the data. Unfortunately, challenges have been encountered when working with the data flow and relations. One of the challenges is that the representation of the data flow between a pair of elements or tasks is insufficiently simplified and formulated, as it considers only a one-to-one data flow relation. In this paper, we discuss how the effectiveness of knowledge representation can be extended in both disciplines. To this end, we introduce a new representation of the data flow and dependency formulation using a flow graph. The flow graph solves the issue of the insufficiency of presenting other relation types, such as many-to-one and one-to-many relations. As an experiment, a new evaluation framework is applied to the Teleclaim process in order to show how this method can provide us with more precise results when compared with other representations.