Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Renormalization group improvement of the effective potential: an EFT approach

We apply effective field theory (EFT) methods to compute the renormalization group improved effective potential for theories with a large mass hierarchy. Our method allows one to compute the effective potential in a systematic expansion in powers of the mass ratio, as well as to sum large logarithms of mass ratios using renormalization group evolution. The effective potential is the sum of one-particle irreducible diagrams (1PI) but information about which diagrams are 1PI is lost after matching to the EFT, since heavy lines get shrunk to a point. We therefore introduce a tadpole condition in place of the 1PI condition, and use the renormalization group improved value of the tadpole in computing the effective potential. We explain why the effective potential computed using an EFT is not the same as the effective potential of the EFT. We illustrate our method using the O(N) model, a theory of two scalars in the unbroken and broken phases, and the Higgs-Yukawa model. Our leading-log result, obtained by integrating the one-loop β-functions, correctly reproduces the log-squared term in explicit two-loop calculations. Our method does not have a Goldstone boson infrared divergence problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A non-intrusive optical (NIO) approach to characterize heliostats in utility-scale power tower plants: Sensitivity study

Power towers, a type of concentrating power tower technology, use a large number of heliostats to concentrate sunlight and produce renewable energy. Optical errors of heliostats can cause drastic losses in the efficiency of power-tower plants. Accurately measuring optical errors is crucial to assessing and improving plant performance. This analysis discusses an innovate non-intrusive optical (NIO) approach to measure slope, canting, and tracking errors by detecting distortions in the reflected tower structure in heliostat images. Further in this work a sensitivity study is carefully conducted to determine the uncertainty requirements to allow the method to calculate slope errors with an accuracy of 0.25 mrad. Measurement uncertainty sources include camera resolution and position uncertainty, tower position uncertainty, and number of collected images. Each uncertainty source is investigated to determine its impact on the accuracy of the slope-error calculation. A combination of theoretical results and experimental results obtained from data collected on a heliostat at Sandia National Laboratories is used to determine and validate uncertainty requirements. The analysis shows that a measurement uncertainty of 0.25 mrad can be realized by realistically controlling uncertainty sources when implementing the NIO method. It demonstrates the superior performance of NIO in performing in-situ optical characterization.

14 SOLAR ENERGY↗
DOE OSTI2020

Soft x-ray spectrum unfold of K-edge filtered x-ray diode arrays using cubic splines

Cubic spline interpolation is able to recover temporally and spectrally resolved soft x-ray flux from an array of K-edge filtered x-ray diodes without the need for a priori assumptions about the spectrum or the geometry of the emitting volume. The mathematics of the cubic spline interpolation is discussed in detail. The analytic nature of the cubic spline solution allows for analytical error propagation, and the method of calculating the error for radiation temperature, spectral power, and confidence intervals of the unfolded spectrally resolved flux is explained. An unfold of a blackbody model demonstrates the accuracy of the cubic spline unfold. Tests of cubic spline performance using spectrally convolved detailed atomic model simulation results have been performed to measure the method's ability to conserve spectral power to within a factor of 2 or better in line-dominated regimes. The unfold is also demonstrated to work when information from the x-ray diode array is limited due to high signal-to-noise ratios or lack of signal due to over-attenuation or over-filtration of the x-ray diode signal. Here, the robustness of the unfold with respect to background subtraction and raw signal processing, signal alignment between diode traces, limited signal information, and initial conditions is discussed. Results from an example analysis of a halfraum drive are presented to demonstrate the capabilities of the unfold in comparison with previously established methods.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

On decoupling the integrals of cosmological perturbation theory

ABSTRACT Perturbation theory (PT) is often used to model statistical observables capturing the translation and rotation-invariant information in cosmological density fields. PT produces higher order corrections by integration over linear statistics of the density fields weighted by kernels resulting from recursive solution of the fluid equations. These integrals quickly become high dimensional and naively require increasing computational resources the higher the order of the corrections. Here, we show how to decouple the integrands that often produce this issue, enabling PT corrections to be computed as a sum of products of independent 1D integrals. Our approach is related to a commonly used method for calculating multiloop Feynman integrals in quantum field theory, the Gegenbauer Polynomial x-Space Technique. We explicitly reduce the three terms entering the 2-loop power spectrum, formally requiring 9D integrations, to sums over successive 1D radial integrals. These 1D integrals can further be performed as convolutions, rendering the scaling of this method Nglog Ng with Ng the number of grid points used for each fast Fourier transform. This method should be highly enabling for upcoming large-scale structure redshift surveys where model predictions at an enormous number of cosmological parameter combinations will be required by Monte Carlo Markov Chain searches for the best-fitting values.

