A Nonstationary Space-Time Gaussian Process Model for Partially Converged Simulations

In the context of expensive numerical experiments, a promising solution for alleviating the computational costs consists of using partially converged simulations instead of exact solutions. The gain in computational time is at the price of precision in the response. This work addresses the issue of fitting a Gaussian process model to partially converged simulation data for further use in prediction. The main challenge consists of the adequate approximation of the error due to partial convergence, which is correlated in both design variables and time directions. Here, we propose fitting a Gaussian process in the joint space of design parameters and computational time. The model is constructed by building a nonstationary covariance kernel that reflects accurately the actual structure of the error. Practical solutions are proposed for solving parameter estimation issues associated with the proposed model. The method is applied to a computational fluid dynamics test case and shows significant improvement in prediction compared to a classical kriging model.

An Optical Survey for Space Debris on Highly Eccentric MEO Orbits

Optical surveys for space debris in high-altitude orbits have been conducted since more than ten years. Originally these efforts concentrated mainly on the geostationary ring (GEO). Corresponding observation strategies, processing techniques and cataloguing approaches have been developed and successfully applied. The ESA GEO surveys, e.g., resulted in the detection of a significant population of small-size debris and later in the discovery of high area-to-mass ratio objects in GEO-like orbits. The observation scenarios were successively adapted to survey the geostationary transfer orbit (GTO) region; and recently surveys to search for debris in the medium Earth orbit (MEO) region of the global navigation satellite constellations were successfully conducted. Comparably less experience (both, in terms of practical observation and strategy definition) is available for eccentric orbits that (at least partly) are in the MEO region, in particular for the Molniya-type orbits. Several breakup events and deliberate fragmentations are known to have taken place in such orbits. Different survey and follow-up strategies for searching space debris objects in highly-eccentric MEO orbits, and to acquire orbits which are sufficiently accurate to catalogue such objects and to maintain their orbits over longer time spans were developed. Simulations were performed to compare the performance of different survey and cataloguing strategies. Eventually, optical observations were conducted in the framework of an ESA study using ESA’s Space Debris Telescope (ESASDT) the 1-m Zeiss telescope located at the Optical Ground Station (OGS) at the Teide Observatory at Tenerife, Spain. A first series of surveys of Molnjya-type orbits was performed between January and April 2013. During these four months survey observations were performed during nine nights. A basic survey consisted of observing a single geocentric field for 10 minutes. If a faint object was found, follow-up observations were performed during the same night to ensure a save rediscovery of the object during the next nights. Additional follow-up observations to maintain the orbits of these newly discovered faint objects were also acquired with AIUB ́s 1m ZIMLAT telescope in Zimmerwald, Switzerland. Eventually 195 basic surveys were performed during these nine nights corresponding to about 32.5 hours of observations. In total 24 uncorrelated faint objects were discovered and all known catalogue objects in the survey fields were detected. On average one uncorrelated object was found every 80 minutes. Some of these objects show a considerable brightness variation and have a high area-to-mass ratio as determined in the orbit estimation process.

Status of the Zimmerwald SLR station

The Zimmerwald SLR station is operated in a monostatic mode with 532nm laser pulses emitted at adjustable frequencies of 90-110Hz with energies slightly less than 10mJ. A rotating shutter protects the CSPAD receiver from the backscatter of the transmit beam. These systems are located below the telescope in an operator room housed within the observatory building with the laser system located in a separated, air-conditioned part of the room. All hardware components may be automatically accessed by the control software and from remote if required. Thanks to the fully automatic and remotely controllable SLR operations, the Zimmerwald station is one of the most productive stations in the ILRS network. Key characteristics of the hardware are shown. Specialities like the tracking of the full GLONASS constellation, one-way ranging to the Lunar Reconnaissance Orbiter, and photon reception from bi-static experiments with the Graz SLR station are highlighted as well.

Temperature Prediction Model for Bone Drilling Based on Density Distribution and In Vivo Experiments for Minimally Invasive Robotic Cochlear Implantation

