Numerical solutions of elliptic inverse problems via the equation error method

To estimate a parameter in an elliptic boundary value problem, the method of equation error chooses the value that minimizes the error in the PDE and boundary condition (the solution of the BVP having been replaced by a measurement). The estimated parameter converges to the exact value as the measured data converge to the exact value, provided Tikhonov regularization is used to control the instability inherent in the problem. The error in the estimated solution can be bounded in an appropriate quotient norm; estimates can be derived for both the underlying (infinite-dimensional) problem and a finite-element discretization that can be implemented in a practical algorithm. Numerical experiments demonstrate the efficacy and limitations of the method.

Inverse thermodilution with conventional pulmonary artery catheters for the assessment of cerebral, hepatic, renal, and femoral blood flow

Assessment of regional blood flow changes is difficult in the clinical setting. We tested whether conventional pulmonary artery catheters (PACs) can be used to measure regional venous blood flows by inverse thermodilution (ITD). Inverse thermodilution was tested in vitro and in vivo using perivascular ultrasound Doppler (USD) flow probes as a reference. In anesthetized pigs, PACs were inserted in jugular, hepatic, renal, and femoral veins, and their measurements were compared with simultaneous USD flow measurements from carotid, hepatic, renal, and femoral arteries and from portal vein. Fluid boluses were injected through the PAC's distal port, and temperature changes were recorded from the proximally located thermistor. Injectates of 2 and 5 mL at 22 degrees C and 4 degrees C were used. Flows were altered by using a roller pump (in vitro), and infusion of dobutamine and induction of cardiac tamponade, respectively. In vitro: At blood flows between 400 mL . min-1 and 700 mL . min-1 (n = 50), ITD and USD correlated well (r = 0.86, P < 0.0001), with bias and limits of agreement of 3 +/- 101 mL . min-1. In vivo: 514 pairs of measurements had to be excluded from analysis for technical reasons, and 976 were analyzed. Best correlations were r = 0.87 (P < 0.0001) for renal flow and r = 0.46 (P < 0.0001) for hepatic flow. No significant correlation was found for cerebral and femoral flows. Inverse thermodilution using conventional PAC compared moderately well with USD for renal but not for other flows despite good in vitro correlation in various conditions. In addition, this method has significant technical limitations.

SPATIAL UNDERSTANDING OF MATRIX INVERSION FOR INVERSE FORCE ESTIMATION

A basic approach to study a NVH problem is to break down the system in three basic elements – source, path and receiver. While the receiver (response) and the transfer path can be measured, it is difficult to measure the source (forces) acting on the system. It becomes necessary to predict these forces to know how they influence the responses. This requires inverting the transfer path. Singular Value Decomposition (SVD) method is used to decompose the transfer path matrix into its principle components which is required for the inversion. The usual approach to force prediction requires rejecting the small singular values obtained during SVD by setting a threshold, as these small values dominate the inverse matrix. This assumption of the threshold may be subjected to rejecting important singular values severely affecting force prediction. The new approach discussed in this report looks at the column space of the transfer path matrix which is the basis for the predicted response. The response participation is an indication of how the small singular values influence the force participation. The ability to accurately reconstruct the response vector is important to establish a confidence in force vector prediction. The goal of this report is to suggest a solution that is mathematically feasible, physically meaningful, and numerically more efficient through examples. This understanding adds new insight to the effects of current code and how to apply algorithms and understanding to new codes.

Characterization of water-soluble organic compounds in ambient aerosol using ultrahigh-resolution elctrospray ionization fourier transform ion cyclotron resonance mass spectrometry.

