Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difference scaling”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Feedstock/pretreatment screening for bioconversion of sugar and lignin streams via deacetylated disc-refining

Recent publications have shown the benefits of deacetylation disc-refining (DDR) as a pretreatment process to deconstruct biomass into sugars and lignin residues. Major advantages of DDR pretreatment over steam and dilute acid pretreatment are the removal of acetyl and lignin during deacetylation. DDR does not generate hydroxymethylfurfural (HMF) and furfural which are commonly produced from steam and dilute acid pretreatments. Acetate, lignin, HMF, and furfural are known inhibitors during enzymatic hydrolysis and fermentation. Another advantage of deacetylation is the production of lignin-rich black liquor, which can be upgraded to other bioproducts. Furthermore, due to the lack of sugar degradation during deacetylation, DDR has significantly less sugar loss than other pretreatment methods. Previous studies for DDR have primarily focused on corn stover, but lacked the investigative studies of other feedstocks. This study was designed to screen various DDR process conditions at pilot scale using three different feedstocks, including corn stover, poplar, and switchgrass. The impact of the pretreatment conditions was evaluated by testing hydrolysates for bioconversion to 2,3-butanediol. Pretreatment of biomass by DDR showed high-conversion-yields and 2,3-BDO fermentation production yields. Techno-economic analysis (TEA) of the pretreatment for biomass to sugar was also developed based on NREL’s Aspen Model. This study shows that the cellulose and hemicellulose in poplar was more recalcitrant than herbaceous feedstocks which ultimately drove up the sugar cost. Switchgrass was also more recalcitrant than corn stover but less than poplar.

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

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗
DOE OSTI2025

Reducing Frequency Bias of Fourier Neural Operators in 3D Seismic Wavefield Simulations Through Multistage Training

The recent development of neural operator (NeurOp) learning for solutions to the elastic wave equation shows promising results and provides the basis for fast large-scale simulations for different seismological applications. In this article, we use the Fourier neural operator (FNO) model to directly solve the 3D Helmholtz wave equation for fast seismic ground-motion simulations on different frequencies and show the frequency bias of the FNO model, that is, it learns the lower frequencies better comparing to the higher frequencies. To reduce the frequency bias, we adopt the multistage FNO training, that is, after training a stage 1 FNO model for estimating the ground motion, we use a second FNO model as the stage 2 to learn from the residual, which greatly reduced the errors on the higher frequencies. By adopting this multistage training, the FNO models show reduced biases on higher frequencies, which enhanced the overall results of the ground-motion simulations. Thus the multistage training FNO improves the accuracy and realism of the ground-motion simulations.

earthquakes↗
DOE OSTI2020

Single T gate in a Clifford circuit drives transition to universal entanglement spectrum statistics

Clifford circuits are insufficient for universal quantum computation or creating t t -designs with t\ge 4 t ≥ 4 . While the entanglement entropy is not a telltale of this insufficiency, the entanglement spectrum of a time evolved random product state is: the entanglement levels are Poisson-distributed for circuits restricted to the Clifford gate-set, while the levels follow Wigner-Dyson statistics when universal gates are used. In this paper we show, using finite-size scaling analysis of different measures of level spacing statistics, that in the thermodynamic limit, inserting a single T (\pi/8) ( π / 8 ) gate in the middle of a random Clifford circuit is sufficient to alter the entanglement spectrum from a Poisson to a Wigner-Dyson distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Cross Correlating Cosmological Probes for Stage-IV Surveys

The next few years will be populated with Stage-IV cosmological surveys. This paper forecasts the cosmological information we expect to constrain using probes from the Rubin Observatory Legacy Survey of Space and Time (LSST) and CMB Stage-4 surveys. We explore the effectiveness of different large-scale-structure probes, including the power spectra of galaxy weak lensing, galaxy positions, Cosmic Microwave Background (CMB) lensing, and their cross-correlation. We use correlated lognormal simulations with the expected redshift distribution, number density of galaxies, and noise levels of the LSST survey, as well as CMB-S4 lensing simulations with anticipated noise levels, to obtain the expected error bars of each of the probes. We investigate the constraining power of the cosmological parameters for the individual probes and their combinations in an idealized scenario, where we do not account for astrophysical and observational parameters such as galaxy bias variations and photometric redshift uncertainties.

