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,315 records · Page 73

Lens models and magnification maps of the six Hubble Frontier Fields clusters

We present strong-lensing models as well as mass and magnification maps for the cores of the six Hubble Space Telescope (HST) Frontier Fields galaxy clusters. Our parametric lens models are constrained by the locations and redshifts of multiple image systems of lensed background galaxies. We use a combination of photometric redshifts and spectroscopic redshifts of the lensed background sources obtained by us (for A2744 and AS1063), collected from the literature, or kindly provided by the lensing community. Using our results, we (1) compare the derived mass distribution of each cluster to its light distribution, (2) quantify the cumulative magnification power of the HST Frontier Fields clusters, (3) describe how our models can be used to estimate the magnification and image multiplicity of lensed background sources at all redshifts and at any position within the cluster cores, and (4) discuss systematic effects and caveats resulting from our modeling methods. We specifically investigate the effect of the use of spectroscopic and photometric redshift constraints on the uncertainties of the resulting models. We find that the photometric redshift estimates of lensed galaxies are generally in excellent agreement with spectroscopic redshifts, where available. However, the flexibility associated with relaxed redshift priors may cause the complexity of large-scale structure that is needed to account for the lensing signal to be underestimated. Our findings thus underline the importance of spectroscopic arc redshifts, or tight photometric redshift constraints, for high precision lens models. All products from our best-fit lens models (magnification, convergence, shear, deflection field) and model simulations for estimating errors are made available via the Mikulski Archive for Space Telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling resistance to the broadly neutralizing antibody PGT121 in people living with HIV-1

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization. However, two individuals experienced long-term ART-free viral suppression following antibody infusion and retained sensitivity to PGT121 upon viral rebound. Here, we develop mathematical models of the HIV-1 dynamics during this phase I clinical trial. We utilize these models to understand the dynamics leading to PGT121 resistance and to identify the mechanisms driving the observed long-term viral control. Our modeling highlights the importance of the relative fitness difference between PGT121 sensitive and resistant subpopulations prior to treatment. Specifically, by fitting our models to data, we identify the treatment-induced competitive advantage of previously existing or newly generated resistant population as a primary driver of resistance. Finally, our modeling emphasizes the high neutralization ability of PGT121 in both participants who exhibited long-term viral control.

60 APPLIED LIFE SCIENCES↗

Forbidden coronal iron-like emission in the Puppis A shock front - The effect of inhomogeneities

CCD images of the shock front at the eastern rim of Pup A, obtained in the forbidden lines of Fe X (6374 A) and Fe XIV (5303 A), are reported and compared to Einstein high-resolution-image soft X-ray data. The observed part of the remnant is complex, containing density irregularities. Optical and X-ray data are consistent in showing a nearly flat gradient of ionization temperature behind the shock. To determine conditions in the shock, scans of surface brightness across it in the optical lines were compared to surface brightnesses predicted by idealized Sedov models. It was not possible to match both the red and green line scans by a simple, single-component model, and the failure is ascribed to the presence of the density inhomogeneities. This result has important implications for the determination of SNR shock-front models by fitting X-ray data with Sedov models.

Teske, Richard G.↗

A Stochastic Model of Space-Time Variability of Mesoscale Rainfall: Statistics of Spatial Averages

A characteristic feature of rainfall statistics is that they depend on the space and time scales over which rain data are averaged. A previously developed spectral model of rain statistics that is designed to capture this property, predicts power law scaling behavior for the second moment statistics of area-averaged rain rate on the averaging length scale L as L right arrow 0. In the present work a more efficient method of estimating the model parameters is presented, and used to fit the model to the statistics of area-averaged rain rate derived from gridded radar precipitation data from TOGA COARE. Statistical properties of the data and the model predictions are compared over a wide range of averaging scales. An extension of the spectral model scaling relations to describe the dependence of the average fraction of grid boxes within an area containing nonzero rain (the "rainy area fraction") on the grid scale L is also explored.

Kundu, Prasun K.↗

Global Variability in Sonic Boom Exposure due to Macroscopic Effects

