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.

Two-Phase Flow In Microchannels: Morphology And Interface Phenomena

The existence and morphology, as well as the dynamics of micro-scale gas-liquid interfaces is investigated numerically and experimentally. These studies can be used to assess liquid management issues in microsystems such as PEMFC gas flow channels, and are meant to open new research perspectives in two-phase flow, particularly in film deposition on non-wetting surfaces. For example the critical plug volume data can be used to deliver desired length plugs, or to determine the plug formation frequency. The dynamics of gas-liquid interfaces, of interest for applications involving small passages (e.g. heat exchangers, phase separators and filtration systems), was investigated using high-speed microscopy - a method that also proved useful for the study of film deposition processes. The existence limit for a liquid plug forming in a mixed wetting channel is determined by numerical simulations using Surface Evolver. The plug model simulate actual conditions in the gas flow channels of PEM fuel cells, the wetting of the gas diffusion layer (GDL) side of the channel being different from the wetting of the bipolar plate walls. The minimum plug volume, denoted as critical volume is computed for a series of GDL and bipolar plate wetting properties. Critical volume data is meant to assist in the water management of PEMFC, when corroborated with experimental data. The effect of cross section geometry is assessed by computing the critical volume in square and trapezoidal channels. Droplet simulations show that water can be passively removed from the GDL surface towards the bipolar plate if we take advantage on differing wetting properties between the two surfaces, to possibly avoid the gas transport blockage through the GDL. High speed microscopy was employed in two-phase and film deposition experiments with water in round and square capillary tubes. Periodic interface destabilization was observed and the existence of compression waves in the gas phase is discussed by taking into consideration a naturally occurring convergent-divergent nozzle formed by the flowing liquid phase. The effect of channel geometry and wetting properties was investigated through two-phase water-air flow in square and round microchannels, having three static contact angles of 20, 80 and 105 degrees. Four different flow regimes are observed for a fixed flow rate, this being thought to be caused by the wetting behavior of liquid flowing in the corners as well as the liquid film stability. Film deposition experiments in wetting and non-wetting round microchannels show that a thicker film is deposited for wetting conditions departing from the ideal 0 degrees contact angle. A film thickness dependence with the contact angle theta as well as the Capillary number, in the form h_R ~ Ca^(2/3)/ cos(theta) is inferred from scaling arguments, for contact angles smaller than 36 degrees. Non-wetting film deposition experiments reveal that a film significantly thicker than the wetting Bretherton film is deposited. A hydraulic jump occurs if critical conditions are met, as given by a proposed nondimensional parameter similar to the Froude number. Film thickness correlations are also found by matching the measured and the proposed velocity derived in the shock theory. The surface wetting as well as the presence of the shock cause morphological changes in the Taylor bubble flow.

Method of Detecting System Function by Measuring Frequency Response

Methods of rapidly measuring an impedance spectrum of an energy storage device in-situ over a limited number of logarithmically distributed frequencies are described. An energy storage device is excited with a known input signal, and aresponse is measured to ascertain the impedance spectrum. An excitation signal is a limited time duration sum-of-sines consisting of a select number offrequencies. In one embodiment, magnitude and phase of each frequency ofinterest within the sum-of-sines is identified when the selected frequencies and sample rate are logarithmic integer steps greater than two. This technique requires a measurement with a duration of one period of the lowest frequency. In another embodiment, where selected frequencies are distributed in octave steps, the impedance spectrum can be determined using a captured time record that is reduced to a half-period of the lowest frequency.

Method of Detecting System Function by Measuring Frequency Response

Real time battery impedance spectrum is acquired using one time record, Compensated Synchronous Detection (CSD). This parallel method enables battery diagnostics. The excitation current to a test battery is a sum of equal amplitude sin waves of a few frequencies spread over range of interest. The time profile of this signal has duration that is a few periods of the lowest frequency. The voltage response of the battery, average deleted, is the impedance of the battery in the time domain. Since the excitation frequencies are known, synchronous detection processes the time record and each component, both magnitude and phase, is obtained. For compensation, the components, except the one of interest, are reassembled in the time domain. The resulting signal is subtracted from the original signal and the component of interest is synchronously detected. This process is repeated for each component.