Gibbins, Grace↗
DOE OSTI2023

FPGA Acceleration of GCN in Light of the Symmetry of Graph Adjacency Matrix

Graph Convolutional Neural Networks (GCNs) are widely used to process large-scale graph data. Different from deep neural networks (DNNs), GCNs are sparse, irregular, and unstructured, posing unique challenges to hardware acceleration with regular processing elements (PEs). In particular, the adjacency matrix of a GCN is extremely sparse, leading to frequent but irregular memory access, low spatial/temporal data locality and poor data reuse. Furthermore, a realistic graph usually consists of unstructured data (e.g., unbalanced distributions), creating significantly different processing times and imbalanced workload for each node in GCN acceleration. To overcome these challenges, we propose an end-to-end hardware-software co-design to accelerate GCNs on resource-constrained FPGAs with the features including: (1) A custom dataflow that leverages symmetry along the diagonal of the adjacency matrix to accelerate feature aggregation for undirected graphs. We utilize either the upper or the lower triangular matrix of the adjacency matrix to perform aggregation in GCN to improve data reuse. (2) Unified compute cores for both aggregation and transform phases, with full support to the symmetry-based dataflow. These cores can be dynamically reconfigured to the systolic mode for transformation or as individual accumulators for aggregation in GCN processing. (3) Preprocessing of the graph in software to rearrange the edges and features to match the custom dataflow. This step improves the regularity in memory access and data reuse in the aggregation phase. Moreover, we quantize the GCN precision from FP32 to INT8 to reduce the memory footprint without losing the inference accuracy. We implement our accelerator design in Intel Stratix10 MX FPGA board with HBM2, and demonstrate 1.3x-110.5x improvement in end-to-end GCN latency as compared to the state-of the-art FPGA implementations, on the graph datasets of Cora, Pubmed, Citeseer and Reddit.

Nair, Gopikrishnan R.↗
DOE OSTI2020

On Simulating the Formation of Structured, Crystalline Systems via Non-classical Pathways

Observations in crystal growth and assembly from recent in situ methods suggest alternative, non-classical crystallization pathways play an important role in the determination of the micro- and meso- structures in crystalline systems. These processes display parallels that cross-cut multiple disciplines investigating crystallization across four orders of magnitude in size scales and widely differing environments, hinting that alternative crystal growth pathways may be a fundamental scheme in natural crystal formation. Using a system of short-range attractive microbeads, we demonstrate that the addition of a small concentration of sub-species incommensurate with the lattice spacing of the dominant species results in a stark change in crystal size and morphology. These changes are attributed to the presence of fleeting, amorphous-like configurations of beads that ultimately change the melting and growth dynamics in preferred directions. From these real-time observations, we hypothesize the amorphous mineral precursors present in biological mineralized tissues undergo similar non-classical crystallization processes resulting in the complex structures found in biomineralization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental Validation of Exact Burst Pressure Solutions for Thick-Walled Cylindrical Pressure Vessels

Burst pressure is one of the critical strength parameters used in the design and operation of pressure vessels because it represents the maximum pressure that a vessel can withstand before failing. Historically, the Barlow formula was used as a design base for estimating burst pressure. However, it does not consider the plastic flow response for ductile steels and is applicable only to thin-walled cylinders (i.e., the diameter to thickness ratio D/t ≥ 20). A new multiaxial plastic yield theory was developed to consider the plastic flow response, and the associated theoretical (i.e., Zhu–Leis) solution of burst pressure was obtained and has gained extensive applications in the pipeline industry because it was validated by different full-scale burst test datasets for large-diameter, thin-walled pipelines in a variety of steel grades from Grade B to X120. The Zhu–Leis flow theory of plasticity was recently extended to thick-walled pressure vessels, and the associated exact flow solution of burst pressure was obtained and is applicable to both thin and thick-walled cylindrical shells. Many full-scale burst tests are available for thin-walled line pipes in the pipeline industry, but limited pressure burst tests exist for thick-walled vessels. To validate the newly developed exact solutions of burst pressure for thick-walled cylinders, this paper conducts a series of burst pressure tests on small-diameter, thick-walled pipes. In particular, six burst tests are carried out for three thick-walled pipes in Grade B carbon steel. These pipes have a nominal diameter of 2.375 inches (60.33 mm) and three nominal wall thicknesses of 0.154, 0.218, and 0.344 inches (3.91, 5.54, and 8.74 mm), leading to D/t = 15.4, 10.9, and 6.9, respectively. With the burst test data, comparisons show that the Zhu–Leis flow solution of burst pressure matches well the burst test data for thick-walled pipes. Thus, these burst tests validate the accuracy of the Zhu–Leis flow solution of burst pressure for thick-walled cylindrical vessels.

