Simultaneous Reduced Basis Approximation of Parameterized Elliptic Eigenvalue Problems

The focus is on a model reduction framework for parameterized elliptic eigenvalue problems by a reduced basis method. In contrast to the standard single output case, one is interested in approximating several outputs simultaneously, namely a certain number of the smallest eigenvalues. For a fast and reliable evaluation of these input-output relations, we analyze a posteriori error estimators for eigenvalues. Moreover, we present different greedy strategies and study systematically their performance. Special attention needs to be paid to multiple eigenvalues whose appearance is parameter-dependent. Our methods are of particular interest for applications in vibro-acoustics

Anisotropic behaviour of human gallbladder walls

Inverse estimation of biomechanical parameters of soft tissues from non-invasive measurements has clinical significance in patient-specific modelling and disease diagnosis. In this paper, we propose a fully nonlinear approach to estimate the mechanical properties of the human gallbladder wall muscles from in vivo ultrasound images. The iteration method consists of a forward approach, in which the constitutive equation is based on a modified Hozapfel–Gasser–Ogden law initially developed for arteries. Five constitutive parameters describing the two orthogonal families of fibres and the matrix material are determined by comparing the computed displacements with medical images. The optimisation process is carried out using the MATLAB toolbox, a Python code, and the ABAQUS solver. The proposed method is validated with published artery data and subsequently applied to ten human gallbladder samples. Results show that the human gallbladder wall is anisotropic during the passive refilling phase, and that the peak stress is 1.6 times greater than that calculated using linear mechanics. This discrepancy arises because the wall thickness reduces by 1.6 times during the deformation, which is not predicted by conventional linear elasticity. If the change of wall thickness is accounted for, then the linear model can used to predict the gallbladder stress and its correlation with pain. This work provides further understanding of the nonlinear characteristics of human gallbladder

A 3D discrete model of the diaphragm and human trunk

In this paper, a 3D discrete model is presented to model the movements of the trunk during breathing. In this model, objects are represented by physical particles on their contours. A simple notion of force generated by a linear actuator allows the model to create forces on each particle by way of a geometrical attractor. Tissue elasticity and contractility are modeled by local shape memory and muscular fibers attractors. A specific dynamic MRI study was used to build a simple trunk model comprised of by three compartments: lungs, diaphragm and abdomen. This model was registered on the real geometry. Simulation results were compared qualitatively as well as quantitatively to the experimental data, in terms of volume and geometry. A good correlation was obtained between the model and the real data. Thanks to this model, pathology such as hemidiaphragm paralysis can also be simulated.Comment: published in: "Lung Modelling", France (2006

Model validation for a noninvasive arterial stenosis detection problem

Copyright @ 2013 American Institute of Mathematical SciencesA current thrust in medical research is the development of a non-invasive method for detection, localization, and characterization of an arterial stenosis (a blockage or partial blockage in an artery). A method has been proposed to detect shear waves in the chest cavity which have been generated by disturbances in the blood flow resulting from a stenosis. In order to develop this methodology further, we use both one-dimensional pressure and shear wave experimental data from novel acoustic phantoms to validate corresponding viscoelastic mathematical models, which were developed in a concept paper [8] and refined herein. We estimate model parameters which give a good fit (in a sense to be precisely defined) to the experimental data, and use asymptotic error theory to provide confidence intervals for parameter estimates. Finally, since a robust error model is necessary for accurate parameter estimates and confidence analysis, we include a comparison of absolute and relative models for measurement error.The National Institute of Allergy and Infectious Diseases, the Air Force Office of Scientific Research, the Deopartment of Education and the Engineering and Physical Sciences Research Council (EPSRC)

Financial innovation and money demand : theory and empirical implementation

Empirically, traditional money demand equations are frequently characterized by periods of"missing money", unstable parameters, and autocorrelated errors. The common practice to solve these problems consists of changing the specification of the regressions once the shifts (which are usually associated to financial innovation) are identified. This paper provides an alternative approach to dealing with the unobservable process of financial innovation. It consists of modelling financial innovation as shocks that have permanent effects on the money demand, analogous to productivity shocks in production functions. This paper describes the theoretical model used and shows the failure of traditional money demand equations using cointegration techiques. It describes a simple GLS-iterative econometric model which allows the authors to recover the path of financial innnovation and obtain sensible estimates of the relevant elasticities. It also shows Monte Carlo simulations to evaluate the behavior of the estimation procedure for particular samples and data generating processes, and to study how robust the procedure to some deviations from the basic assumptions is.Economic Theory&Research,Environmental Economics&Policies,Fiscal&Monetary Policy,Inequality,Financial Intermediation