Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model fitting”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Automatic Segmentation of the Eye in 3D Magnetic Resonance Imaging: A Novel Statistical Shape Model for Treatment Planning of Retinoblastoma

Purpose: Proper delineation of ocular anatomy in 3-dimensional (3D) imaging is a big challenge, particularly when developing treatment plans for ocular diseases. Magnetic resonance imaging (MRI) is presently used in clinical practice for diagnosis confirmation and treatment planning for treatment of retinoblastoma in infants, where it serves as a source of information, complementary to the fundus or ultrasonographic imaging. Here we present a framework to fully automatically segment the eye anatomy for MRI based on 3D active shape models (ASM), and we validate the results and present a proof of concept to automatically segment pathological eyes. Methods and Materials: Manual and automatic segmentation were performed in 24 images of healthy children's eyes (3.29 ± 2.15 years of age). Imaging was performed using a 3-T MRI scanner. The ASM consists of the lens, the vitreous humor, the sclera, and the cornea. The model was fitted by first automatically detecting the position of the eye center, the lens, and the optic nerve, and then aligning the model and fitting it to the patient. We validated our segmentation method by using a leave-one-out cross-validation. The segmentation results were evaluated by measuring the overlap, using the Dice similarity coefficient (DSC) and the mean distance error. Results: We obtained a DSC of 94.90 ± 2.12% for the sclera and the cornea, 94.72 ± 1.89% for the vitreous humor, and 85.16 ± 4.91% for the lens. The mean distance error was 0.26 ± 0.09 mm. The entire process took 14 seconds on average per eye. Conclusion: We provide a reliable and accurate tool that enables clinicians to automatically segment the sclera, the cornea, the vitreous humor, and the lens, using MRI. We additionally present a proof of concept for fully automatically segmenting eye pathology. This tool reduces the time needed for eye shape delineation and thus can help clinicians when planning eye treatment and confirming the extent of the tumor.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Foreground bias from parametric models of far-IR dust emission

We use simple toy models of far-IR dust emission to estimate the accuracy to which the polarization of the cosmic microwave background can be recovered using multi-frequency fits, if the parametric form chosen for the fitted dust model differs from the actual dust emission. Commonly used approximations to the far-IR dust spectrum yield CMB residuals comparable to or larger than the sensitivities expected for the next generation of CMB missions, despite fitting the combined CMB + foreground emission to precision 0.1% or better. The Rayleigh–Jeans approximation to the dust spectrum biases the fitted dust spectral index by Δβ{sub d}=0.2 and the inflationary B-mode amplitude by Δr=0.03. Fitting the dust to a modified blackbody at a single temperature biases the best-fit CMB by Δr>0.003 if the true dust spectrum contains multiple temperature components. A 13-parameter model fitting two temperature components reduces this bias by an order of magnitude if the true dust spectrum is in fact a simple superposition of emission at different temperatures, but fails at the level Δr=0.006 for dust whose spectral index varies with frequency. Restricting the observing frequencies to a narrow region near the foreground minimum reduces these biases for some dust spectra but can increase the bias for others. Data at THz frequencies surrounding the peak of the dust emission can mitigate these biases while providing a direct determination of the dust temperature profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

GLEAM: Galaxy Line Emission & Absorption Modeling

We present Galaxy Line Emission & Absorption Modeling (gleam), a Python tool for fitting Gaussian models to emission and absorption lines in large samples of 1D extragalactic spectra. gleam is tailored to work well in batch mode without much human interaction. With gleam, users can uniformly process a variety of spectra, including galaxies and active galactic nuclei, in a wide range of instrument setups and signal-to-noise regimes. gleam also takes advantage of multiprocessing capabilities to process spectra in parallel. With the goal of enabling reproducible workflows for its users, gleam employs a small number of input files, including a central, user-friendly configuration in which fitting constraints can be defined for groups of spectra and overrides can be specified for edge cases. For each spectrum, gleam produces a table containing measurements and error bars for the detected spectral lines and continuum and upper limits for nondetections. For visual inspection and publishing, gleam can also produce plots of the data with fitted lines overlaid. In the present paper, we describe gleam’s main features, the necessary inputs, expected outputs, and some example applications, including thorough tests on a large sample of optical/infrared multi-object spectroscopic observations and integral field spectroscopic data. gleam is developed as an open-source project hosted at https://github.com/multiwavelength/gleam and welcomes community contributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Free kick instead of cross-validation in maximum-likelihood refinement of macromolecular crystal structures

