A methodological paradigm for patient-specific multi-scale CFD simulations: from clinical measurements to parameter estimates for individual analysis.

Physics-based sound radiation estimation from multiple speakers by combined lumped parameter and reduced-order finite element modeling

Reverse stream flowrouting is a procedure that determines the upstream hydrograph given the downstream hydrograph. This paper presents the development of methodology for Muskingum models parameter estimation for reverse stream flow routing. The standard application of the Muskingum models involves calibration and prediction steps. The calibration step must be performed before the prediction step. The calibration step in a reverse stream flow routing system uses the outflow hydrograph and the inflow at the end period of the inflow hydrograph as the known inputs and Muskingum model parameters are determined by minimizing the error between the remaining portion of the predicted and observed inflow hydrographs. In the present study, methodology for parameter estimation is developed which is based on the concept of minimizing the sum of squares of normalized difference between observed and computed inflows subject to the satisfaction of the routing equation. The parameter estimation problems are formulated as constrained nonlinear optimization problem, and a computational scheme is developed to solve the resulting nonlinear problem. The performance evaluation tests indicate that a fresh calibration is necessary to use the Muskingum models for reverse stream flow routing.

Influence of exhaust gas heating and L/D ratios on the discharge efficiencies for an activated carbon natural gas storage system

Nowadays, experimental techniques allow scientists to have access to large amounts of data. In order to obtain reliable information from the complex systems that produce these data, appropriate analysis tools are needed. The Kalman filter is a frequently used technique to infer, assuming a model of the system, the parameters of the model from uncertain observations. A well-known implementation of the Kalman filter, the unscented Kalman filter (UKF), was recently shown to be able to infer the connectivity of a set of coupled chaotic oscillators. In this work, we test whether the UKF can also reconstruct the connectivity of small groups of coupled neurons when their links are either electrical or chemical synapses. In particular, we consider Izhikevich neurons and aim to infer which neurons influence each other, considering simulated spike trains as the experimental observations used by the UKF. First, we verify that the UKF can recover the parameters of a single neuron, even when the parameters vary in time. Second, we analyze small neural ensembles and demonstrate that the UKF allows inferring the connectivity between the neurons, even for heterogeneous, directed, and temporally evolving networks. Our results show that time-dependent parameter and coupling estimation is possible in this nonlinearly coupled system.

The hydrodynamic cavitation in a Venturi tube is studied both theoretically and experimentally. A lumped parameter model was developed to describe the accumulation and dissipation of energy in the biphasic flow as a function of the bubble population characteristics (mean volume, standard deviation, and void fraction). Resistance, capacitance, inductance, and frequency lumped-parameters were identified applying the fluid conservation equations (mass and momentum) along with electrical/hydraulic analogies. Experiments with 1,2-propanediol were carried out in a hydraulic circuit composed of valves, a pump, and a Venturi nozzle. The acoustic noise generated (at different cavitation regimes) by the passage of the fluid through the tube was acquired with a piezoelectric sensor. After processing the experimental signals, the system frequency at each operation condition was determined. Plausible estimations of the void fraction were obtained at different experimental frequencies by evaluating a theoretical expression of the frequency lumped-parameter. This semi-empirical technique might be a low-cost alternative when the void fraction of a flow needs to be determined and tomography devices are not available.

