Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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 865 records · Page 48

Examining normal modes as fundamental heat carriers in amorphous solids: The case of amorphous silicon

Normal mode decomposition of atomic vibrations has been used to provide microscopic understanding of thermal transport in amorphous solids for decades. In normal mode methods, it is naturally assumed that atoms vibrate around their equilibrium positions, and that individual normal modes are the fundamental vibrational excitations transporting heat. With the abundance of predictions from normal mode methods and experimental measurements now available, we carefully analyze these calculations in amorphous silicon, a model amorphous solid. We find a number of discrepancies, suggesting that treating individual normal modes as fundamental heat carriers may not be accurate in amorphous solids. Furthermore, our classical and ab initio molecular dynamics simulations of amorphous silicon demonstrate a large degree of atomic diffusion, especially at high temperatures, leading to the conclusion that thermal transport in amorphous solids could be better described starting from the perspectives of liquid physics rather than from crystalline solids.

36 MATERIALS SCIENCE↗

Can assembly bias explain the lensing amplitude of the BOSS CMASS sample in a Planck cosmology?

ABSTRACT In this paper, we investigate whether galaxy assembly bias can reconcile the 20–40 ${{\ \rm per\ cent}}$ disagreement between the observed galaxy projected clustering signal and the galaxy–galaxy lensing signal in the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample. We use the suite of abacuscosmos lambda cold dark matter simulations at Planck best-fitting cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy–galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20–40 ${{\ \rm per\ cent}}$ overprediction of the observed galaxy–galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fitting cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a $\sim \! 34{{\ \rm per\ cent}}$ discrepancy between the observed projected clustering and galaxy–galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark Energy Survey Year 3 results: Imprints of cosmic voids and superclusters in the Planck CMB lensing map

ABSTRACT The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set (∼5000 deg2), expanding on previous measurements that used Y1 catalogues (∼1300 deg2). Given the increased statistical power compared to Y1 data, we report a 6.6σ detection of negative CMB convergence (κ) imprints using approximately 3600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the ΛCDM model (parameters Ωm = 0.25, σ8 = 0.8), and the discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_{\kappa }=\kappa _{{\rm DES}}/\kappa _{{\rm MICE}}$ to a simulation-based template with fixed shape and found a moderate 2σ deviation in the signal with Aκ ≈ 0.79 ± 0.12. We also examined the WebSky simulation that is based on a Planck 2018 ΛCDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the 8.4σ level; again with a lower-than-expected Aκ = 0.84 ± 0.10 amplitude. The combination of voids and superclusters yields a 10.3σ detection with an Aκ = 0.82 ± 0.08 constraint on the CMB lensing amplitude, thus the overall signal is 2.3σ weaker than expected from MICE.

79 ASTRONOMY AND ASTROPHYSICS↗

Relativistic R -matrix calculations for the electron-impact excitation of neutral molybdenum

A recent PISCES-B Mod experiment [Nishijima et al., J. Phys. B 43, 225701 (2010)] has revealed up to a factor of 5 discrepancy between measurement and the two existing theoretical models [Badnell et al., J. Phys. B 29, 3683 (1996); Bartschat et al., J. Phys. B 35, 2899 (2002)], providing important diagnostics for M O I . In the following paper we address this issue by employing a relativistic atomic structure and R -matrix scattering calculations to improve upon the available models for future applications and benchmark results against a recent Compact Toroidal Hybrid experiment [Hartwell et al., Fusion Sci. Technol. 72, 76 (2017)]. We determine the atomic structure of M O I using GRASP 0 , which implements the multiconfigurational Dirac-Fock method. Fine structure energies and radiative transition rates are presented and compared to existing experimental and theoretical values. The electron-impact excitation of M O I is investigated using the relativistic R -matrix method and the parallel versions of the Dirac atomic R -matrix codes. Electron-impact excitation cross sections are presented and compared to the few available theoretical cross sections. Throughout, our emphasis is on improving the results for the z 5 P$^{o}_{1, 2, 3}$ → a 5 S 2 , z 7 P$^{o}_{2, 3, 4}$ → a 7 S 3 and y 7 P$^{o}_{2, 3, 4}$ → a 7 S 3 electric dipole transitions of particular relevance for diagnostic work.