The maximum-likelihood free-kick target, which calculates model error estimates from the work set and a randomly displaced model, proved superior in the accuracy and consistency of refinement of crystal structures compared with the maximum-likelihood cross-validation target, which calculates error estimates from the test set and the unperturbed model. The refinement of a molecular model is a computational procedure by which the atomic model is fitted to the diffraction data. The commonly used target in the refinement of macromolecular structures is the maximum-likelihood (ML) function, which relies on the assessment of model errors. The current ML functions rely on cross-validation. They utilize phase-error estimates that are calculated from a small fraction of diffraction data, called the test set, that are not used to fit the model. An approach has been developed that uses the work set to calculate the phase-error estimates in the ML refinement from simulating the model errors via the random displacement of atomic coordinates. It is called ML free-kick refinement as it uses the ML formulation of the target function and is based on the idea of freeing the model from the model bias imposed by the chemical energy restraints used in refinement. This approach for the calculation of error estimates is superior to the cross-validation approach: it reduces the phase error and increases the accuracy of molecular models, is more robust, provides clearer maps and may use a smaller portion of data for the test set for the calculation of R{sub free} or may leave it out completely.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2016

Foreground Bias from Parametric Models of Far-IR Dust Emission

We use simple toy models of far-IR dust emission to estimate the accuracy to which the polarization of the cosmic microwave background can be recovered using multi-frequency fits, if the parametric form chosen for the fitted dust model differs from the actual dust emission. Commonly used approximations to the far-IR dust spectrum yield CMB residuals comparable to or larger than the sensitivities expected for the next generation of CMB missions, despite fitting the combined CMB plus foreground emission to precision 0.1 percent or better. The Rayleigh-Jeans approximation to the dust spectrum biases the fitted dust spectral index by (Delta)(Beta)(sub d) = 0.2 and the inflationary B-mode amplitude by (Delta)(r) = 0.03. Fitting the dust to a modified blackbody at a single temperature biases the best-fit CMB by (Delta)(r) greater than 0.003 if the true dust spectrum contains multiple temperature components. A 13-parameter model fitting two temperature components reduces this bias by an order of magnitude if the true dust spectrum is in fact a simple superposition of emission at different temperatures, but fails at the level (Delta)(r) = 0.006 for dust whose spectral index varies with frequency. Restricting the observing frequencies to a narrow region near the foreground minimum reduces these biases for some dust spectra but can increase the bias for others. Data at THz frequencies surrounding the peak of the dust emission can mitigate these biases while providing a direct determination of the dust temperature profile.

cosmic microwave background↗
DOE OSTI2019

Kinetic Modeling of the Silicothermic Reduction of Manganese Oxide from Slag

