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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗
NASA NTRS2015

Beta Testing of CFD Code for the Analysis of Combustion Systems

A preliminary version of OpenNCC was tested to assess its accuracy in generating steady-state temperature fields for combustion systems at atmospheric conditions using three-dimensional tetrahedral meshes. Meshes were generated from a CAD model of a single-element lean-direct injection combustor, and the latest version of OpenNCC was used to calculate combustor temperature fields. OpenNCC was shown to be capable of generating sustainable reacting flames using a tetrahedral mesh, and the subsequent results were compared to experimental results. While nonreacting flow results closely matched experimental results, a significant discrepancy was present between the code's reacting flow results and experimental results. When wide air circulation regions with high velocities were present in the model, this appeared to create inaccurately high temperature fields. Conversely, low recirculation velocities caused low temperature profiles. These observations will aid in future modification of OpenNCC reacting flow input parameters to improve the accuracy of calculated temperature fields.

Temperature↗
NASA NTRS1981

Influence of gravity field uncertainties on the results from POGO and Magsat geomagnetic surveys

Errors in the gravity models used in satellite position calculations are examined as a possible source of the 0 to 100% variance found between POGO and Magsat magnetic data and the extrapolations of aerial magnetic survey data to satellite heights. For POGO data obtained over the New York Bight region using a relatively poor gravity field (a hybrid spherical harmonic model of degree 7 and order 6 with three higher order resonance terms), the magnitude of the error in the satellite height component is found to be sufficient to account for the amplitude of the discrepancy, however the frequency of the quasi-periodic orbital error is too large to explain the localized nature of the differences. For the case of the Magsat satellite, in which a more accurate gravity model was used, it is found that a 30 mgal gravitational anomaly distributed over a 5 x 5 deg area will produce insufficiently large position errors to account for the variations. The agreement between the two sets of satellite data in the New York Bight region suggests either a consistent error in satellite measurements, or problems with the reduction and processing of the aeromagnetic data.

Taylor, P. T.↗
NASA NTRSDate not supplied

Characterization of Titan Entry Radiative Heating in the Low Density Electric Arc Shock Tube

Prediction of aftbody radiative heating for Dragonfly’s entry at Titan requires accurate models to inform the design of its thermal protection systems. Previously reported investigations of shock layer aftbody radiation in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center focused on the lower altitude/higher density part of the trajectory. This paper presents the first optical emission spectroscopy measurements of an ongoing test series in the new aluminum Low Density Shock Tube (LDST), recently installed at EAST, focusing on the radiative heating characterization through the higher altitude portion of Dragonfly’s entry path. Experimental data collected at a fill pressure of 0.1 Torr and shock velocities between 5.2 and 6.2 km/s show favorable agreement to similar conditions previously tested in the 4" High-Velocity Shock Tube, as well as improvements in the signal to noise ratio and available test time. Simulations employing DPLR and NEQAIR underpredict the experimental radiance, with increasing discrepancy from the NIR towards the UV wavelength range, thus motivating further efforts in improving the current non-Boltzmann model for the CN molecule.

A Fagnani↗
NASA NTRS1991

Non-isotropic turbulence effects on spray combustion

A numerical model for the prediction of local properties of statistically stationary spray-combusting flows is evaluated by comparison with experimental data. To appraise the relative performance of turbulence models, computations were carried out by the k-epsilon model and the algebraic stress model. The present numerical results show the qualitative agreement with experimental data. In terms of overall local flow properties, the algebraic stress model improves a degree of conformity to the experimental data due to its ability to introduce the nonisotropic turbulence effects. Two swirl numbers are considered to investigate the influence of swirl on the droplet evaporation and trajectories, and the effects of droplet/turbulence interactions in flow properties. It is found that the large swirl produces a higher evaporationn rate, and more intensive turbulent mixing and burning. The discrepancies observed in the results are attributed mainly to uncertainties in the initial spray size and velocity distributions, the droplet/wall impingement interaction, the combustion model with the fast chemistry and the turbulence models dealing with the strong streamline curvature and complex interactions between the dispersed droplets and the continuous gas-phase flows.

Kim, Y. M.↗
NASA NTRS1985

Charge and energy dependence of the residence time of cosmic ray nuclei below 15 GeV/nucleon