TracerLPM is an interactive Excel® (2007 or later) workbook program for evaluating groundwater age distributions from environmental tracer data by using lumped parameter models (LPMs). Lumped parameter models are mathematical models of transport based on simplified aquifer geometry and flow configurations that account for effects of hydrodynamic dispersion or mixing within the aquifer, well bore, or discharge area. Five primary LPMs are included in the workbook: piston-flow model (PFM), exponential mixing model (EMM), exponential piston-flow model (EPM), partial exponential model (PEM), and dispersion model (DM). Binary mixing models (BMM) can be created by combining primary LPMs in various combinations. Travel time through the unsaturated zone can be included as an additional parameter. TracerLPM also allows users to enter age distributions determined from other methods, such as particle tracking results from numerical groundwater-flow models or from other LPMs not included in this program. Tracers of both young groundwater (anthropogenic atmospheric gases and isotopic substances indicating post-1940s recharge) and much older groundwater (carbon-14 and helium-4) can be interpreted simultaneously so that estimates of the groundwater age distribution for samples with a wide range of ages can be constrained. TracerLPM is organized to permit a comprehensive interpretive approach consisting of hydrogeologic conceptualization, visual examination of data and models, and best-fit parameter estimation. Groundwater age distributions can be evaluated by comparing measured and modeled tracer concentrations in two ways: (1) multiple tracers analyzed simultaneously can be evaluated against each other for concordance with modeled concentrations (tracer-tracer application) or (2) tracer time-series data can be evaluated for concordance with modeled trends (tracer-time application). Groundwater-age estimates can also be obtained for samples with a single tracer measurement at one point in time; however, prior knowledge of an appropriate LPM is required because the mean age is often non-unique. LPM output concentrations depend on model parameters and sample date. All of the LPMs have a parameter for mean age. The EPM, PEM, and DM have an additional parameter that characterizes the degree of age mixing in the sample. BMMs have a parameter for the fraction of the first component in the mixture. An LPM, together with its parameter values, provides a description of the age distribution or the fractional contribution of water for every age of recharge contained within a sample. For the PFM, the age distribution is a unit pulse at one distinct age. For the other LPMs, the age distribution can be much broader and span decades, centuries, millennia, or more. For a sample with a mixture of groundwater ages, the reported interpretation of tracer data includes the LPM name, the mean age, and the values of any other independent model parameters. TracerLPM also can be used for simulating the responses of wells, springs, streams, or other groundwater discharge receptors to nonpoint-source contaminants that are introduced in recharge, such as nitrate. This is done by combining an LPM or user-defined age distribution with information on contaminant loading at the water table. Information on historic contaminant loading can be used to help evaluate a model's ability to match real world conditions and understand observed contaminant trends, while information on future contaminant loading scenarios can be used to forecast potential contaminant trends.

Children with single ventricle heart disease typically require a series of three operations, (1) Norwood, (2) Glenn, and (3) Fontan, which ultimately results in complete separation of the pulmonary and systemic circuits to improve pulmonary/systemic circulation. In the last stage, the Fontan operation, the inferior vena cava (IVC) is connected to the pulmonary arteries (PAs), allowing the remainder of deoxygenated blood to passively flow to the pulmonary circuit. It is hypothesized that optimizing the Fontan anatomy would lead to decreased power loss and more balanced hepatic flow distribution. One approach to optimizing the geometry is to create a patient-specific digital twin to simulate various configurations of the Fontan conduit, which requires a computational model of the proximal PA anatomy and resistance, as well as the distal Pulmonary Vascular Resistance (PVR), at the Glenn stage. To that end, an optimization pipeline was developed using 3D computational fluid dynamics (CFD) and 0D lumped parameter (LP) simulations to iteratively refine the PVR of each lung by minimizing the simulated flow and pressure error relative to patients' cardiac magnetic resonance (CMR) and catheterization (CATH) data. While the PVR can also be estimated directly by computing the ratio of pressure gradients and flow from CATH and CMR data, the computational approach can separately identify the different components of PVR along the Glenn pathway, allowing for a more detailed depiction of the Glenn vasculature. Results indicate good correlation between the optimized PVR of the CFD and LP models (n = 16), with an intraclass correlation coefficient (ICC) of 0.998 (p = 0.976) and 0.991 (p = 0.943) for the left and right lung, respectively. Furthermore, compared to CMR flow and CATH pressure data, the optimized PVR estimates result in mean outlet flow and pressure errors of less than 5%. The optimized PVR estimates also agree well with the computed PVR estimates from CATH pressure and CMR flow for both lungs, yielding a mean difference of less than 4%.

This paper studies the performance of recursive and batch Bayesian methods for nonlinear model updating. Unscented Kalman filter (UKF) is selected to represent the recursive Bayesian method, and two UKF approaches are investigated and compared, i.e., non-adaptive UKF and adaptive UKF. The proposed new adaptive filter, forgetting factor adaptive UKF, estimates the model parameters and measurement noise covariance in an online manner. The forgetting factor adaptive UKF is based on the principle of matching the covariance of residuals to its theoretical values by updating the measurement noise covariance. The performance of non-adaptive UKF, adaptive UKF and batch Bayesian method are investigated when applied to a numerical nonlinear 3-story 3-bay steel frame structure for parameter estimation of material properties. Different types of modeling errors are considered in the 21 updating models to study the effects of modeling errors on model updating. It is found that adaptive UKF approach provides the most accurate parameter estimations, while batch Bayesian approach gives the smallest errors on response predictions.KeywordsUnscented Kalman filterAdaptive UKFMeasurement noise covarianceBayesian model updatingModeling errors