Slepian, Zachary↗
DOE OSTI2020

Diagnosing the sensitivity of grounding-line flux to changes in sub-ice-shelf melting

Abstract. Using a numerical ice flow model, we study changes in ice shelf buttressing and grounding-line flux due to localized ice thickness perturbations, a proxy for localized changes in sub-ice-shelf melting. From our experiments, applied to idealized (MISMIP+) and realistic (Larsen C) ice shelf domains, we identify a correlation between a locally derived buttressing number on the ice shelf, based on the first principal stress, and changes in the integrated grounding-line flux. The origin of this correlation, however, remains elusive from the perspective of a theoretical or physically based understanding. This and the fact that the correlation is generally much poorer when applied to realistic ice shelf domains motivate us to seek an alternative approach for predicting changes in grounding-line flux. We therefore propose an adjoint-based method for calculating the sensitivity of the integrated grounding-line flux to local changes in ice shelf geometry. We show that the adjoint-based sensitivity is identical to that deduced from pointwise, diagnostic model perturbation experiments. Based on its much wider applicability and the significant computational savings, we propose that the adjoint-based method is ideally suited for assessing grounding-line flux sensitivity to changes in sub-ice-shelf melting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗
DOE OSTI2022

Development of a phonon-based sampling method for thermal neutron scattering data

Simulations of reactor systems require access to accurate nuclear data. For many systems, thermal neutron scattering data can have large effects on the eigenvalue and neutron flux distributions. Inelastic thermal neutron scattering can excite or de-excite vibrational, rotational, and translational modes in a material, so thermal scattering evaluations are often obtained by summing over the number of phonons created/destroyed by a scattering event. In recent years, the thermal scattering cross sections and angular distributions have greatly improved in accuracy, but the format in which this data is delivered to simulation codes has remained virtually unchanged. Thermal scattering data is typically either compiled into large tables and sorted by incoming neutron energy, outgoing neutron energy, scattering angle, and material temperature, or represented as cumulative distribution functions of momentum exchange or energy exchange. Either method can be quite memory intensive when fine bins are used. In an effort to decrease the amount of space that processed thermal scattering data requires, an alternate format is proposed. The phonon-based sampling method introduced here can sample the number of phonons excited for each collision, the change in neutron energy, and the scattering angle while avoiding pre-computed angular bins and limiting the amount of data that is dependent on incoming energy. Through this method, the generation and storage of large interpolation tables is avoided, which could have benefits in both memory storage and accuracy. While the initial implementation of this method is slower than current alternatives, it is significantly more resistant to grid coarseness errors and has good potential for improvement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Specifications of FIPD Fission Gas Release Data

All fission gas release data stored in the Fuels Irradiation & Physics Database (FIPD) was originally measured using the Gas Assay, Sample and Recharge (GASR) System in the Hot Fuel Examination Facility (HFEF). It is therefore called GASR data in FIPD. During the measurement of a sample, such as an irradiated EBR-II fuel element/capsule, a pinhole-sized region near the top of the element plenum was melted by a laser. Plenum gas then expanded into a calibrated volume (note: in this document, “sample” and “capsule/element” are used interchangeably consistent with GASR documents in FIPD). The pressure rise in the volume was recorded. Helium backfilling and expansion was then performed to determine the sample (e.g., fuel element plenum) volume using Boyle’s Law and assuming ideal gas behavior at constant temperature. With the plenum volume and the recorded pressure rise, the sample (e.g., fuel element plenum) pressure was derived with assumption of ideal gas law behavior. The plenum volume and pressure as well as the cladding temperature during the measurement were collected (GASR data in FIPD). Other records associated with the fission gas release data include: raw GASR data records including volumes and post-puncture pressures of seal head/sealing head and manifold, calibration data, backfilling gas pressure data, and the data analysis records. A sample(s) of the fission gas released from the plenum was collected by the GASR system into sample bottles. The chemical and isotopic composition of the gas sample was analyzed separately from GASR data, and will be discussed in a separate specification. The plenum volume, pressure, and cladding temperature during the measurement are typically utilized to determine the number of moles of gas in the plenum. This quantity is often compared to the number of moles of gas generated by fission events. However, calculating these values and their associated uncertainty is beyond the scope of this document, as it necessitates additional assumptions. The most important document to understand the FIPD fission gas data is the GASR operational manual (title: Gas Assay, Sample and Recharge System (GASR) operation and maintenance manual, HFEF/N OMM 4381, DOC. NO. W0018-0032-ES-00). This manual provides: (1) description of the GASR and the functions of each component (laser drilling, welding, seal head/sealing head, manifold, vacuum system, sample system, purge and gas tag system, etc.); (2) step-by-step guidance on calibrations, operations, and measurements; and (3) maintenance procedures and other details relating to the structure and operation of the GASR. Note that the original GASR operated until 2020. A new GASR with the same design and measurement methodology was installed in 2021. The specifications of the GASR presented on the HFEF website at this time are consistent with the ones given in the operational manual. The methods to calculate the plenum volume and pressures were not included in the operational manual, but were recorded in the legacy data analysis files. Details of the methods are given in Chapter 3.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Open momentum space method for the Hofstadter butterfly and the quantized Lorentz susceptibility