Smyth, R. T.↗

Neutron measurements from antineutrino hydrocarbon reactions

Charged-current antineutrino interactions on a hydrocarbon scintillator in the MINERvA detector are used to study activity from their final-state neutrons. To ensure that most of the neutrons are from the primary interaction, rather than hadronic reinteractions in the detector, the sample is limited to momentum transfers below 0.8 GeV / c . From 16 129 interactions, 15 246 neutral particle candidates are observed. The reference simulation predicts 64% of these candidates are due to neutrons from the antineutrino interaction directly but also overpredicts the number of candidates by 15% overall. This discrepancy is beyond the standard uncertainty estimates for models of neutrino interactions and neutron propagation in the detector. We explore these two aspects of the models using the measured distributions for energy deposition, time of flight, position, and speed. We also use multiplicity distributions to evaluate the presence of a two-nucleon knockout process. These results provide critical new information toward a complete description of the hadronic final state of neutrino interactions, which is vital to neutrino oscillation experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Critical role of isopotential surfaces for magnetostatic ponderomotive forces

By producing localized wave regions at the ends of an open-field-line magnetic confinement system, ponderomotive walls can be used to differentially confine different species in the plasma. Furthermore, if the plasma is rotating, this wall can be magnetostatic in the laboratory frame, resulting in simpler engineering and better power flow. However, recent work on such magnetostatic walls has shown qualitatively different potentials than those found in the earlier, nonrotating theory. Here, using a simple slab model of a ponderomotive wall, we resolve this discrepancy. We show that the form of the ponderomotive potential in the comoving plasma frame depends on the assumption made about the electrostatic potential in the laboratory frame. If the laboratory-frame potential is unperturbed by the magnetic oscillation, one finds a parallel-polarized wave in the comoving frame, while if each field line remains equipotential throughout the perturbation region, one finds a perpendicularly polarized wave. As a result, this in turn dramatically changes the averaged ponderomotive force experienced by a charged particle along the field line, not only its scaling, but also its direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomistic Simulations of the Elastic Compression of Platinum Nanoparticles

Abstract The elastic behavior of nanoparticles depends strongly on particle shape, size, and crystallographic orientation. Many prior investigations have characterized the elastic modulus of nanoscale particles using experiments or simulations; however their reported values vary widely depending on the methods for measurement and calculation. To understand these discrepancies, we used classical molecular dynamics simulation to model the compression of platinum nanoparticles with two different polyhedral shapes and a range of sizes from 4 to 20 nm, loaded in two different crystal orientations. Multiple standard methods were used to calculate the elastic modulus from stress-vs-strain data for each nanoparticle. The magnitudes and particle-size dependence of the resulting moduli varied with calculation method and, even for larger nanoparticles where bulk-like behavior may be expected, the effective elastic modulus depended strongly on shape and orientation. Analysis of per-atom stress distributions indicated that the shape- and orientation-dependence arise due to stress triaxiality and inhomogeneity across the particle. When the effective elastic modulus was recalculated using a representative volume element in the center of a large nanoparticle, the elastic modulus had the expected value for each orientation and was shape independent. It is only for single-digit nanoparticles that meaningful differences emerged, where even the very center of the particle had a lower modulus due to the effect of the surface. These findings provide better understanding of the elastic properties of nanoparticles and disentangle geometric contributions (such as stress triaxiality and spatial inhomogeneity) from true changes in elastic properties of the nanoscale material.

36 MATERIALS SCIENCE↗

Micromechanics of Composite Tubes by Curved and Flat RVEs

A novel method for analyzing tubular composites was developed by obtaining a curved representative volume element (RVE) from the transformation of a flat RVE. There was found to be a discrepancy in the stress between the two RVE models of approximately 15%. The curved RVE model was validated using an isotropic material.