42 ENGINEERING↗
DOE OSTI2023

Strongly Heated Turbulent Flow in a Channel with Pin Fins

Large-eddy simulations (LES) were performed to study the turbulent flow in a channel of height H with a staggered array of pin fins with diameter D = H/2 as a function of heating loads that are relevant to the cooling of turbine blades and vanes. The following three heating loads were investigated—wall-to-coolant temperatures of Tw/Tc = 1.01, 2.0, and 4.0 - where the Reynolds number at the channel inlet was 10,000 and the back pressure at the channel outlet was 1 bar. For the LES, two different subgrid-scale models—the dynamic kinetic energy model (DKEM) and the wall-adapting local eddy-viscosity model (WALE)—were examined and compared. This study was validated by comparing with data from direct numerical simulation and experimental measurements. The results obtained show high heating loads to create wall jets next to all heated surfaces that significantly alter the structure of the turbulent flow. Results generated on effects of heat loads on the mean and fluctuating components of velocity and temperature, turbulent kinetic energy, the anisotropy of the Reynolds stresses, and velocity-temperature correlations can be used to improve existing RANS models.

42 ENGINEERING↗
DOE OSTI2020

Subspace Methods in Multi-Parameter Seismic Full Waveform Inversion

In full waveform inversion (FWI) high-resolution subsurface model parameters are sought. FWI is normally treated as a nonlinear least-squares inverse problem, in which the minimum of the corresponding misfit function is found by updating the model parameters. When multiple elastic or acoustic properties are solved for, simple gradient methods tend to confuse parameter classes. This is referred to as parameter cross-talk; it leads to incorrect model solutions, poor convergence and strong dependence on the scaling of the different parameter types. Determining step lengths in a subspace domain, rather than directly in terms of gradients of different parameters, is a potentially valuable approach to address this problem. The particular subspace used can be defined over a span of different sets of data or different parameter classes, provided it involves a small number of vectors compared to those contained in the whole model space. Additionally, in a subspace method, the basis vectors are defined first, and a local minimum is found in the space spanned by these. We examine the application of the subspace method within acoustic FWI in determining simultaneously updates for velocity and density. We first discuss the choice of basis vectors to construct the spanned space, from linear updates by distinguishing only the contributions of different parameter classes towards nonlinear updates by adding the contributions of higher-order perturbations of each parameter class. The numerical character of FWI solutions generated via subspace methods involving different basis vectors is then analyzed and compared with traditional FWI methods. The subspace methods can provide better reconstructions of the model, especially for the velocity, as well as improved convergence rates, while the computational costs are still comparable with the traditional FWI methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Burning plasma achieved in inertial fusion

The achievement of obtaining a burning plasma is a critical step toward self-sustaining fusion energy. A burning plasma is a fusion plasma where the alpha-particles created by the deuterium-tritium (DT) fusion reactions are the primary source of heating in the plasma, which is necessary to sustain and propagate the fusion reaction to enable high energy gain. After decades of fusion research, a burning plasma state has finally been achieved. Herein, we report upon the first burning-plasma experiments; this state was achieved using a strategy to increase the capsule spatial scale via two different implosion concepts, on the US National Ignition Facility. These experiments show energies from self-heating in excess of the mechanical work injected into the implosions satisfying several burning plasma metrics, the last experiment additionally shows that the fusion self-heating is greater than losses from radiation and heat conduction. These experiments triple the fusion yield performance and show significantly higher yield amplification from self-heating than prior results; remaining degradations can be reduced for even higher fusion performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

A theory of sea scatter at large incident angles

