A new adaptive multiple modelling approach for non-linear and non-stationary systems

This paper proposes a novel adaptive multiple modelling algorithm for non-linear and non-stationary systems. This simple modelling paradigm comprises K candidate sub-models which are all linear. With data available in an online fashion, the performance of all candidate sub-models are monitored based on the most recent data window, and M best sub-models are selected from the K candidates. The weight coefficients of the selected sub-model are adapted via the recursive least square (RLS) algorithm, while the coefficients of the remaining sub-models are unchanged. These M model predictions are then optimally combined to produce the multi-model output. We propose to minimise the mean square error based on a recent data window, and apply the sum to one constraint to the combination parameters, leading to a closed-form solution, so that maximal computational efficiency can be achieved. In addition, at each time step, the model prediction is chosen from either the resultant multiple model or the best sub-model, whichever is the best. Simulation results are given in comparison with some typical alternatives, including the linear RLS algorithm and a number of online non-linear approaches, in terms of modelling performance and time consumption.

PMIP4-CMIP6: the contribution of the Paleoclimate Modelling Intercomparison Project to CMIP6

The goal of the Palaeoclimate Modelling Intercomparison Project (PMIP) is to understand the response of the climate system to changes in different climate forcings and to feedbacks. Through comparison with observations of the environmental impacts of these climate changes, or with climate reconstructions based on physical, chemical or biological records, PMIP also addresses the issue of how well state-of-the-art models simulate climate changes. Palaeoclimate states are radically different from those of the recent past documented by the instrumental record and thus provide an out-of-sample test of the models used for future climate projections and a way to assess whether they have the correct sensitivity to forcings and feedbacks. Five distinctly different periods have been selected as focus for the core palaeoclimate experiments that are designed to contribute to the objectives of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). This manuscript describes the motivation for the choice of these periods and the design of the numerical experiments, with a focus upon their novel features compared to the experiments performed in previous phases of PMIP and CMIP as well as the benefits of common analyses of the models across multiple climate states. It also describes the information needed to document each experiment and the model outputs required for analysis and benchmarking.

Modelling the line variations from the wind-wind shock emissions of WR 30a

The study of Wolf-Rayet stars plays an important role in evolutionary theories of massive stars. Among these objects, similar to 20 per cent are known to be in binary systems and can therefore be used for the mass determination of these stars. Most of these systems are not spatially resolved and spectral lines can be used to constrain the orbital parameters. However, part of the emission may originate in the interaction zone between the stellar winds, modifying the line profiles and thus challenging us to use different models to interpret them. In this work, we analysed the He II lambda 4686 angstrom + C IV lambda 4658 angstrom blended lines of WR 30a (WO4+O5) assuming that part of the emission originate in the wind-wind interaction zone. In fact, this line presents a quiescent base profile, attributed to the WO wind, and a superposed excess, which varies with the orbital phase along the 4.6-d period. Under these assumptions, we were able to fit the excess spectral line profile and central velocity for all phases, except for the longest wavelengths, where a spectral line with constant velocity seems to be present. The fit parameters provide the eccentricity and inclination of the binary orbit, from which it is possible to constrain the stellar masses.

Ab initio calculation of the BCC Fe-Al-Mo (Iron-Aluminum-Molybdenum) phase diagram: Implications for the nature of the tau(2) phase

The metastable phase diagram of the BCC-based ordering equilibria in the Fe-Al-Mo system has been calculated via a truncated cluster expansion, through the combination of Full-Potential-Linear augmented Plane Wave (FP-LAPW) electronic structure calculations and of Cluster Variation Method (CVM) thermodynamic calculations in the irregular tetrahedron approximation. Four isothermal sections at 1750 K, 2000 K, 2250 K and 2500 K are calculated and correlated with recently published experimental data on the system. The results confirm that the critical temperature for the order-disorder equilibrium between Fe(3)Al-D0(3) and FeAl-B2 is increased by Mo additions, while the critical temperature for the FeAl-B2/A2 equilibrium is kept approximately invariant with increasing Mo contents. The stabilization of the Al-rich A2 phase in equilibrium with overstoichiometric B2-(Fe,Mo)Al is also consistent with the attribution of the A2 structure to the tau(2) phase, stable at high temperatures in overstoichiometric B2-FeAl. (C) 2009 Elsevier Ltd. All rights reserved.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

