Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Improvement of a mixture experiment model relating the component proportions to the size of nanonized itraconazole particles in extemporary suspensions

A paper by Foglio Bonda et al. published previously in this journal (2016, Vol. 83, pp. 175–183) discussed the use of mixture experiment design and modeling methods to study how the proportions of three components in an extemporaneous oral suspension affected the mean diameter of drug particles (Z ave ). The three components were itraconazole (ITZ), Tween 20 (TW20), and Methocel® E5 (E5). This commentary addresses some errors and other issues in the previous paper, and also discusses an improved model relating proportions of ITZ, TW20, and E5 to Z ave . The improved model contains six of the 10 terms in the full-cubic mixture model, which were selected using a different cross-validation procedure than used in the previous paper. In conclusion, compared to the four-term model presented in the previous paper, the improved model fit the data better, had excellent cross-validation performance, and the predicted Z ave of a validation point was within model uncertainty of the measured value.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2016

WE-H-BRA-03: Development of a Model to Include the Evolution of Resistant Tumor Subpopulations Into the Treatment Optimization Process for Schedules Involving Targeted Agents in Chemoradiation Therapy

Purpose: To develop a model that includes the process of resistance development into the treatment optimization process for schedules that include targeted therapies. Further, to validate the approach using clinical data and to apply the model to assess the optimal induction period with targeted agents before curative treatment with chemo-radiation in stage III lung cancer. Methods: Growth of the tumor and its subpopulations is modeled by Gompertzian growth dynamics, resistance induction as a stochastic process. Chemotherapy induced cell kill is modeled by log-cell kill dynamics, targeted agents similarly but restricted to the sensitive population. Radiation induced cell kill is assumed to follow the linear-quadratic model. The validation patient data consist of a cohort of lung cancer patients treated with tyrosine kinase inhibitors that had longitudinal imaging data available. Results: The resistance induction model was successfully validated using clinical trial data from 49 patients treated with targeted agents. The observed recurrence kinetics, with tumors progressing from 1.4–63 months, result in tumor growth equaling a median volume doubling time of 92 days [34–248] and a median fraction of pre-existing resistance of 0.035 [0–0.22], in agreement with previous clinical studies. The model revealed widely varying optimal time points for the use of curative therapy, reaching from ∼1m to >6m depending on the patient’s growth rate and amount of pre-existing resistance. This demonstrates the importance of patient-specific treatment schedules when targeted agents are incorporated into the treatment. Conclusion: We developed a model including evolutionary dynamics of resistant sub-populations with traditional chemotherapy and radiation cell kill models. Fitting to clinical data yielded patient specific growth rates and resistant fraction in agreement with previous studies. Further application of the model demonstrated how proper timing of chemo-radiation could minimize the probability of resistance, increasing tumor control significantly.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A PRELIMINARY JUPITER MODEL

In anticipation of new observational results for Jupiter's axial moment of inertia and gravitational zonal harmonic coefficients from the forthcoming Juno orbiter, we present a number of preliminary Jupiter interior models. We combine results from ab initio computer simulations of hydrogen–helium mixtures, including immiscibility calculations, with a new nonperturbative calculation of Jupiter's zonal harmonic coefficients, to derive a self-consistent model for the planet's external gravity and moment of inertia. We assume helium rain modified the interior temperature and composition profiles. Our calculation predicts zonal harmonic values to which measurements can be compared. Although some models fit the observed (pre-Juno) second- and fourth-order zonal harmonics to within their error bars, our preferred reference model predicts a fourth-order zonal harmonic whose absolute value lies above the pre-Juno error bars. This model has a dense core of about 12 Earth masses and a hydrogen–helium-rich envelope with approximately three times solar metallicity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Spectral properties of cool stars: Extended abundance analysis of 1,617 planet-search stars

