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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Accurate Time-Dependent Traveling-Wave Tube Model Developed for Computational Bit-Error-Rate Testing

The phenomenal growth of the satellite communications industry has created a large demand for traveling-wave tubes (TWT's) operating with unprecedented specifications requiring the design and production of many novel devices in record time. To achieve this, the TWT industry heavily relies on computational modeling. However, the TWT industry's computational modeling capabilities need to be improved because there are often discrepancies between measured TWT data and that predicted by conventional two-dimensional helical TWT interaction codes. This limits the analysis and design of novel devices or TWT's with parameters differing from what is conventionally manufactured. In addition, the inaccuracy of current computational tools limits achievable TWT performance because optimized designs require highly accurate models. To address these concerns, a fully three-dimensional, time-dependent, helical TWT interaction model was developed using the electromagnetic particle-in-cell code MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm). The model includes a short section of helical slow-wave circuit with excitation fed by radiofrequency input/output couplers, and an electron beam contained by periodic permanent magnet focusing. A cutaway view of several turns of the three-dimensional helical slow-wave circuit with input/output couplers is shown. This has been shown to be more accurate than conventionally used two-dimensional models. The growth of the communications industry has also imposed a demand for increased data rates for the transmission of large volumes of data. To achieve increased data rates, complex modulation and multiple access techniques are employed requiring minimum distortion of the signal as it is passed through the TWT. Thus, intersymbol interference (ISI) becomes a major consideration, as well as suspected causes such as reflections within the TWT. To experimentally investigate effects of the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers, in addition to acquiring the required digital hardware. As an alternative, the time-domain TWT interaction model developed here provides the capability to establish a computational test bench where ISI or bit error rate can be simulated as a function of TWT operating parameters and component geometries. Intermodulation products, harmonic generation, and backward waves can also be monitored with the model for similar correlations. The advancements in computational capabilities and corresponding potential improvements in TWT performance may prove to be the enabling technologies for realizing unprecedented data rates for near real time transmission of the increasingly larger volumes of data demanded by planned commercial and Government satellite communications applications. This work is in support of the Cross Enterprise Technology Development Program in Headquarters' Advanced Technology & Mission Studies Division and the Air Force Office of Scientific Research Small Business Technology Transfer programs.

Kory, Carol L.↗
DOE OSTI2018

Evaluation of models determined by neutron diffraction and proposed improvements to their validation and deposition

The Protein Data Bank (PDB) contains a growing number of models that have been determined using neutron diffraction or a hybrid method that combines X-ray and neutron diffraction. The advantage of neutron diffraction experiments is that the positions of all atoms can be determined, including H atoms, which are hardly detectable by X-ray diffraction. This allows the determination of protonation states and the assignment of H atoms to water molecules. Because neutrons are scattered differently by hydrogen and its isotope deuterium, neutron diffraction in combination with H/D exchange can provide information on accessibility, dynamics and chemical lability. In this study, the deposited data, models and model-to-data fit for all PDB entries that used neutron diffraction as the source of experimental data have been analysed. In many cases, the reported R work and R free values were not reproducible. In such cases, the model and data files were analysed to identify the reasons for this mismatch. The issues responsible for the discrepancies are summarized and explained. The analysis unveiled limitations to the annotation, deposition and validation of models and data, and a lack of community-wide accepted standards for the description of neutron models and data, as well as deficiencies in current model refinement tools. Most of the issues identified concern the handling of H atoms. Since the primary use of neutron macromolecular crystallography is to locate and directly visualize H atoms, it is important to address these issues, so that the deposited neutron models allow the retrieval of the maximum amount of information with the smallest effort of manual intervention. A path forward to improving the annotation, validation and deposition of neutron models and hybrid X-ray and neutron models is suggested.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

CFD-DEM study of bubble properties in a cylindrical fluidized bed of Geldart Group D particles and comparison with prior MRI data

Numerical simulations of a freely bubbling cylindrical fluidized bed are carried out using a coupled computational fluid dynamics and discrete element method (CFD-DEM) model and compared to recent experimental data. The experiments were conducted using high-resolution and high-frequency magnetic resonance imaging providing high-fidelity data of the bubbling within a central 10 mm slice of the bed. Qualitatively, we nd more regular (i.e., less chaotic) structures in the simulated beds than observed experimentally. Quantitatively, however, the bubble diameter and number of bubbles as a function of height within the bed is predicted well by the base model. Unfortunately, the regularity in the simulations manifests as a considerable discrepancy in the speed of the (dense) emulsion phase. The simulated velocity probability distribution functions show an accumulation of low-speed regions and deficiency of high-speed regions. A simple parametric study of the base model is also carried out considering many of the most common CFD-DEM modeling parameters. It is found that the fluid grid size, geometry resolution, transfer kernel and drag law did not have a significant effect on bubble or particle dynamics.

42 ENGINEERING↗
DOE OSTI2017

Shock Acceleration Model for the Toothbrush Radio Relic

Although many of the observed properties of giant radio relics detected in the outskirts of galaxy clusters can be explained by relativistic electrons accelerated at merger-driven shocks, significant puzzles remain. In the case of the so-called Toothbrush relic, the shock Mach number estimated from X-ray observations (M{sub X}≈1.2--1.5) is substantially weaker than that inferred from the radio spectral index (M{sub rad}≈2.8). Toward understanding such a discrepancy, we here consider the following diffusive shock acceleration (DSA) models: (1) weak-shock models with M{sub s}≲2 and a preexisting population of cosmic-ray electrons (CRe) with a flat energy spectrum, and (2) strong-shock models with M{sub s}≈3 and either shock-generated suprathermal electrons or preexisting fossil CRe. We calculate the synchrotron emission from the accelerated CRe, following the time evolution of the electron DSA, and the subsequent radiative cooling and postshock turbulent acceleration (TA). We find that both models could reproduce reasonably well the observed integrated radio spectrum of the Toothbrush relic, but the observed broad transverse profile requires the stochastic acceleration by downstream turbulence, which we label “turbulent acceleration” or TA to distinguish it from DSA. Moreover, to account for the almost uniform radio spectral index profile along the length of the relic, the weak-shock models require a preshock region over 400 kpc with a uniform population of preexisting CRe with a high cutoff energy (≳40 GeV). Due to the short cooling time, it is challenging to explain the origin of such energetic electrons. Therefore, we suggest the strong-shock models with low-energy seed CRe (≲150 MeV) are preferred for the radio observations of this relic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

An improved model for the accurate calculation of parallel heat fluxes at the JET bulk tungsten outer divertor

Parallel heat flux calculations at the JET divertor have been based on the assumption that all incoming heat is due to the projection of the heat flux parallel to the magnetic line, q ∥ , plus a constant background. This simplification led to inconsistencies during the analysis of a series of dedicated tungsten melting experiments performed in 2013, for which infrared (IR) thermography surface measurements could not be recreated through simulations unless the parallel heat flux was reduced by 80% for L-mode and 60% for H-mode. We give an explanation for these differences using a new IR inverse analysis code, a set of geometrical corrections, and most importantly an additional term for the divertor heat flux accounting for non-parallel effects such as cross-field transport, recycled neutrals or charge exchange. This component has been evaluated comparing four different geometries with impinging angles varying from 2 to 90°. Its magnitude corresponds to 1.2%–1.9% of q ∥ , but because it is not affected by the magnetic projection, it accounts for up to 20%–30% of the tile surface heat flux. The geometrical corrections imply a further reduction of 24% of the measured heat flux. In addition, the application of the new inverse code increases the accuracy of the tile heat flux calculation, eliminating any previous discrepancy. The parallel heat flux computed with this new model is actually much lower than previously deduced by inverse analysis of IR temperatures—40% for L-mode and 50% for H-mode—while being independent of the geometry on which it is measured. This main result confirms the validity of the optical projection as long as a non-constant and non-parallel component is considered. For a given total heating power, the model predicts over 10% reduction of the maximum tile surface heat flux compared to strict optical modelling, as well as a 30% reduced sensitivity to manufacturing and assembling tolerances. Furthermore these conclusions, along with the improvement in the predictability of the divertor thermal behaviour, are critical for JET future DT operations, and are also directly applicable to the design of the ITER divertor monoblocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Surprising similarities in model and observational aerosol radiative forcing estimates

Abstract. The radiative forcing from aerosols (particularly through their interaction with clouds) remains one of the most uncertain components of the human forcing of the climate. Observation-based studies have typically found a smaller aerosol effective radiative forcing than in model simulations and were given preferential weighting in the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). With their own sources of uncertainty, it is not clear that observation-based estimates are more reliable. Understanding the source of the model and observational differences is thus vital to reduce uncertainty in the impact of aerosols on the climate. These reported discrepancies arise from the different methods of separating the components of aerosol forcing used in model and observational studies. Applying the observational decomposition to global climate model (GCM) output, the two different lines of evidence are surprisingly similar, with a much better agreement on the magnitude of aerosol impacts on cloud properties. Cloud adjustments remain a significant source of uncertainty, particularly for ice clouds. However, they are consistent with the uncertainty from observation-based methods, with the liquid water path adjustment usually enhancing the Twomey effect by less than 50 %. Depending on different sets of assumptions, this work suggests that model and observation-based estimates could be more equally weighted in future synthesis studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Thermodynamic analysis of the stability of planar interfaces between coexisting phases and its application to supercooled water