The relative abundance of nuclear species measured in cosmic rays at Earth has often been interpreted with the simple leaky box model. For this model to be consistent an essential requirement is that the escape length does not depend on the nuclear species. The discrepancy between escape length values derived from iron secondaries and from the B/C ratio was identified by Garcia-Munoz and his co-workers using a large amount of experimental data. Ormes and Protheroe found a similar trend in the HEAO data although they questioned its significance against uncertainties. They also showed that the change in the B/C ratio values implies a decrease of the residence time of cosmic rays at low energies in conflict with the diffusive convective picture. These conclusions crucially depend on the partial cross section values and their uncertainties. Recently new accurate cross sections of key importance for propagation calculations have been measured. Their statistical uncertainties are often better than 4% and their values significantly different from those previously accepted. Here, these new cross sections are used to compare the observed B/C+O and (Sc to Cr)/Fe ratio to those predicted with the simple leaky box model.

Soutoul, A.↗
NASA NTRS2016

Data-Driven Radiative Hydrodynamic Modeling of the 2014 March 29 X1.0 Solar Flare

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in Mg II hk, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DSTIBIS) in H(alpha) 6563A and Ca II 8542A, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multithreaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated Mg II line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of approx. 800 km below the transition region can turn the simulated Mg II line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

Da Costa, Fatima Rubio↗
NASA NTRS1997

Evaluation of Tropospheric Water Vapor Simulations from the Atmospheric Model Intercomparison Project

Simulations of humidity from 28 general circulation models for the period 1979-88 from the Atmospheric Model Intercomparison Project are compared with observations from radiosondes over North America and the globe and with satellite microwave observations over the Pacific basin. The simulations of decadal mean values of precipitable water (W) integrated over each of these regions tend to be less moist than the real atmosphere in all three cases; the median model values are approximately 5% less than the observed values. The spread among the simulations is larger over regions of high terrain, which suggests that differences in methods of resolving topographic features are important. The mean elevation of the North American continent is substantially higher in the models than is observed, which may contribute to the overall dry bias of the models over that area. The authors do not find a clear association between the mean topography of a model and its mean W simulation, however, which suggests that the bias over land is not purely a matter of orography. The seasonal cycle of W is reasonably well simulated by the models, although over North America they have a tendency to become moister more quickly in the spring than is observed. The interannual component of the variability of W is not well captured by the models over North America. Globally, the simulated W values show a signal correlated with the Southern Oscillation index but the observations do not. This discrepancy may be related to deficiencies in the radiosonde network, which does not sample the tropical ocean regions well. Overall, the interannual variability of W, as well as its climatology and mean seasonal cycle, are better described by the median of the 28 simulations than by individual members of the ensemble. Tests to learn whether simulated precipitable water, evaporation, and precipitation values may be related to aspects of model formulation yield few clear signals, although the authors find, for example, a tendency for the few models that predict boundary layer depth to have large values of evaporation and precipitation. Controlled experiments, in which aspects of model architecture are systematically varied within individual models, may be necessary to elucidate whether and how model characteristics influence simulations.

Gaffen, Dian J.↗
DOE OSTI2021

High frequency impedance characteristics of a tunable microplasma device

Computational studies on high frequency impedance characteristics of a microplasma device are reported. While microplasma is ignited using a primary excitation signal, frequency response of plasma impedance is determined by a secondary high frequency probe signal with significantly lower voltage amplitude such that it does not influence the plasma parameters. Here, the computational model utilizing the drift–diffusion approximation is first validated by comparing with experimental data for microplasmas ignited at pressures ranging from 1 to 5 Torr. In spite of quantitative discrepancies, good overall agreement is obtained between the measured frequency response of impedance of the discharge. Comparisons are also presented for various plasma parameters including mean electron number density, sheath thickness, mean electron temperature, and collision frequency that were inferred from the impedance measurements. The computational model is then used to perform simulations of near-atmospheric pressure microplasmas with the probe signal frequency ranging from 3 to 20 GHz. The simulations demonstrate the presence of a resonance frequency at which the impedance vanishes. More importantly, it is shown that this resonant frequency can be tuned effectively by suitably modifying the operating parameters (gap size, pressure, and excitation voltage). The simulated impedance characteristics are used to determine the effective plasma inductance and capacitance using a non-linear fitting approach, thereby showing the dependence of these electrical parameters on the plasma operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

Chandra Contaminant Migration Model

High volatility cleans OBFs and low volatility produces a high build-up at OBF centers; only a narrow (factor of 2 or less) volatility range produces the observed spatial pattern. Simulations predict less accumulation above outer S-array CCDs; this may explain, in part, gratings/imaging C/MnL discrepancies. Simulations produce a change in center accumulation due solely to DH heater ON/OFF temperature change; but a 2nd contaminant and perhaps a change in source rate is also required. Emissivity E may depend on thickness; another model parameter. Additional physics, e.g., surface migration, is not warranted at this time. At t approx. 14 yrs, model produced 0.22 grams of contaminant, 0.085 grams remaining within ACIS cavity; 7 percent (6mg) on OBFs.