A model was developed to describe the mixed rate control kinetics in systems with MnO-SiO{sub 2}-CaO-Al{sub 2}O{sub 3} slags and Fe-Si metal droplets. During the reaction of Fe-Si droplets and slag, Mn{sup 2+} transport was found to be part of controlling the system but could not sufficiently describe all initial conditions (chiefly, changes to initial silicon content). The current model describes the kinetics of the stated system and offers answers to the question of rate control; the model has been fitted to nine datasets of varying initial conditions including initial [Si], initial (MnO), initial droplet size, and reaction temperature. The fitted mass transfer coefficients for metal and slag were 2.3 × 10{sup −4} and 6.7 × 10{sup −4} m/s, respectively; these values are constant across the nine datasets. Previous claims about the efficiency of silicon usage in reducing manganese have been modified; it appears that the formation of silica is favored throughout the reaction, but that the formation of a silicon monoxide gas layer on the metal surface dramatically impacts the rate of reaction. As a measure of overall fit of the model, the average of the root-mean-square errors for all datasets is 14 pct. Mass transport in slag is twice as influential on rate control as the metal phase. The simultaneous transport of both Mn{sup 2+} and silicate dimers controls mass transport in the slag. Both [Si] and [Mn] can control from the metal side but are dependent on the initial conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A thermodynamic model to predict the aqueous solubility of hydrocarbon mixtures at two-phase hydrate-liquid water equilibrium

In order to understand the fate and transport of hydrocarbons and hydrocarbon mixtures in the deep sea and underlying sediments it requires an accurate determination of the two-phase hydrate (H)-liquid water (L w ) thermodynamic equilibrium in the absence of a free gas phase. In addition to controlling hydrate formation directly from the aqueous phase, the H-L w equilibrium provides the aqueous solubility of dissolving hydrate. The aqueous solubility of hydrocarbon mixture at H-L w equilibrium was calculated based on the van der Waals and Platteeuw statistical thermodynamic model and the Holder model. Langmuir constants were calculated from cell potential parameters obtained from ab initio intermolecular potentials and thus the model contains no fitting parameters. The model accurately predicts pure methane, ethane, and propane experimental H-L w equilibrium solubility data, including new data provided here for methane. Here, we present hydrocarbon solubilities in water and seawater at H-L w equilibrium at deep sea conditions for both methane and a thermogenic mixture based on the methane–ethane–propane ratio from the 2010 Macondo oil spill. We also present model predictions of the proportion of occupied large and small cages at H-L w equilibrium for both systems, as well as hydrocarbon ratios in the hydrate phase for the thermogenic mixture.

58 GEOSCIENCES↗
DOE OSTI2014

Random Forests to Predict Rectal Toxicity Following Prostate Cancer Radiation Therapy

Purpose: To propose a random forest normal tissue complication probability (RF-NTCP) model to predict late rectal toxicity following prostate cancer radiation therapy, and to compare its performance to that of classic NTCP models. Methods and Materials: Clinical data and dose-volume histograms (DVH) were collected from 261 patients who received 3-dimensional conformal radiation therapy for prostate cancer with at least 5 years of follow-up. The series was split 1000 times into training and validation cohorts. A RF was trained to predict the risk of 5-year overall rectal toxicity and bleeding. Parameters of the Lyman-Kutcher-Burman (LKB) model were identified and a logistic regression model was fit. The performance of all the models was assessed by computing the area under the receiving operating characteristic curve (AUC). Results: The 5-year grade ≥2 overall rectal toxicity and grade ≥1 and grade ≥2 rectal bleeding rates were 16%, 25%, and 10%, respectively. Predictive capabilities were obtained using the RF-NTCP model for all 3 toxicity endpoints, including both the training and validation cohorts. The age and use of anticoagulants were found to be predictors of rectal bleeding. The AUC for RF-NTCP ranged from 0.66 to 0.76, depending on the toxicity endpoint. The AUC values for the LKB-NTCP were statistically significantly inferior, ranging from 0.62 to 0.69. Conclusions: The RF-NTCP model may be a useful new tool in predicting late rectal toxicity, including variables other than DVH, and thus appears as a strong competitor to classic NTCP models.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Impact of hydrate saturation on water permeability in hydrate-bearing sediments

