Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MODELING”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Modeling surface motion effects in N{sub 2} dissociation on W(110): Ab initio molecular dynamics calculations and generalized Langevin oscillator model

Accurately modeling surface temperature and surface motion effects is necessary to study molecule-surface reactions in which the energy dissipation to surface phonons can largely affect the observables of interest. We present here a critical comparison of two methods that allow to model such effects, namely, the ab initio molecular dynamics (AIMD) method and the generalized Langevin oscillator (GLO) model, using the dissociation of N{sub 2} on W(110) as a benchmark. AIMD is highly accurate as the surface atoms are explicitly part of the dynamics, but this advantage comes with a large computational cost. The GLO model is much more computationally convenient, but accounts for lattice motion effects in a very approximate way. Results show that, despite its simplicity, the GLO model is able to capture the physics of the system to a large extent, returning dissociation probabilities which are in better agreement with AIMD than static-surface results. Furthermore, the GLO model and the AIMD method predict very similar energy transfer to the lattice degrees of freedom in the non-reactive events, and similar dissociation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Measuring and modeling diffuse scattering in protein X-ray crystallography

X-ray diffraction has the potential to provide rich information about the structural dynamics of macromolecules. To realize this potential, both Bragg scattering, which is currently used to derive macromolecular structures, and diffuse scattering, which reports on correlations in charge density variations, must be measured. Until now, measurement of diffuse scattering from protein crystals has been scarce because of the extra effort of collecting diffuse data. Here, we present 3D measurements of diffuse intensity collected from crystals of the enzymes cyclophilin A and trypsin. The measurements were obtained from the same X-ray diffraction images as the Bragg data, using best practices for standard data collection. To model the underlying dynamics in a practical way that could be used during structure refinement, we tested translation–libration–screw (TLS), liquid-like motions (LLM), and coarse-grained normal-modes (NM) models of protein motions. The LLM model provides a global picture of motions and was refined against the diffuse data, whereas the TLS and NM models provide more detailed and distinct descriptions of atom displacements, and only used information from the Bragg data. Whereas different TLS groupings yielded similar Bragg intensities, they yielded different diffuse intensities, none of which agreed well with the data. In contrast, both the LLM and NM models agreed substantially with the diffuse data. In conclusion, these results demonstrate a realistic path to increase the number of diffuse datasets available to the wider biosciences community and indicate that dynamics-inspired NM structural models can simultaneously agree with both Bragg and diffuse scattering.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

PHOTOMETRIC REDSHIFTS AND QUASAR PROBABILITIES FROM A SINGLE, DATA-DRIVEN GENERATIVE MODEL

We describe a technique for simultaneously classifying and estimating the redshift of quasars. It can separate quasars from stars in arbitrary redshift ranges, estimate full posterior distribution functions for the redshift, and naturally incorporate flux uncertainties, missing data, and multi-wavelength photometry. We build models of quasars in flux-redshift space by applying the extreme deconvolution technique to estimate the underlying density. By integrating this density over redshift, one can obtain quasar flux densities in different redshift ranges. This approach allows for efficient, consistent, and fast classification and photometric redshift estimation. This is achieved by combining the speed obtained by choosing simple analytical forms as the basis of our density model with the flexibility of non-parametric models through the use of many simple components with many parameters. We show that this technique is competitive with the best photometric quasar classification techniques-which are limited to fixed, broad redshift ranges and high signal-to-noise ratio data-and with the best photometric redshift techniques when applied to broadband optical data. We demonstrate that the inclusion of UV and NIR data significantly improves photometric quasar-star separation and essentially resolves all of the redshift degeneracies for quasars inherent to the ugriz filter system, even when included data have a low signal-to-noise ratio. For quasars spectroscopically confirmed by the SDSS 84% and 97% of the objects with Galaxy Evolution Explorer UV and UKIDSS NIR data have photometric redshifts within 0.1 and 0.3, respectively, of the spectroscopic redshift; this amounts to about a factor of three improvement over ugriz-only photometric redshifts. Our code to calculate quasar probabilities and redshift probability distributions is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