Here we develop a generic k ∙ p open momentum space method for calculating the Hofstadter butterfly of both continuum (moiré) models and tight-binding models, where the quasimomentum is directly substituted by the Landau level (LL) operators. By taking a LL cutoff (and a reciprocal lattice cutoff for continuum models), one obtains the Hofstadter butterfly with in-gap spectral flows. For continuum models such as the moiré model for twisted bilayer graphene, our method gives a sparse Hamiltonian, making it much more efficient than existing methods. The spectral flows in the Hofstadter gaps can be understood as edge states on a momentum space boundary, from which one can determine the two integers (t ν , s ν ) of a gap ν satisfying the Diophantine equation. The spectral flows can also be removed to obtain a clear Hofstadter butterfly. While t ν is known as the Chern number, our theory identifies s ν as a dual Chern number for the momentum space, which corresponds to a quantized Lorentz susceptibility γ xy = eBs ν .

2-dimensional systems↗
DOE OSTI2020

Model of Operation-and-Maintenance Costs for Photovoltaic Systems

This article presents a method for calculating costs associated with operation and maintenance (O&M) of photovoltaic (PV) systems. It compiles details regarding the cost and frequency of multiple O&M services to estimate annual O&M costs ($\$$/year) for each year of an analysis period, the net present value ($\$$) of life cycle costs accumulated over the analysis period, and the reserve account amount ($\$$). Here we show that this method is an improvement over the previous averaged or levelized per-unit ($\$$/kW/year) valuations for estimating PV O&M costs, because it allows a detailed selection of services to perform based on system size, market served (e.g., residential, commercial, or utility), type and configuration of system components (e.g., micro-, string or central inverter) and site and environmental conditions (e.g., snow, pollen, bird populations). This model also distinguishes costs that vary from year to year and increase at different rates over time because of heuristic failure distributions (e.g. Weibull or Lognormal distribution) based on actuarial data for many of the services. This cost model was created by the PV O&M Working Group of researchers and industry, sponsored by DOE Solar Energy Technologies Office, and has been published in an on-line version hosted by SunSpec Alliance at apsuite.sunspec.org. A spreadsheet version is included with this paper as supplementary material.

14 SOLAR ENERGY↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
DOE OSTI2022

Neutron Leakage Spectra Sensitivities for ICSBEP Benchmarks

Neutron leakage spectra have been measured, simulated, and investigated by many groups. These spectra have many uses, including determining shielding requirements and calculating dose for radiation protection purposes, validating nuclear data, and determining material composition. It has even been proposed to use measurements of neutron leakage spectra to determine the soil composition of Mars. As part of the Los Alamos National Laboratory (LANL) Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) Laboratory Directed Research Development (LDRD) project, the authors are working on developing methods to calculate sensitivities to neutron leakage spectra. This may allow nuclear data evaluators to better use neutron leakage spectra data to constrain or adjust nuclear data. This may also be particularly useful as many neutron leakage spectra measurements do not require fissile material and can be performed with well characterized neutron sources. Neutron leakage spectra measurements are usually performed by placing a detector at the outside of a nuclear system. This may be a reactor, a neutron generator, or a neutron source. Certain detection systems can use pulse height data to infer neutron energies, other detector systems rely on other ways of determining neutron energy (for example, a Bonner sphere with multiple moderator thicknesses can be used to measure neutron spectrum). These measurements typically rely on some unfolding of the measured results. For pulsed systems like the "Livermore Pulsed Spheres," time-of-flight information can also help to constrain the neutron energy spectra data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Many-body theory of phonon-induced spin relaxation and decoherence