We present a catalog of uniformly determined stellar properties and abundances for 1,617 F, G, and K stars using an automated spectral synthesis modeling procedure. All stars were observed using the HIRES spectrograph at Keck Observatory. Our procedure used a single line list to fit model spectra to observations of all stars to determine effective temperature, surface gravity, metallicity, projected rotational velocity, and the abundances of 15 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, and Y). Sixty percent of the sample had Hipparcos parallaxes and V-band photometry, which we combined with the spectroscopic results to obtain mass, radius, and luminosity. Additionally, we used the luminosity, effective temperature, metallicity and α-element enhancement to interpolate in the Yonsei–Yale isochrones to derive mass, radius, gravity, and age ranges for those stars. Finally, we determined new relations between effective temperature and macroturbulence for dwarfs and subgiants. Our analysis achieved precisions of 25 K in T{sub eff}, 0.01 dex in [M/H], 0.028 dex for logg, and 0.5 km s{sup −1} in vsin i based on multiple observations of the same stars. The abundance results were similarly precise, between ∼0.01 and ∼0.04 dex, though trends with respect to T{sub eff} remained for which we derived empirical corrections. The trends, though small, were much larger than our uncertainties and are shared with published abundances. We show that changing our model atmosphere grid accounts for most of the trend in [M/H] between 5000 and 5500 K, indicating a possible problem with the atmosphere models or opacities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

AGN Coronae through a Jet Perspective

This paper presents an in-depth look at the jet and coronal properties of 41 active galactic nuclei (AGNs). Utilizing the highest quality NuSTAR, XMM-Newton, and NRAO VLA Sky Survey 1.4 GHz data, we find that the radio Eddington luminosity inversely scales with X-ray reflection fraction, and positively scales with the distance between the corona and the reflected regions in the disk. We next investigate a model fit to the data that predicts the corona is outflowing and propagates into the large-scale jet. We find this model describes the data well and predicts that the corona has mildly relativistic velocities, 0.04<β<0.40. We discuss our results in the context of disk–jet connections in AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Central Engine and Host Galaxy of RXJ 1301.9+2747: A Multiwavelength View of a Low-mass Black Hole Active Galactic Nuclei with Ultra-soft X-Ray Emission

RXJ 1301.9+2747 is an optically identified very-low-mass AGN candidate with M{sub BH}∼1×10{sup 6} M{sub ⊙}, which shows extremely soft X-ray emission and unusual X-ray variability in the form of short-lived flares. We present an analysis of multiwavelength observations of RXJ 1301.9+2747 in order to study the properties of the active nucleus and its host galaxy. The UV-to-X-ray spectrum in the quiescent state can be well and self-consistently described by a thermal and a Comptonized emission from the accretion disk, with the black body dominating ∼70% of the X-rays in the 0.2–2 keV. The same model can describe the X-ray spectrum in the flare state, but the Comptonized component becomes dominant (∼80%). The best fit implies an Eddington ratio of ∼0.14 and a black-hole mass of (1.7−2.8)×10{sup 6} M {sub ⊙}, in agreement with the estimation from the optical data within errors. However, the best-fitting model under predicts the optical flux for the HST point source by a factor of ∼2. The excess of nuclear optical emission could be attributed to a nuclear stellar cluster, which is frequently seen in low-mass AGNs. The X-ray to optical spectral slope (α{sub ox}) is lower than in most other active galaxies, which may be attributed to intrinsically X-ray weakness due to very little hot and optically thin coronal emission. We performed a pilot search for weak or hidden broad emission lines using optical spectropolarimetry observations, but no polarized broad lines are detected. The host galaxy appears to be a disk galaxy with a boxy pseudobulge or nuclear bar accounting for ∼15% of the total starlight, which is consistent with the general characteristics of the host of low-mass AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Observations of the J = 10 manifold of the pure rotational band of phosphine on Saturn

Saturn was observed in the vicinity of the J = 10 manifold of the pure rotational band of phosphine on 1984 July 10 and 12 from NASA's Kuiper Airborne Observatory with the facility far-infrared cooled grating spectrometer. On each night observations of the full disk plus rings were made at 4 to 6 discrete wavelengths which selectively sampled the manifold and the adjacent continuum. The previously reported detection of this manifold is confirmed. After subtraction of the flux due to the rings, the data are compared with disk-averaged models of Saturn. It is found that PH3 must be strongly depleted above the thermal inversion (approx. 70 mbar). The best fitting models consistent with other observational constaints indicate that PH3 is significantly depleted at even deeper atmospheric levels ( or = 500 mbar), implying an eddy diffusion coefficient for Saturn of 10 to the 4 cm sq/sec.

