Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional regression”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-J-59: Assessment of Dose Response Distribution in Individual Human Tumor

Purpose: To fulfill precision radiotherapy via adaptive dose painting by number, voxel-by-voxel dose response or radio-sensitivity in individual human tumor needs to be determined in early treatment to guide treatment adaptation. In this study, multiple FDG PET images obtained pre- and weekly during the treatment course were utilized to determine the distribution/spectrum of dose response parameters in individual human tumors. Methods: FDG PET/CT images of 18 HN cancer patients were used in the study. Spatial parametric image of tumor metabolic ratio (dSUV) was created following voxel by voxel deformable image registration. Each voxel value in dSUV was a function of pre-treatment baseline SUV and treatment delivered dose, and used as a surrogate of tumor survival fraction (SF). Regression fitting with break points was performed using the LQ-model with tumor proliferation for the control and failure group of tumors separately. The distribution and spectrum of radiation sensitivity and growth in individual tumors were determined and evaluated. Results: Spectrum of tumor dose-sensitivity and proliferation in the controlled group was broad with α in tumor survival LQ-model from 0.17 to 0.8. It was proportional to the baseline SUV. Tlag was about 21∼25 days, and Tpot about 0.56∼1.67 days respectively. Commonly tumor voxels with high radio-sensitivity or larger α had small Tlag and Tpot. For the failure group, the radio-sensitivity α was low within 0.05 to 0.3, but did not show clear Tlag. In addition, tumor voxel radio-sensitivity could be estimated during the early treatment weeks. Conclusion: Dose response distribution with respect to radio-sensitivity and growth in individual human tumor can be determined using FDG PET imaging based tumor metabolic ratio measured in early treatment course. The discover is critical and provides a potential quantitative objective to implement tumor specific precision radiotherapy via adaptive dose painting by number.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-DE-201-07: Measurement of Real-Time Dose for Tandem and Ovoid Brachytherapy Procedures Using a High Precision Optical Fiber Radiation Detector

Purpose: Development of a novel on-line dosimetry tool is needed to move toward patient-specific quality assurance measurements for Ir-192 HDR brachytherapy to verify accurate dose delivery to the intended location. This work describes the development and use of a nano-crystalline yttrium oxide inorganic scintillator based optical-fiber detector capable of acquiring real-time high-precision dose measurements during tandem and ovoid (T&O) gynecological (GYN) applicator Ir-192 HDR brachytherapy procedures. Methods: An optical-fiber detector was calibrated by acquiring light output measurements in liquid water at 3, 5, 7, and 9cm radial source-detector-distances from an Ir-192 HDR source. A regression model was fit to the data to describe the relative light output per unit dose (TG-43 derived) as a function of source-detector-distance. Next, the optical-fiber detector was attached to a vaginal balloon fixed to a Varian Fletcher-Suit-Delclos-style applicator (to mimic clinical setup), and localized by acquiring high-resolution computed tomography (CT) images. To compare the physical point dose to the TPS calculated values (TG-43 and Acuros-BV), a phantom measurement was performed, by submerging the T&O applicator in a liquid water bath and delivering a treatment template representative of a clinical T&O procedure. The fiber detector collected scintillation signal as a function of time, and the calibration data was applied to calculate both real-time dose rate, and cumulative dose. Results: Fiber cumulative dose values were 100.0cGy, 94.3cGy, and 348.9cGy from the tandem, left ovoid, and right ovoid dwells, respectively (total of 443.2cGy). A plot of real time dose rate during the treatment was also acquired. The TPS values at the fiber location were 458.4cGy using TG-43, and 437.6cGy using Acuros-BV calculated as Dm,m (per TG-186). Conclusion: The fiber measured dose value agreement was 3% vs TG-43 and −1% vs Acuros-BV. This fiber detector opens up new possibilities for performing patient-specific quality assurance for Ir-192 HDR GYN procedures. Funding from Coulter Foundation, Duke Bio-medical Engineering. Company is being created around the detector technology. Duke holds patents on the technology.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1982