Supersonic flight over land has been prohibited since 1973 due to the loudness of sonic booms. NASA is building the X-59 aircraft as part of its Quesst mission to demonstrate low-loudness shaped sonic booms, or “sonic thumps.” The Quesst mission will gather human perception data via a series of community noise surveys across the USA. The noise dose and perceptual response data will be provided to the International Civil Aviation Organization (ICAO) and the Federal Aviation Administration for use in determining potential future supersonic aircraft noise certification standards, effectively changing the prohibition from a speed limit to a noise limit. These noise regulations must be globally effective, as long travel distances see the largest benefit to supersonic flight. The state of the atmosphere through which a sonic boom travels affects the size of the region exposed to sound, the “carpet width” (CW), as well as the loudness. The focus of this dissertation is to understand and quantify the expected loudness and CW of sonic booms due to the macroscopic atmospheric effects around the world. A pair of large-scale propagation simulation studies were conducted using the NASA PCBoom code to compare predicted sonic boom loudness and CW statistics first across the USA and then across the world. For the USA study, near-field data of the X-59 in steady cruise was propagated at 4 cardinal headings at 138 locations through 5 years of Climate Forecast System Version 2 (CFSv2) atmospheric profiles. Results of a bootstrap forest predictor screening model indicated the importance of climate zone, latitude, ground elevation, season, and heading. It also noted the unimportance of time of day for predicting loudness and CW. The data is visualized in aggregate, and then broken out geographically, by season and heading, and by climate zone. Multiple linear regression models were fit to the data from the 138 locations so that estimates of the loudness and CW can be produced anywhere in the US. The results can aid in planning when and where to fly the X-59. For the global study, near-field data from three aircraft, the X-59 in a quiet and loud configuration, B-58, and Concorde, were propagated at four cardinal headings through data from three atmospheric models, the CFSv2, the Global Forecast System (GFS), and the ECMWF Reanalysis Version 5 (ERA5), at 100 global locations over 1 year. Results of a bootstrap forest predictor screening model indicated the importance of climate zone, ground elevation, season, and heading. Similar to the US study, the model indicated time of day was not an important predictor. The model also indicated that choice of weather model was not important, so the atmospheric model data are effectively interchangeable. The ERA5 model was chosen for use in an extension of the study to include 18 additional locations to ensure sampling of every climate zone. Loudness and CW results are shown in aggregate, and split geographically and by heading, season, and climate. Multiple linear regression models were fit to the data from the 118 locations so that estimates of loudness and CW can be produced around the world. N-waves and shaped booms did not have the same global variability. Koppen-Geiger climate zones were used as the climate zone definition for the global study. These are available as present-day and future climate projections. Making use of the multiple linear regression models, the future climate zones were input to estimate the effect of the changing climate on sonic boom loudness and CW. Results indicate that a changing climate would have little impact on the effectiveness of noise regulations.

X-59↗

Curve fitting methods for solar radiation data modeling

This paper studies the use of several type of curve fitting method to smooth the global solar radiation data. After the data have been fitted by using curve fitting method, the mathematical model of global solar radiation will be developed. The error measurement was calculated by using goodness-fit statistics such as root mean square error (RMSE) and the value of R{sup 2}. The best fitting methods will be used as a starting point for the construction of mathematical modeling of solar radiation received in Universiti Teknologi PETRONAS (UTP) Malaysia. Numerical results indicated that Gaussian fitting and sine fitting (both with two terms) gives better results as compare with the other fitting methods.

14 SOLAR ENERGY↗

New method for fitting coefficients in standard model effective theory

We present an alternative method for carrying out a principal-component analysis of Wilson coefficients in standard model effective field theory (SMEFT). The method is based on singular-value decomposition (SVD). The SVD method provides information about the sensitivity of experimental observables to physics beyond the standard model that is not accessible in the Fisher-information method. In principle, the SVD method can also have computational advantages over diagonalization of the Fisher information matrix. We demonstrate the SVD method by applying it to the dimension-6 coefficients for the process of top-quark decay to a b quark and a W boson and use this example to illustrate some pitfalls in widely used fitting procedures. We also outline an iterative procedure for applying the SVD method to dimension-8 SMEFT coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Combined selection of the dynamic model and modeling error in nonlinear aeroelastic systems using Bayesian Inference