Thandaga Nagaraju, Hemanth↗

Evaluation and Validation of the n+ 63,65 Cu Cross Sections [Slides]

This presentation discusses the objective of the research which is to update n+ 63,65 Cu cross section evaluations with recently measured data to resolve discrepancies in benchmark performance. It also discusses the models used which are the R -matrix analysis and the analysis of angular distribution coefficients. Additionally, the validation methods that were used are discussed, including the Rez shielding benchmark and the ICSBEP criticality benchmarks. In conclusion, he n+ 63,65 Cu cross sections have been updated via R-matrix analysis up to 100 keV. An increased average capture cross section and the adoption of experimentally based Legendre coefficients lead to improved performance in reactivity benchmarks. In the fast region, the adoption of the JENDL-4.0 cross sections above 4.0 MeV improves the performance in shielding benchmarks. Ultimately, the n 63,65 Cu ENDF files will be submitted to ENDF/B-VIII.1.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic In-Vacuum Field Mapping System for the Muon g-2 Experiment

The Muon g-2 experiment at Fermilab will measure the anomalous magnetic moment of the muon, $a^{}_\mu$ with a precision of 140ppb, a four-fold improvement over the former BNL E821 result. If the current discrepancy of about 3.5$\sigma$ between the Standard Model value for $a^{}_\mu$ and the former experiments was real, then the new experiment would verify the deviation with >7$\sigma$ significance. The experiment uses muons with the magic momentum of 3.09GeV/c stored in a very homogeneous magnetic field of a 45m long storage ring. The determination of $a^{}_\mu$ requires the precise knowledge of both the muon spin precession, $\omega^{}_a$, and the magnetic field, $\omega^{}_p$, expressed in terms of the free-proton Larmor frequency. The field measurement uses the high precision tool of pulsed proton Nuclear Magnetic Resonance and requires cross-calibrating a set of probes: a very accurate NMR probe calibrates the 17 probes on the in-vacuum field mapping system (the so-called trolley). This trolley in turn is used to map the field over the full azimuth of the storage ring. The trolley measurement also provides the calibration of the 378 fixed probes around the ring that constantly monitor the field drift. The trolley is particularly important as it measures the relevant field that is seen by the stored muons and this article presents the details of this field mapping system. This work was supported in part by the US DOE, Fermilab and US DOE OHEP under contract No. KA2201020.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Investigation of Flux Spreading in a Light-Trapping, Planar-Cavity Receiver for Enclosed Solar Particle Heating

