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,603 records · Page 89

Open-source release of CGMF 1.1 and Integration into the MCNP6.3 ® Code [Slides]

As a result of a multi-year NA-22 project, CGMF was integrated into MCNP6.2 and publicly released. CGMF was open-sourced and publicly released and MCNP6.3 was updated to include the latest version and is in the process of being publicly released. Current and future plans include global optimization and uncertainty quantification within CGMF, model parameter fitting such that CGMF may be used in ENDF/B evaluations, and improving both standalone and MCNP-integrated CGM (non-fission) simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Carrier Dynamics of Polar, Semipolar, and Nonpolar InGaN/GaN LEDs Measured by Small-Signal Electroluminescence

The carrier dynamics in InGaN/GaN LEDs are directly tied to their efficiency and maximum modulation speed, which are important metrics for solid-state lighting, displays, and optical communication. We measure the carrier dynamics of nonpolar, semipolar, and commercial c-plane InGaN/GaN LEDs using small-signal electroluminescence methods. Rate equations and a small-signal circuit are used to model the carrier dynamics. The model is fit to the optical frequency response and input impedance of the LEDs to extract the carrier lifetimes and the recombination rates. The results offer insight into the underlying causes of efficiency droop and the green gap, and inform device design strategies.

42 ENGINEERING↗

Scaling hyporheic nitrogen cycling in large river alluvial aquifers

Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.

54 ENVIRONMENTAL SCIENCES↗

Weak Gravitational Lensing of Low Surface Brightness Galaxies in the Dark Energy Year 3 Catalog

We present galaxy-galaxy lensing measurements of a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. Given their faintness, the use of standard observational techniques proves challenging. Weak gravitational lensing probes both the baryonic and dark matter content of galaxies, rendering it a powerful technique to estimate LSBG masses. The LSBG lens sample consists of 23, 790 total extended galaxies separated into red and blue color types at \(g-i\ge 0.60\) and \(g-i\le 0.60\), respectively. We use the Y3 \sc{metacalibration} shape catalog as the source sample, with a number density of \(5.59 \mathrm{gal}/\mathrm{arcmin}^2\). We measure the tangential shear around the lens galaxies across angular scales of \(0.25-400\) arcmin and find a signal-to-noise of 6.67 for red galaxies, 2.17 for blue galaxies, and 5.30 for the total sample. We fit a model built from two NFW profiles corresponding to the LSBG dark matter subhalo and host halo to the red LSBG shear measurements with an MCMC. We estimate the host halo mass at \(7.3 ^{+2.0}_{-1.6}\times 10^{12} M_\mathrm{\odot}\). We place a 95\% upper bound on the subhalo mass at \(3\times 10^{11} M_\mathrm{\odot}\). We utilize the lens sample photometry to obtain an estimate of the red LSBG stellar mass distribution. We compare the ratio between the stellar mass and the subhalo mass to the parameterized, satellite-specific SHMR. This work represents the first example of an attempted constraint on the masses of LSBGs using weak gravitational lensing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

CryoPDK Development for 22nm FDSOI CryoCMOS

Cryogenic Process Design Kits (PDKs) are an indispensable tool in the design of complex integrated circuits across a wide spectrum of applications, from noble element detectors to Quantum Information Science, Superconducting Nanowire Single Photon Detectors (SNSPDs), and precision atomic clocks. The development of PDK-compatible SPICE models is a complex endeavor requiring test structures, measurements, model extraction and fitting. We will present the cryogenic modeling and development of a cryo-PDK for a 22nm FDSOI CMOS process for operation at 3.8 Kelvin.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Risk Ratio and Risk Difference Estimation in Case-cohort Studies