Radar sea scatter at large incident angles (30 deg or greater) is investigated using a model which takes into account theoretically-calculated small- and large-scale wave perturbations as well as measured sea surface slope densities and sea spectra obtained from empirical observations. Effects of incident angle, polarization, frequency, wind speed and anisotropic characteristics of the sea surface may be studied by application of the model. In particular, it is found that the increase in backscattering coefficients with wind velocity is due primarily to the growth of the sea spectrum and secondarily to the interaction between small- and large-scale perturbations. Furthermore, differences in upwind and downwind observations may be attributed to anisotropic properties of the sea spectrum. Theoretical findings derived from the model are compared with experimental data, and good agreement is found.

Chan, H. L.↗
NASA NTRS1977

In-flight gyro drift rate calibration on the Viking Orbiters

The drift rates of the attitude control gyros onboard the Viking Orbiters were calibrated several times during flight. The calibration was performed by engaging the gyro control of attitude for a six-hour period, recording the Viking Orbiter coordinates of the sun and the roll reference star as a function of time, and processing these data by a computer program. The computer program calculated drift rates for each data increment, and then smoothed the calculated rate versus the time function for all the data on the basis of a Gauss-Markov gyro-drift model. The effects of data noise were minimized by using a relatively long time constant in the Gauss-Markov model and by empirically determining a set of corrections for celestial inertial scale factor error differences. In this way the rms error of the data fit was lessened.

Breckenridge, W. G.↗
NASA NTRS1980

Time-average fluxes of heavy ions at the geostationary orbit

Measurements of average proton, helium, carbon, and oxygen fluxes at 6.6 earth radii are reported. The data represent averages obtained on ATS 6 between June 15 and Oct. 3, 1974. The energy range covered was 0.36-1.1 MeV for protons and approximately 1-4 MeV for the heavier ions. The results indicate that above about 1 MeV the heavier ion fluxes dominate over the proton flux on the energy/ion scale. Using two different spectral dependencies to fit the data, the carbon to oxygen concentration ratio for energies above 1 MeV was found to be 0.43 for the power law spectrum and 0.44 for the exponential spectrum. Thus in either case the abundance ratio is consistent with the solar origin of the particles. Similarly, the helium to oxygen concentration ratio is found to be 0.17 for the power law spectrum and 0.28 for the exponential spectrum. This is inconsistent with the solar wind ratio, which is about 10. The results quoted above are based on a very small portion of the distribution function above 1 MeV.

Konradi, A.↗
NASA NTRS1981

Metal-halide mixtures for latent heat energy storage

Alkali metal and alkali halide mixtures are identified which may be suitable for thermal energy storage at temperatures above 600 C. The use of metal-halides is appropriate because of their tendency to form two immiscible melts with a density difference, which reduces scale formation and solidification on heat transfer surfaces. Also, the accumulation of phase change material along the melt interface is avoided by the self-dispersing characteristic of some metal-halides, in particular Sr-SrCl2, Ba-BaCl2, and Ba-BaBr2 mixtures. Further advantages lie in their high thermal conductivities, ability to cope with thermal shock, corrosion inhibition, and possibly higher energy densities.

Chen, K.↗
NASA NTRS1983

Charge-exchange in the magnetosheaths of Venus and Mars - A comparison

The amount of solar wind absorption due to charge-exchange in the Martian magnetosheath is evaluated and found to be about an order of magnitude less than that in the Venus magnetosheath. This difference might explain the observed difference in the scaled position and shape between the shocks at Venus and Mars. The lower solar wind absorption for Mars is attributable to the less dense hot oxygen corona of Mars compared to Venus.

Russell, C. T.↗
NASA NTRS1983

Theoretical quasar emission-line profiles. II - Effects of source function vibration and total thickness

By calculating the actual limb brightening functions for the emission-line clouds, rather than using the approximation that all observed radiation originates from the same part of the growth curve, the present consideration of radiation transfer effects in the broad emission line regions of quasars and Seyfert 1 galaxies allows the study of clouds whose total optical thickness is intermediate between the limiting cases of an optically thin Doppler core and Voigt wing emissions. The ability to predict minor differences in observed line profiles in a given object in terms of radiative transfer effects leads to the proposal of a test which may be able to distinguish between radial and orbital geometries. The test involves the difference in width between emission lines from the same family having different optical depth scales.

Puetter, R. C.↗