Haas, M. R.↗
NASA NTRS1985

Observations of the J = 10 manifold of the pure rotational band of phosphine on Saturn

Saturn was observed in the vicinity of the J = 10 manifold of the pure rotational band of phosphine on 1984 July 10 and 12 from NASA's Kuiper Airborne Observatory with the facility far-infrared cooled grating spectrometer. On each night observations of the full disk plus rings were made at 4 to 6 discrete wavelengths which selectively sampled the manifold and the adjacent continuum. The previously reported detection of this manifold is confirmed. After subtraction of the flux due to the rings, the data are compared with disk-averaged models of Saturn. It is found that PH3 must be strongly depleted above the thermal inversion (approx. 70 mbar). The best fitting models consistent with other observational constraints indicate that PH3 is significantly depleted at even deeper atmospheric levels (or = 500 mbar), implying an eddy diffusion coefficient for Saturn of 10 to the 4 cm sq/sec.

Haas, M. R.↗
NASA NTRS2004

Mass Anomalies on Ganymede

Radio Doppler data from two Ganymede encounters (G1 and G2) on the first two orbits in the Galileo mission have been analyzed previously for gravity information . For a satellite in hydrostatic equilibrium, its gravitational field can be modeled adequately by a truncated spherical harmonic series of degree two. However, a fourth degree field is required in order to fit the second Galileo flyby (G2). This need for a higher degree field strongly suggests that Ganymede s gravitational field is perturbed by a gravity anomaly near the G2 closest approach point (79.29 latitude, 123.68 west longitude). In fact, a plot of the Doppler residuals , after removal of the best-fit model for the zero degree term (GM) and the second degree moments (J2 and C22), suggests that if an anomaly exists, it is located downtrack of the closest approach point, closer to the equator.

Schubert, G.↗
NASA NTRS2011

Modelling the Huge, Herschel-Resolved Debris Ring around HD 207129

Debris disks, which are inferred from the observed infrared excess to be ensembles of dust, rocks. and probably planetesimals, arc common features of stellar systems. As the mechanisms of their fonnation and evolution are linked to those of planetary bodies, they provide valuable infonnation. The few well-resolved debris disks are even more valuable because they can serve as modelling benchmarks and help resolve degeneracies in modelling aspects such as typical grain sizes and distances. Here, we present an analysis of the HD 207129 debris disk, based on its well-covered spectral energy distribution and Herschel/PACS images obtained in the framework of the DUNES (DUst around NEarby Stars) programme. We use an empirical power-law approach to the distribution of dust and we then model the production and removal of dust by means of collisions, direct radiation pressure, and drag forces. The resulting best-fit model contains a total of nearly 10(sup -2) Earth masses in dust, with typical grain sizes in the planetesimal beh ranging from 4 to 7 micrometers. We constrain the dynamical excitation to be low, which results in very long collisional lifetimes and a drag that notably fills the inner gap, especially at 70 micrometers. The radial distribution stretches from well within 100 AU in an unusual, outward-riSing slope towards a rather sharp outer edge at about 170-190 AU. The inner edge is therefore smoother than that reported for Fomalhaut, but the contribution from the extended halo of barely bound grains is similarly small. Both slowly self-stirring and planetary perturbations could potentially have formed and shaped this disk.

Lohne, T.↗
DOE OSTI2023

Pore-resolved investigation of turbulent open channel flow over a randomly packed permeable sediment bed

Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers, Re K = 2.56 , 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). Here, the mean flow and shear penetration depths increase with Re K and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.

turbulence simulation↗
DOE OSTI2023

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

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

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗
DOE OSTI2015

Automated contouring error detection based on supervised geometric attribute distribution models for radiation therapy: A general strategy