Modelling water adsorption on Au(210) surfaces: II. Monte Carlo simulations

Canonical Monte Carlo simulations for the Au(210)/H(2)O interface, using a force field recently proposed by us, are reported. The results exhibit the main features normally observed in simulations of water molecules in contact with different noble metal surfaces. The calculations also assess the influence of the surface topography on the structural aspects of the adsorbed water and on the distribution of the water molecules in the direction normal to the metal surface plane. The adsorption process is preferential at sites in the first layer of the metal. The analysis of the density profiles and dipole moment distributions points to two predominant orientations. Most of the molecules are adsorbed with the molecular plane parallel to surface, while others adsorb with one of the O-H bonds parallel to the surface and the other bond pointing towards the bulk liquid phase. There is also evidence of hydrogen bond formation between the first and second solvent layers at the interface. (c) 2007 Elsevier B.V. All rights reserved.

Face Detection and Facial Feature Localization for multi-pose faces and complex backgroundimages

The objective of this thesis work, is to propose an algorithm to detect the faces in a digital image with complex background. A lot of work has already been done in the area of face detection, but drawback of some face detection algorithms is the lack of ability to detect faces with closed eyes and open mouth. Thus facial features form an important basis for detection. The current thesis work focuses on detection of faces based on facial objects. The procedure is composed of three different phases: segmentation phase, filtering phase and localization phase. In segmentation phase, the algorithm utilizes color segmentation to isolate human skin color based on its chrominance properties. In filtering phase, Minkowski addition based object removal (Morphological operations) has been used to remove the non-skin regions. In the last phase, Image Processing and Computer Vision methods have been used to find the existence of facial components in the skin regions.This method is effective on detecting a face region with closed eyes, open mouth and a half profile face. The experiment’s results demonstrated that the detection accuracy is around 85.4% and the detection speed is faster when compared to neural network method and other techniques.

TCA Evaluation : Lab Measurements, Modelling and System Simulations

The study reported here is part of a large project for evaluation of the Thermo-Chemical Accumulator (TCA), a technology under development by the Swedish company ClimateWell AB. The studies concentrate on the use of the technology for comfort cooling. This report concentrates on measurements in the laboratory, modelling and system simulation. The TCA is a three-phase absorption heat pump that stores energy in the form of crystallised salt, in this case Lithium Chloride (LiCl) with water being the other substance. The process requires vacuum conditions as with standard absorption chillers using LiBr/water. Measurements were carried out in the laboratories at the Solar Energy Research Center SERC, at Högskolan Dalarna as well as at ClimateWell AB. The measurements at SERC were performed on a prototype version 7:1 and showed that this prototype had several problems resulting in poor and unreliable performance. The main results were that: there was significant corrosion leading to non-condensable gases that in turn caused very poor performance; unwanted crystallisation caused blockages as well as inconsistent behaviour; poor wetting of the heat exchangers resulted in relatively high temperature drops there. A measured thermal COP for cooling of 0.46 was found, which is significantly lower than the theoretical value. These findings resulted in a thorough redesign for the new prototype, called ClimateWell 10 (CW10), which was tested briefly by the authors at ClimateWell. The data collected here was not large, but enough to show that the machine worked consistently with no noticeable vacuum problems. It was also sufficient for identifying the main parameters in a simulation model developed for the TRNSYS simulation environment, but not enough to verify the model properly. This model was shown to be able to simulate the dynamic as well as static performance of the CW10, and was then used in a series of system simulations. A single system model was developed as the basis of the system simulations, consisting of a CW10 machine, 30 m2 flat plate solar collectors with backup boiler and an office with a design cooling load in Stockholm of 50 W/m2, resulting in a 7.5 kW design load for the 150 m2 floor area. Two base cases were defined based on this: one for Stockholm using a dry cooler with design cooling rate of 30 kW; one for Madrid with a cooling tower with design cooling rate of 34 kW. A number of parametric studies were performed based on these two base cases. These showed that the temperature lift is a limiting factor for cooling for higher ambient temperatures and for charging with fixed temperature source such as district heating. The simulated evacuated tube collector performs only marginally better than a good flat plate collector if considering the gross area, the margin being greater for larger solar fractions. For 30 m2 collector a solar faction of 49% and 67% were achieved for the Stockholm and Madrid base cases respectively. The average annual efficiency of the collector in Stockholm (12%) was much lower than that in Madrid (19%). The thermal COP was simulated to be approximately 0.70, but has not been possible to verify with measured data. The annual electrical COP was shown to be very dependent on the cooling load as a large proportion of electrical use is for components that are permanently on. For the cooling loads studied, the annual electrical COP ranged from 2.2 for a 2000 kWh cooling load to 18.0 for a 21000 kWh cooling load. There is however a potential to reduce the electricity consumption in the machine, which would improve these figures significantly. It was shown that a cooling tower is necessary for the Madrid climate, whereas a dry cooler is sufficient for Stockholm although a cooling tower does improve performance. The simulation study was very shallow and has shown a number of areas that are important to study in more depth. One such area is advanced control strategy, which is necessary to mitigate the weakness of the technology (low temperature lift for cooling) and to optimally use its strength (storage).