Permeability of the hydrate-bearing sediments vitally affects the hydrate dissociation process as well as the rate and efficiency of gas production. Reported permeability values are observed to be widely are very scattered owing to the dependence on multiple factors such as experimental conditions and test procedures. It is critical that the permeability is measured accurately to enable prediction of long-term gas production using numerical simulation for an economic development of hydrate-bearing reservoirs. In this study, the tetrahydrofuran (THF) hydrates that exhibit pore habits of the pore-filling pattern are formed in sediments, and the water permeability is quantified as a function of hydrate saturation. The results reveal that a decrease in the permeability for an increase in the hydrate saturation. Based on the measured values that are bounded by Kozeny grain-coating and pore-filling models, the fitting parameters of the empirical permeability models are indicated for use in numerical simulation of long-term hydrate dissociation. Wave velocity measurements reveal pore habits of pore-filling pattern of THF hydrates in the sediments with S h < 0.5, and a gradual transition to patchy and load-bearing pore-habits when S h > 0.6. The numerical simulation results obtained using the complementary pore-network model suggest that the THF hydrates in the sediments may form in small clusters with an average patch size of ~4 pores.

58 GEOSCIENCES↗
DOE OSTI2013

Towards rapid nanoscale measurement of strain in III-nitride heterostructures

We report the structural and compositional nanoscale characterization of InAlN/GaN nanoribbon-structured high electron mobility transistors (HEMTs) through the use of geometric phase analysis (GPA) and nanobeam electron diffraction (NBED). The strain distribution in the HEMT layer is quantified and compared to the expected strain profile for the nominal structure predicted by finite element analysis (FEA). Using the experimental strain results, the actual structure is determined and used to modify the FEA model. The improved fit of the model demonstrates that GPA and NBED provide a powerful platform for routine and rapid characterization of strain in III-V semiconducting device systems leading to insights into device evolution during processing and future device optimization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Comment: Spurious Correlation and Other Observations on Experimental Design for Engineering Dimensional Analysis

This article discusses the paper "Experimental Design for Engineering Dimensional Analysis" by Albrecht et al. (2013, Technometrics). That paper provides and overview of engineering dimensional analysis (DA) for use in developing DA models. The paper proposes methods for generating model-robust experimental designs to supporting fitting DA models. The specific approach is to develop a design that maximizes the efficiency of a specified empirical model (EM) in the original independent variables, subject to a minimum efficiency for a DA model expressed in terms of dimensionless groups (DGs). This discussion article raises several issues and makes recommendations regarding the proposed approach. Also, the concept of spurious correlation is raised and discussed. Spurious correlation results from the response DG being calculated using several independent variables that are also used to calculate predictor DGs in the DA model.

Dimensional analysis↗
DOE OSTI2017

Oxygen evolution from olivine Mn 1– x M x PO 4 ( M =Fe,Ni,Al,Mg) delithiated cathode materials

Olivine LiMnPO 4 is a promising cathode material for Li-ion batteries. One drawback of this material is the propensity of its delithiated phase, MnPO 4 , to evolve oxygen gas above approximately 200 °C. During thermal runaway of cells, this oxygen gas can burn the electrolyte and other cell components and thereby jeopardize safety. Partial substitution of Mn with M = Fe, Ni, Al, or Mg has been used to improve the lithium intercalation kinetics of Li x MnPO 4 ; however, the effect of these substitutions on oxygen evolution is not fully documented. In this paper, we calculate phase diagrams and oxygen evolution diagrams for these Mn 1–x M x PO 4 delithiated cathode materials. To generate the phase diagrams, we use subregular solid-solution models and fit the energetic parameters of these models to density functional theory calculations of special quasirandom structures. The resulting thermodynamic models describe the effect of mixing on the initial temperature of oxygen evolution and on the cumulative amount of oxygen evolution at elevated temperatures. We find that addition of Fe increases the initial temperature and decreases the cumulative amount of oxygen evolution. Mn 0.5 Fe 0.5 PO 4 exhibits an initial temperature 50 °C higher than MnPO 4 and releases 70% less oxygen gas at 300 °C. Al is insoluble in MnPO 4 , so addition of Al has no affect on the initial temperature. However, Al addition does slightly decrease the amount of oxygen evolution due to an inactive AlPO 4 component. Mg and Ni both decrease the initial temperature of oxygen evolution, and therefore may worsen the safety of MnPO 4 .