Swartz, Douglas A.↗
NASA NTRS1999

Aerothermodynamic Calculations on X-34 at Mach 6 Wind Tunnel Conditions

The effects of Reynolds number and turbulence on surface heat-transfer rates are numerically investigated for a 0.015 scale X-34 vehicle at wind tunnel conditions. Laminar heating rates, non-dimensionalized by Fay-Riddell stagnation heating, do not change appreciably with an order of magnitude variation in Reynolds number. Modeling a turbulent versus laminar boundary layer at the same Reynolds number increases the windside heating by a factor of four, portions on the leeside by a factor of two, and causes a 30 percent increase in wing leading edge heating. A discrepancy between laminar and turbulent heating trends on the windside centerline is explained by the presence of attached windside vortices in the laminar solutions, structures that are inhibited by the turbulence modeling.

Wood, William A.↗
NASA NTRS1988

Zodiacal emission. I - Dust near the earth's orbit

The infrared emission of interplanetary dust near the earth's orbit is derived from IRAS observations of the gradient of in-ecliptic brightness tangent to the earth's orbit, and the annual variation of the ecliptic polar brightness. Models with five grain constituents and three size distributions are compared with the observations. The observed emission is twice as bright as predicted; this discrepancy is due either to calibration errors or to enhanced radiative efficiency of 'fluffy' particles. Graphite and magnetite particles are ruled out because they are too hot. The size distribution is constrained to be less steep than that derived from lunar microcrater studies, and is consistent with that obtained by earth-orbiting satellites. Nonhomogeneous silicate grains with 3 percent graphite impurity produce the best fit to the spectrum. The model predictions extend from 3 microns to 1 mm, for use in analyzing future infrared background observations.

Reach, William T.↗
NASA NTRS1995

Parametrics on 2D Navier-Stokes analysis of a Mach 2.68 bifurcated rectangular mixed-compression inlet

The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a two-dimensional (2D) Navier-Stokes flow solver. Parametric studies were performed on turbulence models, computational grids and bleed models. The computer flowfield was substantially different from the original inviscid design, due to interactions of shocks, boundary layers, and bleed. Good agreement with experimental data was obtained in many aspects. Many of the discrepancies were thought to originate primarily from 3D effects. Therefore, a balance should be struck between expending resources on a high fidelity 2D simulation, and the inherent limitations of 2D analysis. The solutions were fairly insensitive to turbulence models, grids and bleed models. Overall, the k-e turbulence model, and the bleed models based on unchoked bleed hole discharge coefficients or uniform velocity are recommended. The 2D Navier-Stokes methods appear to be a useful tool for the design and analysis of supersonic inlets, by providing a higher fidelity simulation of the inlet flowfield than inviscid methods, in a reasonable turnaround time.

Mizukami, M.↗
DOE OSTI2022

Fast Reactor Physics Model Verification Studies using ARC and PyARC Workflows

PyARC was recently developed at Argonne National Laboratory to automate many of the tasks required in the ARC (Argonne Reactor Computation) fast reactor simulation workflow, from input file generation, code execution, data transfer between ARC codes, and output postprocessing. PyARC will likely be the path forward to train new users of the ARC codes with the goal of wide adoption by the national laboratories, academia, and industry. In particular, for the ANL-JAEA collaboration under the Civil Nuclear Working Group (CNWG) project agreement NE-01, PyARC will be used to model the Joyo and EBR-II reactors for comparisons with measured data and calculated results from JAEA (Task 3: Fast Reactor Fuel and Core). As an additional avenue for verification and validation, this report investigates the use of PyARC towards a variety of existing ARC-based reactor models, in order to understand its efficacy in replicating the behavior of base ARC codes and better understand any limitations within modeling realistic fast reactor problems. To this end, PyARC was used to model the Joyo MKI, RBEC Benchmark-M, PRISM Mod-B, and EBR-II Run 138B cores, and its results were compared to those from existing ARC-based models. It was found that for hexagonal-based geometries PyARC was able to replicate the behavior of ARC codes to within 10 pcm for small reactor cores, and ~150pcm difference in eigenvalue for larger cores. These discrepancies are attributed primarily to differences in local mesh refinement options between ARC and PyARC, which currently cannot be resolved with PyARC’s latest version (1.6.0). In some of these cases, PyARC was used to model steady-state problems with initial core compositions originating from a prior REBUS depletion calculation. While PyARC was not designed to support such steady-state calculations, workarounds were applied to replicate the behavior of ARC-based calculations as closely as possible. Thus, these results demonstrate the wide extent to which they can be applied to fast reactor problems while still providing immense benefit to the user in terms of automating and standardizing common routines within the fast reactor analysis workflow. This study concluded that PyARC will be suitable for modeling the steady-state conditions of the EBR-II and Joyo fast reactors as part of the CNWG project agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗
DOE OSTI2019