Atmospheric aerosol water-soluble organic compounds (WSOC) exist in a complex mixture of thousands of organic compounds which may have a significant influence on the climate-relevant properties of the atmospheric aerosol. To understand the potential influences, the ambient aerosol was collected at a nonurban mountainous site near Steamboat Springs, CO. The WSOC fraction was analyzed using positive and negative electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry. Approximately 2400 and 4000 molecular formulas were identified from the detected positive and negative ions, respectively. The formulas contained carbon (C), hydrogen (H), oxygen (O), nitrogen (N), and sulfur (S) atoms over the mass range of 100-800 Da in both ionization modes. The number range of double bond equivalents (DBE), the mean O:C, H:C, and oxidation state of carbon for the positive ions were 0 – 18, 0.25 ± 0.15, 1.39 ± 0.29, and -0.89 ± 0.23, respectively. Comparatively, the negative ion values were 0 – 14, 0.53 ± 0.20, 1.48 ± 0.30, and -0.41 ± 0.45, respectively. Overall, the positive ion molecular formulas were less oxygenated than negative ions as seen with the lower O:C and OSc values. Molecular formulas of the positive ions classified as aliphatic, olefinic, and aromatic compound classes based on the aromaticity index values. Aliphatic compounds were the CHNO and CHO formulas that had mean DBE values of about 5 and 3, respectively. However, a majority of the CHOS, CHNOS, and CHS formulas were defined as olefinic compounds and had mean DBE values of about 12, 13, and 10, respectively. Overall, more than half of the assigned molecular formulas contained sulfur and were olefinic to aromatic compounds with a DBE range of 7-18. Source of the unsaturated sulfur containing compounds is currently unknown. Several nitrogen containing compounds were in common with the field and laboratory studies of the biomass burning aerosol and aged secondary organic aerosol products of the limonene ozonolysis.

Real-Time Joint Coupling of the Spine for Inverse Kinematics

In this paper we propose a simple model for the coupling behavior of the human spine for an inverse kinematics framework. Our spine model exhibits anatomically correct motions of the vertebrae of virtual mannequins by coupling standard swing and revolute joint models. The adjustement of the joints is made with several simple (in)equality constraints, resulting in a reduction of the solution space dimensionality for the inverse kinematics solver. By reducing the solution space dimensionality to feasible spine shapes, we prevent the inverse kinematics algorithm from providing infeasible postures for the spine.In this paper, we exploit how to apply these simple constraints to the human spine by a strict decoupling of the swing and torsion motion of the vertebrae. We demonstrate the validity of our approach on various experiments.

The von Mises-Fisher distribution of the first exit point from the hypersphere of the drifted Brownian motion and the density of the first exit time

A characterization is provided for the von Mises–Fisher random variable, in terms of first exit point from the unit hypersphere of the drifted Wiener process. Laplace transform formulae for the first exit time from the unit hypersphere of the drifted Wiener process are provided. Post representations in terms of Bell polynomials are provided for the densities of the first exit times from the circle and from the sphere.

Distance-based kriging relying on proxy simulations for inverse conditioning

Let us consider a large set of candidate parameter fields, such as hydraulic conductivity maps, on which we can run an accurate forward flow and transport simulation. We address the issue of rapidly identifying a subset of candidates whose response best match a reference response curve. In order to keep the number of calls to the accurate flow simulator computationally tractable, a recent distance-based approach relying on fast proxy simulations is revisited, and turned into a non-stationary kriging method where the covariance kernel is obtained by combining a classical kernel with the proxy. Once the accurate simulator has been run for an initial subset of parameter fields and a kriging metamodel has been inferred, the predictive distributions of misfits for the remaining parameter fields can be used as a guide to select candidate parameter fields in a sequential way. The proposed algorithm, Proxy-based Kriging for Sequential Inversion (ProKSI), relies on a variant of the Expected Improvement, a popular criterion for kriging-based global optimization. A statistical benchmark of ProKSI’s performances illustrates the efficiency and the robustness of the approach when using different kinds of proxies.

Inverse relationship between fractionated electrograms and atrial fibrosis in persistent atrial fibrillation: combined magnetic resonance imaging and high-density mapping

