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 1,297 records · Page 72

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation and benchmarking of two particle-in-cell codes for a glow discharge

The two particle-in-cell codes EDIPIC and LSP are benchmarked and validated for a parallel-plate glow discharge in helium, in which the axial electric field had been carefully measured, primarily to investigate and improve the fidelity of their collision models. The scattering anisotropy of electron-impact ionization, as well as the value of the secondary-electron emission yield, are not well known in this case. The experimental uncertainty for the emission yield corresponds to a factor of two variation in the cathode current. If the emission yield is tuned to make the cathode current computed by each code match the experiment, the computed electric fields are in excellent agreement with each other, and within about 10% of the experimental value. The non-monotonic variation of the width of the cathode fa ll with the applied voltage seen in the experiment is reproduced by both codes. The electron temperature in the negative glow is within experimental error bars for both codes, but the density of slow trapped electrons is underestimated. A more detailed code comparison don e for several synthetic cases of electron-beam injection into helium gas shows that the codes are in excellent agreement for ionization rate, as well as for elastic and excitation collisions with isotropic scattering pattern. The remaining significant discrepancies between the two codes are due to differences in their electron binary-collision models, and for anisotropic scattering due to elastic and excitation collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Solar Wind

Parker's model of a spherically expanding corona, the "solar wind," is compared with D. E. Blackwell's observations of the 1954 minimum equatorial corona. A significant discrepancy is found between the predicted and the observed electron densities at distances from the sun greater than 20 solar radii. Blackwell's data are found to be consistent with a model in which the corona expands mostly within a disk less than 25 solar radii thick, lying within the sun's equatorial plane. The thickness of the disk as a function of distance from the sun is qualitatively explained in terms of magnetic pressure. The solar wind is found to have a considerable effect on the lunar atmosphere. First, the calculated density of the lunar atmosphere is greatly reduced by collisions with protons in the solar wind. If the flux of particles in this wind has the conventional values ranging between 10(exp 10) to 10(exp 11) per sq cm-sec, the calculations yield a lunar pressure of 10(exp -13) atmosphere of argon, in agreement with the value predicted by Elsmore and Whitfield on the basis of observations on the occultation of radio stars. Second, following a suggestion by Gold, it was found that the collisions of solar-wind protons with the lunar surface produce an atmosphere of cold neutral hydrogen with a density of 10(exp 5) per cu cm at the lunar surface. The density falls off at greater distances in accordance with the inverse-square law. Estimates indicate that the interaction of solar particles with the neutral hydrogen will produce an extended lunar ionosphere with a density of the order of 400 protons/cu cm in the vicinity of the moon.

Herring, J. R.↗

Quantitative Studies of the Optical and UV Spectra of Galactic Early B Supergiants