LOW MACH NUMBER MODELING OF CONVECTION IN HELIUM SHELLS ON SUB-CHANDRASEKHAR WHITE DWARFS. I. METHODOLOGY

We assess the robustness of a low Mach number hydrodynamics algorithm for modeling helium shell convection on the surface of a white dwarf in the context of the sub-Chandrasekhar model for Type Ia supernovae. We use the low Mach number stellar hydrodynamics code, MAESTRO, to perform three-dimensional, spatially adaptive simulations of convection leading up to the point of the ignition of a burning front. We show that the low Mach number hydrodynamics model provides a robust description of the system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

CARBON DEFLAGRATION IN TYPE Ia SUPERNOVA. I. CENTRALLY IGNITED MODELS

A leading model for Type Ia supernovae (SNe Ia) begins with a white dwarf near the Chandrasekhar mass that ignites a degenerate thermonuclear runaway close to its center and explodes. In a series of papers, we shall explore the consequences of ignition at several locations within such dwarfs. Here we assume central ignition, which has been explored before, but is worth revisiting, if only to validate those previous studies and to further elucidate the relevant physics for future work. A perturbed sphere of hot iron ash with a radius of {approx}100 km is initialized at the middle of the star. The subsequent explosion is followed in several simulations using a thickened flame model in which the flame speed is either fixed-within the range expected from turbulent combustion-or based on the local turbulent intensity. Global results, including the explosion energy and bulk nucleosynthesis (e.g., {sup 56}Ni of 0.48-0.56 M{sub Sun }) turn out to be insensitive to this speed. In all completed runs, the energy released by the nuclear burning is adequate to unbind the star, but not enough to give the energy and brightness of typical SNe Ia. As found previously, the chemical stratification observed in typical events is not reproduced. These models produce a large amount of unburned carbon and oxygen in central low velocity regions, which is inconsistent with spectroscopic observations, and the intermediate mass elements and iron group elements are strongly mixed during the explosion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Time-resolved analysis of Fermi gamma-ray bursts with fast- and slow-cooled synchrotron photon models

Time-resolved spectroscopy is performed on eight bright, long gamma-ray bursts (GRBs) dominated by single emission pulses that were observed with the Fermi Gamma-Ray Space Telescope. Fitting the prompt radiation of GRBs by empirical spectral forms such as the Band function leads to ambiguous conclusions about the physical model for the prompt radiation. Moreover, the Band function is often inadequate to fit the data. The GRB spectrum is therefore modeled with two emission components consisting of optically thin non-thermal synchrotron radiation from relativistic electrons and, when significant, thermal emission from a jet photosphere, which is represented by a blackbody spectrum. To produce an acceptable fit, the addition of a blackbody component is required in five out of the eight cases. We also find that the low-energy spectral index α is consistent with a synchrotron component with α = –0.81 ± 0.1. This value lies between the limiting values of α = –2/3 and α = –3/2 for electrons in the slow- and fast-cooling regimes, respectively, suggesting ongoing acceleration at the emission site. The blackbody component can be more significant when using a physical synchrotron model instead of the Band function, illustrating that the Band function does not serve as a good proxy for a non-thermal synchrotron emission component. The temperature and characteristic emission-region size of the blackbody component are found to, respectively, decrease and increase as power laws with time during the prompt phase. In addition, we find that the blackbody and non-thermal components have separate temporal behaviors as far as their respective flux and spectral evolutions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Toward a better understanding of the GRB phenomenon: A new model for GRB PROMPT emission and its effects on the new L{sub i}{sup NT}– E{sub peak,i}{sup rest,NT} relation