Here, we report a Bayesian framework for concurrent selection of physics-based models and (modeling) error models. We investigate the use of colored noise to capture the mismatch between the predictions of calibrated models and observational data that cannot be explained by measurement error alone within the context of Bayesian estimation for stochastic ordinary differential equations. Proposed models are characterized by the average data-fit, a measure of how well a model fits the measurements, and the model complexity measured using the Kullback–Leibler divergence. The use of a more complex error models increases the average data-fit but also increases the complexity of the combined model, possibly over-fitting the data. Bayesian model selection is used to find the optimal physical model as well as the optimal error model. The optimal model is defined using the evidence, where the average data-fit is balanced by the complexity of the model. The effect of colored noise process is illustrated using a nonlinear aeroelastic oscillator representing a rigid NACA0012 airfoil undergoing limit cycle oscillations due to complex fluid–structure interactions. Several quasi-steady and unsteady aerodynamic models are proposed with colored noise or white noise for the model error. The use of colored noise improves the predictive capabilities of simpler models.

42 ENGINEERING↗

Pharmacokinetic Modeling of Intranasal Scopolamine in Plasma Saliva and Urine

An intranasal gel formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness. The bioavailability and pharmacokinetics (PK) were evaluated under the Food and Drug Administration guidelines for clinical trials for an Investigative New Drug (IND). The aim of this project was to develop a PK model that can predict the relationship between plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trial with INSCOP. METHODS: Twelve healthy human subjects were administered three dose levels (0.1, 0.2 and 0.4 mg) of INSCOP. Serial blood, saliva and urine samples were collected between 5 min to 24 h after dosing and scopolamine concentrations measured by using a validated LC-MS-MS assay. Pharmacokinetic Compartmental models, using actual dosing and sampling times, were built using Phoenix (version 1.2). Model discrimination was performed, by minimizing the Akaike Information Criteria (AIC), maximizing the coefficient of determination (r²) and by comparison of the quality of fit plots. RESULTS: The best structural model to describe scopolamine disposition after INSCOP administration (minimal AIC =907.2) consisted of one compartment for plasma, saliva and urine respectively that were inter-connected with different rate constants. The estimated values of PK parameters were compiled in Table 1. The model fitting exercises revealed a nonlinear PK for scopolamine between plasma and saliva compartments for K21, Vmax and Km. CONCLUSION: PK model for INSCOP was developed and for the first time it satisfactorily predicted the PK of scopolamine in plasma, saliva and urine after INSCOP administration. Using non-linear PK yielded the best structural model to describe scopolamine disposition between plasma and saliva compartments, and inclusion of non-linear PK resulted in a significant improved model fitting. The model can be utilized to predict scopolamine plasma concentration using saliva and/or urine data that allows non-invasive assessment of pharmacotherapeutics of scopolamine in space and other remote environments without requiring blood sampling.

Wu, L.↗

Spectral variations in the low-mass X-ray binary 4U 0614+09

We report the X-ray spectra of a binary 4U 0614+09 observed with EXOSAT on five different occasions in 1984 and 1986. In four observations in 1984, the source intensity was 'low' compared to its intensity in the 1986 observation. The source was steady on timescales of 10 to 1000 s in each of the 1984 observations but highly variable during the 1986 variations. The best-fit spectral models for the continuum are a blackbody spectrum plus an emission from a Comptonized region. In addition, a component due to reflection from a cold disk around the source was needed as an absorbiton edge due to cold or ionized iron was seen near 7 keV with good significance in most of the spectra, particularly in the 'low' state. The relation of the hardness ratios with intensity follow the behavior seen in the 'Atoll' sources. We have studied the behavior of the hardness ratios with the variations in the parameters of the best-fit spectral models. Most of the variability is associated with the blackbody component thought to arise from the boundary layer on the surface of the neutron star. During the 'high' state in 1986, the blackbody temperature is higher, its size smaller, and its contribution significantly greater than in the 'low' state. The other notable difference between the two states is the lower opacity in the Comptonized region, and a significantly higher absorption in the 'high' state. A weak line emission (equivalent width = 10 to 40 eV) is detected during the 'high' state. The line energy lies between 6 and 7 keV wh ich implies its origin from ionized iron.

Singh, K. P.↗

Maximum Likelihood Analysis of Low Energy CDMS II Germanium Data

