Safety-critical medical device development using the UPP2SF model translation tool

Software-based control of life-critical embedded systems has become increasingly complex, and to a large extent has come to determine the safety of the human being. For example, implantable cardiac pacemakers have over 80,000 lines of code which are responsible for maintaining the heart within safe operating limits. As firmware-related recalls accounted for over 41% of the 600,000 devices recalled in the last decade, there is a need for rigorous model-driven design tools to generate verified code from verified software models. To this effect, we have developed the UPP2SF model-translation tool, which facilitates automatic conversion of verified models (in UPPAAL) to models that may be simulated and tested (in Simulink/Stateflow). We describe the translation rules that ensure correct model conversion, applicable to a large class of models. We demonstrate how UPP2SF is used in themodel-driven design of a pacemaker whosemodel is (a) designed and verified in UPPAAL (using timed automata), (b) automatically translated to Stateflow for simulation-based testing, and then (c) automatically generated into modular code for hardware-level integration testing of timing-related errors. In addition, we show how UPP2SF may be used for worst-case execution time estimation early in the design stage. Using UPP2SF, we demonstrate the value of integrated end-to-end modeling, verification, code-generation and testing process for complex software-controlled embedded systems. © 2014 ACM.

Aplicación de cadenas absorbentes de Markov a un modelo de autofecundación en alohexaploides

p.21-27

A case based model to assist in the design process in the manufacture of furniture products

This paper describes progress on a project to utilise case based reasoning methods in the design and manufacture of furniture products. The novel feature of this research is that cases are represented as structures in a relational database of products, components and materials. The paper proposes a method for extending the usual "weighted sum" over attribute similarities for a ·single table to encompass relational structures over several tables. The capabilities of the system are discussed, particularly with respect to differing user objectives, such as cost estimation, CAD, cutting scheme re-use, and initial design. It is shown that specification of a target case as a relational structure combined with suitable weights can fulfil several user functions. However, it is also shown that some user functions cannot satisfactorily be specified via a single target case. For these functions it is proposed to allow the specification of a set of target cases. A derived similarity measure between individuals and sets of cases is proposed.

Generalized Markov branching processes

A generalized Markov Brnching Process (GMBP) is a Markov branching model where the infinitesimal branching rates are modified with an interaction index. It is proved that there always exists only one GMBP. An associated differential-integral equation is derived. The extinction probalility and the mean and conditional mean extinction times are obtained. Ergodicity and stability of GMBP with resurrection are also considered. Easy checking criteria are established for ordinary and strong ergodicty. The equilibrium distribution is given in an elegant closed form. The probability meaning of our results is clear and thus explained.

Human respiration rate estimation using ultra-wideband distributed cognitive radar system

It has been shown that remote monitoring of pulmonary activity can be achieved using ultra-wideband (UWB) systems, which shows promise in home healthcare, rescue, and security applications. In this paper, we first present a multi-ray propagation model for UWB signal, which is traveling through the human thorax and is reflected on the air/dry-skin/fat/muscle interfaces. A geometry-based statistical channel model is then developed for simulating the reception of UWB signals in the indoor propagation environment. This model enables replication of time-varying multipath profiles due to the displacement of a human chest. Subsequently, a UWB distributed cognitive radar system (UWB-DCRS) is developed for the robust detection of chest cavity motion and the accurate estimation of respiration rate. The analytical framework can serve as a basis in the planning and evaluation of future measurement programs. We also provide a case study on how the antenna beamwidth affects the estimation of respiration rate based on the proposed propagation models and system architecture

Communication system model for information rate evaluation of differential detection over time-varying channels

A communication system model for mutual information performance analysis of multiple-symbol differential M-phase shift keying over time-correlated, time-varying flat-fading communication channels is developed. This model is a finite-state Markov (FSM) equivalent channel representing the cascade of the differential encoder, FSM channel model and differential decoder. A state-space approach is used to model channel phase time correlations. The equivalent model falls in a class that facilitates the use of the forward backward algorithm, enabling the important information theoretic results to be evaluated. Using such a model, one is able to calculate mutual information for differential detection over time-varying fading channels with an essentially finite time set of correlations, including the Clarke fading channel. Using the equivalent channel, it is proved and corroborated by simulations that multiple-symbol differential detection preserves the channel information capacity when the observation interval approaches infinity.

Sun glint estimation in marine satellite images: a comparison of results from calculation and radiative transfer modeling