Determination of airplane aerodynamic parameters from flight data at high angles of attack

The problem of determining airplane aerodynamic model equations and estimating the associated parameters from flight data taken at high angles of attack is addressed. Two representations of the aerodynamic function based on the polynomial and spline representations are given. Then the technique of building an adequate model using a stepwise regression is presented with examples demonstrating the construction of the model and various approaches to model verification.

Klein, V.↗
NASA NTRS1999

"Novel Techniques in Non-Stationary Analysis of Rotorcraft Vibration Signitures"

This research effort produced new methods to analyze the performance of linear predictors that track non-stationary processes. Specifically, prediction methods have been applied to the vibration pattern of rotorcraft drivetrains. This analysis is part or a larger rotorcraft Health and Usage Monitoring System (HUMS) that can diagnose immediate failures of the subsystems, as indicated by abrupt change in the vibration signature, and prognosticate future health, by examining the vibration patterns against long-term trends. This problem is described by a earlier joint paper co-authored by members of the funding agency and the recipient institutions prior to this grant effort. Specific accomplishments under this grant include the following: (1) Definition of a framework for analysis of non-stationary time-series estimation using the coefficients of an adaptive filter. (2) Description of a novel method of combining short-term predictor error and long-term regression error to analyze the performance of a non-stationary predictor. (3) Formulation of a multi-variate probability density function that quantifies the performance of a adaptive predictor by using the short- and long-term error variables in a Gamma function distribution. and (4) Validation of the mathematical formulations with empirical data from NASA flight tests and simulated data to illustrate the utility beyond the domain of vibrating machinery.

Meng, Teresa↗
NASA NTRS2003

Pulsed Phase Lock Loop Device for Monitoring Intracranial Pressure During Space Flight

We have developed an ultrasonic device to monitor ICP waveforms non-invasively from cranial diameter oscillations using a NASA-developed pulsed phase lock loop (PPLL) technique. The purpose of this study was to attempt to validate the PPLL device for reliable recordings of ICP waveforms and analysis of ICP dynamics in vivo. METHODS: PPLL outputs were recorded in patients during invasive ICP monitoring at UCSD Medical Center (n=10). RESULTS: An averaged linear regression coefficient between ICP and PPLL waveform data during one cardiac cycle in all patients is 0.88 +/- 0.02 (mean +/- SE). Coherence function analysis indicated that ICP and PPLL waveforms have high correlation in the lst, 2nd, and 3rd harmonic waves associated with a cardiac cycle. CONCLUSIONS: PPLL outputs represent ICP waveforms in both frequency and time domains. PPLL technology enables in vivo evaluation of ICP dynamics non-invasively, and can acquire continuous ICP waveforms during spaceflight because of compactness and non-invasive nature.

Ueno, Toshiaki↗
NASA NTRS2004

NASA Experimental Program to Stimulate Competitive Research: South Carolina

The use of an appropriate relationship model is critical for reliable prediction of future urban growth. Identification of proper variables and mathematic functions and determination of the weights or coefficients are the key tasks for building such a model. Although the conventional logistic regression model is appropriate for handing land use problems, it appears insufficient to address the issue of interdependency of the predictor variables. This study used an alternative approach to simulation and modeling urban growth using artificial neural networks. It developed an operational neural network model trained using a robust backpropagation method. The model was applied in the Myrtle Beach region of South Carolina, and tested with both global datasets and areal datasets to examine the strength of both regional models and areal models. The results indicate that the neural network model not only has many theoretic advantages over other conventional mathematic models in representing the complex urban systems, but also is practically superior to the logistic model in its capability to predict urban growth with better - accuracy and less variation. The neural network model is particularly effective in terms of successfully identifying urban patterns in the rural areas where the logistic model often falls short. It was also found from the area-based tests that there are significant intra-regional differentiations in urban growth with different rules and rates. This suggests that the global modeling approach, or one model for the entire region, may not be adequate for simulation of a urban growth at the regional scale. Future research should develop methods for identification and subdivision of these areas and use a set of area-based models to address the issues of multi-centered, intra- regionally differentiated urban growth.