Purpose: One of the most critical steps in radiation therapy treatment is accurate tumor and critical organ-at-risk (OAR) contouring. Both manual and automated contouring processes are prone to errors and to a large degree of inter- and intraobserver variability. These are often due to the limitations of imaging techniques in visualizing human anatomy as well as to inherent anatomical variability among individuals. Physicians/physicists have to reverify all the radiation therapy contours of every patient before using them for treatment planning, which is tedious, laborious, and still not an error-free process. In this study, the authors developed a general strategy based on novel geometric attribute distribution (GAD) models to automatically detect radiation therapy OAR contouring errors and facilitate the current clinical workflow. Methods: Considering the radiation therapy structures’ geometric attributes (centroid, volume, and shape), the spatial relationship of neighboring structures, as well as anatomical similarity of individual contours among patients, the authors established GAD models to characterize the interstructural centroid and volume variations, and the intrastructural shape variations of each individual structure. The GAD models are scalable and deformable, and constrained by their respective principal attribute variations calculated from training sets with verified OAR contours. A new iterative weighted GAD model-fitting algorithm was developed for contouring error detection. Receiver operating characteristic (ROC) analysis was employed in a unique way to optimize the model parameters to satisfy clinical requirements. A total of forty-four head-and-neck patient cases, each of which includes nine critical OAR contours, were utilized to demonstrate the proposed strategy. Twenty-nine out of these forty-four patient cases were utilized to train the inter- and intrastructural GAD models. These training data and the remaining fifteen testing data sets were separately employed to test the effectiveness of the proposed contouring error detection strategy. Results: An evaluation tool was implemented to illustrate how the proposed strategy automatically detects the radiation therapy contouring errors for a given patient and provides 3D graphical visualization of error detection results as well. The contouring error detection results were achieved with an average sensitivity of 0.954/0.906 and an average specificity of 0.901/0.909 on the centroid/volume related contouring errors of all the tested samples. As for the detection results on structural shape related contouring errors, an average sensitivity of 0.816 and an average specificity of 0.94 on all the tested samples were obtained. The promising results indicated the feasibility of the proposed strategy for the detection of contouring errors with low false detection rate. Conclusions: The proposed strategy can reliably identify contouring errors based upon inter- and intrastructural constraints derived from clinically approved contours. It holds great potential for improving the radiation therapy workflow. ROC and box plot analyses allow for analytically tuning of the system parameters to satisfy clinical requirements. Future work will focus on the improvement of strategy reliability by utilizing more training sets and additional geometric attribute constraints.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

THE SPECTRAL ENERGY DISTRIBUTION OF THE COLDEST KNOWN BROWN DWARF

WISE J085510.83–071442.5 (hereafter WISE 0855–0714) is the coldest known brown dwarf (∼250 K) and the fourth-closest known system to the Sun (2.2 pc). It has been previously detected only in the J band and two mid-IR bands. To better measure its spectral energy distribution (SED), we have performed deep imaging of WISE 0855–0714 in six optical and near-IR bands with Gemini Observatory, the Very Large Telescope, and the Hubble Space Telescope . Five of the bands show detections, although one detection is marginal (S/N ∼ 3). We also have obtained two epochs of images with the Spitzer Space Telescope for use in refining the parallax of the brown dwarf. By combining astrometry from this work and previous studies, we have derived a parallax of 0.449 ± 0.008″ (2.23 ± 0.04 pc). We have compared our photometry for WISE 0855–0714 to data for known Y dwarfs and to the predictions of three suites of models by Saumon et al. and Morley et al. that are defined by the presence or absence of clouds and nonequilibrium chemistry. Our estimates of Y − J and J − H for WISE 0855–0714 are redder than colors of other Y dwarfs, confirming a predicted reversal of near-IR colors to redder values at temperatures below 300–400 K. In color–magnitude diagrams, no single suite of models provides a clearly superior match to the sequence formed by WISE 0855–0714 and other Y dwarfs. Instead, the best-fitting model changes from one diagram to the next. Similarly, all of the models have substantial differences from the SED of WISE 0855–0714. As a result, we are currently unable to constrain the presence of clouds or nonequilibrium chemistry in its atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyα Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyα forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit–accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible IGM thermal histories by varying the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen at redshifts 2.2 ≤ z ≤ 5.0 and with the helium Lyα opacity at redshifts 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyα forest data to break model degeneracies, we find that the IGM experienced two main reheating events over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z R ≈ 6.0 with a first IGM temperature peak of T 0 ≃ 1.3 × 10 4 K and is followed by the reionization of He ii that completes by z R ≈ 3.0 and yields a second temperature peak of T 0 ≃ 1.4 × 10 4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