Approximate Model Checking of Real-Time Systems for Linear Duration Invariants

Real-time systems are usually modelled with timed automata and real-time requirements relating to the state durations of the system are often specifiable using Linear Duration Invariants, which is a decidable subclass of Duration Calculus formulas. Various algorithms have been developed to check timed automata or real-time automata for linear duration invariants, but each needs complicated preprocessing and exponential calculation. To the best of our knowledge, these algorithms have not been implemented. In this paper, we present an approximate model checking technique based on a genetic algorithm to check real-time automata for linear durration invariants in reasonable times. Genetic algorithm is a good optimization method when a problem needs massive computation and it works particularly well in our case because the fitness function which is derived from the linear duration invariant is linear. ACM Computing Classification System (1998): D.2.4, C.3.

Practical stability with respect to model mismatch of approximate discrete-time output feedback control

This paper establishes a practical stability result for discrete-time output feedback control involving mismatch between the exact system to be stabilised and the approximating system used to design the controller. The practical stability is in the sense of an asymptotic bound on the amount of error bias introduced by the model approximation, and is established using local consistency properties of the systems. Importantly, the practical stability established here does not require the approximating system to be of the same model type as the exact system. Examples are presented to illustrate the nature of our practical stability result.

Fractal geometry model for effective thermal conductivity of three-phase porous media

An approximate model, a fractal geometry model, for the effective thermal conductivity of three-phase/unsaturated porous media is proposed based on the thermal-electrical analogy technique and on statistical self-similarity of porous media. The proposed thermal conductivity model is expressed as a function of porosity (related to stage n of Sierpinski carpet), ratio of areas, ratio of component thermal conductivities, and saturation. The recursive algorithm for the thermal conductivity by the proposed model is presented and found to be quite simple. The model predictions are compared with the existing measurements. Good agreement is found between the present model predictions and the existing experimental data. This verifies the validity of the proposed model. (C) 2004 American Institute of Physics.

Simulation model for driving dynamics, energy use and power supply

Dissertação para obtenção do grau de Mestre em Engenharia Electrotécnica Ramo de Energia

Consistent approximate models of the global atmosphere: shallow, deep, hydrostatic, quasi-hydrostatic and non-hydrostatic

We study global atmosphere models that are at least as accurate as the hydrostatic primitive equations (HPEs), reviewing known results and reporting some new ones. The HPEs make spherical geopotential and shallow atmosphere approximations in addition to the hydrostatic approximation. As is well known, a consistent application of the shallow atmosphere approximation requires omission of those Coriolis terms that vary as the cosine of latitude and of certain other terms in the components of the momentum equation. An approximate model is here regarded as consistent if it formally preserves conservation principles for axial angular momentum, energy and potential vorticity, and (following R. Müller) if its momentum component equations have Lagrange's form. Within these criteria, four consistent approximate global models, including the HPEs themselves, are identified in a height-coordinate framework. The four models, each of which includes the spherical geopotential approximation, correspond to whether the shallow atmosphere and hydrostatic (or quasi-hydrostatic) approximations are individually made or not made. Restrictions on representing the spatial variation of apparent gravity occur. Solution methods and the situation in a pressure-coordinate framework are discussed. © Crown copyright 2005.

Data Flow Analysis and the Linear Programming Model

* The research is supported partly by INTAS: 04-77-7173 project, http://www.intas.be

Dynamics of loop formation in a semiflexible polymer

The dynamics of loop formation by linear polymer chains has been a topic of several theoretical and experimental studies. Formation of loops and their opening are key processes in many important biological processes. Loop formation in flexible chains has been extensively studied by many groups. However, in the more realistic case of semiflexible polymers, not much results are available. In a recent study [K. P. Santo and K. L. Sebastian, Phys. Rev. E 73, 031923 (2006)], we investigated opening dynamics of semiflexible loops in the short chain limit and presented results for opening rates as a function of the length of the chain. We presented an approximate model for a semiflexible polymer in the rod limit based on a semiclassical expansion of the bending energy of the chain. The model provided an easy way to describe the dynamics. In this paper, using this model, we investigate the reverse process, i.e., the loop formation dynamics of a semiflexible polymer chain by describing the process as a diffusion-controlled reaction. We make use of the ``closure approximation'' of Wilemski and Fixman [G. Wilemski and M. Fixman, J. Chem. Phys. 60, 878 (1974)], in which a sink function is used to represent the reaction. We perform a detailed multidimensional analysis of the problem and calculate closing times for a semiflexible chain. We show that for short chains, the loop formation time tau decreases with the contour length of the polymer. But for longer chains, it increases with length obeying a power law and so it has a minimum at an intermediate length. In terms of dimensionless variables, the closing time is found to be given by tau similar to L-n exp(const/L), where n=4.5-6. The minimum loop formation time occurs at a length L-m of about 2.2-2.4. These are, indeed, the results that are physically expected, but a multidimensional analysis leading to these results does not seem to exist in the literature so far.

