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 937 records · Page 52

Tape-in-tape-out (TITO): a new approach for in-situ contact resistance measurements in high temperature superconducting CORC® cables

One of the ongoing development challenges with ReBCO high-temperature superconducting (HTS) cables is normal zone initiation and local heating, which is associated with over-critical current flowing through and around local performance reductions in individual tapes. Although inter-tape contact resistances are well-reported for individual tapes and HTS cables, these measurements are a challenge in CORC® cables as current percolates through much of the cable. In this work, developments in tape-in-tape-out (TITO) automated individual tape powering experiments are presented, and a simplified modeling approach for current percolation in CORC® cables is developed. An optimization is formulated to fit the model parameters to a large set of TITO experiments on a single cable, allowing the layer-dependent inter-tape contact resistance to be extracted. Measurements are presented and discussed for a straight CORC® cable and a cable bent to a 152 mm and 76 mm radius. The approach provides both quantitative insight and facilitates better understanding of current sharing phenomena in HTS cables.

Teyber, Reed↗

GPS measurements of crustal deformation associated with the 22 April 1991, Valle de la Estrella , Costa Rica earthquake

Results are presented for analysis of Global Positioning System (GPS) measurements made at sites in Costa Rica in February and July 1991. Significant horizontal and vertical displacements relative to February positions were observed. Differences were found in GPS derived vertical and horizontal displacements as compared to other types of geodetic measurements of uplift in the coastal regions. A slip dislocation model which fits the GPS measured displacement was computed. Differences between these data sets and their associated models were investigated in light of unmodeled slip heterogeneity on the fault and post-seismic displacements.

Lundgren, Paul R.↗

System Identification for eVTOL Aircraft Using Simulated Flight Data

This paper describes a system identification method for electric vertical takeoff and landing (eVTOL) aircraft. The approach merges fixed-wing and rotary-wing modeling techniques with new strategies to develop a modeling method for eVTOL vehicles using flight test data. The eVTOL aircraft system identification approach is demonstrated through application to the NASA LA-8 tandem tilt-wing, distributed electric propulsion aircraft using a high-fidelity flight dynamics simulation. Orthogonal phase-optimized multisine inputs are applied to each control surface and propulsor at numerous flight conditions throughout the flight envelope to collect informative flight data. An aero-propulsive model is identified at each flight condition using the equation-error method in the frequency domain. The local model parameters are then blended to create a global model across the nominal flight envelope. Parameter estimation results are shown to provide a good fit to modeling data and have good prediction capability. The methodology is developed with a discussion of unique eVTOL vehicle aerodynamic characteristics and practical strategies intended to inform future flight-based system identification efforts for eVTOL aircraft.

system identification↗

Development and Validation of an Initial Electric Motor Rotor Vibration Model

An overview of the development of a low fidelity prediction tool for electric motor rotor resonant modes and frequencies. These predictions are useful for the prediction of UAM motor noise during the design phase. Finite Element Analysis of motor rotors has led to the development of a parameterized model to predict these modes shapes. This has in turn allowed to a curve-fit based model from parametric sweeps. These curves can predict resonant frequencies and shapes for 1-5kW class motors without running a simulation.

acoustics↗

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗

Reconstructing the Universe: Testing the Mutual Consistency of the Pantheon and SDSS/eBOSS BAO Data Sets with Gaussian Processes

We test the mutual consistency between the baryon acoustic oscillation measurements from the eBOSS SDSS final release and the Pantheon supernova compilation in a model-independent fashion using Gaussian process regression. We also test their joint consistency with the ΛCDM model in a model-independent fashion. We also use Gaussian process regression to reconstruct the expansion history that is preferred by these two data sets. While this methodology finds no significant preference for model flexibility beyond ΛCDM, we are able to generate a number of reconstructed expansion histories that fit the data better than the best-fit ΛCDM model. These example expansion histories may point the way toward modifications to ΛCDM. We also constrain the parameters Ω{sub k} and H {sub 0} r {sub d} both with ΛCDM and with Gaussian process regression. We find that H {sub 0} r {sub d} = 10,030 ± 130 km s{sup −1} and Ω{sub k} = 0.05 ± 0.10 for ΛCDM and that H {sub 0} r {sub d} = 10,040 ± 140 km s{sup −1} and Ω{sub k} = 0.02 ± 0.20 for the Gaussian process case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Impact of SN1987A with its Circumstellar Ring