This chapter deals with the problem of non-linear and non-Gaussian states and parameters estimation using Bayesian methods. The performances of various conventional and state-of-the-art state estimation techniques are compared when they are utilized to achieve this objective. These techniques include the Extended Kalman Filter (EKF), Unscented Kalman Filter (UKF), and Particle Filter (PF). In the current work, the authors consider two systems (biological model and power system) to perform evaluation of estimation algorithms. The results of the comparative studies show that the UKF provides a higher accuracy than the EKF due to the limited ability of EKF to accurately estimate the mean and covariance matrix of the estimated states through lineralization of the nonlinear process model. The results also show that the PF provides a significant improvement over the UKF because, unlike UKF, PF is not restricted by linear-Gaussian assumptions which greatly extends the range of problems that can be tackled.

In the first part of this work [1], we presented the methodology for simultaneous soil moisture and hydraulic parameter estimation. In [1], sensitivity analysis and orthogonalization projection were used to select a subset of the most estimable states and hydraulic parameters at every sampling time, an extended Kalman filtering approach that can take into account the variable selection results was also presented. The methodology was demonstrated to be effective using a simulated case study.In the second part of this work, we further apply the simultaneous estimation method to a real large-scale agriculture field in Lethbridge, Alberta, Canada. The field is equipped with a centre pivot irrigation system and microwave soil moisture sensors are mounted on the centre pivot. Cross-validation results obtained for two specific days revealed that the state and parameter estimation procedure improved the accuracy of the soil moisture estimates by 24% and 43%.

We used ultrasonic shear waves for nondestructive defect detection in a billet using transmitted waves. We utilized the deviation the time-of-flight (TOF) obtained by cross-correlation of transmitted waves of a defect-free reference plane and that of a measurement plane containing a defect. We compared the performance of longitudinal waves and shear waves at different wavelengths in detecting the diameter of a circular defect in two-dimensional (2D) simulation and the TOF for a cylindrical defect while changing the vibration direction of shear waves in three-dimensional (3D) simulation. Shear waves detected defects better than longitudinal waves in the 2D simulation, especially at wavelengths of 1.4–2.4 mm. In the 3D simulation, the maximum TOF was larger when the vibration direction was perpendicular to the defect’s major axis than when it was parallel in the measurement using shear waves. This suggests a defect’s shape can be estimated by measurement using shear waves.

SummaryThis paper addresses the problem of nonlinear time‐varying state and parameter estimation of induction machines (IMs) on the basis of a third‐order electrical model. The objectives of this paper are threefold. The first objective is to propose the use of an improved particle filter (IPF) with better proposal distribution for nonlinear and non‐Gaussian state and parameter estimation. The second objective is to extend the state and parameter estimation techniques (i.e., extended Kalman filter (EKF), unscented Kalman filter (UKF), particle filter (PF), and IPF) to better handle nonlinear and non‐Gaussian processes withouta prioristate information, by utilizing a time‐varying assumption of statistical parameters. In this case, the state vector to be estimated at any instant is assumed to follow a Gaussian model, where the expectation and the covariance matrix are both random. The third objective is to compare the performances of EKF, UKF, PF, and IPF in estimating the states of the power process model representing the IM (i.e, the rotor speed, the rotor flux, the stator flux, the rotor resistance, and the magnetizing inductance) and their abilities to estimate some of the key system parameters, which are needed to define the IM process model. The results show that the IPF provides a significant improvement over the PF because, unlike the PF, which depends on the choice of sampling distribution used to estimate the posterior distribution, the IPF yields an optimum choice of the sampling distribution, which also accounts for the observed data. This conclusion is also supported by the experimental results. Copyright © 2014 John Wiley & Sons, Ltd.