Genetic programming metamodel for rotating beams

This paper investigates the use of Genetic Programming (GP) to create an approximate model for the non-linear relationship between flexural stiffness, length, mass per unit length and rotation speed associated with rotating beams and their natural frequencies. GP, a relatively new form of artificial intelligence, is derived from the Darwinian concept of evolution and genetics and it creates computer programs to solve problems by manipulating their tree structures. GP predicts the size and structural complexity of the empirical model by minimizing the mean square error at the specified points of input-output relationship dataset. This dataset is generated using a finite element model. The validity of the GP-generated model is tested by comparing the natural frequencies at training and at additional input data points. It is found that by using a non-dimensional stiffness, it is possible to get simple and accurate function approximation for the natural frequency. This function approximation model is then used to study the relationships between natural frequency and various influencing parameters for uniform and tapered beams. The relations obtained with GP model agree well with FEM results and can be used for preliminary design and structural optimization studies.

Asymptotic Bounds on the Power-Delay Tradeoff for Fading Point-to-Point Links From Geometric Bounds on the Stationary Distribution of the Queue Length

The optimal power-delay tradeoff is studied for a time-slotted independently and identically distributed fading point-to-point link, with perfect channel state information at both transmitter and receiver, and with random packet arrivals to the transmitter queue. It is assumed that the transmitter can control the number of packets served by controlling the transmit power in the slot. The optimal tradeoff between average power and average delay is analyzed for stationary and monotone transmitter policies. For such policies, an asymptotic lower bound on the minimum average delay of the packets is obtained, when average transmitter power approaches the minimum average power required for transmitter queue stability. The asymptotic lower bound on the minimum average delay is obtained from geometric upper bounds on the stationary distribution of the queue length. This approach, which uses geometric upper bounds, also leads to an intuitive explanation of the asymptotic behavior of average delay. The asymptotic lower bounds, along with previously known asymptotic upper bounds, are used to identify three new cases where the order of the asymptotic behavior differs from that obtained from a previously considered approximate model, in which the transmit power is a strictly convex function of real valued service batch size for every fade state.

Effects of Noise and Detection Error in a Dynamic Adaptive Optics System

It is well known that noise and detection error can affect the performances of an adaptive optics (AO) system. Effects of noise and detection error on the phase compensation effectiveness in a dynamic AO system are investigated by means of a pure numerical simulation in this paper. A theoretical model for numerically simulating effects of noise and detection error in a static AO system and a corresponding computer program were presented in a previous article. A numerical simulation of effects of noise and detection error is combined with our previous numeral simulation of a dynamic AO system in this paper and a corresponding computer program has been compiled. Effects of detection error, readout noise and photon noise are included and investigated by a numerical simulation for finding the preferred working conditions and the best performances in a practical dynamic AO system. An approximate model is presented as well. Under many practical conditions such approximate model is a good alternative to the more accurate one. A simple algorithm which can be used for reducing the effect of noise is presented as well. When signal to noise ratio is very low, such method can be used to improve the performances of a dynamic AO system.

Optimal guidance of low-thrust interplanetary space vehicles

The low-thrust guidance problem is defined as the minimum terminal variance (MTV) control of a space vehicle subjected to random perturbations of its trajectory. To accomplish this control task, only bounded thrust level and thrust angle deviations are allowed, and these must be calculated based solely on the information gained from noisy, partial observations of the state. In order to establish the validity of various approximations, the problem is first investigated under the idealized conditions of perfect state information and negligible dynamic errors. To check each approximate model, an algorithm is developed to facilitate the computation of the open loop trajectories for the nonlinear bang-bang system. Using the results of this phase in conjunction with the Ornstein-Uhlenbeck process as a model for the random inputs to the system, the MTV guidance problem is reformulated as a stochastic, bang-bang, optimal control problem. Since a complete analytic solution seems to be unattainable, asymptotic solutions are developed by numerical methods. However, it is shown analytically that a Kalman filter in cascade with an appropriate nonlinear MTV controller is an optimal configuration. The resulting system is simulated using the Monte Carlo technique and is compared to other guidance schemes of current interest.

Nonlinear drillstring dynamics analysis