Following is final report on the study "The Impact of SN1987A with its Circumstellar Ring", which is now complete. In 1994, it was predicted that the blast wave from SN1987A would strike the circumstellar ring in AD1999+/-3, and that the ring would brighten by several hundreds optically as the shock entered the ring. It was also predicted that the emission lines from the shocked ring would have linewidths of a few hundred km/s. Today, we see a "hot spot" on the ring that first appeared in Hubble Space Telescope (HST) images in 1995 and has doubled in brightness between August 1997 and February 1998. Moreover, spectra from STIS show that the emission lines from the hot spot have widths and blue shifts of order 200 km/s, just as we predicted in 1994. Our guess of 1999+/-3 for the impact time was lucky, because we assumed that the hot gas in the bubble between the supernova and the ring had a low density (approximately 10 cm(exp -3)). But it was pointed out that the ROSAT observations of soft X-rays from SN 1987A implied that the intervening gas had a higher density, (approximately 100 cm(exp -3)), which would delay the impact until ca. 2007. Others developed hydrodynamic models to fit the ROSAT X-ray emission spectrum and came to the same conclusion. But these models were oversimplified in that they assumed that the ring was perfectly round. Now we see clearly that the hot spot is a peninsula that protrudes inward from the ring -- the first spot on the ring to be struck by the blast wave. No doubt there are other protrusions on the ring, which we may expect to light up in the next few years until they finally merge to set the entire ring ablaze, probably within 5 - 10 years. It was also predicted that the X-ray emission must be accompanied by optical and ultraviolet emission from atoms in the supernova debris and in the circumstellar gas that cross the reverse shock and the blast wave, respectively, and that the Lyman-alpha and NV emission lines should be bright enough to see with the STIS.

McCray, Richard↗

Modelling above-ground biomass stock over Norway using national forest inventory data with ArcticDEM and Sentinel-2 data

Boreal forests constitute a large portion of the global forest area, yet they are undersampled through field surveys, and only a few remotely sensed data sources provide structural information wall-to-wall throughout the boreal domain. ArcticDEM is a collection of high-resolution (2 m) space-borne stereogrammetric digital surface models (DSM) covering the entire land area north of 60° of latitude. The free-availability of ArcticDEM data offers new possibilities for aboveground biomass mapping (AGB) across boreal forests, and thus it is necessary to evaluate the potential for these data to map AGB over alternative open-data sources (i.e., Sentinel-2). This study was performed over the entire land area of Norway north of 60° of latitude, and the Norwegian national forest inventory (NFI) was used as a source of field data composed of accurately geolocated field plots (n=7710) systematically distributed across the study area. Separate random forest models were fitted using NFI data, and corresponding remotely sensed data consisting of either: i) a canopy height model (ArcticCHM) obtained by subtracting a high-quality digital terrain model (DTM) from the ArcticDEM DSM height values, ii) Sentinel-2 (S2), or iii) a combination of the two (ArcticCHM+S2). Furthermore, we assessed the effect of the forest- and terrain-specific factors on the models’ predictive accuracy. The best model (,i.e., ArcticCHM+S2) explained nearly 60% of the variance of the training set, which translated in the largest accuracy in terms of root mean square error (RMSE=41.4 t/ha). This result highlights the synergy between 3D and multispectral data in AGB modelling. Furthermore, this study showed that despite the importance of ArcticCHM variables, the S2 model performed slightly better than ArcticCHM model. This finding highlights some of the limitations of ArcticDEM, which, despite the unprecedented spatial resolution, is highly heterogeneous due to the blending of multiple acquisitions across different years and seasons. We found that both forest- and terrain-specific characteristics affected the uncertainty of the ArcticCHM+S2 model and concluded that the combined use of ArcticCHM and Sentinel-2 represents a viable solution for AGB mapping across boreal forests. The synergy between the two data sources allowed for a reduction of the saturation effects typical of multispectral data while ensuring the spatial consistency in the output predictions due to the removal of artifacts and data voids present in ArcticCHM data. While the main contribution of this study is to provide the first evidence of the best-case-scenario (i.e., availability of accurate terrain models) that ArcticDEM data can provide for large-scale AGB modelling, it remains critically important for other studies to investigate how ArcticDEM may be used in areas where no DTMs are available as is the case for large portions of the boreal zone.