25 ENERGY STORAGE↗
DOE OSTI2014

HERMES impact for the access of Compton form factors

We utilize the DVCS asymmetry measurements of the HERMES collaboration for access to Compton form factors in the deeply virtual regime and to generalized parton distributions. In particular, the (almost) complete measurement of DVCS observables allows us to map various asymmetries into the space of Compton form factors, where we still rely in this analysis on dominance of twist-two associated Compton form factors. We compare this one-to-one map with local Compton form factor fits and a model dependent global fit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

FuelsGen

This toolset generates heterogeneous spatial layouts, aimed at representing forest fuels applicable to wildland fire modeling. Methods are included for calibrating (fitting) this model to data extracted from plot survey information.

Gattiker, James↗
DOE OSTI2021

Using heterogeneous data to identify signatures of dengue outbreaks at fine spatio-temporal scales across Brazil

Dengue virus remains a significant public health challenge in Brazil, and seasonal preparation efforts are hindered by variable intra- and interseasonal dynamics. Here, we present a framework for characterizing weekly dengue activity at the Brazilian mesoregion level from 2010–2016 as time series properties that are relevant to forecasting efforts, focusing on outbreak shape, seasonal timing, and pairwise correlations in magnitude and onset. In addition, we use a combination of 18 satellite remote sensing imagery, weather, clinical, mobility, and census data streams and regression methods to identify a parsimonious set of covariates that explain each time series property. The models explained 54% of the variation in outbreak shape, 38% of seasonal onset, 34% of pairwise correlation in outbreak timing, and 11% of pairwise correlation in outbreak magnitude. Regions that have experienced longer periods of drought sensitivity, as captured by the “normalized burn ratio,” experienced less intense outbreaks, while regions with regular fluctuations in relative humidity had less regular seasonal outbreaks. Both the pairwise correlations in outbreak timing and outbreak trend between mesoresgions were best predicted by distance. Our analysis also revealed the presence of distinct geographic clusters where dengue properties tend to be spatially correlated. Forecasting models aimed at predicting the dynamics of dengue activity need to identify the most salient variables capable of contributing to accurate predictions. Our findings show that successful models may need to leverage distinct variables in different locations and be catered to a specific task, such as predicting outbreak magnitude or timing characteristics, to be useful. This advocates in favor of “adaptive models” rather than “one-size-fits-all” models. The results of this study can be applied to improving spatial hierarchical or target-focused forecasting models of dengue activity across Brazil.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Overview of Algorithms for Using Particle Morphology in Pre-Detonation Nuclear Forensics

A major goal in pre-detonation nuclear forensics is to infer the processing conditions and/or facility type that produced radiological material. This review paper focuses on analyses of particle size, shape, texture (“morphology”) signatures that could provide information on the provenance of interdicted materials. For example, uranium ore concentrates (UOC or yellowcake) include ammonium diuranate (ADU), ammonium uranyl carbonate (AUC), sodium diuranate (SDU), magnesium diuranate (MDU), and others, each prepared using different salts to precipitate U from solution. Once precipitated, UOCs are often dried and calcined to remove adsorbed water. The products can be allowed to react further, forming uranium oxides UO3, U3O8, or UO2 powders, whose surface morphology can be indicative of precipitation and/or calcination conditions used in their production. This review paper describes statistical issues and approaches in using quantitative analyses of measurements such as particle size and shape to infer production conditions. Statistical topics include multivariate t tests (Hotelling’s T 2 ), design of experiments, and several machine learning (ML) options including decision trees, learning vector quantization neural networks, mixture discriminant analysis, and approximate Bayesian computation (ABC). ABC is emphasized as an attractive option to include the effects of model uncertainty in the selected and fitted forward model used for inferring processing conditions.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