Gamma-ray burst (GRB) prompt emission spectra in the keV–MeV energy range are usually considered to be adequately fitted with the empirical Band function. Recent observations with the Fermi Gamma-ray Space Telescope (Fermi) revealed deviations from the Band function, sometimes in the form of an additional blackbody (BB) component, while on other occasions in the form of an additional power law (PL) component extending to high energies. In this article we investigate the possibility that the three components may be present simultaneously in the prompt emission spectra of two very bright GRBs (080916C and 090926A) observed with Fermi, and how the three components may affect the overall shape of the spectra. While the two GRBs are very different when fitted to a single Band function, they look like “twins” in the three-component scenario. Through fine-time spectroscopy down to the 100 ms timescale, we follow the evolution of the various components. We succeed in reducing the number of free parameters in the three-component model, which results in a new semi-empirical model—but with physical motivations—to be competitive with the Band function in terms of number of degrees of freedom. From this analysis using multiple components, the Band function is globally the most intense component, although the additional PL can overpower the others in sharp time structures. The Band function and the BB component are the most intense at early times and globally fade across the burst duration. The additional PL is the most intense component at late time and may be correlated with the extended high-energy emission observed thousands of seconds after the burst with Fermi/Large Area Telescope. Unexpectedly, this analysis also shows that the additional PL may be present from the very beginning of the burst, where it may even overpower the other components at low energy. We investigate the effect of the three components on the new time-resolved luminosity–hardness relation in both the observer and rest frames and show that a strong correlation exists between the flux of the non-thermal Band function and its E{sub peak} only when the three components are fitted simultaneously to the data (i.e., F{sub i}{sup NT}–E{sub peak,i}{sup NT} relation). In addition, this result points toward a universal relation between those two quantities when transposed to the central engine rest frame for all GRBs (i.e., L{sub i}{sup NT}–E{sub peak,i}{sup rest,NT} relation). We discuss a possible theoretical interpretation of the three spectral components within this new empirical model. We suggest that (i) the BB component may be interpreted as the photosphere emission of a magnetized relativistic outflow, (ii) the Band component has synchrotron radiation in an optically thin region above the photosphere, either from internal shocks or magnetic field dissipation, and (iii) the extra PL component extending to high energies likely has an inverse Compton origin of some sort, even though its extension to a much lower energy remains a mystery.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

QUASI-STATIC MODEL OF MAGNETICALLY COLLIMATED JETS AND RADIO LOBES. II. JET STRUCTURE AND STABILITY

This is the second in a series of companion papers showing that when an efficient dynamo can be maintained by accretion disks around supermassive black holes in active galactic nuclei, it can lead to the formation of a powerful, magnetically driven, and mediated helix that could explain both the observed radio jet/lobe structures and ultimately the enormous power inferred from the observed ultrahigh-energy cosmic rays. In the first paper, we showed self-consistently that minimizing viscous dissipation in the disk naturally leads to jets of maximum power with boundary conditions known to yield jets as a low-density, magnetically collimated tower, consistent with observational constraints of wire-like currents at distances far from the black hole. In this paper we show that these magnetic towers remain collimated as they grow in length at nonrelativistic velocities. Differences with relativistic jet models are explained by three-dimensional magnetic structures derived from a detailed examination of stability properties of the tower model, including a broad diffuse pinch with current profiles predicted by a detailed jet solution outside the collimated central column treated as an electric circuit. We justify our model in part by the derived jet dimensions in reasonable agreement with observations. Using these jet properties, we also discuss the implications for relativistic particle acceleration in nonrelativistically moving jets. The appendices justify the low jet densities yielding our results and speculate how to reconcile our nonrelativistic treatment with general relativistic MHD simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Background model systematics for the Fermi GeV excess

The possible gamma-ray excess in the inner Galaxy and the Galactic center (GC) suggested by Fermi-LAT observations has triggered a large number of studies. It has been interpreted as a variety of different phenomena such as a signal from WIMP dark matter annihilation, gamma-ray emission from a population of millisecond pulsars, or emission from cosmic rays injected in a sequence of burst-like events or continuously at the GC. We present the first comprehensive study of model systematics coming from the Galactic diffuse emission in the inner part of our Galaxy and their impact on the inferred properties of the excess emission at Galactic latitudes 2° < |b| < 20° and 300 MeV to 500 GeV. We study both theoretical and empirical model systematics, which we deduce from a large range of Galactic diffuse emission models and a principal component analysis of residuals in numerous test regions along the Galactic plane. We show that the hypothesis of an extended spherical excess emission with a uniform energy spectrum is compatible with the Fermi-LAT data in our region of interest at 95% CL. Assuming that this excess is the extended counterpart of the one seen in the inner few degrees of the Galaxy, we derive a lower limit of 10.0° (95% CL) on its extension away from the GC. We show that, in light of the large correlated uncertainties that affect the subtraction of the Galactic diffuse emission in the relevant regions, the energy spectrum of the excess is equally compatible with both a simple broken power-law of break energy E{sub break} = 2.1 ± 0.2 GeV, and with spectra predicted by the self-annihilation of dark matter, implying in the case of b-bar b final states a dark matter mass of m{sub χ}=49{sup +6.4}{sub −5.4} GeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Complementary test of the dark matter self-interaction in dark U(1) model by direct and indirect dark matter detection

