Frequency-domain identification of continuous-time ARMA models from sampled data

The subject of this paper is the direct identification of continuous-time autoregressive moving average (CARMA) models. The topic is viewed from the frequency domain perspective which then turns the reconstruction of the continuous-time power spectral density (CT-PSD) into a key issue. The first part of the paper therefore concerns the approximate estimation of the CT-PSD from uniformly sampled data under the assumption that the model has a certain relative degree. The approach has its point of origin in the frequency domain Whittle likelihood estimator. The discrete- or continuous-time spectral densities are estimated from equidistant samples of the output. For low sampling rates the discrete-time spectral density is modeled directly by its continuous-time spectral density using the Poisson summation formula. In the case of rapid sampling the continuous-time spectral density is estimated directly by modifying its discrete-time counterpart.     

The effect of the microfluidic diodicity on the efficiency of valve-less rectification micropumps using Lattice Boltzmann Method

The efficiency of the valve-less rectification micropump depends primarily on the microfluidic diodicity (the ratio of the backward pressure drop to the forward pressure drop). In this study, different rectifying structures, including the conventional structures (nozzle/diffuser and Tesla structures), were investigated at very low Reynolds numbers (between 0.2 and 60). The rectifying structures were characterized with respect to their design, and a numerical approach was illustrated to calculate the diodicity for the rectifying structures. In this study, the microfluidic diodicity was evaluated numerically for different rectifying structures including half circle, semicircle, heart, triangle, bifurcation, nozzle/diffuser, and Tesla structures. The Lattice Boltzmann Method (LBM) was utilized as a numerical method to simulate the fluid flow in the microscale. The results suggest that at very low Reynolds number flow, rectification and multifunction micropumping may be achievable by using a number of the presented structures. The results for the conventional structures agree with the reported results.     

Managing dual warehouses with an incentive policy for deteriorating items

Distributors in a supply chain usually limit their own warehouse in finite capacity for cost reduction and excess stock is held in a rent warehouse. In this study, we examine inventory control for deteriorating items in a two-warehouse setting. Assuming that there is an incentive offered by a rent warehouse that allows the rental fee to decrease over time, the objective of this study is to maximise the joint profit of the manufacturer and the distributor. An optimisation procedure is developed to derive the optimal joint economic lot size policy. Several criteria are identified to select the most appropriate warehouse configuration and inventory policy on the basis of storage duration of materials in a rent warehouse. Sensitivity analysis is done to examine the results of model robustness. The proposed model enables a manufacturer with a channel distributor to coordinate the use of alternative warehouses, and to maximise the joint profit of the manufacturer and the distributor.     

A Modeling Framework for Troubleshooting Automotive Systems

This article presents a novel framework for modeling the troubleshooting process for automotive systems such as trucks and buses. We describe how a diagnostic model of the troubleshooting process can be created using event-driven, nonstationary, dynamic Bayesian networks. Exact inference in such a model is in general not practically possible. Therefore, we evaluate different approximate methods for inference based on the Boyen–Koller algorithm. We identify relevant model classes that have particular structure such that inference can be made with linear time complexity. We also show how models created using expert knowledge can be tuned using statistical data. The proposed learning mechanism can use data that is collected from a heterogeneous fleet of modular vehicles that can consist of different components. The proposed framework is evaluated both theoretically and experimentally on an application example of a fuel injection system.     

Computational performance of a parallelized three-dimensional high-order spectral element toolbox

In this paper, a comprehensive performance review of an MPI-based high-order three-dimensional spectral element method C++ toolbox is presented. The focus is put on the performance evaluation of several aspects with a particular emphasis on the parallel efficiency. The performance evaluation is analyzed with the help of a time prediction model based on a parameterization of the application and the hardware resources. Two tailor-made benchmark cases in computational fluid dynamics (CFD) are introduced and used to carry out this review, stressing the particular interest for clusters with up to thousands of cores. Some problems in the parallel implementation have been detected and corrected. The theoretical complexities with respect to the number of elements, to the polynomial degree, and to communication needs are correctly reproduced. It is concluded that this type of code has a nearly perfect speedup on machines with thousands of cores, and is ready to make the step to next-generation petaFLOP machines.     

Robust interactive cutting based on an adaptive octree simulation mesh