Determination of menthol in plasma and urine of rats and humans by headspace solid phase microextraction and gas chromatography--mass spectrometry

A method for the determination of menthol and menthol glucuronide (M-G) after enzymatic hydrolysis in plasma and urine of rats and humans was developed using headspace solid phase microextraction and gas chromatography-mass spectrometry in the selected ion monitoring mode (HS-SPME/GC-MS). The assay linearity for plasma ranged from 5 to 1000 ng/ml. The limit of quantification (LOQ) in plasma was 5 ng/ml. The intra- and inter-day precision for menthol and M-G were < or = 18.1% R.S.D. at the LOQ and < or = 4.0% at higher concentrations. Menthol and M-G were determined in rat and human plasma and urine after administration of menthol.

Collision Detection: Broad Phase Adaptation from Multi-Core to Multi-GPU Architecture

We present in this paper several contributions on the collision detection optimization centered on hardware performance. We focus on the broad phase which is the first step of the collision detection process and propose three new ways of parallelization of the well-known Sweep and Prune algorithm. We first developed a multi-core model takes into account the number of available cores. Multi-core architecture enables us to distribute geometric computations with use of multi-threading. Critical writing section and threads idling have been minimized by introducing new data structures for each thread. Programming with directives, like OpenMP, appears to be a good compromise for code portability. We then proposed a new GPU-based algorithm also based on the "Sweep and Prune" that has been adapted to multi-GPU architectures. Our technique is based on a spatial subdivision method used to distribute computations among GPUs. Results show that significant speed-up can be obtained by passing from 1 to 4 GPUs in a large-scale environment.

Net CO 2 surface emissions at Bern, Switzerland inferred from ambient observations of CO 2 , δ(O 2 /N 2 ), and 222 Rn using a customized radon tracer inversion

The 222Radon tracer method is a powerful tool to estimate local and regional surface emissions of, e.g., greenhouse gases. In this paper we demonstrate that in practice, the method as it is commonly used, produces inaccurate results in case of nonhomogeneously spread emission sources, and we propose a different approach to account for this. We have applied the new methodology to ambient observations of CO2 and 222Radon to estimate CO2 surface emissions for the city of Bern, Switzerland. Furthermore, by utilizing combined measurements of CO2 and δ(O2/N2) we obtain valuable information about the spatial and temporal variability of the main emission sources. Mean net CO2 emissions based on 2 years of observations are estimated at (11.2 ± 2.9) kt km−2 a−1. Oxidative ratios indicate a significant influence from the regional biosphere in summer/spring and fossil fuel combustion processes in winter/autumn. Our data indicate that the emissions from fossil fuels are, to a large degree, related to the combustion of natural gas which is used for heating purposes.

Precipitation over the past four centuries in the Dieshan Mountains as inferred from tree rings: An introduction to an HHT-based method

To improve our understanding of the Asian monsoon system, we developed a hydroclimate reconstruction in a marginal monsoon shoulder region for the period prior to the industrial era. Here, we present the first moisture sensitive tree-ring chronology, spanning 501 years for the Dieshan Mountain area, a boundary region of the Asian summer monsoon in the northeastern Tibetan Plateau. This reconstruction was derived from 101 cores of 68 old-growth Chinese pine (Pinus tabulaeformis) trees. We introduce a Hilbert–Huang Transform (HHT) based standardization method to develop the tree-ring chronology, which has the advantages of excluding non-climatic disturbances in individual tree-ring series. Based on the reliable portion of the chronology, we reconstructed the annual (prior July to current June) precipitation history since 1637 for the Dieshan Mountain area and were able to explain 41.3% of the variance. The extremely dry years in this reconstruction were also found in historical documents and are also associated with El Niño episodes. Dry periods were reconstructed for 1718–1725, 1766–1770 and 1920–1933, whereas 1782–1788 and 1979–1985 were wet periods. The spatial signatures of these events were supported by data from other marginal regions of the Asian summer monsoon. Over the past four centuries, out-of-phase relationships between hydroclimate variations in the Dieshan Mountain area and far western Mongolia were observed during the 1718–1725 and 1766–1770 dry periods and the 1979–1985 wet period.