Energy Efficiency through Thermal Energy Storage - Evaluation of the Possibilities for the Swedish Building Stock, Phase 1

As a first step in assessing the potential of thermal energy storage in Swedish buildings, the current situation of the Swedish building stock and different storage methods are discussed in this paper. Overall, many buildings are from the 1960’s or earlier having a relatively high energy demand, creating opportunities for large energy savings. The major means of heating are electricity for detached houses and district heating for multi dwelling houses and premises. Cooling needs are relatively low but steadily increasing, emphasizing the need to consider energy storage for both heat and cold. The thermal mass of a building is important for passive storage of thermal energy but this has not been considered much when constructing buildings in Sweden. Instead, common ways of storing thermal energy in Swedish buildings today is in water storage tanks or in the ground using boreholes, while latent thermal energy storage is still very uncommon.

Multi-agent simulation of sawmill yard operations

This paper reports the findings of using multi-agent based simulation model to evaluate the sawmill yard operations within a large privately owned sawmill in Sweden, Bergkvist Insjön AB in the current case. Conventional working routines within sawmill yard threaten the overall efficiency and thereby limit the profit margin of sawmill. Deploying dynamic work routines within the sawmill yard is not readily feasible in real time, so discrete event simulation model has been investigated to be able to report optimal work order depending on the situations. Preliminary investigations indicate that the results achieved by simulation model are promising. It is expected that the results achieved in the current case will support Bergkvist-Insjön AB in making optimal decisions by deploying efficient work order in sawmill yard.

Multivariate modelling of landuse on in-stream salinity over multiple spatial scales

An impediment to sustainable dryland salinity management is the lack of information on contributing factors. GIS and satellite imagery now offer a cost-effective means of generating relevant land and water resource information for integrated regional management of salinity. In this paper the relationships between patterns in land uselcover distribution and base flow salt concentration in streams (indicated by EC) are investigated and modelled. The Glenelg-Hopkins area is a large regional watershed in southwest Victoria, Australia, covering approximately 2.6 million ha. It is currently estimated that 27,400 ha of land is affected by dryland salinity and this is predicted to rapidly increase in the next decade' if current conditions prevail. Salt concentration data from five gauging stations were analysed with multi-temporal land use maps obtained from satellite imagery. Multiple regression analyses demonstrated that the variables Native Vegetation and Dry/and Grain Cropping were the most significant influences on in~stream salinity in the whole catchment (1=88.9%) and 500 m V=88.3%) and 100 m riparian buffers (1=86.9%) during times of base flow. The implications for future land use planning, effectiveness of riparian zones and revegetation programmes is discussed. This work also demonstrates the utility of applying nmltivariate statistical analyses, spatial statistics, and remote sensing with data integrated in a GIS framework for the purpose of predicting and managing the regional salinity threat.

Landscape characteristics associated with species richness and occurrence of small native mammals inhabiting a coastal heathland: a spatial modelling approach