We present an adaptive octree based approach for interactive cutting of deformable objects. Our technique relies on efficient refine- and node split-operations. These are sufficient to robustly represent cuts in the mechanical simulation mesh. A high-resolution surface embedded into the octree is employed to represent a cut visually. Model modification is performed in the rest state of the object, which is accomplished by back-transformation of the blade geometry. This results in an improved robustness of our approach. Further, an efficient update of the correspondences between simulation elements and surface vertices is proposed. The robustness and efficiency of our approach is underlined in test examples as well as by integrating it into a prototype surgical simulator.     

A quantum chemical study of the mechanism of action of Vitamin K epoxide reductase (VKOR) II. Transition states

A reaction path including transition states is generated for the Silverman mechanism [R.B. Silverman, Chemical model studies for the mechanism of Vitamin K epoxide reductase, J. Am. Chem. Soc. 103 (1981) 5939-5941] of action for Vitamin K epoxide reductase (VKOR) using quantum mechanical methods (B3LYP/6-311G**). VKOR, an essential enzyme in mammalian systems, acts to convert Vitamin K epoxide, formed by Vitamin K carboxylase, to its (initial) quinone form for cellular reuse. This study elaborates on a prior work that focused on the thermodynamics of VKOR [D.W. Deerfield II, C.H. Davis, T. Wymore, D.W. Stafford, L.G. Pedersen, Int. J. Quant. Chem. 106 (2006) 2944-2952]. The geometries of proposed model intermediates and transition states in the mechanism are energy optimized. We find that once a key disulfide bond is broken, the reaction proceeds largely downhill. An important step in the conversion of the epoxide back to the quinone form involves initial protonation of the epoxide oxygen. We find that the source of this proton is likely a free mercapto group rather than a water molecule. The results are consistent with the current view that the widely used drug Warfarin likely acts by blocking binding of Vitamin K at the VKOR active site and thereby effectively blocking the initiating step. These results will be useful for designing more complete QM/MM studies of the enzymatic pathway once three-dimensional structural data is determined and available for VKOR.     

Efficient angular spectrum decomposition of acoustic sources. I. Theory

The angular spectrum decomposition is evaluated in terms of plane wave angular range, angular resolution, and spatial aliasing error using two-dimensional FFT (2-D FFT). The algorithm makes possible the source plane decomposition of normal velocity and pressure fields radiated by transducers of arbitrary shape, with significantly faster results achievable for planar sources. Although the angular spectrum is equally applicable to fields far from the transducer, the efficient calculation is derived specifically for fields in or very close to the source plane. An antialiasing algorithm is proposed that allows the source to be discretized with fewer sample points for a given accuracy than required with simple discretization techniques. Guidelines for the selection of sampling interval, discretization size, etc. are developed on an application-specific basis and indicate the best ratio of numerical accuracy to computational cost.     

Maximization of eigenvalues using topology optimization

Topology optimization is used to optimize the eigenvalues of plates. The results are intended especially for MicroElectroMechanical Systems (MEMS) but can be seen as more general. The problem is not formulated as a case of reinforcement of an existing structure, so there is a problem related to localized modes in low density areas. The topology optimization problem is formulated using the SIMP method. Special attention is paid to a numerical method for removing localized eigenmodes in low density areas. The method is applied to numerical examples of maximizing the first eigenfrequency. One example is a practical MEMS application; a probe used in an Atomic Force Microscope (AFM). For the AFM probe the optimization is complicated by a constraint on the stiffness and constraints on higher order eigenvalues. Received June 10, 1999     

Yield and transient effects during the plastic deformation of solid polymers

Tensile tests were performed on seven commercial polymers at 22° C and at constant true strain rates of 10^–4 to 10^–1 sec^–1. The constant strain rates were imposed on the minimum section of each sample with the aid of a diametral transducer, an exponential function generator and a closed-loop hydraulic testing machine. The polymers investigated were: high and low density polyethylene, polytetrafluoroethylene, polypropylene. polyvinylchloride and polyamide 6 and 66. True yield drops were observed in the rigid glassy polymers, whereas yielding was more gradual in the semi-crystalline or plasticized polymers. Strain rate change tests were also performed, during which one order of magnitude increases and decreases were imposed on the specimens. Normal transients were observed at small strains in the samples containing a rubbery phase, while the transients were of an inverse nature in the samples containing a glassy phase. With an increase in the strain at which the change was initiated, the normal transients changed in character to inverse. Transient tests were also performed in which straining was interrupted to permit a period of stress relaxation or of holding in the unloaded condition prior to the resumption of straining. A quantitative model is proposed, based on the dynamics of plastic waves which accounts for the transition from normal to inverse transient behaviour with increasing strain, and also explains the opposite effects of stress relaxation and of specimen unloading on the restraining transients.