We report on the results of a search for a Weakly Interacting Massive Particle (WIMP) signal in low-energy data of the Cryogenic Dark Matter Search experiment using a maximum likelihood analysis. A background model is constructed using GEANT4 to simulate the surface-event background from Pb210decay-chain events, while using independent calibration data to model the gamma background. Fitting this background model to the data results in no statistically significant WIMP component. In addition, we also perform fits using an analytic ad hoc background model proposed by Collar and Fields, who claimed to find a large excess of signal-like events in our data. Finally, we confirm the strong preference for a signal hypothesis in their analysis under these assumptions, but excesses are observed in both single- and multiple-scatter events, which implies the signal is not caused by WIMPs, but rather reflects the inadequacy of their background model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Temperature dependent simulation of diamond depleted Schottky PIN diodes

Diamond is considered as an ideal material for high field and high power devices due to its high breakdown field, high lightly doped carrier mobility, and high thermal conductivity. The modeling and simulation of diamond devices are therefore important to predict the performances of diamond based devices. In this context, we use Silvaco{sup ®} Atlas, a drift-diffusion based commercial software, to model diamond based power devices. The models used in Atlas were modified to account for both variable range and nearest neighbor hopping transport in the impurity bands associated with high activation energies for boron doped and phosphorus doped diamond. The models were fit to experimentally reported resistivity data over a wide range of doping concentrations and temperatures. We compare to recent data on depleted diamond Schottky PIN diodes demonstrating low turn-on voltages and high reverse breakdown voltages, which could be useful for high power rectifying applications due to the low turn-on voltage enabling high forward current densities. Three dimensional simulations of the depleted Schottky PIN diamond devices were performed and the results are verified with experimental data at different operating temperatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of band offsets and shunt resistance in CdTe solar cells through temperature and intensity dependence of open circuit voltage and photoluminescence

Band offsets at the back contact and front window layer in CdTe-based solar cells affect photovoltaic performance and challenge standard characterization methods. By analyzing the temperature and excitation dependence of both open circuit voltage and absolute photoluminescence intensity, we show that the effects band offsets can be separated from the effects of recombination and shunting. Solar cells were grown with MgZnO window layers and compared to cells with CdS window layers containing varying amounts of oxygen. It was demonstrated that band alignment rather than reduced recombination velocity is the reason for the success of MgZnO as a front interface contact. An assortment of thin back contact interlayers were also deposited, and a PbTe interlayer showed some promise as an Ohmic contact to the CdTe, though it appears to induce a photoconductive shunt. Finally, we show that the shunting resistance given by a standard current-voltage curve technique generally does not represent a physically meaningful quantity unless it is well below one kiloOhm square cm.

14 SOLAR ENERGY↗

Extended magnetic exchange interactions in the high-temperature ferromagnet MnBi

The high-temperature ferromagnet MnBi continues to receive attention as a candidate to replace rare-earth-containing permanent magnets in applications above room temperature. This is due to a high Curie temperature, large magnetic moments, and a coercivity that increases with temperature. The synthesis of MnBi also allows for crystals that are free of interstitial Mn, enabling more direct access to the key interactions underlying the physical properties of binary Mn-based ferromagnets. In this work, we use inelastic neutron scattering to measure the spin waves of MnBi in order to characterize the magnetic exchange at low temperature. Consistent with the spin reorientation that occurs below 140 K, we do not observe a spin gap in this system above our experimental resolution. A Heisenberg model was fit to the spin wave data in order to characterize the long-range nature of the exchange. It was found that interactions up to sixth nearest neighbor are required to fully parametrize the spin waves. Surprisingly, the nearest-neighbor term is antiferromagnetic, and the realization of a ferromagnetic ground state relies on the more numerous ferromagnetic terms beyond nearest neighbor, suggesting that the ferromagnetic ground state arises as a consequence of the long-ranged interactions in the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Temperature-dependent optical measurements of the dominant recombination mechanisms in InAs/InAsSb type-2 superlattices

Temperature-dependent measurements of carrier recombination rates using a time-resolved optical pump-probe technique are reported for mid-wave infrared InAs/InAs{sub 1−x}Sb{sub x} type-2 superlattices (T2SLs). By engineering the layer widths and alloy compositions, a 16 K band-gap of ∼235 ± 10 meV was achieved for five unintentionally and four intentionally doped T2SLs. Carrier lifetimes were determined by fitting lifetime models based on Shockley-Read-Hall (SRH), radiative, and Auger recombination processes to the temperature and excess carrier density dependent data. The minority carrier (MC), radiative, and Auger lifetimes were observed to generally increase with increasing antimony content and decreasing layer thickness for the unintentionally doped T2SLs. The MC lifetime is limited by SRH processes at temperatures below 200 K in the unintentionally doped T2SLs. The extracted SRH defect energy levels were found to be near mid-bandgap. Also, it is observed that the MC lifetime is limited by Auger recombination in the intentionally doped T2SLs with doping levels greater than n ∼ 10{sup 16} cm{sup −3}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimal Methods for Estimating Cactus Pear Biomass Using Cladode Dimensions of Morphologically Diverse Accessions