Estimation of aortic pressure waveforms from 4D phase-contrast MRI

Several approaches for the non-invasive MRI-based measurement of the aortic pressure waveform over the heart cycle have been proposed in the last years. These methods are normally based on time-resolved, two-dimensional phase-contrast sequences with uni-directionally encoded velocities (2D PC-MRI). In contrast, three-dimensional acquisitions with tridirectional velocity encoding (4D PC-MRI) have been shown to be a suitable data source for detailed investigations of blood flow and spatial blood pressure maps. In order to avoid additional MR acquisitions, it would be advantageous if the aortic pressure waveform could also be computed from this particular form of MRI. Therefore, we propose an approach for the computation of the aortic pressure waveform which can be completely performed using 4D PC-MRI. After the application of a segmentation algorithm, the approach automatically computes the aortic pressure waveform without any manual steps. We show that our method agrees well with catheter measurements in an experimental phantom setup and produces physiologically realistic results in three healthy volunteers.

Electrochemical simulation of phase I metabolism for 21 drugs using four different working electrodes in an automated screening setup with MS detection.

BACKGROUND Electrochemical conversion of xenobiotics has been shown to mimic human phase I metabolism for a few compounds. MATERIALS & METHODS Twenty-one compounds were analyzed with a semiautomated electrochemical setup and mass spectrometry detection. RESULTS The system was able to mimic some metabolic pathways, such as oxygen gain, dealkylation and deiodination, but many of the expected and known metabolites were not produced. CONCLUSION Electrochemical conversion is a useful approach for the preparative synthesis of some types of metabolites, but as a screening method for unknown phase I metabolites, the method is, in our opinion, inferior to incubation with human liver microsomes and in vivo experiments with laboratory animals, for example.

Combining Monte Carlo methods with coherent wave optics for the simulation of phase-sensitive X-ray imaging

Phase-sensitive X-ray imaging shows a high sensitivity towards electron density variations, making it well suited for imaging of soft tissue matter. However, there are still open questions about the details of the image formation process. Here, a framework for numerical simulations of phase-sensitive X-ray imaging is presented, which takes both particle- and wave-like properties of X-rays into consideration. A split approach is presented where we combine a Monte Carlo method (MC) based sample part with a wave optics simulation based propagation part, leading to a framework that takes both particle- and wave-like properties into account. The framework can be adapted to different phase-sensitive imaging methods and has been validated through comparisons with experiments for grating interferometry and propagation-based imaging. The validation of the framework shows that the combination of wave optics and MC has been successfully implemented and yields good agreement between measurements and simulations. This demonstrates that the physical processes relevant for developing a deeper understanding of scattering in the context of phase-sensitive imaging are modelled in a sufficiently accurate manner. The framework can be used for the simulation of phase-sensitive X-ray imaging, for instance for the simulation of grating interferometry or propagation-based imaging.

Multimodal treatment for high-risk prostate cancer with high-dose intensity-modulated radiation therapy preceded or not by radical prostatectomy, concurrent intensified-dose docetaxel and long-term androgen deprivation therapy: results of a prospective phase II trial.

