Numerical solution of optimal control problems with constant control delays

We investigate a class of optimal control problems that exhibit constant exogenously given delays in the control in the equation of motion of the differential states. Therefore, we formulate an exemplary optimal control problem with one stock and one control variable and review some analytic properties of an optimal solution. However, analytical considerations are quite limited in case of delayed optimal control problems. In order to overcome these limits, we reformulate the problem and apply direct numerical methods to calculate approximate solutions that give a better understanding of this class of optimization problems. In particular, we present two possibilities to reformulate the delayed optimal control problem into an instantaneous optimal control problem and show how these can be solved numerically with a stateof- the-art direct method by applying Bock’s direct multiple shooting algorithm. We further demonstrate the strength of our approach by two economic examples.

Distributed sensor failure detection in sensor networks

We investigate the problem of distributed sensors' failure detection in networks with a small number of defective sensors, whose measurements differ significantly from the neighbor measurements. We build on the sparse nature of the binary sensor failure signals to propose a novel distributed detection algorithm based on gossip mechanisms and on Group Testing (GT), where the latter has been used so far in centralized detection problems. The new distributed GT algorithm estimates the set of scattered defective sensors with a low complexity distance decoder from a small number of linearly independent binary messages exchanged by the sensors. We first consider networks with one defective sensor and determine the minimal number of linearly independent messages needed for its detection with high probability. We then extend our study to the multiple defective sensors detection by modifying appropriately the message exchange protocol and the decoding procedure. We show that, for small and medium sized networks, the number of messages required for successful detection is actually smaller than the minimal number computed theoretically. Finally, simulations demonstrate that the proposed method outperforms methods based on random walks in terms of both detection performance and convergence rate.

Comparative Evaluation of Pelvic Allograft Selection Methods

This paper presents a firsthand comparative evaluation of three different existing methods for selecting a suitable allograft from a bone storage bank. The three examined methods are manual selection, automatic volume-based registration, and automatic surface-based registration. Although the methods were originally published for different bones, they were adapted to be systematically applied on the same data set of hemi-pelvises. A thorough experiment was designed and applied in order to highlight the advantages and disadvantages of each method. The methods were applied on the whole pelvis and on smaller fragments, thus producing a realistic set of clinical scenarios. Clinically relevant criteria are used for the assessment such as surface distances and the quality of the junctions between the donor and the receptor. The obtained results showed that both automatic methods outperform the manual counterpart. Additional advantages of the surface-based method are in the lower computational time requirements and the greater contact surfaces where the donor meets the recipient.

Performance of a new, 3D-monitor based random-dot stereotest for children under 4 years of age

BACKGROUND: The aim of this study was to determine the performance of a new, 3D-monitor based, objective stereotest in children under the age of four. METHODS: Random-dot circles (diameter 10 cm, crossed, disparity of 0.34 degrees) randomly changing their position were presented on an 3D-monitor while eye movements were monitored by infrared photo-oculography. If > or = 3 consecutive stimuli were seen, a positive response was assumed. One hundred thirty-four normal children aged 2 months to 4 years (average 17+/-15.3 months) were examined. RESULTS: Below the age of 12 months, we were not able to obtain a response to the 3D stimulus. For older children the following rates of positive responses were found: 12-18 months 25%, 18-24 months 10%, 24-30 months 16%, 30-36 months 57%, 36-42 months 100%, and 42-48 months 91%. Multiple linear logistic regression showed a significant influence on stimulus recognition of the explanatory variables age (p<0.00001) and child cooperation (p<0.001), but not of gender (p>0.1). CONCLUSIONS: This 3D-monitor based stereotest allows an objective measurement of random-dot stereopsis in younger children. It might open new ways to screen children for visual abnormalities and to study the development of stereovision. However, the current experimental setting does not allow determining random-dot stereopsis in children younger than 12 months.

Statistical Analysis of the heavy neutral Atoms measured by IBEX