Two-phase simulations are commonly used to evaluate coexistence conditions, interfacial tensions, and other thermodynamic properties associated with first-order phase transitions. Calculation of these properties is often simplified when the interfaces between the two phases are flat or planar. In this paper, we derive a general thermodynamic criterion for selecting simulation cell dimensions to stabilize planar interfaces in phase-separated fluid-fluid systems with respect to homogeneous, single-phase states. The resulting expression is validated by analyzing the effects of simulation cell dimensions on the formation of planar liquid-vapor interfaces in the Lennard-Jones fluid and in the TIP4P/2005 model of water. We also perform large scale molecular dynamics simulations to study metastable liquid-liquid phase separation in the ST2 and TIP4P/2005 models of water under deeply supercooled conditions. Our simulations confirm the stability of a liquid-liquid interface in ST2, and they demonstrate that the corresponding interface for TIP4P/2005 can be stabilized by judiciously choosing the simulation cell aspect ratio in a manner consistent with the thermodynamic criterion. We posit that this sensitivity to the simulation cell aspect ratio may explain discrepancies between previous studies examining liquid-liquid separation in models of supercooled water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Ultraviolet emission of the Vela supernova remnant and the Cygnus Loop

Ultraviolet emission spectra of several filaments in the Vela supernova remnant and the Cygnus Loop are compared with theoretical models of radiative shocks. Shock velocities range from 85 to 130 km/s. Little or no grain depletion is apparent. Several of the observed filaments deviate from the steady-flow models in the sense that much of the photoionization-recombination zone is absent. It is shown that resonance line scattering within the emitting region can account for the weakness of the C II 1335 and C IV 1550 lines. Evaporation fronts, which are basic to the McKee and Ostriker picture of supernova remnant evolution, do not contribute significantly to the observed emission. Some discrepancies in UV emission predicted by different sets of shock models are traced to differing atomic rates used in the calculations. A fluorescent emission line of H2 may be observable in some supernova remnants.