BACKGROUND The optimal management of high-risk prostate cancer remains uncertain. In this study we assessed the safety and efficacy of a novel multimodal treatment paradigm for high-risk prostate cancer. METHODS This was a prospective phase II trial including 35 patients with newly diagnosed high-risk localized or locally advanced prostate cancer treated with high-dose intensity-modulated radiation therapy preceded or not by radical prostatectomy, concurrent intensified-dose docetaxel-based chemotherapy and long-term androgen deprivation therapy. Primary endpoint was acute and late toxicity evaluated with the Common Terminology Criteria for Adverse Events version 3.0. Secondary endpoint was biochemical and clinical recurrence-free survival explored with the Kaplan-Meier method. RESULTS Acute gastro-intestinal and genito-urinary toxicity was grade 2 in 23% and 20% of patients, and grade 3 in 9% and 3% of patients, respectively. Acute blood/bone marrow toxicity was grade 2 in 20% of patients. No acute grade ≥ 4 toxicity was observed. Late gastro-intestinal and genito-urinary toxicity was grade 2 in 9% of patients each. No late grade ≥ 3 toxicity was observed. Median follow-up was 63 months (interquartile range 31-79). Actuarial 5-year biochemical and clinical recurrence-free survival rate was 55% (95% confidence interval, 35-75%) and 70% (95% confidence interval, 52-88%), respectively. CONCLUSIONS In our phase II trial testing a novel multimodal treatment paradigm for high-risk prostate cancer, toxicity was acceptably low and mid-term oncological outcome was good. This treatment paradigm, thus, may warrant further evaluation in phase III randomized trials.

A quantitiative LC-MS/MS method for the measurement of arachidonic acid, prostanoids, endocannabinoids, N-acylethanolamines and steroids in human plasma.

Free arachidonic acid is functionally interlinked with different lipid signaling networks including those involving prostanoid pathways, the endocannabinoid system, N-acylethanolamines, as well as steroids. A sensitive and specific LC-MS/MS method for the quantification of arachidonic acid, prostaglandin E2, thromboxane B2, anandamide, 2-arachidonoylglycerol, noladin ether, lineoyl ethanolamide, oleoyl ethanolamide, palmitoyl ethanolamide, steroyl ethanolamide, aldosterone, cortisol, dehydroepiandrosterone, progesterone, and testosterone in human plasma was developed and validated. Analytes were extracted using acetonitrile precipitation followed by solid phase extraction. Separations were performed by UFLC using a C18 column and analyzed on a triple quadrupole MS with electron spray ionization. Analytes were run first in negative mode and, subsequently, in positive mode in two independent LC-MS/MS runs. For each analyte, two MRM transitions were collected in order to confirm identity. All analytes showed good linearity over the investigated concentration range (r>0.98). Validated LLOQs ranged from 0.1 to 190ng/mL and LODs ranged from 0.04 to 12.3ng/mL. Our data show that this LC-MS/MS method is suitable for the quantification of a diverse set of bioactive lipids in plasma from human donors (n=32). The determined plasma levels are in agreement with the literature, thus providing a versatile method to explore pathophysiological processes in which changes of these lipids are implicated.

Comparison of (31)P saturation and inversion magnetization transfer in human liver and skeletal muscle using a clinical MR system and surface coils.

(31)P MRS magnetization transfer ((31)P-MT) experiments allow the estimation of exchange rates of biochemical reactions, such as the creatine kinase equilibrium and adenosine triphosphate (ATP) synthesis. Although various (31)P-MT methods have been successfully used on isolated organs or animals, their application on humans in clinical scanners poses specific challenges. This study compared two major (31)P-MT methods on a clinical MR system using heteronuclear surface coils. Although saturation transfer (ST) is the most commonly used (31)P-MT method, sequences such as inversion transfer (IT) with short pulses might be better suited for the specific hardware and software limitations of a clinical scanner. In addition, small NMR-undetectable metabolite pools can transfer MT to NMR-visible pools during long saturation pulses, which is prevented with short pulses. (31)P-MT sequences were adapted for limited pulse length, for heteronuclear transmit-receive surface coils with inhomogeneous B1 , for the need for volume selection and for the inherently low signal-to-noise ratio (SNR) on a clinical 3-T MR system. The ST and IT sequences were applied to skeletal muscle and liver in 10 healthy volunteers. Monte-Carlo simulations were used to evaluate the behavior of the IT measurements with increasing imperfections. In skeletal muscle of the thigh, ATP synthesis resulted in forward reaction constants (k) of 0.074 ± 0.022 s(-1) (ST) and 0.137 ± 0.042 s(-1) (IT), whereas the creatine kinase reaction yielded 0.459 ± 0.089 s(-1) (IT). In the liver, ATP synthesis resulted in k = 0.267 ± 0.106 s(-1) (ST), whereas the IT experiment yielded no consistent results. ST results were close to literature values; however, the IT results were either much larger than the corresponding ST values and/or were widely scattered. To summarize, ST and IT experiments can both be implemented on a clinical body scanner with heteronuclear transmit-receive surface coils; however, ST results are much more robust against experimental imperfections than the current implementation of IT.