space-borne imagery↗

Using Random Forest Models to Predict Organizational Violence

We present a methodology to access the proclivity of an organization to commit violence against nongovernment personnel. We fitted a Random Forest model using the Minority at Risk Organizational Behavior (MAROS) dataset. The MAROS data is longitudinal; so, individual observations are not independent. We propose a modification to the standard Random Forest methodology to account for the violation of the independence assumption. We present the results of the model fit, an example of predicting violence for an organization; and finally, we present a summary of the forest in a "meta-tree,"

Levine, Burton↗

On the reconciliation of simultaneous microwave imaging and hard X-ray observations of a solar flare

Microwave imaging data for a small flare with simultaneous hard X-ray spectral observations are compared. The X-ray data suggest that the power-law index delta of the energy distribution of the radiating electrons is 5.3 (thick-target) which differs significantly from the estimate (delta = 1.4) from a homogeneous optically-thin gyrosynchrotron model which fits the radio observations well. In order to reconcile these results, a double power-law energy spectrum is investigated for the energetic electrons in the flare, as assumed by other authors: the power law is steep at low energies and much flatter at the higher energies which produce the bulk of the microwaves. The emission of soft photons by the flat tail strongly contributes to the observed hard X-ray range and would flatten the spectrum there. A thin-target model for the X-ray emission is also inconsistent with radio data. An inhomogenous gyrosychrotron model with a number of free parameters and containing an electron distribution given by the thick-target X-ray model could be made to fit the radio data.

Nitta, N.↗

Gemini Planet Imager Spectroscopy of the Dusty Substellar Companion HD 206893 B

We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The J, H, K1, and K2 spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 to 1800 K for effective temperature and from 3.0 to 5.0 for log(g), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of submicron dust particles, similar to dereddening approaches applied to the peculiar red field L dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from 4 yr of astrometric monitoring is consistent with a ∼30 yr period, an orbital inclination of 147°, and a semimajor axis of 10 au, well within the estimated disk inner radius of ∼50 au. As one of a very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion–disk interaction and provide a benchmark for substellar and planetary atmospheric study.

47 OTHER INSTRUMENTATION↗

Kepler and TESS Observations of PG 1159-035

Abstract PG 1159-035 is the prototype of the PG 1159 hot (pre-)white dwarf pulsators. This important object was observed during the Kepler satellite K2 mission for 69 days in 59 s cadence mode and by the TESS satellite for 25 days in 20 s cadence mode. We present a detailed asteroseismic analysis of those data. We identify a total of 107 frequencies representing 32 ℓ = 1 modes, 27 frequencies representing 12 ℓ = 2 modes, and eight combination frequencies. The combination frequencies and the modes with very high k values represent new detections. The multiplet structure reveals an average splitting of 4.0 ± 0.4 μ Hz for ℓ = 1 and 6.8 ± 0.2 μ Hz for ℓ = 2, indicating a rotation period of 1.4 ± 0.1 days in the region of period formation. In the Fourier transform of the light curve, we find a significant peak at 8.904 ± 0.003 μ Hz suggesting a surface rotation period of 1.299 ± 0.002 days. We also present evidence that the observed periods change on timescales shorter than those predicted by current evolutionary models. Our asteroseismic analysis finds an average period spacing for ℓ = 1 of 21.28 ± 0.02 s. The ℓ = 2 modes have a mean spacing of 12.97 ± 0.4 s. We performed a detailed asteroseismic fit by comparing the observed periods with those of evolutionary models. The best-fit model has T eff = 129, 600 ± 11 100 K, M * = 0.565 ± 0.024 M ⊙ , and log g = 7.41 − 0.54 + 0.38 , within the uncertainties of the spectroscopic determinations. We argue for future improvements in the current models, e.g., on the overshooting in the He-burning stage, as the best-fit model does not predict excitation for all of the pulsations detected in PG 1159-035.