This paper studies the dynamical response of a rotary drilling system with a drag bit, using a lumped parameter model that takes into consideration the axial and torsional vibration modes of the bit. These vibrations are coupled through a bit-rock interaction law. At the bit-rock interface, the cutting process introduces a state-dependent delay, while the frictional process is responsible for discontinuous right-hand sides in the equations governing the motion of the bit. This complex system is characterized by a fast axial dynamics compared to the slow torsional dynamics. A dimensionless formulation exhibits a large parameter in the axial equation, enabling a two-time-scales analysis that uses a combination of averaging methods and a singular perturbation approach. An approximate model of the decoupled axial dynamics permits us to derive a pseudoanalytical expression of the solution of the axial equation. Its averaged behavior influences the slow torsional dynamics by generating an apparent velocity weakening friction law that has been proposed empirically in earlier work. The analytical expression of the solution of the axial dynamics is used to derive an approximate analytical expression of the velocity weakening friction law related to the physical parameters of the system. This expression can be used to provide recommendations on the operating parameters and the drillstring or the bit design in order to reduce the amplitude of the torsional vibrations. Moreover, it is an appropriate candidate model to replace empirical friction laws encountered in torsional models used for control. © 2009 Society for Industrial and Applied Mathematics.

C–H local modes in cyclobutene. II. Laser photoacoustic studies 10 000–17 000 cm−1. Vibrational structure and C–H local mode dynamics

In part I of this study [Baggott, Clase, and Mills, Spectrochim. Acta Part A 42, 319 (1986)] we presented FTIR spectra of gas phase cyclobutene and modeled the v=1–3 stretching states of both olefinic and methylenic C–H bonds in terms of a local mode model. In this paper we present some improvements to our original model and make use of recently derived ‘‘x,K relations’’ to find the equivalent normal mode descriptions. The use of both the local mode and normal mode approaches to modeling the vibrational structure is described in some detail. We present evidence for Fermi resonance interactions between the methylenic C–H stretch overtones and ring C–C stretch vibrations, revealed in laser photoacoustic spectra in the v=4–6 region. An approximate model vibrational Hamiltonian is proposed to explain the observed structure and is used to calculate the dynamics of the C–H stretch local mode decay resulting from interaction with lower frequency ring modes. The implications of our experimental and theoretical studies for mode‐selective photochemistry are discussed briefly.

Measuring toughness and the cohesive stress-displacement relationship by the essential work of fracture concept

The complete fracture behaviour of ductile double edge notched tension (DENT) specimen is analysed with an approximate model, which is then used to discuss the essential work of fracture (EWF) concept. The model results are compared with the experimental results for an aluminium alloy 6082-O. The restrictions on the ligament size for valid application of the EWF method are discussed with the aid of the model. The model is used to suggest an improved method of obtaining the cohesive stress-displacement relationship for the fracture process zone (FPZ).

Control design and genetic algorithm optimization for electrostatic MEMS

This dissertation discusses structural-electrostatic modeling techniques, genetic algorithm based optimization and control design for electrostatic micro devices. First, an alternative modeling technique, the interpolated force model, for electrostatic micro devices is discussed. The method provides improved computational efficiency relative to a benchmark model, as well as improved accuracy for irregular electrode configurations relative to a common approximate model, the parallel plate approximation model. For the configuration most similar to two parallel plates, expected to be the best case scenario for the approximate model, both the parallel plate approximation model and the interpolated force model maintained less than 2.2% error in static deflection compared to the benchmark model. For the configuration expected to be the worst case scenario for the parallel plate approximation model, the interpolated force model maintained less than 2.9% error in static deflection while the parallel plate approximation model is incapable of handling the configuration. Second, genetic algorithm based optimization is shown to improve the design of an electrostatic micro sensor. The design space is enlarged from published design spaces to include the configuration of both sensing and actuation electrodes, material distribution, actuation voltage and other geometric dimensions. For a small population, the design was improved by approximately a factor of 6 over 15 generations to a fitness value of 3.2 fF. For a larger population seeded with the best configurations of the previous optimization, the design was improved by another 7% in 5 generations to a fitness value of 3.0 fF. Third, a learning control algorithm is presented that reduces the closing time of a radiofrequency microelectromechanical systems switch by minimizing bounce while maintaining robustness to fabrication variability. Electrostatic actuation of the plate causes pull-in with high impact velocities, which are difficult to control due to parameter variations from part to part. A single degree-of-freedom model was utilized to design a learning control algorithm that shapes the actuation voltage based on the open/closed state of the switch. Experiments on 3 test switches show that after 5-10 iterations, the learning algorithm lands the switch with an impact velocity not exceeding 0.2 m/s, eliminating bounce.