The halo dark matter (DM) can be captured by the Sun if its final velocity after the collision with a nucleus in the Sun is less than the escape velocity. We consider a selfinteracting dark matter (SIDM) model where U(1) gauge symmetry is introduced to account for the DM self-interaction. Such a model naturally leads to isospin violating DM-nucleon interaction, although isospin symmetric interaction is still allowed as a special case. We present the IceCube-PINGU 2σ sensitivity to the parameter range of the above model with 5 years of search for neutrino signature from DM annihilation in the Sun. This indirect detection complements the direct detection by probing those SIDM parameter ranges which are either the region for very small m{sub χ} or the region opened up due to isospin violations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Modeling CMB lensing cross correlations with CLEFT

A new generation of surveys will soon map large fractions of sky to ever greater depths and their science goals can be enhanced by exploiting cross correlations between them. In this paper we study cross correlations between the lensing of the CMB and biased tracers of large-scale structure at high z. We motivate the need for more sophisticated bias models for modeling increasingly biased tracers at these redshifts and propose the use of perturbation theories, specifically Convolution Lagrangian Effective Field Theory (CLEFT). Since such signals reside at large scales and redshifts, they can be well described by perturbative approaches. We compare our model with the current approach of using scale independent bias coupled with fitting functions for non-linear matter power spectra, showing that the latter will not be sufficient for upcoming surveys. Furthermore, we illustrate our ideas by estimating σ 8 from the auto- and cross-spectra of mock surveys, finding that CLEFT returns accurate and unbiased results at high z. We report on uncertainties due to the redshift distribution of the tracers, and several avenues for future development.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Quenching measurements and modeling of a boron-loaded organic liquid scintillator

Organic liquid scintillators are used in a wide variety of applications in experimental nuclear and particle physics. Boron-loaded scintillators are particularly useful for detecting neutron captures, due to the high thermal neutron capture cross section of 10 B. These scintillators are commonly used in neutron detectors, including the DarkSide-50 neutron veto, where the neutron may produce a signal when it scatters o protons in the scintillator or when it captures on 10 B. Reconstructing the energy of these recoils is complicated by scintillation quenching. Understanding how nuclear recoils are quenched in these scintillators is an important and dicult problem. In this article, we present a set of measurements of neutron-induced proton recoils in a boron-loaded organic liquid scintillator at recoil energies ranging from 57-467 keV, and we compare these measurements to predictions from different quenching models. We and that a modified Birks' model whose denominator is quadratic in dE=dx best describes the measurements, with χ 2 /NDF = 1:6. This result will help model nuclear recoil scintillation in similar detectors and can be used to improve their neutron tagging efficiency.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Improving streamflow predictions across CONUS by integrating advanced machine learning models and diverse data

Accurate streamflow prediction is crucial to understand climate impacts on water resources and develop effective adaption strategies. A global long short-term memory (LSTM) model, using data from multiple basins, can enhance streamflow prediction, yet acquiring detailed basin attributes remains a challenge. To overcome this, we introduce the Geo-vision transformer (ViT)-LSTM model, a novel approach that enriches LSTM predictions by integrating basin attributes derived from remote sensing with a ViT architecture. Applied to 531 basins across the Contiguous United States, our method demonstrated superior prediction accuracy in both temporal and spatiotemporal extrapolation scenarios. Geo-ViT-LSTM marks a significant advancement in land surface modeling, providing a more comprehensive and effective tool for better understanding the environment responses to climate change.

Tayal, Kshitij↗
DOE OSTI2015