The normal unscented Kalman filter (UKF) suffers from performance degradation and even divergence while mismatch between the noise distribution assumed to be known as a priori by UKF and the true ones in a real system. In this paper, a novel adaptive UKF (AUKF) is developed and applied to nonlinear joint estimation of both time-varying states and parameters. A cost function is built based on the error between the covariance matrices of innovation and their corresponding estimations. An adaptive algorithm is then designed to online update the covariance of the process noise by minimizing the cost function. The updated covariance is further fed back into the normal UKF. As a result of using such an adaptive mechanism, the robustness of conventional UKF is substantially improved with respect to the uncertain and time-varying noise covariance in the real system. To illustrate this mechanism, simulations are conducted on the dynamics of an unmanned helicopter by jointly estimating both the states and model errors. The improvements of the proposed AUKF are demonstrated by comparing the results with and without the adaptive mechanism.

Requirements of highly precise pointing performance have been imposed on recently developed spacecrafts for a variety of missions. The stringent requirements have called on on-orbit estimation of spacecraft dynamics parameters and calibration of on-board sensors as indispensible practices. Consequently, on-orbit estimation of the mass moment of inertia of spacecraft has been a major issue mostly due to the changes by solar panel deployment and a large portion of fuel consumption (Creamer et al., 1996; Ahmed et al., 1998; Bordany et al., 2000; VanDyke et al., 2004; Myung et al., 2007; Myung & Bang, 2008; Sekhavat et al., 2009). As for measurement sensors, on-board calibration of alignment and bias errors of attitude and rate sensors is one of the main concerns of attitude sensor calibration researches (Pittelkau, 2001 & 2002, Lai et al., 2003). Pittelkau (2002) proposed an attitude estimator based on the Kalman filter (Kalman, 1960), in which spacecraft attitude quaternion, rate sensor misalignment and bias, and star tracker misalignments are taken into consideration as states, whereas the body rate is dealt as a synthesized signal by the estimates. Lai at al. (2003) derived a method for alignment estimation of attitude and rate sensors based on the unscented Kalman filter (UKF) (Julier and Uhlmann, 1997). Ma and Jiang (2005) presented spacecraft attitude estimation and calibration based only magnetometer measurements using an UKF. An interesting point is that we need predesigned 3-axis excitation manoeuvres of spacecraft for both dynamics parameter estimation and sensor calibration. Therefore, this study is motivated to merge above estimation and calibration processes into a single filtering problem. It is noteworthy that poor information of moments of inertia is to be treated as a system uncertainty while the rate sensor model errors are to be incorporated into the measurement process. As a filtering algorithm, this study employs a UKF. Extended Kalman filters (EKFs) have been successfully applied to the nonlinear attitude estimation problem (Crassidis et al., 2007). Hybrid estimation using the EKF has been reported by Myung at al. (2007). However, the EKF estimates using the first order linearization, which may lead to instability of the filter (ValDyke et al., 2004). The UKF approximates the nonlinear model to the second order by spreading points 1 sigma apart from the a priori mean. Performing nonlinear

Many interplanetary mission concepts can benefit from autonomous orbit estimation, particularly during critical mission phases. Previous studies have examined the feasibility of optical navigation using nanosatellite class instruments. While promising, these techniques are not without drawbacks. Convergence of the navigation estimates are often sensitive to errors in initial state estimates. This thesis compares various methods to perform nonlinear estimation for autonomous optical navigation. These methods include an extended Kalman filter (EKF), an unscented Kalman filter (UKF), a particle filter (PF), a fixed-lag smoother (FLS), and moving horizon estimation (MHE). The EKF, UKF, and PF can be implemented in real time, while the FLS and MHE implement a delay into the estimation process. To compare the performance of each state estimator three initial reference scenarios around Mars were considered: a hyperbolic flyby, an elliptic orbit and a orbital maneuver using observations of Mars and its moons. Parameter estimation was also explored, where the mass of Mars was to be estimated as a reference parameter in both the hyperbolic and elliptical trajectories. One last reference scenario included a low Earth orbit (LEO) using observations of satellites in a geosynchronous equatorial orbit. In each case, the FLS and MHE showed similar or better performance over each state estimator but at the cost of an increased computation time with respect to the reference EKF. Similarly the UKF was able to provide improved results withe respect to the EKF. While, the PF provided poor estimates in the Mars trajectories but improvements were seen from the UKF and EKF in the LEO scenario.