OBJECTIVES This study sought to evaluate the relationship between fibrosis imaged by delayed-enhancement (DE) magnetic resonance imaging (MRI) and atrial electrograms (Egms) in persistent atrial fibrillation (AF). BACKGROUND Atrial fractionated Egms are strongly related to slow anisotropic conduction. Their relationship to atrial fibrosis has not yet been investigated. METHODS Atrial high-resolution MRI of 18 patients with persistent AF (11 long-lasting persistent AF) was registered with mapping geometry (NavX electro-anatomical system (version 8.0, St. Jude Medical, St. Paul, Minnesota)). DE areas were categorized as dense or patchy, depending on their DE content. Left atrial Egms during AF were acquired using a high-density, 20-pole catheter (514 ± 77 sites/map). Fractionation, organization/regularity, local mean cycle length (CL), and voltage were analyzed with regard to DE. RESULTS Patients with long-lasting persistent versus persistent AF had larger left atrial (LA) surface area (134 ± 38 cm(2) vs. 98 ± 9 cm(2), p = 0.02), a higher amount of atrial DE (70 ± 16 cm(2) vs. 49 ± 10 cm(2), p = 0.01), more complex fractionated atrial Egm (CFAE) extent (54 ± 16 cm(2) vs. 28 ± 15 cm(2), p = 0.02), and a shorter baseline AF CL (147 ± 10 ms vs. 182 ± 14 ms, p = 0.01). Continuous CFAE (CFEmean [NavX algorithm that quantifies Egm fractionation] <80 ms) occupied 38 ± 19% of total LA surface area. Dense DE was detected at the left posterior left atrium. In contrast, the right posterior left atrium contained predominantly patchy DE. Most CFAE (48 ± 14%) occurred at non-DE LA sites, followed by 41 ± 12% CFAE at patchy DE and 11 ± 6% at dense DE regions (p = 0.005 and p = 0.008, respectively); 19 ± 6% CFAE sites occurred at border zones of dense DE. Egms were less fractionated, with longer CL and lower voltage at dense DE versus non-DE regions: CFEmean: 97 ms versus 76 ms, p < 0.0001; local CL: 153 ms versus 143 ms, p < 0.0001; mean voltage: 0.63 mV versus 0.86 mV, p < 0.0001. CONCLUSIONS Atrial fibrosis as defined by DE MRI is associated with slower and more organized electrical activity but with lower voltage than healthy atrial areas. Ninety percent of continuous CFAE sites occur at non-DE and patchy DE LA sites. These findings are important when choosing the ablation strategy in persistent AF.

Universally Applicable Model for the Quantitative Determination of Lake Sediment Composition Using Fourier Transform Infrared Spectroscopy

ABSTRACT: Fourier transform infrared spectroscopy (FTIRS) can provide detailed information on organic and minerogenic constituents of sediment records. Based on a large number of sediment samples of varying age (0�340 000 yrs) and from very diverse lake settings in Antarctica, Argentina, Canada, Macedonia/Albania, Siberia, and Sweden, we have developed universally applicable calibration models for the quantitative determination of biogenic silica (BSi; n = 816), total inorganic carbon (TIC; n = 879), and total organic carbon (TOC; n = 3164) using FTIRS. These models are based on the differential absorbance of infrared radiation at specific wavelengths with varying concentrations of individual parameters, due to molecular vibrations associated with each parameter. The calibration models have low prediction errors and the predicted values are highly correlated with conventionally measured values (R = 0.94�0.99). Robustness tests indicate the accuracy of the newly developed FTIRS calibration models is similar to that of conventional geochemical analyses. Consequently FTIRS offers a useful and rapid alternative to conventional analyses for the quantitative determination of BSi, TIC, and TOC. The rapidity, cost-effectiveness, and small sample size required enables FTIRS determination of geochemical properties to be undertaken at higher resolutions than would otherwise be possible with the same resource allocation, thus providing crucial sedimentological information for climatic and environmental reconstructions.

Fourier transform infrared spectroscopy, a new method for rapid determination of total organic and inorganic carbon and biogenic silica concentration in lake sediments

Abstract We demonstrate the use of Fourier transform infrared spectroscopy (FTIRS) to make quantitative measures of total organic carbon (TOC), total inorganic carbon (TIC) and biogenic silica (BSi) concentrations in sediment. FTIRS is a fast and costeffective technique and only small sediment samples are needed (0.01 g). Statistically significant models were developed using sediment samples from northern Sweden and were applied to sediment records from Sweden, northeast Siberia and Macedonia. The correlation between FTIRS-inferred values and amounts of biogeochemical constituents assessed conventionally varied between r = 0.84–0.99 for TOC, r = 0.85– 0.99 for TIC, and r = 0.68–0.94 for BSi. Because FTIR spectra contain information on a large number of both inorganic and organic components, there is great potential for FTIRS to become an important tool in paleolimnology.