Raymond, J. C.↗
NASA NTRS2003

Toroidal wave frequency at L=6-10: AMPTE/CCE Observations and Comparison with Theoretical Model

Magnetospheric standing Alfven waves are guided along the ambient magnetic field and their frequency depends on the mass density of the plasma distributed along the field lines. These properties allow us to use Alfen waves to map time-dependent phenomena between space and ground and to estimate the mass density. In this paper we present a statistical study of the spatial variation of the fundamental frequency f(sub T1) of toroidal-mode standing Alfen waves in the L range from 6 to 10, where L indicates the maximum geocentric distance on the field line. The data used for this analysis are energetic particle flux anisotropy (proxy of transverse electric field) and magnetic field measurements from the Active Magnetospheric Particle Tracer Explorers Charge Composition Explorer (CCE) spacecraft. Using CCE data covering 4 years we obtained approximately 5000 20-min intervals containing a clear signature of toroidal waves. The median f(sub T1), is 6-10 mHz at L=7 and decreases to 4-8 mHz at L = . The frequency tends to be lower at noon than at midnight. The observed frequencies are compared with numerically derived frequencies using an empirical mass density model and a magnetic field model. We found a good agreement for 12-24 MLT but a large discrepancy near MLT= 3. This may indicate that the flux tubes at this local time are more heavily loaded than specified in the Gallagher et al. model.