Sutton, Michael A.↗
NASA NTRS1997

Assessment of a remote sensing-based model for predicting malaria transmission risk in villages of Chiapas, Mexico

A blind test of two remote sensing-based models for predicting adult populations of Anopheles albimanus in villages, an indicator of malaria transmission risk, was conducted in southern Chiapas, Mexico. One model was developed using a discriminant analysis approach, while the other was based on regression analysis. The models were developed in 1992 for an area around Tapachula, Chiapas, using Landsat Thematic Mapper (TM) satellite data and geographic information system functions. Using two remotely sensed landscape elements, the discriminant model was able to successfully distinguish between villages with high and low An. albimanus abundance with an overall accuracy of 90%. To test the predictive capability of the models, multitemporal TM data were used to generate a landscape map of the Huixtla area, northwest of Tapachula, where the models were used to predict risk for 40 villages. The resulting predictions were not disclosed until the end of the test. Independently, An. albimanus abundance data were collected in the 40 randomly selected villages for which the predictions had been made. These data were subsequently used to assess the models' accuracies. The discriminant model accurately predicted 79% of the high-abundance villages and 50% of the low-abundance villages, for an overall accuracy of 70%. The regression model correctly identified seven of the 10 villages with the highest mosquito abundance. This test demonstrated that remote sensing-based models generated for one area can be used successfully in another, comparable area.

Insect Vectors/growth & development↗
DOE OSTI2018

Prediction of Stable and Active (Oxy-Hydro) Oxide Nanoislands on Noble-Metal Supports for Electrochemical Oxygen Reduction Reaction

Developing cost-effective oxygen electrocatalysts with high activity and stability is key to their commercialization. However, economical earth-abundant catalysts based on first-row transition-metal oxides suffer from low electrochemical stability, which is difficult to improve without compromising their activity. In this paper, using density functional theory calculations, we demonstrate that noble-metal supports lead to bifunctional enhancement of both the stability and the oxygen reduction reaction (ORR) activity of metal (oxy-hydro) oxide nanoislands. We observe a significant stabilization of supported nanoislands beyond the intrinsic stability limits of bulk phases, which originates from a favorable lattice mismatch and reductive charge transfer from oxophilic supports. Additionally, we discover that interfacial active sites (located between the nanoisland and the support) reinforce the binding strength of reaction intermediates, hence boosting ORR activity. Considering that both stability and activity lead to discovery of CoOOH|Pt, NiOOH|Ag, and FeO 2 |Ag as viable systems for alkaline ORR, we then use a multivariant linear regression method to identify elementary descriptors for efficient screening of promising cost-effective nanoisland|support catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Robust Machine Learning Inference from X-ray Absorption Near Edge Spectra through Featurization

X-ray absorption spectroscopy (XAS) is a commonly employed technique for characterizing functional materials. In particular, X-ray absorption near edge spectra (XANES) encode local coordination and electronic information, and machine learning approaches to extract this information are of significant interest. To date, most ML approaches for XANES have primarily focused on using the raw spectral intensities as input, overlooking the potential benefits of incorporating spectral transformations and dimensionality reduction techniques into ML predictions. Here, in this work, we focused on systematically comparing the impact of different featurization methods on the performance of ML models for XAS analysis. We evaluated the classification and regression capabilities of these models on computed data sets and validated their performance on previously unseen experimental data sets. Our analysis revealed an intriguing discovery: the cumulative distribution function feature achieves both high prediction accuracy and exceptional transferability. This remarkably robust performance can be attributed to its tolerance to horizontal shifts in the spectra, which is crucial when validating models using experimental data. While this work exclusively focuses on XANES analysis, we anticipate that the methodology presented here will hold promise as a versatile asset to the broader spectroscopy community.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Dose-Escalation Study for Cardiac Radiosurgery in a Porcine Model