Ramucirumab versus placebo as second-line treatment in patients with advanced hepatocellular carcinoma following first-line therapy with sorafenib (REACH): a randomised, double-blind, multicentre, phase 3 trial

BACKGROUND VEGF and VEGF receptor-2-mediated angiogenesis contribute to hepatocellular carcinoma pathogenesis. Ramucirumab is a recombinant IgG1 monoclonal antibody and VEGF receptor-2 antagonist. We aimed to assess the safety and efficacy of ramucirumab in advanced hepatocellular carcinoma following first-line therapy with sorafenib. METHODS In this randomised, placebo-controlled, double-blind, multicentre, phase 3 trial (REACH), patients were enrolled from 154 centres in 27 countries. Eligible patients were aged 18 years or older, had hepatocellular carcinoma with Barcelona Clinic Liver Cancer stage C disease or stage B disease that was refractory or not amenable to locoregional therapy, had Child-Pugh A liver disease, an Eastern Cooperative Oncology Group performance status of 0 or 1, had previously received sorafenib (stopped because of progression or intolerance), and had adequate haematological and biochemical parameters. Patients were randomly assigned (1:1) to receive intravenous ramucirumab (8 mg/kg) or placebo every 2 weeks, plus best supportive care, until disease progression, unacceptable toxicity, or death. Randomisation was stratified by geographic region and cause of liver disease with a stratified permuted block method. Patients, medical staff, investigators, and the funder were masked to treatment assignment. The primary endpoint was overall survival in the intention-to-treat population. This study is registered with ClinicalTrials.gov, number NCT01140347. FINDINGS Between Nov 4, 2010, and April 18, 2013, 565 patients were enrolled, of whom 283 were assigned to ramucirumab and 282 were assigned to placebo. Median overall survival for the ramucirumab group was 9·2 months (95% CI 8·0-10·6) versus 7·6 months (6·0-9·3) for the placebo group (HR 0·87 [95% CI 0·72-1·05]; p=0·14). Grade 3 or greater adverse events occurring in 5% or more of patients in either treatment group were ascites (13 [5%] of 277 patients treated with ramucirumab vs 11 [4%] of 276 patients treated with placebo), hypertension (34 [12%] vs ten [4%]), asthenia (14 [5%] vs five [2%]), malignant neoplasm progression (18 [6%] vs 11 [4%]), increased aspartate aminotransferase concentration (15 [5%] vs 23 [8%]), thrombocytopenia (13 [5%] vs one [<1%]), hyperbilirubinaemia (three [1%] vs 13 [5%]), and increased blood bilirubin (five [2%] vs 14 [5%]). The most frequently reported (≥1%) treatment-emergent serious adverse event of any grade or grade 3 or more was malignant neoplasm progression. INTERPRETATION Second-line treatment with ramucirumab did not significantly improve survival over placebo in patients with advanced hepatocellular carcinoma. No new safety signals were noted in eligible patients and the safety profile is manageable. FUNDING Eli Lilly and Co.