Concentrating solar thermal power (CSP) technology development has recently focused on increasing the operating temperatures to accommodate high efficiency power cycles and thermochemical processes. Inert solid particles as heat transfer media enable solar receivers to operate above 700 degrees Celsius resulting in increased system thermal efficiency compared to the conventional molten salt based CSP system. An open-cavity falling-particle solar receiver that can efficiently heat particles by direct heating from concentrated solar radiation faces challenges with large particle losses from wind and unable to support thermochemical reactions. A light-trapping, planar cavity reiver (LTPCR) where particles are indirectly heated can significantly minimize the particle losses during the operation, support thermochemical reactions, and offer scalability potential. The LTPCR features an array of vertical planar receiver/absorber panels arranged within a cavity configuration. Concentrated solar radiation from heliostats is focused onto the receiver walls, where heat is indirectly transferred to solid particles flowing inside the receiver channels. Heat transfer occurs through direct contact between the receiver panel walls and particles, and can be enhanced by fluidizing particles with air. This fluidization increases particle-wall contact and extends particle residence time, maximizing heat transfer efficiency. The unique vertical planar receiver structure originated from a near-blackbody tubular light absorber, effectively distributing the incoming solar beam spread across the panel walls and trapping light. This flux spreading effect, driven by cosine projection, converts high incident solar flux into a lower, more uniform heat flux on the panel walls. This redistribution enhances heat transfer efficiency between particle-wall or reaction gases-wall, while preventing localized overheating of the receiver panel. Indirect planar cavity solar receivers completely separate solid particles from the ambient environment that can greatly reduce the thermal losses in heated particles resulting in high efficiency at high temperatures above 700 degrees Celsius. This design ensures no particle losses to the environment during the operation while open-cavity designs can experience significant particle losses from wind. An experimental investigation was conducted to observe flux spreading on the receiver panel wall. A lab-scale prototype planar receiver, fabricated using Haynes 230 alloy, was tested under direct concentrated solar radiation using the high-flux solar furnace (HFSF) facility at NREL. The experiment was performed under normal peak radiative heat fluxes ranging from 800 to 1900 kW/m2. A temperature distribution on the panel wall was measured using a thermal imaging camera (FLIR A 6600). To prevent overheating at the receiver front tip, prism-shaped heat shields (Zircar UNIFROM C1) were placed in front of the receiver, and their influence on flux spreading was also studied. Absorbed flux distribution on the panel wall was modeled using SolTrace. The total solar power and flux distributions delivered from HFSF were determined based on the heliostat mirror optical properties, direct normal irradiance (DNI) on the on-sun testing days, peak flux measurement during the on-sun testing, and shutter/attenuator settings Due to the large incident angles of the solar beam on the panel wall, the angular optical properties of Haynes 230 alloy and Zircar heat shields were incorporated into the model. This flux distribution model was then integrated into a computational fluid dynamics (CFD) simulation to predict the receiver panel wall temperature, which was compared with the experimental measurements. Both prediction and measurements identified a temperature hotspot at the backside of the panel, indicating that the incident solar beam can fully reach to the rear of the receiver. The heat shields positioned at the front of the receiver effectively reduced the excessive temperature rise at the receiver front tip. Overall, the temperature was well distributed over the panel wall, with a minor hotspot at the back of the receiver. The model slightly overpredicted the temperature, possibly due to discrepancies in optical properties of the panel and an underprediction of thermal loss in the receiver. The advancement of the particle LTPCR offers a viable alternative to open-cavity receivers by addressing particle loss issues. Additionally, it presents a pathway for enabling solar thermochemical processes, extending CSP technology beyond power generation to fuel and chemical production.

14 SOLAR ENERGY↗

An explication of the radiometric method for size and albedo determination

A radiometric model for determination of diameters, albedos, and temperatures of asteroids and other airless planets is presented. Integrated photometry of planets is reviewed, and attention is given to problems of stellar calibration, monochromatic corrections, and corrections for the temperature variations across the planetary sphere. A previous model, which assumes that an asteroid emits infrared radiation as a smooth sphere, is modified to account for observed sunward beaming of infrared emissions, apparently due to craters. In contrast to some studies of the lunar emission, the observable flux ratio between a cratered sphere and a smooth sphere is calculated for large ranges of wavelength, temperature, and phase angle. Revised diameters based on the crater model are given for 84 asteroids. The systematic discrepancy between radiometrically- and polarimetrically-obtained albedos is accounted for; further areas of disagreement between radiometry and polarimetry are discussed, and research to resolve them is suggested.

Hansen, O. L.↗

The composition of ultra-heavy cosmic rays

Preliminary results on the composition of ultraheavy nuclei in the cosmic rays are now available from experiments on HEAO-3 and Ariel-VI satellites. These results are compared with the composition predicted for reasonable propagation models using source abundances similar to those of solar system material. The effects of varying model parameters and of preferential acceleration are discussed. Discrepancies between predicted and observed abundances are larger than can be accounted for solely by first ionization potential effects and are interpreted in terms of a nucleosynthesis process.

Protheroe, R. J.↗

Model-based estimation and prediction of task-imposed mental workload

Mental workload has been an area of intensive research for better than a decade. One specific area of interest in aircrew related workload research is concerned with the development of quantitative indices of workload in aircraft piloting tasks. This paper presents a model-based approach for quantifying mental workload in operational terms. The suggested modeling framework is based on an interpreted Petri net characterization of a task in which 'places' are equated to specific task-related activities and 'transitions' are viewed as internal or external forcing events. It is shown that within this framework quantitative assessments can be made of both cumulative and instantaneous workload associated with the performance of a task and its individual component subtasks. It is suggested that insights gained from analyzing task-specific workload within this modeling paradigm can suggest plausible explanations for reconciling discrepancies between subjectively elicited workload estimates and behavioral/performance measures.