Background: In case-cohort studies with binary outcomes, ordinary logistic regression analyses have been widely used because of their computational simplicity. However, the resultant odds ratio estimates cannot be interpreted as relative risk measures unless the event rate is low. The risk ratio and risk difference are more favorable outcome measures that are directly interpreted as effect measures without the rare disease assumption. Methods: We provide pseudo-Poisson and pseudo-normal linear regression methods for estimating risk ratios and risk differences in analyses of case-cohort studies. These multivariate regression models are fitted by weighting the inverses of sampling probabilities. Also, the precisions of the risk ratio and risk difference estimators can be improved using auxiliary variable information, specifically by adapting the calibrated or estimated weights, which are readily measured on all samples from the whole cohort. Finally, we provide computational code in R (R Foundation for Statistical Computing, Vienna, Austria) that can easily perform these methods. Results: Through numerical analyses of artificially simulated data and the National Wilms Tumor Study data, accurate risk ratio and risk difference estimates were obtained using the pseudo-Poisson and pseudo-normal linear regression methods. Also, using the auxiliary variable information from the whole cohort, precisions of these estimators were markedly improved. Conclusion: The ordinary logistic regression analyses may provide uninterpretable effect measure estimates, and the risk ratio and risk difference estimation methods are effective alternative approaches for case-cohort studies. These methods are especially recommended under situations in which the event rate is not low.

60 APPLIED LIFE SCIENCES↗

Basis light-front quantization approach to nucleon

We obtain the light-front wavefunctions for the nucleon in the valence quark Fock space from an effective Hamiltonian, which includes the transverse and longitudinal confinement and the one gluon exchange interaction with fixed coupling. The wavefunctions are generated by solving the eigenvalue equation in a basis light front quantization. Fitting the model parameters, the wavefunctions lead to good simultaneous description of electromagnetic form factors, radii, and parton distribution functions for the proton.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measuring Intelligence with the Sandia Matrices: Psychometric Review and Recommendations for Free Raven-Like Item Sets

The Sandia Matrices are a free alternative to the Raven’s Progressive Matrices (RPMs). This study offers a psychometric review of Sandia Matrices items focused on two of the most commonly investigated issues regarding the RPMs: (a) dimensionality and (b) sex differences. Model-data fit of three alternative factor structures are compared using confirmatory multidimensional item response theory (IRT) analyses, and measurement equivalence analyses are conducted to evaluate potential sex bias. Although results are somewhat inconclusive regarding factor structure, results do not show evidence of bias or mean differences by sex. Finally, although the Sandia Matrices software can generate infinite items, editing and validating items may be infeasible for many researchers. Further, to aide implementation of the Sandia Matrices, we provide scoring materials for two brief static tests and a computer adaptive test. Implications and suggestions for future research using the Sandia Matrices are discussed.

60 APPLIED LIFE SCIENCES↗

A mathematical approach to survey electrochemical impedance spectroscopy for aging in lithium-ion batteries

Electrochemical impedance spectroscopy (EIS) is a valuable technique to detect the health status and aging phenomena in lithium-ion batteries (LiB). Equivalent circuit modeling (ECM) is conventionally used when interpreting EIS data and gaining physical insights into the aging mechanisms. However, performing ECM is resource intensive and expert-level of knowledge is usually required to select suitable models and fitting parameters. This article presents a quick and user-friendly data analysis algorithm as an alternative to ECM by mathematically fitting geometric features in Nyquist plots and obtaining the growth trends of the features. The evolving trends in the Nyquist plots, such as chord lengths of the arcs and interception points, are consistent with the growth of resistance components obtained using ECM with R 2 values from 0.67 to 0.99, and therefore can be used as indicators of battery aging. Our results show that the quick-fitting approach is suitable for analyzing a series of EIS data acquired during battery cycling and identifying the underlying aging mechanisms.

Chen, Bor-Rong↗

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗

Urban Green Spaces, Greenness Exposure and Species Richness in Residential Environments and Relations with Physical Activity and BMI in Portuguese Adolescents