We investigate the directional distribution of heavy neutral atoms in the heliosphere by using heavy neutral maps generated with the IBEX-Lo instrument over three years from 2009 to 2011. The interstellar neutral (ISN) O&Ne gas flow was found in the first-year heavy neutral map at 601 keV and its flow direction and temperature were studied. However, due to the low counting statistics, researchers have not treated the full sky maps in detail. The main goal of this study is to evaluate the statistical significance of each pixel in the heavy neutral maps to get a better understanding of the directional distribution of heavy neutral atoms in the heliosphere. Here, we examine three statistical analysis methods: the signal-to-noise filter, the confidence limit method, and the cluster analysis method. These methods allow us to exclude background from areas where the heavy neutral signal is statistically significant. These methods also allow the consistent detection of heavy neutral atom structures. The main emission feature expands toward lower longitude and higher latitude from the observational peak of the ISN O&Ne gas flow. We call this emission the extended tail. It may be an imprint of the secondary oxygen atoms generated by charge exchange between ISN hydrogen atoms and oxygen ions in the outer heliosheath.

Polygonal cracks in the seasonal semi-translucent CO2 ice layer in Martian polar areas

In this paper, we use morphological and numerical methods to test the hypothesis that seasonally formed fracture patterns in the Martian polar regions result from the brittle failure of seasonal CO2 slab ice. The observations by the High Resolution Imaging Science Experiment (HiRISE) of polar regions of Mars show very narrow dark elongated linear patterns that are observed during some periods of time in spring, disappear in summer and re-appear again in the following spring. They are repeatedly formed in the same areas but they do not repeat the exact pattern from year to year. This leads to the conclusion that they are cracks formed in the seasonal ice layer. Some of models of seasonal surface processes rely on the existence of a transparent form of CO2 ice, so-called slab ice. For the creation of the observed cracks the ice is required to be a continuous media, not an agglomeration of relatively separate particles like a firn. The best explanation for our observations is a slab ice with relatively high transparency in the visible wavelength range. This transparency allows a solid state green-house effect to act underneath the ice sheet raising the pressure by sublimation from below. The trapped gas creates overpressure and the ice sheet breaks at some point creating the observed cracks. We show that the times when the cracks appear are in agreement with the model calculation, providing one more piece of evidence that CO2 slab ice covers polar areas in spring.

Focal laminar cortical infarcts following aneurysmal subarachnoid haemorrhage

INTRODUCTION: The aim of this prospective study was to analyse small band-like cortical infarcts after subarachnoid haemorrhage (SAH) using magnetic resonance imaging (MRI) with reference to additional digital subtraction angiography (DSA). METHODS: In a 5-year period between January 2002 and January 2007 10 out of 188 patients with aneurysmal SAH were evaluated (one patient Hunt and Hess grade I, one patient grade II, four patients grade III, two patients grade IV, and two patients grade V). The imaging protocol included serially performed MRI with diffusion- and perfusion-weighted images (DWI/PWI) at three time points after aneurysm treatment, and cerebral vasospasm (CVS) was analysed on follow-up DSA on day 7+/-3 after SAH. RESULTS: The lesions were located in the frontal lobe (n=10), in the insular cortex (n=3) and in the parietal lobe (n=1). The band-like infarcts occurred after a mean time interval of 5.8 days (range 3-10 days) and showed unexceptional adjacent thick sulcal clots. Seven out of ten patients with cortical infarcts had no or mild CVS, and in the remaining three patients DSA disclosed moderate (n=2) or severe (n=1) CVS. CONCLUSION: The infarct pattern after aneurysmal SAH includes cortical band-like lesions. In contrast to territorial infarcts or lacunar infarcts in the white matter which develop as a result of moderate or severe proximal and/or distal vasospasm visible on angiography, the cortical band-like lesions adjacent to sulcal clots may also develop without evidence of macroscopic vasospasm, implying a vasospastic reaction of the most distal superficial and intraparenchymal vessels.

Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Modeling asteroid collisions and impact processes

As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.