The intensity and location of Sun glint in two Medium Resolution Imaging Spectrometer (MERIS) images was modeled using a radiative transfer model that includes elevation features as well as the slope of the sea surface. The results are compared to estimates made using glint flagging and correction approaches used within standard atmospheric correction processing code. The model estimate gives a glint pattern with a similar width but lower peak level than any current method, or than that estimated by a radiative transfer model with surfaces that include slope but not height. The MERIS third reprocessing recently adopted a new slope statistics model for Sun glint correction; the results show that this model is an outlier with respect to both the elevation model and other slope statistics models and we recommend that its adoption should be reviewed.

Mechanisms of long-term decline in size of lesser sandeels in the North Sea explored using a growth and phenology model

The lesser sandeel Ammodytes marinus is a key species in the North Sea ecosystem, transferring energy from planktonic producers to top predators. Previous studies have shown a long-term decline in the size of 0-group sandeels in the western North Sea, but they were unable to pinpoint the mechanism (later hatching, slower growth or changes in size-dependent mortality) or cause. To investigate the first 2 possibilities we combined 2 independent time series of sandeel size, namely data from chick-feeding Atlantic puffins Fratercula arctica and from the Continuous Plankton Recorder (CPR), in a novel statistical model implemented using Markov Chain Monte Carlo (MCMC). The model estimated annual mean length on 1 July, as well as hatching date and growth rate for sandeels from 1973 to 2006. Mean length-at-date declined by 22% over this period, corresponding to a 60% decrease in energy content, with a sharper decline since 2002. Up to the mid-1990s, the decline was associated with a trend towards later hatching. Subsequently, hatching became earlier again, and the continued trend towards smaller size appears to have been driven by lower growth rates, particularly in the most recent years, although we could not rule out changes in size-dependent mortality. Our findings point to major changes in key aspects of sandeel life history, which we consider are most likely due to direct and indirect temperature-related changes over a range of biotic factors, including the seasonal distribution of copepods and intra- and inter-specific competition with planktivorous fish. The results have implications both for the many predators of sandeels and for age and size of maturation in this aggregation of North Sea sandeels.

Air-sea exchange of methanol and acetone during HiWinGS: Estimation of air phase, water phase gas transfer velocities

The air-sea fluxes of methanol and acetone were measured concurrently using a proton-transfer-reaction mass spectrometer (PTR-MS) with the eddy covariance (EC) technique during the High Wind Gas Exchange Study (HiWinGS) in 2013. The seawater concentrations of these compounds were also measured twice daily with the same PTR-MS coupled to a membrane inlet. Dissolved concentrations near the surface ranged from 7 to 28 nM for methanol and from 3 to 9 nM for acetone. Both gases were consistently transported from the atmosphere to the ocean as a result of their low sea surface saturations. The largest influxes were observed in regions of high atmospheric concentrations and strong winds (up to 25 m s(-1)). Comparison of the total air-sea transfer velocity of these two gases (K-a), along with the in situ sensible heat transfer rate, allows us to constrain the individual gas transfer velocity in the air phase (k(a)) and water phase (k(w)). Among existing parameterizations, the scaling of k(a) from the COARE model is the most consistent with our observations. The k(w) we estimated is comparable to the tangential (shear driven) transfer velocity previously determined from measurements of dimethyl sulfide. Lastly, we estimate the wet deposition of methanol and acetone in our study region and evaluate the lifetimes of these compounds in the surface ocean and lower atmosphere with respect to total (dry plus wet) atmospheric deposition.

Judging Time Intevals using a model of perceptuo-motor control

Estimating a time interval and temporally coordinating movements in space are fundamental skills, but the relationships between these different forms of timing, and the neural processes that they incur, are not well understood. While different theories have been proposed to account for time perception, time estimation, and the temporal patterns of coordination, there are no general mechanisms which unify these various timing skills. This study considers whether a model of perceptuo-motor timing, the tau(GUIDE), can also describe how certain judgements of elapsed time are made. To evaluate this, an equation for determining interval estimates was derived from the tau(GUIDE) model and tested in a task where participants had to throw a ball and estimate when it would hit the floor. The results showed that in accordance with the model, very accurate judgements could be made without vision (mean timing error -19.24 msec), and the model was a good predictor of skilled participants' estimate timing. It was concluded that since the tau(GUIDE) principle provides temporal information in a generic form, it could be a unitary process that links different forms of timing.

Stochastic protein folding simulation in the three-dimensional HP-model