We undertake an optical and ultraviolet spectroscopic analysis of a sample of 20 Galactic B0-B5 supergiants of luminosity classes Ia, Ib, Iab, and II. Fundamental stellar parameters are obtained from optical diagnostics and a critical comparison of the model predictions to observed UV spectral features is made. Methods. Fundamental parameters (e.g., T(sub eff), log L(sub *), mass-loss rates and CNO abundances) are derived for individual stars using CMFGEN, a nLTE, line-blanketed model atmosphere code. The impact of these newly derived parameters on the Galactic B supergiant Ten scale, mass discrepancy, and wind-momentum luminosity relation is examined. Results. The B supergiant temperature scale derived here shows a reduction of about 1000-3000 K compared to previous results using unblanketed codes. Mass-loss rate estimates are in good agreement with predicted theoretical values, and all of the 20 BO-B5 supergiants analysed show evidence of CNO processing. A mass discrepancy still exists between spectroscopic and evolutionary masses, with the largest discrepancy occuring at log (L/(solar)L approx. 5.4. The observed WLR values calculated for B0-B0.7 supergiants are higher than predicted values, whereas the reverse is true for B1-B5 supergiants. This means that the discrepancy between observed and theoretical values cannot be resolved by adopting clumped (i.e., lower) mass-loss rates as for O stars. The most surprising result is that, although CMFGEN succeeds in reproducing the optical stellar spectrum accurately, it fails to precisely reproduce key UV diagnostics, such as the N v and C IV P Cygni profiles. This problem arises because the models are not ionised enough and fail to reproduce the full extent of the observed absorption trough of the P Cygni profiles. Conclusions. Newly-derived fundamental parameters for early B supergiants are in good agreement with similar work in the field. The most significant discovery, however, is the failure of CMFGEN to predict the correct ionisation fraction for some ions. Such findings add further support to revising the current standard model of massive star winds, as our understanding of these winds is incomplete without a precise knowledge of the ionisation structure and distribution of clumping in the wind. Key words. techniques: spectroscopic - stars: mass-loss - stars: supergiants - stars: abundances - stars: atmospheres - stars: fundamental parameters

Searle, S. C.↗

IAEA Photonuclear Data Library 2019

Here, photo-induced reaction cross section data are of importance for a variety of current or emerging applications, such as radiation shielding design and radiation transport analyses, calculations of absorbed dose in the human body during radiotherapy, physics and technology of fission reactors (influence of photo-reactions on neutron balance) and fusion reactors (plasma diagnostics and shielding), activation analyses, safeguards and inspection technologies, nuclear waste transmutation, medical isotope production and astrophysical applications. To address these data needs the IAEA Photonuclear Data library was produced in 1999, containing evaluated photo-induced cross sections and neutron spectra for 164 nuclides which were deemed relevant for the applications. Since the release of the IAEA Photonuclear Data Library however, new experimental data as well as new methods to assess the reliability of experimental cross sections have become available. Theoretical models and input parameters used to evaluate photo-induced reactions have improved significantly over the years. In addition, new measurements of partial photoneutron cross sections using mono-energetic photon beams and advanced neutron detection systems have been performed allowing for the validation of the evaluations and assessments of the experimental data. Furthermore, technological advances have led to the construction of new and more powerful gamma-beam facilities, therefore new data needs are emerging. We report our coordinated efforts to address these data needs and present the results of the new up-to-date evaluations included in the new updated IAEA Photonuclear Data Library consisting of 219 nuclides. The new library includes 188 new evaluations produced by the CRP evaluators, and one evaluation taken from the JENDL/PD-2016 library, while 20 evaluations were retained from the previous 1999 IAEA Photonuclear Data Library. In most of the cases, the photon energy goes up to 200 MeV. A total of 55 nuclides are new in this library reflecting the progress in measurements but also the developing data needs. In this paper we discuss the new assessment method and make recommendations to the user community in cases where the experimental data are discrepant and the assessments disagree. In addition, in the absence of experimental data, we present model predictions for photo-induced reaction cross section on nuclides of potential interest to medical radioisotope production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An Assessment of the Drag Models in the Case of a Shock Interacting With a Fixed Bed of Point Particles

Abstract In this work, three-dimensional Euler–Lagrange (EL) point-particle simulations of a shock wave interacting with a fixed bed of particles are carried out. The results from the particle-resolved (PR) simulations are used to assess the performance of the point-particle drag models during short time scales. We demonstrate that in a one-way coupled regime, the point-particle simulations recover the dominant gas dynamic features of the flow and are in a good agreement with the exact Riemann solution of a shock traveling through a sudden area contraction. Although the PR simulations are inviscid, we show that a dissipative drag is necessary to predict the mean behavior of the gas. As a model for the inviscid shock-induced (SI) drag two different models are presented in lieu of the quasi-steady drag. Finally, two-way coupled simulations are performed at four different particle volume fractions {0.10, 0.15, 0.20, 0.25} and three different incident shock Mach numbers {1.22, 1.66, 3.0} and compared against the data from PR inviscid simulations. At a lower Mach number (1.22), averaged flow quantities from the two-way coupled simulations agree well with the PR simulations. As the Mach number increases, we observe that the discrepancies between the point-particle and the PR simulations grow. A sensitivity analysis of the drag models involved reveals a strong influence of the inviscid-unsteady force on the gas quantities especially in the case of a strong shock interacting with a dense bed of particles. The use of Mach correlation beyond the subcritical regime coupled with the model for volume fraction correction is identified as a probable cause for the additional drag.

Engineering↗

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗

Global long-range transport and lung cancer risk from polycyclic aromatic hydrocarbons shielded by coatings of organic aerosol

Polycyclic aromatic hydrocarbons (PAHs) have toxic impacts on ecosystems and human health. Laboratory measurements show that one of the most carcinogenic PAHs, benzo(a)pyrene, which is adsorbed on surfaces of soot particles, reacts very quickly with atmospheric oxidants like ozone within ~2 hours. Yet, field observations indicate that it actually persists for much longer periods in the atmosphere, and this large discrepancy is not well understood. Driven by novel experimental understanding, we develop a new modelling approach, whereby particle-bound BaP is shielded from oxidation by a coating of viscous organic aerosol (OA). We show that application of this new approach in a global climate model leads to higher atmospheric BaP concentrations that agree much better with measurements, compared to the default model, as well as stronger long-range transport and greater deposition fluxes. Here, this new approach also predicts elevated lung-cancer risk from PAHs. Predicted oxidation of BaP is highest over a tropical belt where OA is liquid-like.

54 ENVIRONMENTAL SCIENCES↗

Speculations on the W-mass measurement at CDF*

Abstract The W mass determination at the Tevatron CDF experiment reported a deviation from the SM expectation at the 7 σ level. We discuss a few possible interpretations and their collider implications. We perform electroweak global fits under various frameworks and assumptions. We consider three types of electroweak global fits in the effective-field-theory framework: the S - T , S - T - , and eight-parameter flavor-universal one. We discuss the amounts of tensions between different measurements reflected in these fits and the corresponding shifts in central values of these parameters. With these electroweak fit pictures in hand, we present a few different classes of models and discuss their compatibility with these results. We find that while explaining the discrepancy, the single gauge boson extensions face strong LHC direct search constraints unless the is fermiophobic (leptophobic), which can be realized if extra vector fermions (leptons) mix with the SM fermions (leptons). Vector-like top partners can partially generate the needed shift to the electroweak observables. The compatibility with the top squark is also studied in detail. We find that the non-degenerate top squark soft masses enhance the needed operator coefficients, enabling an allowed explanation compatible with current LHC measurements. Overall, more theoretical and experimental developments are highly in demand to reveal the physics behind this discrepancy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Validation of Tendril TrueHome Using Software-to-Software Comparison

This study performed comparative evaluation of EnergyPlus version 8.6 and Tendril TrueHome, two physics-based home energy simulation models, to identify differences in energy consumption predictions between the two programs and resolve discrepancies between them. EnergyPlus is considered a benchmark, best-in-class software tool for building energy simulation. This exercise sought to improve both software tools through additional evaluation/scrutiny.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Satellite measurements of ionic concentrations applied to low altitude incoherent scatter interpretations

The determination of ion and electron temperatures and electron density from incoherent scatter radar data has in the past involved assumptions concerning the ionospheric composition below 250 km. Using a very large data base of measurements of O(+), O2(+), NO(+), and N2(+) made by the Atmosphere Explorer C satellite, a model has been developed of the ionic concentration between 130 and 300 km as a function of solar zenith angle. As the effect of the new model is to increase the radar determinations of the temperature, this may explain previous discrepancies between Langmuir probe and incoherent scatter temperatures, where the comparisons have been made at these low altitudes

Torr, M. R.↗

Chlorine chemistry in diffuse interstellar clouds

The results of recent observations of the column densities of chlorine species in diffuse interstellar clouds are inconsistent with the column densities predicted by model calculations. The use of new atomic and molecular data fails to decrease the discrepancy. It is speculated that the exothermic reaction Cl(+) yields HCl(+) + H may not be rapid at the low temperatures prevailing in diffuse interstellar clouds.

Black, J. H.↗

On the importance of interstellar helium for the propagation of heavy cosmic rays

The influence of interstellar He on the fragmentation of heavy cosmic rays in the interstellar medium (ISM) has long been a controversial subject. While H-induced cross section data are now avialable over broad mass and energy ranges, little data for He-induced fragmentation exists. With the recent reports of accurate measurements of the secondary/primary ratios in cosmic rays and of H-induced cross sections the problem of including interstellar He in propagation calculations becomes even more critical. As is argued the escape lengths lambda e deduced from the B/C+) and Sc-Cr/Fe ratios cannot be reconciled within the frame of a simple leaky box model assuming the ISM composed of pure H. It is quite remarkable that the discrepancy is especially large in the GeV region where (1) secondary/primary ratios measured by several groups agree fairly well and (2) fragmentation cross sections have been recently measured with good accuracy.