Takahashi, Kazue↗
NASA NTRS1985

Effects of soil tillage on the microwave emission of soils

In order to understand the interactions of soil properties and microwave emission better, a series of field experiments were conducted in 1984. Small plots were measured with a truck-mounted passive microwave radiometer operating at 1.4 GHz. These data were collected concurrent with ground observations of soil moisture and bulk density. Treatment effects studied included different soil moisture contents and bulk densities. Evaluations of the data have shown that commonly used models of the dielectric properties of wet soils do not explain the observations obtained in these experiments. This conclusion was based on the fact that the roughness parameters determined through optimization were significantly larger than those observed in similar investigations. These discrepancies are most likely due to the soil structure. Commonly used models assume a homogeneous three phase mixture of soil solids, air and water. Under tilled conditions the soil is actually a two phase mixture of aggregates and voids. Appropriate dielectric models for this tilled condition were evaluated and found to explain the observations. These results indicate that previous conclusions concerning the effects of surface roughness in tilled fields may be incorrect, and they may explain some of the inconsistencies encountered in roughness modeling.

Jackson, T. J.↗
DOE OSTI2025

Benchmarking a Defeatured Geant4 NIF Model with HTOAD and HNED Neutron Data

This study benchmarks a simplified Geant4 NIF Target Chamber (TC) model against foil measurements in two instruments: the HTOAD and HNED Snout. The number of product atoms per source neutron per gram ( N 0 /n/g ) is compared across multiple locations and material configurations. The model reproduces high energy threshold reactions dominated by 14 MeV neutrons to within a few percent, depending on location. Discrepancies occur where scattering and moderation are significant for low energy threshold reactions. The defeatured Geant4 model achieves comparable statistical precision with ∼400 times less CPU time than a full fidelity TC model in MCNP. The results can be produced locally on a laptop without the need for high performance computing. Low energy biasing can be improved by preserving room return pathways and implementing IRDFF cross section evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Call for Participation: Collaborative Benchmarking of Functional-Structural Root Architecture Models. The Case of Root Water Uptake

Three-dimensional models of root growth, architecture and function are becoming important tools that aid the design of agricultural management schemes and the selection of beneficial root traits. However, while benchmarking is common in many disciplines that use numerical models, such as natural and engineering sciences, functional-structural root architecture models have never been systematically compared. The following reasons might induce disagreement between the simulation results of different models: different representation of root growth, sink term of root water and solute uptake and representation of the rhizosphere. Presently, the extent of discrepancies is unknown, and a framework for quantitatively comparing functional-structural root architecture models is required. We propose, in a first step, to define benchmarking scenarios that test individual components of complex models: root architecture, water flow in soil and water flow in roots. While the latter two will focus mainly on comparing numerical aspects, the root architectural models have to be compared at a conceptual level as they generally differ in process representation. Therefore, defining common inputs that allow recreating reference root systems in all models will be a key challenge. In a second step, benchmarking scenarios for the coupled problems are defined. We expect that the results of step 1 will enable us to better interpret differences found in step 2. This benchmarking will result in a better understanding of the different models and contribute toward improving them. Improved models will allow us to simulate various scenarios with greater confidence and avoid bugs, numerical errors or conceptual misunderstandings. This work will set a standard for future model development.

Plant Sciences↗
NASA NTRS2001

Determination of Top-of-Atmosphere Longwave Radiative Fluxes: A Comparison Between Two Approaches Using ScaRaB Data