Environmental factors play an important role in obesity-related behaviors. Evidence indicates significant associations between weight and urban green spaces in adults, but it is not clear whether this relationship applies to adolescents. Therefore, our aim was to determine the associations between urban green spaces, greenness exposure and species richness in residential environments with physical activity and body mass index. Sixty-two adolescents between 12 and 18 years of age answered a self-administered questionnaire, providing information on height, weight, age, sex and home address. Data on socioeconomic deprivation were obtained from the European Index of Deprivation for Small Portuguese Areas. Physical activity levels were assessed using accelerometers. Urban green space counts and the normalized difference vegetation index values were measured using buffers along the roads with distances of 300, 500, 1000 and 1500 m from each participant’s residence. To quantify the species richness, the species richness index was used. Linear regression models were fitted to analyze whether urban green spaces, exposure to green spaces and species richness counts for each distance were associated with physical activity and self-reported body mass index. We did not find significant associations between the independent variables and the probability of overweight or obesity. The relationship between environmental variables, adolescents’ physical activity and body weight seems to be complex and further studies may contribute to better understanding of the topic.

Melo, Juliana (ORCID:0000000301275340)↗

Measurement of shunt resistance and conduction band offset in Cu(In,Ga)Se 2 solar cells through joint analysis of temperature and intensity dependence of open-circuit voltage and photoluminescence

Joint analysis studies of open-circuit voltage and photoluminescence intensity (PL-I) are reported for CuIn 1-x GaxSe 2 (CIGSe) solar cells. A range of compositions are investigated, including constant x=0.35 and x=0.55 as well as a graded composition profile having a minimum of x=0.25. Both the open circuit voltage and PL-I are measured as functions of temperature and illumination intensity. With these two measurements, a full model-based fitting of the temperature and illumination dependence allows extraction of the effects of window layer band offset and shunt resistance, in addition to bulk and interface recombination parameters. Finally, to quantitatively analyze the two distinct measurements jointly, the absolute PL-I is measured to obtain quasi-Fermi-level splitting.

14 SOLAR ENERGY↗

A soft X-ray spectral episode for the clocked burster, GS 1826–24 as measured by Swift and NuSTAR

We report on NuSTAR and Swift observations of a soft state of the neutron star low-mass X-ray binary GS 1826–24, commonly known as the “clocked” burster. The transition to the soft state was recorded in 2014 June through an increase of the 2–20 keV source intensity measured by MAXI, simultaneous with a decrease of the 15–50 keV intensity measured by Swift/BAT. The episode lasted approximately two months, after which the source returned to its usual hard state. We analyze the broadband spectrum measured by Swift/XRT and NuSTAR and estimate the accretion rate during the soft episode to be ≈13% m-dot {sub Edd}, within the range of previous observations. However, the best-fit spectral model, adopting the double Comptonization used previously, exhibits significantly softer components. We detect seven type-I X-ray bursts, all significantly weaker (and with shorter rise and decay times) than observed previously. The burst profiles and recurrence times vary significantly, ruling out the regular bursts that are typical for this source. One burst exhibited photospheric radius expansion and we estimate the source distance as (5.7±0.2) ξ{sub b}{sup −1/2} kpc, where ξ{sub b} parameterizes the possible anisotropy of the burst emission. The observed soft state may most likely be interpreted as a change in accretion geometry at about similar bolometric luminosity as in the hard state. The different burst behavior can therefore be attributed to this change in accretion flow geometry, but the fundamental cause and process for this effect remain unclear.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE CONTRIBUTION OF FLUORESCENCE TO LYα HALOS AROUND STAR-FORMING GALAXIES

We quantify the contribution of Lyα fluorescence to observed spatially extended Lyα halos around Lyα emitters at redshift z = 3.1. The key physical quantities that describe the fluorescent signal include (i) the distribution of cold gas in the circumgalactic medium (CGM); we explore simple analytic models and fitting functions to recent hydrodynamical simulations; and (ii) local variations in the ionizing background due to ionizing sources that cluster around the central galaxy. We account for clustering by boosting the observationally inferred volumetric production rate of ionizing photons, ϵ{sub LyC}, by a factor of 1+ξ{sub LyC}(r), in which ξ{sub LyC}(r) quantifies the clustering of ionizing sources around the central galaxy. We compute ξ{sub LyC}(r) by assigning an “effective” bias parameter to the ionizing sources. This novel approach allows us to quantify our ignorance of the population of ionizing sources in a simple parametrized form. We find a maximum enhancement in the local ionizing background in the range 50–200 at r ∼ 10 physical kpc. For spatially uncorrelated ionizing sources and fluorescing clouds we find that fluorescence can contribute up to ∼ 50%–60% of the observed spatially extended Lyα emission. We briefly discuss how future observations can shed light on the nature of Lyα halos around star-forming galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