Surgical robots have been proposed ex vivo to drill precise holes in the temporal bone for minimally invasive cochlear implantation. The main risk of the procedure is damage of the facial nerve due to mechanical interaction or due to temperature elevation during the drilling process. To evaluate the thermal risk of the drilling process, a simplified model is proposed which aims to enable an assessment of risk posed to the facial nerve for a given set of constant process parameters for different mastoid bone densities. The model uses the bone density distribution along the drilling trajectory in the mastoid bone to calculate a time dependent heat production function at the tip of the drill bit. Using a time dependent moving point source Green's function, the heat equation can be solved at a certain point in space so that the resulting temperatures can be calculated over time. The model was calibrated and initially verified with in vivo temperature data. The data was collected in minimally invasive robotic drilling of 12 holes in four different sheep. The sheep were anesthetized and the temperature elevations were measured with a thermocouple which was inserted in a previously drilled hole next to the planned drilling trajectory. Bone density distributions were extracted from pre-operative CT data by averaging Hounsfield values over the drill bit diameter. Post-operative [Formula: see text]CT data was used to verify the drilling accuracy of the trajectories. The comparison of measured and calculated temperatures shows a very good match for both heating and cooling phases. The average prediction error of the maximum temperature was less than 0.7 °C and the average root mean square error was approximately 0.5 °C. To analyze potential thermal damage, the model was used to calculate temperature profiles and cumulative equivalent minutes at 43 °C at a minimal distance to the facial nerve. For the selected drilling parameters, temperature elevation profiles and cumulative equivalent minutes suggest that thermal elevation of this minimally invasive cochlear implantation surgery may pose a risk to the facial nerve, especially in sclerotic or high density mastoid bones. Optimized drilling parameters need to be evaluated and the model could be used for future risk evaluation.

Dissolution-Precipitation in Porous Media: Experiments and Modelling

The evolution of porosity due to dissolution/precipitation processes of minerals and the associated change of transport parameters are of major interest for natural geological environments and engineered underground structures. We designed a reproducible and fast to conduct 2D experiment, which is flexible enough to investigate several process couplings implemented in the numerical code OpenGeosys-GEM (OGS-GEM). We investigated advective-diffusive transport of solutes, effect of liquid phase density on advective transport, and kinetically controlled dissolution/precipitation reactions causing porosity changes. In addition, the system allowed to investigate the influence of microscopic (pore scale) processes on macroscopic (continuum scale) transport. A Plexiglas tank of dimension 10 × 10 cm was filled with a 1 cm thick reactive layer consisting of a bimodal grain size distribution of celestite (SrSO4) crystals, sandwiched between two layers of sand. A barium chloride solution was injected into the tank causing an asymmetric flow field to develop. As the barium chloride reached the celestite region, dissolution of celestite was initiated and barite precipitated. Due to the higher molar volume of barite, its precipitation caused a porosity decrease and thus also a decrease in the permeability of the porous medium. The change of flow in space and time was observed via injection of conservative tracers and analysis of effluents. In addition, an extensive post-mortem analysis of the reacted medium was conducted. We could successfully model the flow (with and without fluid density effects) and the transport of conservative tracers with a (continuum scale) reactive transport model. The prediction of the reactive experiments initially failed. Only the inclusion of information from post-mortem analysis gave a satisfactory match for the case where the flow field changed due to dissolution/precipitation reactions. We concentrated on the refinement of post-mortem analysis and the investigation of the dissolution/precipitation mechanisms at the pore scale. Our analytical techniques combined scanning electron microscopy (SEM) and synchrotron X-ray micro-diffraction/micro-fluorescence performed at the XAS beamline (Swiss Light Source). The newly formed phases include an epitaxial growth of barite micro-crystals on large celestite crystals (epitaxial growth) and a nano-crystalline barite phase (resulting from the dissolution of small celestite crystals) with residues of celestite crystals in the pore interstices. Classical nucleation theory, using well-established and estimated parameters describing barite precipitation, was applied to explain the mineralogical changes occurring in our system. Our pore scale investigation showed limits of the continuum scale reactive transport model. Although kinetic effects were implemented by fixing two distinct rates for the dissolution of large and small celestite crystals, instantaneous precipitation of barite was assumed as soon as oversaturation occurred. Precipitation kinetics, passivation of large celestite crystals and metastability of supersaturated solutions, i.e. the conditions under which nucleation cannot occur despite high supersaturation, were neglected. These results will be used to develop a pore scale model that describes precipitation and dissolution of crystals at the pore scale for various transport and chemical conditions. Pore scale modelling can be used to parameterize constitutive equations to introduce pore-scale corrections into macroscopic (continuum) reactive transport models. Microscopic understanding of the system is fundamental for modelling from the pore to the continuum scale.