Bioenergy Supply and Environmental Impacts on Cropland: Insights from Multi‐market Forecasts in a Great Lakes Subregional Bioeconomic Model

Abstract Using subregional models of crop production choices in central Wisconsin and southwest Michigan, we predict biomass production, land use, and environmental impacts with details that are unavailable from national scale models. When biomass prices are raised exogenously, we find that the subregional models overestimate the supply, the land use, and the beneficial environmental aspects of perennial biomass crops. Multi‐market price feedbacks tied to realistic policy parameters predict high threshold absolute prices for biomass to enter production, resulting in intensified production of biomass from annual grain crops with damaging environmental impacts. Multi‐market feedbacks also predict regional specialization in energy biomass production in areas with lower yields of food crops. Policies promoting biofuels will not necessarily generate environmental benefits in the absence of environmental regulations.

Egbendewe‐Mondzozo, Aklesso↗
DOE OSTI2018

Accelerating flux balance calculations in genome-scale metabolic models by localizing the application of loopless constraints

Abstract Background Genome-scale metabolic network models and constraint-based modeling techniques have become important tools for analyzing cellular metabolism. Thermodynamically infeasible cycles (TICs) causing unbounded metabolic flux ranges are often encountered. TICs satisfy the mass balance and directionality constraints but violate the second law of thermodynamics. Current practices involve implementing additional constraints to ensure not only optimal but also loopless flux distributions. However, the mixed integer linear programming problems required to solve become computationally intractable for genome-scale metabolic models. Results We aimed to identify the fewest needed constraints sufficient for optimality under the loopless requirement. We found that loopless constraints are required only for the reactions that share elementary flux modes representing TICs with reactions that are part of the objective function. We put forth the concept of localized loopless constraints (LLCs) to enforce this minimal required set of loopless constraints. By combining with a novel procedure for minimal null-space calculation, the computational time for loopless flux variability analysis (ll-FVA) is reduced by a factor of 10–150 compared to the original loopless constraints and by 4–20 times compared to the current fastest method Fast-SNP with the percent improvement increasing with model size. Importantly, LLCs offer a scalable strategy for loopless flux calculations for multi-compartment/multi-organism models of large sizes, for example, shortening the CPU time for ll-FVA from 35 h to less than 2 h for a model with more than104 reactions. Availability and implementation Matlab functions are available in the Supplementary Material or at https://github.com/maranasgroup/lll-FVA Supplementary information Supplementary data are available at Bioinformatics online.

Chan, Siu H. J. (ORCID:000000027707656X)↗
DOE OSTI2023

Systematics in asteroseismic modelling: application of a correlated noise model for oscillation frequencies

ABSTRACT The detailed modelling of stellar oscillations is a powerful approach to characterizing stars. However, poor treatment of systematics in theoretical models leads to misinterpretations of stars. Here, we propose a more principled statistical treatment for the systematics to be applied to fitting individual mode frequencies with a typical stellar model grid. We introduce a correlated noise model based on a Gaussian process (GP) kernel to describe the systematics given that mode frequency systematics are expected to be highly correlated. We show that tuning the GP kernel can reproduce general features of frequency variations for changing model input physics and fundamental parameters. Fits with the correlated noise model better recover stellar parameters than traditional methods that either ignore the systematics or treat them as uncorrelated noise.

Li, Tanda (ORCID:0000000163962563)↗
DOE OSTI2021

Universality class for loopless invasion percolation models and a percolation avalanche burst model for hydraulic fracturing

Invasion percolation is a model that was originally proposed to describe growing networks of fractures. In this work, we describe a loopless algorithm on random lattices, coupled with an avalanche-based model for bursts. The model reproduces the characteristic b-value seismicity and spatial distribution of bursts consistent with earthquakes resulting from hydraulic fracturing (“fracking”). We test models for both site invasion percolation and bond invasion percolation. These have differences on the scale of site and bond lengths l. But since the networks are characterized by their large-scale behavior, l << L, we find small differences between scaling exponents. Though data may not differentiate between models, our results suggest that both models belong to different universality classes.

58 GEOSCIENCES↗