Development of the multigroup cross section library for the CASL neutronics simulator MPACT: Method and procedure

The MPACT neutronics module of the Consortium for Advanced Simulation of Light Water Reactors (CASL) core simulator is a 3D whole core transport code being developed for the CASL toolset known as the Virtual Environment for Reactor Analysis (VERA). Key characteristics of the MPACT code include a subgroup method for resonance self-shielding and a whole-core transport solver with a 2D/1D synthesis method. The AMPX and SCALE code packages have been significantly improved to include intermediate resonance parameters, as well as new modules to enhance the accuracy of resonance data and scattering matrices for the Bondarenko approach and transport cross sections of hydrogen. In addition, CASL XSTools have been developed to generate the MPACT multigroup cross section library. A simple super-homogenization method for 238U has been developed to resolve reaction rate discrepancy issues caused by angle-dependent total cross sections, poor scattering matrices, and a poor resonance interference model at coarse group structure. Here, due to the limitation of computing capacity that can occur even in high performance computing, a new 51-group structure has been developed for efficient simulation that is applicable for both pressurized water reactor (PWR) and boiling water reactor (BWR) simulations. A new ENDF/B-VII.1 MPACT 51-group library has been developed using the AMPX/SCALE and XSTools code packages, and it has been verified and validated by performing various benchmark calculations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Quantitative Evaluation of Nanosecond Pulsed Laser-Induced Photomodification of Plasmonic Gold Nanoparticles

The rapid growth of gold nanoparticle applications in laser therapeutics and diagnostics has brought about the need for establishing innovative standardized test methods for evaluation of safety and performance of these technologies and related medical products. Furthermore, given the incomplete and inconsistent data on nanoparticle photo-modification thresholds provided in the literature, further elucidation of processes that impact the safety and effectiveness of laser-nanoparticle combination products is warranted. Therefore, we present a proof-of-concept study on an analytical experimental test methodology including three approaches (transmission electron microscopy, dynamic light scattering, and spectrophotometry) for experimental evaluation of damage thresholds in nanosecond pulsed laser-irradiated gold nanospheres, and compared our results with a theoretical model and prior studies. This thorough evaluation of damage threshold was performed based on irradiation with a 532nm nanosecond-pulsed laser over a range of nanoparticle diameters from 20 to 100nm. Experimentally determined damage thresholds were compared to a theoretical heat transfer model of pulsed laser-irradiated nanoparticles and found to be in reasonably good agreement, although some significant discrepancies with prior experimental studies were found. This study and resultant dataset represent an important foundation for developing a standardized test methodology for determination of laser-induced nanoparticle damage thresholds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Possible octupole deformation of 208 Pb and the ultracentral v 2 to v 3 puzzle

Recent measurements have established the sensitivity of ultracentral heavy-ion collisions to the deformation parameters of nonspherical nuclei. In the case of 129 Xe collisions, a quadrupole deformation of the nuclear profile led to an enhancement of elliptic flow in the most central collisions. In 208 Pb collisions a discrepancy exists in similar centralities, where either elliptic flow is overpredicted or triangular flow is underpredicted by hydrodynamic models; this is known as the v 2 -to-v 3 puzzle in ultracentral collisions. Motivated by low-energy nuclear structure calculations, we consider the possibility that 208Pb nuclei could have a pear-shape deformation (octupole), which has the effect of increasing triangular flow in central PbPb collisions. Using the recent data from ALICE and ATLAS, we reexamine the v 2 -to-v 3 puzzle in ultracentral collisions, including new constraints from recent measurements of the triangular cumulant ratio v 3 {4}/v 3 {2} and comparing two different hydrodynamic models. We find that while an octupole deformation would slightly improve the ratio between v 2 and v 3 , it is at the expense of a significantly worse triangular flow cumulant ratio. In fact, the latter observable prefers no octupole deformation, with β 3 ≲ 0.0375 for 208 Pb, and is therefore consistent with the expectation for a doubly-magic nucleus even at top collider energies. Furthermore, the v 2 -to-v 3 puzzle remains a challenge for hydrodynamic models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