Purpose: To perform a proof-of-principle dose-escalation study to radiosurgically induce scarring in cardiac muscle tissue to block veno-atrial electrical connections at the pulmonary vein antrum, similar to catheter ablation. Methods and Materials: Nine mini-pigs underwent pretreatment magnetic resonance imaging (MRI) evaluation of heart function and electrophysiology assessment by catheter measurements in the right superior pulmonary vein (RSPV). Immediately after examination, radiosurgery with randomized single-fraction doses of 0 and 17.5-35 Gy in 2.5-Gy steps were delivered to the RSPV antrum (target volume 5-8 cm{sup 3}). MRI and electrophysiology were repeated 6 months after therapy, followed by histopathologic examination. Results: Transmural scarring of cardiac muscle tissue was noted with doses ≥32.5 Gy. However, complete circumferential scarring of the RSPV was not achieved. Logistic regressions showed that extent and intensity of fibrosis significantly increased with dose. The 50% effective dose for intense fibrosis was 31.3 Gy (odds ratio 2.47/Gy, P<.01). Heart function was not affected, as verified by MRI and electrocardiogram evaluation. Adjacent critical structures were not damaged, as verified by pathology, demonstrating the short-term safety of small-volume cardiac radiosurgery with doses up to 35 Gy. Conclusions: Radiosurgery with doses >32.5 Gy in the healthy pig heart can induce circumscribed scars at the RSPV antrum noninvasively, mimicking the effect of catheter ablation. In our study we established a significant dose-response relationship for cardiac radiosurgery. The long-term effects and toxicity of such high radiation doses need further investigation in the pursuit of cardiac radiosurgery for noninvasive treatment of atrial fibrillation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Single toxin dose-response models revisited

The goal of this paper is to offer a rigorous analysis of the sigmoid shape single toxin dose-response relationship. The toxin efficacy function is introduced and four special points, including maximum toxin efficacy and inflection points, on the dose-response curve are defined. The special points define three phases of the toxin effect on mortality: (1) toxin concentrations smaller than the first inflection point or (2) larger then the second inflection point imply low mortality rate, and (3) concentrations between the first and the second inflection points imply high mortality rate. Probabilistic interpretation and mathematical analysis for each of the four models, Hill, logit, probit, and Weibull is provided. Two general model extensions are introduced: (1) the multi-target hit model that accounts for the existence of several vital receptors affected by the toxin, and (2) model with a nonzero mortality at zero concentration to account for natural mortality. Special attention is given to statistical estimation in the framework of the generalized linear model with the binomial dependent variable as the mortality count in each experiment, contrary to the widespread nonlinear regression treating the mortality rate as continuous variable. The models are illustrated using standard EPA Daphnia acute (48 h) toxicity tests with mortality as a function of NiCl or CuSO{sub 4} toxin. - Highlights: • The paper offers a rigorous study of a sigmoid dose-response relationship. • The concentration with highest mortality rate is rigorously defined. • A table with four special points for five morality curves is presented. • Two new sigmoid dose-response models have been introduced. • The generalized linear model is advocated for estimation of sigmoid dose-response relationship.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1985

Automatic rendezvous and docking systems functional and performance requirements

A generalized mission design scheme which utilizes a standard mission profile for all OMV rendezvous operations, recognizes typical operational constraints, and minimizes propellant penalties due to nodal regression effects was developed. This scheme has been used to demonstrate a unified guidance and navigation maneuver processor (the UMP), which supports all mission phases through station-keeping. The initial demonstration version of the Orbital Rendezvous Mission Planner (ORMP) was provided for evaluation purposes, and program operation was discussed.

Source record↗
DOE OSTI2022

ASP (Adaptive Splines for Prediction) [SWR-22-66]