THE DYNAMICAL RELATIONSHIP BETWEEN THE BAR AND SPIRAL PATTERNS OF NGC 1365

Theories that attempt to explain the dynamical relationship between bar and spiral patterns in galactic disks make different predictions about the radial profile of the pattern speed. These are tested for the H-alpha bar and spiral patterns of NGC 1365. The radial profile of the pattern speed is measured by fitting mathematical models that are based on the Tremaine–Weinberg method. The results show convincing evidence for the bar rotating at a faster rate than the spiral pattern, inconsistent with a global wave mode or a manifold. There is evidence for mode coupling of the bar and spiral patterns at the overlap of corotation and inner Lindblad resonances (ILRs), but the evidence is unreliable and inconsistent. The results are the most consistent with the bar and spiral patterns being dynamically distinct features. The pattern speed of the bar begins near an ILR and ends near the corotation resonance (CR). The radial profile of the pattern speed beyond the bar most closely resembles what is expected for coupled spiral modes and tidal interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

THE DIMENSIONLESS AGE OF THE UNIVERSE: A RIDDLE FOR OUR TIME

We present the interesting coincidence of cosmology and astrophysics that points toward a dimensionless age of the Universe H {sub 0} t {sub 0} that is close to one. Despite cosmic deceleration for 9 Gyr and acceleration since then, we find H {sub 0} t {sub 0} = 0.96 ± 0.01 for the ΛCDM model that fits SN Ia data from Pan-STARRS, CMB power spectra, and baryon acoustic oscillations. Similarly, astrophysical measures of stellar ages and the Hubble constant derived from redshifts and distances point to H {sub 0} t ∼ 1.0 ± 0.1. The wide range of possible values for H {sub 0} t {sub 0} realized during cosmic evolution means that we live at what appears to be a special time. This “synchronicity problem” is not precisely the same as the usual coincidence problem , because there are combinations of Ω{sub M} and Ω{sub Λ} for which the usual coincidence problem holds but for which H {sub 0} t {sub 0} is not close to 1.

79 ASTRONOMY AND ASTROPHYSICS↗

EXOPLANETARY DETECTION BY MULTIFRACTAL SPECTRAL ANALYSIS

Owing to technological advances, the number of exoplanets discovered has risen dramatically in the last few years. However, when trying to observe Earth analogs, it is often difficult to test the veracity of detection. We have developed a new approach to the analysis of exoplanetary spectral observations based on temporal multifractality, which identifies timescales that characterize planetary orbital motion around the host star and those that arise from stellar features such as spots. Without fitting stellar models to spectral data, we show how the planetary signal can be robustly detected from noisy data using noise amplitude as a source of information. For observation of transiting planets, combining this method with simple geometry allows us to relate the timescales obtained to primary and secondary eclipse of the exoplanets. Making use of data obtained with ground-based and space-based observations we have tested our approach on HD 189733b. Moreover, we have investigated the use of this technique in measuring planetary orbital motion via Doppler shift detection. Finally, we have analyzed synthetic spectra obtained using the SOAP 2.0 tool, which simulates a stellar spectrum and the influence of the presence of a planet or a spot on that spectrum over one orbital period. We have demonstrated that, so long as the signal-to-noise-ratio ≥ 75, our approach reconstructs the planetary orbital period, as well as the rotation period of a spot on the stellar surface.

79 ASTRONOMY AND ASTROPHYSICS↗