Current allometric methods for photosynthetic-stem (cladode) plants, such as cactus pear (Opuntia spp.), require refinement to be used in field settings in which diverse accessions are grown. We analysed cladode dimensional data using 14 accessions representing four species and two hybrids to quantify statistically significant morphological differences among accessions and derived cross-accession models to approximate cladode fresh weight. A Box model using cladode dimensions (e.g., length, width, thickness and diameter) and factorial combinations of these measures (e.g., length*width*thickness*diameter vs. fresh weight) resulted in the highest coefficient of determination (R 2 = 0.95 general fit) across all accessions for estimating fresh weight along with parsimony estimates using the Schwarz–Bayes Criterion (SBC), which assesses the most consistent performance on individual accessions. A Fitting-box modelling approach used the measured cladode area captured using ImageJ (R 2 = 0.93 general fit). Lastly, an Elliptical model used an elliptical approximation for the measured area and performed well over all accessions (R 2 = 0.94 general fit) while avoiding extensive manual measurements. These models meet or exceed the performance of previously published approaches when applied across morphologically diverse accessions, providing efficient tools for nondestructive estimation of cactus pear biomass under the conditions tested.

Opuntia↗

Multivariate analysis of gamma spectra to characterize used nuclear fuel

Abstract—The Multi-Isotope Process (MIP) Monitor provides an efficient approach to monitoring the process conditions in used nuclear fuel reprocessing facilities to support process verification and validation. The MIP Monitor applies multivariate analysis to gamma spectroscopy of reprocessing streams in order to detect small changes in the gamma spectrum, which may indicate changes in process conditions. This research extends the MIP Monitor by characterizing a used fuel sample after initial dissolution according to the type of reactor of origin (pressurized or boiling water reactor), initial enrichment, burn up, and cooling time. Simulated gamma spectra were used to develop and test three fuel characterization algorithms. The classification and estimation models employed are based on the partial least squares regression (PLS) algorithm. A PLS discriminate analysis model was developed which perfectly classified reactor type. Locally weighted PLS models were fitted on-the-fly to estimate continuous fuel characteristics. Burn up was predicted within 0.1% root mean squared percent error (RMSPE) and both cooling time and initial enrichment within approximately 2% RMSPE. This automated fuel characterization can be used to independently verify operator declarations of used fuel characteristics and inform the MIP Monitor anomaly detection routines at later stages of the fuel reprocessing stream to improve sensitivity to changes in operational parameters and material diversions.

Multivariate Analysis for Characterization of Used↗

Liquidus temperature in the spinel primary phase field: A comparison between optical and crystal fraction methods

Liquidus temperature (TL) was measured for simulated high-level waste borosilicate glasses covering a Hanford composition region, using an optical method and a crystal-fraction extrapolation method with X-ray diffraction data from isothermal heat treatments. The 38 glasses encompassed a one-component-at-a-time variation of a 16-component matrix. The TL values ranged from 1006°C to 1603°C. First-order polynomial models were fitted to data to obtain coefficients in terms of the effect of 1 mass% component addition on the TL: TL-increasing oxides Cr2O3 (264°C), “Others” (minor components, 163°C), oxides of noble metals (137°C), NiO (91°C), Al2O3 and Fe2O3 (~19–21°C); TL-decreasing oxides K2O (-26°C), Na2O (-41°C), and Li2O (-68°C); oxides of little effect MnO, P2O5, ZrO2, F, Bi2O3, SiO2, B2O3, and CaO (9 to -12°C). Also presented are temperatures (T1%) at which glasses contain 1 vol% of spinel as these values are considered relevant to the Hanford Tank Waste Treatment and Immobilization Plant. The measured and estimated values are compared and contrasted and the effect of TL and T1% on glass formulation is discussed.

crystal fraction↗