In the coastal region of south-western Victoria, Australia, populations of native small mammal species are restricted to patches of suitable habitat in a highly fragmented landscape. The size and spatial arrangement of these patches is likely to influence both the occupancy and richness of species at a location. Geographic Information System (GIS)-based habitat models of the species richness of native small mammals, and individual species  occurrences, were developed to produce maps displaying the spatial  configuration of suitable habitat. Models were generated using either generalised linear Poisson regression (for species richness) or logistic regression (for species occurrences) with species richness or  presence/absence as the dependent variable and landscape variables, extracted from both GIS data layers and multi-spectral digital imagery, as the predictor variables. A multi-model inference approach based on the Akaike Information Criterion was used and the resulting model was applied in a GIS framework to extrapolate predicted richness/likelihood of occurrence across the entire area of the study. A negative association between species  richness and elevation, habitat complexity and sun index indicated that richness within the study area decreases with increasing altitude, vertical vegetation structure and exposure to solar radiation. Landform  characteristics were important (to varying degrees) in determining habitat occupancy for all of the species examined, while the influence of habitat complexity was important for only one of the species. Performance of all but one of the models generated using presence/absence data was high, as indicated by the area under the curve of a receiver-operating characteristic plot. The effective conservation of the small mammal species in the area of concern is likely to depend on management actions that promote the protection of the critical habitats identified in the models.

Multi-temporal land use mapping using remotely sensed techniques and the integration of a pollutant load model in a GIS

Estimation of nutrient load production based on multi-temporal remotely sensed land-use data for the Glenelg-Hopkins region in southwest Victoria, Australia, is discussed. Changes in land use were mapped using archived Landsat data and computerized classification techniques. Land-use
change is unparalleled in recent history, with 13% of the region transformed in the last decade. Total nitrogen and phosphorus loading were estimated using an export coefficient model. The analysis demonstrated a disturbing increase in nitrogen and phosphorus loadings from 1995 to 2002. Whilst such increases were suspected from past anecdotal and ad-hoc evidence, our modelling quantitatively estimated such increases and thus demonstrated the enormous potential of using remote sensing and GIS for monitoring land-use change and hence improve land-use management.

Modelling habitat suitability of the swamp antechinus (Antechinus minimus maritimus) in the coastal heathlands of southern Victoria, Australia

In recent years, predictive habitat distribution models, derived by combining multivariate statistical analyses with Geographic Information System (GIS) technology, have been recognised for their utility in conservation planning. The size and spatial arrangement of suitable habitat can influence the long-term persistence of some faunal species. In southwestern Victoria, Australia, populations of the rare swamp antechinus (Antechinus minimus maritimus) are threatened by further fragmentation of suitable habitat. In the current study, a spatially explicit habitat suitability model was developed for A. minimus that incorporated a measure of vegetation structure. Models were generated using logistic regression with species presence or absence as the dependent variable and landscape variables, extracted from both GIS data layers and multi-spectral digital imagery, as the predictors. The most parsimonious model, based on the Akaike Information Criterion, was spatially extrapolated in the GIS. Probability of species presence was used as an index of habitat suitability. A negative association between A. minimus presence and both elevation and habitat complexity was evidenced, suggesting a preference for relatively low altitudes and a vegetation structure of low vertical complexity. The predictive performance of the selected model was shown to be high (91%), indicating a good fit of the model to the data. The proportion of the study area predicted as suitable habitat for A. minimus (Probability of occurrence greater-or-equal, slanted0.5) was 11.7%. Habitat suitability maps not only provide baseline information about the spatial arrangement of potentially suitable habitat for a species, but they also help to refine the search for other populations, making them an important conservation tool.

Modelling the Provenance of Data in Autonomous Systems

Determining the provenance of data, i.e. the process that led to that data, is vital in many disciplines. For example, in science, the process that produced a given result must be demonstrably rigorous for the result to be deemed reliable. A provenance system supports applications in recording adequate documentation about process executions to answer queries regarding provenance, and provides functionality to perform those queries. Several provenance systems are being developed, but all focus on systems in which the components are textitreactive, for example Web Services that act on the basis of a request, job submission system, etc. This limitation means that questions regarding the motives of autonomous actors, or textitagents, in such systems remain unanswerable in the general case. Such questions include: who was ultimately responsible for a given effect, what was their reason for initiating the process and does the effect of a process match what was intended to occur by those initiating the process? In this paper, we address this limitation by integrating two solutions: a generic, re-usable framework for representing the provenance of data in service-oriented architectures and a model for describing the goal-oriented delegation and engagement of agents in multi-agent systems. Using these solutions, we present algorithms to answer common questions regarding responsibility and success of a process and evaluate the approach with a simulated healthcare example.