79 ASTRONOMY AND ASTROPHYSICS↗

Dust scattering in the high-state eclipse of Cen X-3

Exosat data from two high-state eclipses of Cen X-3 are analyzed for the effects of dust scattering. The light curves show bright (2 per cent of the out-of-eclipse flux), smoothly varying emission which reaches its minimum after mideclipse. Of this emission, the soft flux falls monotonically with orbital phase. This soft component with dust scattering is identified and spectral models are fitted to the asymmetric emission, that is, to spectra obtained by subtracting from a spectrum at a given orbital phase before phase zero the corresponding spectrum from the equal and opposite phase after phase zero. It is found that the model of Molnar and Mauche, used to fit dust scattering in Cyg X-3, can explain the observed spectra and yields values for the distance to Cen X-3 and for the size of the dust grains of 5.4 kpc and 0.33 micron, consistent with previous observations in other wavebands.

Day, C. S. R.↗

CF + excitation in the interstellar medium

The detection of CF + in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF + abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF + in collisions with para- and ortho-H 2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF + emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF + as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗

The kinetics of SARS-CoV-2 infection based on a human challenge study

Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection.

59 BASIC BIOLOGICAL SCIENCES↗

Battery Performance, Thermal, and Life Modeling for Southern California Edison (Cooperative Research and Development Final Report)

The objective of this project is to develop battery performance models, thermal models and life models for SCE to use in simulation and analysis of grid services involving energy storage. NREL will fit the models to three separate chemistries using test data provided by SCE. The dataset for each chemistry will include electrical, thermal and aging response to different temperatures and cycling conditions, measured under a variety of cell-level and module-level experiments described in SCE test protocol documents. Once tuned to a test dataset, the battery system model software will provide predictions of battery energy and power loss, and thus cycle and calendar life, for any energy storage grid service power profile of interest to SCE. The automated battery life modeling and simulation tool can also be applied to any other grid battery systems as long as necessary input data is available to accelerate battery model development and battery lifetime analysis.

25 ENERGY STORAGE↗

Analysis of MIR-18 results for physical and biological dosimetry: radiation shielding effectiveness in LEO

We compare models of radiation transport and biological response to physical and biological dosimetry results from astronauts on the Mir space station. Transport models are shown to be in good agreement with physical measurements and indicate that the ratio of equivalent dose from the Galactic Cosmic Rays (GCR) to protons is about 3/2:1 and that this ratio will increase for exposures to internal organs. Two biological response models are used to compare to the Mir biodosimetry for chromosome aberration in lymphocyte cells; a track-structure model and the linear-quadratic model with linear energy transfer (LET) dependent weighting coefficients. These models are fit to in vitro data for aberration formation in human lymphocytes by photons and charged particles. Both models are found to be in reasonable agreement with data for aberrations in lymphocytes of Mir crew members: however there are differences between the use of LET dependent weighting factors and track structure models for assigning radiation quality factors. The major difference in the models is the increased effectiveness predicted by the track model for low charge and energy ions with LET near 10 keV/micrometers. The results of our calculations indicate that aluminum shielding, although providing important mitigation of the effects of trapped radiation, provides no protective effect from the galactic cosmic rays (GCR) in low-earth orbit (LEO) using either equivalent dose or the number of chromosome aberrations as a measure until about 100 g/cm 2 of material is used.

Flight Experiment↗

Single Event Test Methodologies and System Error Rate Analysis for Triple Modular Redundant Field Programmable Gate Arrays

We present a test methodology for estimating system error rates of Field Programmable Gate Arrays (FPGAs) mitigated with Triple Modular Redundancy (TMR). The test methodology is founded in a mathematical model, which is also presented. Accelerator data from 90 nm Xilins Military/Aerospace grade FPGA are shown to fit the model. Fault injection (FI) results are discussed and related to the test data. Design implementation and the corresponding impact of multiple bit upset (MBU) are also discussed.

error rate calculations↗