Ferrando, P.↗

Comparisons of SiO maser and long-period variable positions in the R Aquarii and Omicron Ceti binary systems

It has been determined that the absolute position of the centroid of SiO maser-emitting spots toward both R Aqr and Omicron Cet are coincident with the position of the long-period variables (LPVs) in these binary systems to within the errors of measurement. The SiO positions were determined with the Hat Creek interferometer, while the LPV positions were determined with the 8-in transit circle of the US Naval Observatory. These results contradict an earlier report of a circumbinary SiO maser far removed from the LPV in the R Aqr binary system; statistical reasons for the discrepancy are suggested. High-resolution spectra of both sources are presented and possible models are discussed.

Hollis, J. M.↗

Eclipse maps of the accretion disk in the classical nova V Persei

The surface brightness distribution in the accretion disk of V Per is reconstructed using a maximum entropy technique. The resulting brightness temperature gradient is much less steep than that from a standard steady state, optically thick disk for all reasonable mass accretion rates. A model including a rim which obscures part of the disk does not resolve this discrepancy. However, if the inner disk is disrupted out to about 0.15-0.25R sub L1 (for example by a magnetic field), the brightness temperature gradient can be very close to T sub b varies as R exp -3/4 (1-b(Rin/R)exp 1/2) law. The value of Rin needed to match this steady state distribution depends on the distance to the system, its size and brightness, and the rotation rate of the white dwarf.

Wood, Janet H.↗

Modeling of turbulent supersonic H2-air combustion with a multivariate beta PDF

Recent calculations of turbulent supersonic reacting shear flows using an assumed multivariate beta PDF (probability density function) resulted in reduced production rates and a delay in the onset of combustion. This result is not consistent with available measurements. The present research explores two possible reasons for this behavior: use of PDF's that do not yield Favre averaged quantities, and the gradient diffusion assumption. A new multivariate beta PDF involving species densities is introduced which makes it possible to compute Favre averaged mass fractions. However, using this PDF did not improve comparisons with experiment. A countergradient diffusion model is then introduced. Preliminary calculations suggest this to be the cause of the discrepancy.

Baurle, R. A.↗

Extraterrestrial accretion and glacial cycles

We propose that the approx. 100-k.y. cycle seen in terrestrial glaciation is due to changes in meteor flux that come from changes in the Earth's orbit. This model can explain a 70-k.y. 'anomalous' period in climate data and the apparent discrepancy between present extraterrestrial fluxes and those in oceanic sediments. It can be tested by measuring Ir densities in sediments and ice during glacials and interglacials.

Muller, R. A.↗