Madni, A. M.↗

A numerical investigation of a subsonic jet in a crossflow

The flowfield induced by a single, subsonic jet exhausting perpendicularly from a flat plate into a subsonic crossflow has been numerically investigated. The test case was chosen to match available experimental data where the jet Mach number was 0.78, and the freestream Mach number was 0.13. Time-averaged solutions were obtained using the thin-layer Navier-Stokes equations and two overlapping grids. The solutions were sensitive to the radial grid clustering near the edge of the jet and to the far-field boundary conditions. Experimental data comparisons were required to determine the most appropriate jet grid and satisfactory boundary conditions. Globally, the solutions converged in about 6000 iterations. The computational results accurately showed the deflected jet and associated contrarotating vortices. The fine clustered grid in the region upstream of the jet exit allowed the horseshoe vortex in the boundary layer near the jet exit to be captured. Most importantly for aircraft applications, the computed plate pressure distributions compared favorably with the experimental data over most of the surface. However, in the wake region immediately downstream of the jet exit, where there is extensive flow separation, some discrepancies with experimental data were observed. Two turbulence models were used in this study: (1) the zero-equation, two layer Baldwin-Lomax turbulence model; and (2) one-equation Baldwin-Barth turbulence model. The turbulence models gave results which generally compared no better with experimental data than the laminar computation results.

Chiu, Stephen H.↗

Reduction of parasitic lasing

A technique was developed which carefully retro-reflects precisely controlled amounts of light back into a laser system thereby intentionally forcing the laser system components to oscillate in a new resonator called the parasitic oscillator. The parasitic oscillator uses the laser system to provide the gain and an external mirror is used to provide the output coupling of the new resonator. Any change of gain or loss inside the new resonator will directly change the lasing threshold of the parasitic oscillator. This change in threshold can be experimentally measured as a change in the absolute value of reflectivity, provided by the external mirror, necessary to achieve lasing in the parasitic oscillator. Discrepancies between experimental data and a parasitic oscillator model are direct evidence of optical misalignment or component performance problems. Any changes in the optical system can instantly be measured as a change in threshold for the parasitic oscillator. This technique also enables aligning the system for maximum parasitic suppression with the system fully operational.

Storm, Mark E.↗

Magnetic Field Suppression of Flow in Semiconductor Melt

One of the most promising approaches for the reduction of convection during the crystal growth of conductive melts (semiconductor crystals) is the application of magnetic fields. Current technology allows the experimentation with very intense static fields (up to 80 KGauss) for which nearly convection free results are expected from simple scaling analysis in stabilized systems (vertical Bridgman method with axial magnetic field). However, controversial experimental results were obtained. The computational methods are, therefore, a fundamental tool in the understanding of the phenomena accounting during the solidification of semiconductor materials. Moreover, effects like the bending of the isomagnetic lines, different aspect ratios and misalignments between the direction of the gravity and magnetic field vectors can not be analyzed with analytical methods. The earliest numerical results showed controversial conclusions and are not able to explain the experimental results. Although the generated flows are extremely low, the computational task is a complicated because of the thin boundary layers. That is one of the reasons for the discrepancy in the results that numerical studies reported. Modeling of these magnetically damped crystal growth experiments requires advanced numerical methods. We used, for comparison, three different approaches to obtain the solution of the problem of thermal convection flows: (1) Spectral method in spectral superelement implementation, (2) Finite element method with regularization for boundary layers, (3) Multiquadric method, a novel method with global radial basis functions, that is proven to have exponential convergence. The results obtained by these three methods are presented for a wide region of Rayleigh and Hartman numbers. Comparison and discussion of accuracy, efficiency, reliability and agreement with experimental results will be presented as well.

Fedoseyev, A. I.↗

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