Two conceptually different approaches (broadband-based ERBE and narrowband-based ISCCP approaches), used to derive the TOA longwave radiative fluxes, are compared using the ScaRaB simultaneous narrowband and broadband measurements. This study directly shows that the ERBE MLE-derived cloud covers implicitly contain some information on the cloud optical properties. A spurious view-zenith-angle dependence of the MLE scene identification scheme is confirmed by this study. Except for very thin cirrus clouds, differences between the ERBE and ISCCP approaches are in general < 10 W/sq m for the TOA LW radiative fluxes. For clear pixels, the model calculated (ISCCP approach) TOA LW radiances are systematically smaller than the observations. Though the bias is found to be correlated on the column precipitable water amount, the exact source of this discrepancy remains undetermined and merits further study. Compared with the radiative transfer model used in this study, the ERBE LW ADMs are too weakly limb-darkened for optically thin clouds, but too strongly limb-darkened for optically thick clouds, indicating that more accurate instantaneous TOA LW flux estimations from the ERBE approach would require additional cloud classes based on cloud height and optical thickness.

Chen, Ting↗
NASA NTRS1987

On the interpretation of the Sandage period-shift effect among globular-cluster RR Lyrae variables

A new grid of canonical zero-age horizontal branch (ZAHB) models is constructed to study the Sandage period-shift effect, and a number of possible explanations for the discrepancy between the observed and predicted period shifts are explored. Various factors which determine the period shift in ZAHB models are discussed in detail, showing that possible uncertainties in the input physics used for computing the preceding red-giant evolution are unlikely to resolve the discrepancy. Some test calculations for the ZAHB phase are presented which suggest than an increase in the metal opacity at temperatures around a million K might substantially increase the predicted period shift. It is found that a nonsolar CNO to Fe ratio is not a plausible explanation of the period-shift discrepancy, and the effects of evolution off the ZAHB in the period-effective temperature diagram are discussed. It is shown that evolution can increase the predicted period shift but not to the extent required by the Sandage effect. The possibility that rotation might explain the discrepancy is considered and rejected.

Sweigart, Allen V.↗
DOE OSTI2023

Comparing Measured Driver Behavior Distributions to Results from Car-Following Models using SUMO and Real-World Vehicle Trajectories from Radar: SUMO Default vs. Radar-Measured CF model Parameters

In this study, the physical principles governing car-following (CF) behavior and their impact on traffic flow at signalized intersections are investigated. High temporal-resolution radar data is used to provide valuable insights into actual CF behavior, including acceleration, deceleration, and time headway distribution. Demand-calibrated SUMO simulations are run using empirical CF parameter distributions, and three CF models are evaluated: IDM, EIDM, and Krauss. By emulating radar data in SUMO and processing simulated vehicle traces, discrepancies between empirical and simulated parameter distributions are identified. Further analysis includes comparisons with default SUMO CF model parameters. The findings reveal that measured accelerations differ from CF model parameter accelerations and using the empirical value ($\mu = 0.89m/s^2$) leads to unrealistic simulations that fail volume-based calibration. Default parameters for all three models reasonably approximate the mean and median of measured parameters, but fail to capture the true distribution shape, partly due to homogeneity when using default parameters. The results show that it is more effective to simulate with the default parameters provided by SUMO rather than using measurements of real-world distributions without additional calibration. Future work will investigate closing the loop between the measured real-world and SUMO distributions using traditional calibration tactics, as well as assess the impact of calibrated vs. default CF parameters on simulation outputs like fuel consumption.

Schrader, Max (ORCID:0000000336737672)↗
NASA NTRS1983

Dust particles in the anomalous tail of Comet Kohoutek (1973 XII)

Intermediate results from theoretical models of the anomalous spikelike dust distribution observed in the tail of Comet Kohoutek using the Skylab white-light coronograph on December 28.996, 1973, are reported. The time evolution is studied using the model of Finson and Probstein (1968) and assuming power-law variation of the dust-production rate and the particle-size distribution. Brightness profiles are calculated and presented graphically for comparison with those derived from the observation photograph: the model predicts a longer duration of the spike phenomenon than was observed. It is shown qualitatively that this discrepancy can be removed by introducing a finite dust-emission velocity (e.g., 1 m/s) into the model.

Keller, H. U.↗