First-principles calculations enable accurate predictions of electronic interactions and dynamics. However, computing the electron spin dynamics remains challenging. The spin-orbit interaction causes various dynamical phenomena that couple with phonons, such as spin precession and spin-flip e-ph scattering, which are difficult to describe with current first-principles calculations. Here, in this work, we show a rigorous framework to study phonon-induced spin relaxation and decoherence, by computing the spin-spin correlation function and its vertex corrections due to e-ph interactions. We apply this approach to a model system and develop corresponding first-principles calculations of spin relaxation in GaAs. Our vertex-correction formalism is shown to capture the Elliott-Yafet, Dyakonov-Perel, and strong-precession mechanisms—three independent spin decoherence regimes with distinct physical origins—thereby unifying their theoretical treatment and calculation. Our method is general and enables quantitative studies of spin relaxation, decoherence, and transport in a wide range of materials and devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Analysis of Data Fusion Between Waveform Events and Radionuclide Detections Reported in 2021 by the International Data Centre

During this project, we assessed the percentage of events reported by the International Data Centre (IDC) that are candidates for evaluation using data fusion methods. We calculated the percentage of waveform (seismic, hydroacoustic, infrasound, SHI) events that may be associated with radionuclide (RN) detections and vice versa. For this task, data fusion means identifying RN detections that may be associated with SHI events. The task of producing a fused SHI and RN bulletin is included in the IDC operational manual. As of 2022, the production of a fused event bulletin remains incomplete, largely because methods for producing the bulletin are under development. Nonetheless, simple methods for data fusion have been developed and preliminary software tools are available at the IDC to associate SHI events and RN detections. This report is a first look at the percentage of SHI events that may be candidates for fusion with RN detections or bulletins. One year of IDC bulletin data (2021) is considered in this report. It is important to note that there were no known nuclear explosions during this time period. Hence, we expect that events during this time period would generate SHI and RN signals. If many SHI events are candidates for fusion with RN data, greater priority should be placed on establishing screening procedures before and/or during the fusion process. By contrast, if the number of candidate events is diminishingly small, then data fusion is not a driver for further development of event or detection screening at the IDC. LLNL reviewed the events reported in the Standard Event Bulletin (SEB) and determined which events are candidates for fusion with the Reviewed Radionuclide Reports (RRR). Much of the work involves establishing a metric or metrics for determining whether a Reviewed Event Bulletin (REB) event is a candidate for data fusion.

58 GEOSCIENCES↗
DOE OSTI2021

Revisiting adiabatic fraction estimations in cumulus clouds: high-resolution simulations with a passive tracer

The process of mixing in warm convective clouds and its effects on microphysics are crucial for an accurate description of cloud fields, weather, and climate. Still, they remain open questions in the field of cloud physics. Adiabatic regions in the cloud could be considered non-mixed areas and therefore serve as an important reference to mixing. For this reason, the adiabatic fraction (AF) is an important parameter that estimates the mixing level in the cloud in a simple way. Here, we test different methods of AF calculations using high-resolution (10 m) simulations of isolated warm cumulus clouds. The calculated AFs are compared with a normalized concentration of a passive tracer, which is a measure of dilution by mixing. This comparison enables the examination of how well the AF parameter can determine mixing effects and the estimation of the accuracy of different approaches used to calculate it. Comparison of three different methods to derive AF, with the passive tracer, shows that one method is much more robust than the others. Moreover, this method's equation structure also allows for the isolation of different assumptions that are often practiced when calculating AF such as vertical profiles, cloud-base height, and the linearity of AF with height. The use of a detailed spectral bin microphysics scheme allows an accurate description of the supersaturation field and demonstrates that the accuracy of the saturation adjustment assumption depends on aerosol concentration, leading to an underestimation of AF in pristine environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Many-Body Majorana Braiding without an Exponential Hilbert Space

Qubits built out of Majorana zero modes constitute the primary path toward topologically protected quantum computing. Simulating the braiding process of multiple Majorana zero modes corresponds to the quantum dynamics of a superconducting many-body system. It is crucial to study the Majorana dynamics both in the presence of all other quasiparticles and for reasonably large system sizes. We present a method to calculate arbitrary many-body wave functions as well as their expectation values, correlators, and overlaps from time evolved single-particle states of a superconductor, allowing for significantly larger system sizes. We calculate the fidelity, transition probabilities, and joint parities of Majorana pairs to track the quality of the braiding process. Here, we show how the braiding success depends on the speed of the braid. Moreover, we demonstrate the topological CNOT two-qubit gate as an example of two-qubit entanglement. Our Letter opens the path to test and analyze the many theoretical implementations of Majorana qubits. Moreover, this method can be used to study the dynamics of any noninteracting superconductor.

Majorana bound states↗
DOE OSTI2026

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