We present results from three-dimensional protein folding simulations in the HP-model on ten benchmark problems. The simulations are executed by a simulated annealing-based algorithm with a time-dependent cooling schedule. The neighbourhood relation is determined by the pull-move set. The results provide experimental evidence that the maximum depth D of local minima of the underlying energy landscape can be upper bounded by D < n(2/3). The local search procedure employs the stopping criterion (In/delta)(D/gamma) where m is an estimation of the average number of neighbouring conformations, gamma relates to the mean of non-zero differences of the objective function for neighbouring conformations, and 1-delta is the confidence that a minimum conformation has been found. The bound complies with the results obtained for the ten benchmark problems. (c) 2008 Elsevier Ltd. All rights reserved.

Plant-wide predictive control for a thermal power plant based on a physical plant model

A constrained non-linear, physical model-based, predictive control (NPMPC) strategy is developed for improved plant-wide control of a thermal power plant. The strategy makes use of successive linearisation and recursive state estimation using extended Kalman filtering to obtain a linear state-space model. The linear model and a quadratic programming routine are used to design a constrained long-range predictive controller One special feature is the careful selection of a specific set of plant model parameters for online estimation, to account for time-varying system characteristics resulting from major system disturbances and ageing. These parameters act as nonstationary stochastic states and help to provide sufficient degrees-of-freedom to obtain unbiased estimates of controlled outputs. A 14th order non-linear plant model, simulating the dominant characteristics of a 200 MW oil-fired pou er plant has been used to test the NPMPC algorithm. The control strategy gives impressive simulation results, during large system disturbances and extremely high rate of load changes, right across the operating range. These results compare favourably to those obtained with the state-space GPC method designed under similar conditions.

Contribution of 30 Biomarkers to 10-year cardiovascular risk estimation in 2 population cohorts: The MONICA, risk, genetics, archiving and monograph (MORGAM) biomarker project

Background— Cardiovascular risk estimation by novel biomarkers needs assessment in disease-free population cohorts, followed up for incident cardiovascular events, assaying the serum and plasma archived at baseline. We report results from 2 cohorts in such a continuing study.
Methods and Results— Thirty novel biomarkers from different pathophysiological pathways were evaluated in 7915 men and women of the FINRISK97 population cohort with 538 incident cardiovascular events at 10 years (fatal or nonfatal coronary or stroke events), from which a biomarker score was developed and then validated in the 2551 men of the Belfast Prospective Epidemiological Study of Myocardial Infarction (PRIME) cohort (260 events). No single biomarker consistently improved risk estimation in FINRISK97 men and FINRISK97 women and the Belfast PRIME Men cohort after allowing for confounding factors; however, the strongest associations (with hazard ratio per SD in FINRISK97 men) were found for N-terminal pro-brain natriuretic peptide (1.23), C-reactive protein (1.23), B-type natriuretic peptide (1.19), and sensitive troponin I (1.18). A biomarker score was developed from the FINRISK97 cohort with the use of regression coefficients and lasso methods, with selection of troponin I, C-reactive protein, and N-terminal pro-brain natriuretic peptide. Adding this score to a conventional risk factor model in the Belfast PRIME Men cohort validated it by improved c-statistics (P=0.004) and integrated discrimination (P<0.0001) and led to significant reclassification of individuals into risk categories (P=0.0008).
Conclusions— The addition of a biomarker score including N-terminal pro-brain natriuretic peptide, C-reactive protein, and sensitive troponin I to a conventional risk model improved 10-year risk estimation for cardiovascular events in 2 middle-aged European populations. Further validation is needed in other populations and age groups.

Estimation of in situ stresses using anisotropic elasticity and suction measurements

Sampling and specimen preparation produce changes in mean effective stresses and pore water pressures, even with ‘perfect sampling’. The paper takes an existing simplified three-parameter cross-anisotropic elastic model and uses it to model these changes. The required ratio of cross-anisotropic parameters J/3G* can be obtained from standard CIU triaxial tests. If measurements are also made of suctions in unloaded specimens in the laboratory, then a combination of J/3G*, the measured suction, and the effective overburden pressure permits an estimation of the horizontal effective pressure and the K 0 ‘at rest’ coefficient. This can be helpful in numerical modelling that needs to start from in situ conditions, and in planning pressure levels for reconsolidation of clay specimens in the laboratory. Tests were done on Belfast Upper Boulder Clay from a depth of 28 m. Values of horizontal in situ effective stress estimated from these measurements compare favorably with conventional estimates of the ‘at rest’ coefficient K 0 and the overconsolidation ratio. Estimates of horizontal stress in London Clay were made using published data and the results compared with actual measurements. Again reasonable agreement was obtained.