This software package provides functionality to fit piecewise cubic splines to model the relationship between a univariate X (independent variable) and a univariate Y (dependent variable). The output is akin to linear regression but has higher "capacity" in that it can learn more complicated relationships despite leveraging only a single predictor. It thus can learn very reasonable relationships between X and Y with minimal effort. The canonical intended use case is for predicting electric demand (Y) from temperature (X). There can be complicated, non-linear relationships between these two variables, but these can be reasonably approximated by the sum of: 1) an "aggregate" relationship between daily average temperature and daily average demand, and 2) intra-day patterns represented as hourly offsets from the daily average load. The spline functionality is used for learning a robust "aggregate" relationship without the inconvenience of manual feature engineering while still ensuring robust results. The model fitting workhorse is R's built-in smooth.spline. This package provides a suite of options for controlling how smooth.spline fits a spline to the data. In particular, it is designed to return a model with a single critical point (i.e., a single location along the domain/support where the first derivative is zero). This is to ensure that predicted changes in electric demand are always positive as temperatures become more extreme. This helps to prevent overfitting, particularly when the input data has few observations or has significant uncertainty from other sources.

Murphy, Sinnott↗
NASA NTRS2013

Estimated Probability of a Cervical Spine Injury During an ISS Mission

Introduction: The Integrated Medical Model (IMM) utilizes historical data, cohort data, and external simulations as input factors to provide estimates of crew health, resource utilization and mission outcomes. The Cervical Spine Injury Module (CSIM) is an external simulation designed to provide the IMM with parameter estimates for 1) a probability distribution function (PDF) of the incidence rate, 2) the mean incidence rate, and 3) the standard deviation associated with the mean resulting from injury/trauma of the neck. Methods: An injury mechanism based on an idealized low-velocity blunt impact to the superior posterior thorax of an ISS crewmember was used as the simulated mission environment. As a result of this impact, the cervical spine is inertially loaded from the mass of the head producing an extension-flexion motion deforming the soft tissues of the neck. A multibody biomechanical model was developed to estimate the kinematic and dynamic response of the head-neck system from a prescribed acceleration profile. Logistic regression was performed on a dataset containing AIS1 soft tissue neck injuries from rear-end automobile collisions with published Neck Injury Criterion values producing an injury transfer function (ITF). An injury event scenario (IES) was constructed such that crew 1 is moving through a primary or standard translation path transferring large volume equipment impacting stationary crew 2. The incidence rate for this IES was estimated from in-flight data and used to calculate the probability of occurrence. The uncertainty in the model input factors were estimated from representative datasets and expressed in terms of probability distributions. A Monte Carlo Method utilizing simple random sampling was employed to propagate both aleatory and epistemic uncertain factors. Scatterplots and partial correlation coefficients (PCC) were generated to determine input factor sensitivity. CSIM was developed in the SimMechanics/Simulink environment with a Monte Carlo wrapper (MATLAB) used to integrate the components of the module. Results: The probability of generating an AIS1 soft tissue neck injury from the extension/flexion motion induced by a low-velocity blunt impact to the superior posterior thorax was fitted with a lognormal PDF with mean 0.26409, standard deviation 0.11353, standard error of mean 0.00114, and 95% confidence interval [0.26186, 0.26631]. Combining the probability of an AIS1 injury with the probability of IES occurrence was fitted with a Johnson SI PDF with mean 0.02772, standard deviation 0.02012, standard error of mean 0.00020, and 95% confidence interval [0.02733, 0.02812]. The input factor sensitivity analysis in descending order was IES incidence rate, ITF regression coefficient 1, impactor initial velocity, ITF regression coefficient 2, and all others (equipment mass, crew 1 body mass, crew 2 body mass) insignificant. Verification and Validation (V&V): The IMM V&V, based upon NASA STD 7009, was implemented which included an assessment of the data sets used to build CSIM. The documentation maintained includes source code comments and a technical report. The software code and documentation is under Subversion configuration management. Kinematic validation was performed by comparing the biomechanical model output to established corridors.

Brooker, John E.↗