Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Step”

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

Beyond Taylor’s hypothesis: a novel volumetric reconstruction of velocity and density fields for variable‑density and shear flows

Abstract This work presents a novel numerical procedure for reconstructing volumetric density and velocity fields from planar laser-induced fluorescence (PLIF) and stereoscopic particle image velocimetry (SPIV) data. This new method is theoretically and practically demonstrated to provide more accurate 3D vortical structures and density fields in high shear flows than reconstruction methods based on the mean convective velocity. While Taylor’s hypothesis of frozen turbulence is commonly applied by using the local mean streamwise velocity, the proposed algorithm uses the measured local instantaneous velocity for data convection. It consists of a step-by-step reconstruction based on a mixed Lagrangian–Eulerian solver that includes the 3D interpolation of scattered flow data and that relaxes the Taylor’s hypothesis by iterative enforcement of the incompressibility constraint on the velocity field. This methodology provides 3D fields with temporal resolution, spatial resolution, and accuracy comparable to that of real 3D snapshots, thus providing a practical alternative to tomographic measurements. The procedure is validated using numerical data of the constant-density channel flow available on the Johns Hopkins University Turbulence Database (JHTDB), showing the accurate reconstruction of the 3D velocity field. The algorithm is applied to an experimental dataset of PLIF and SPIV measurements of a variable-density jet flow, demonstrating its capability to provide 3D velocity and density fields that are more consistent with the Navier–Stokes equations compared to the mean flow convective method. Graphic abstract

42 ENGINEERING↗
DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Accretion disks, quasars and cosmology: meandering towards understanding

Abstract As Setti and Woltjer noted back in 1973, one can use quasars to construct the Hubble diagram; however, the actual application of the idea was not that straightforward. It took years to implement the proposition successfully. Most ways to employ quasars for cosmology now require an advanced understanding of their structure, step by step. We briefly review this progress, with unavoidable personal biases, and concentrate on bright unobscured sources. We will mention the problem of the gas flow character close to the innermost stable circular orbit near the black hole, as discussed five decades ago. This problem later led to the development of the slim disk scenario and is recently revived in the context of Magnetically Arrested Disks (MAD) and Standard and Normal Evolution (SANE) models. We also discuss the hot or warm corona issue, which is still under debate and complicates the analysis of X-ray reflection. We present the scenario of the formation of the low ionization part of the Broad Line Region as a failed wind powered by radiation pressure acting on dust (Failed Radiatively Driven Dusty Outflow – FRADO). Next, we examine the cosmological constraints currently achievable with quasars, primarily concentrating on light echo methods (continuum time delays and spectral-line time delays to the continuum) that are (or should be) incorporating the progress mentioned above. Finally, we briefly discuss prospects in this lively subject area.

Astronomy & Astrophysics↗
DOE OSTI2020

High precision control and deep learning-based corn stand counting algorithms for agricultural robot

This paper presents high precision control and deep learning-based corn stand counting algorithms for a low-cost, ultra-compact 3D printed and autonomous field robot for agricultural operations. Currently, plant traits, such as emergence rate, biomass, vigor, and stand counting, are measured manually. This is highly labor-intensive and prone to errors. The robot, termed TerraSentia, is designed to automate the measurement of plant traits for efficient phenotyping as an alternative to manual measurements. In this paper, we formulate a Nonlinear Moving Horizon Estimator that identifies key terrain parameters using onboard robot sensors and a learning-based Nonlinear Model Predictive Control that ensures high precision path tracking in the presence of unknown wheel-terrain interaction. Moreover, we develop a machine vision algorithm designed to enable an ultra-compact ground robot to count corn stands by driving through the fields autonomously. The algorithm leverages a deep network to detect corn plants in images, and a visual tracking model to re-identify detected objects at different time steps. We collected data from 53 corn plots in various fields for corn plants around 14 days after emergence (stage V3 - V4). The robot predictions have agreed well with the ground truth with C robot =1.02×C human -0.86 and a correlation coefficient R=0.96. The mean relative error given by the algorithm is -3.78%, and the standard deviation is 6.76%. These results indicate a first and significant step towards autonomous robot-based real-time phenotyping using low-cost, ultra-compact ground robots for corn and potentially other crops.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A holistic framework integrating plant-microbe-mineral regulation of soil bioavailable nitrogen

Abstract Soil organic nitrogen (N) is a critical resource for plants and microbes, but the processes that govern its cycle are not well-described. To promote a holistic understanding of soil N dynamics, we need an integrated model that links soil organic matter (SOM) cycling to bioavailable N in both unmanaged and managed landscapes, including agroecosystems. We present a framework that unifies recent conceptual advances in our understanding of three critical steps in bioavailable N cycling: organic N (ON) depolymerization and solubilization; bioavailable N sorption and desorption on mineral surfaces; and microbial ON turnover including assimilation, mineralization, and the recycling of microbial products. Consideration of the balance between these processes provides insight into the sources, sinks, and flux rates of bioavailable N. By accounting for interactions among the biological, physical, and chemical controls over ON and its availability to plants and microbes, our conceptual model unifies complex mechanisms of ON transformation in a concrete conceptual framework that is amenable to experimental testing and translates into ideas for new management practices. This framework will allow researchers and practitioners to use common measurements of particulate organic matter (POM) and mineral-associated organic matter (MAOM) to design strategic organic N-cycle interventions that optimize ecosystem productivity and minimize environmental N loss.

Environmental Sciences & Ecology↗
DOE OSTI2019

Developing AlMn Films for Argonne TES Fabrication

The reference design for the next-generation cosmic microwave background (CMB) experiment, CMB-S4, relies on large arrays of transition-edge sensor (TES) bolometers coupled to Superconducting Quantum Interference Device (SQUID)-based readout systems. Mapping the CMB to near cosmic variance limits will enable the search for signatures of inflation and constrain dark energy and neutrino physics. AlMn TESes provide simple film manufacturing and highly uniform arrays over large areas to meet the requirements of the CMB-S 4 experiment. TES parameters such as critical temperature and normal resistance must be tuned to experiment specifications and can be varied based on geometry and steps in the fabrication process such as deposition layering, geometry, and baking time and temperature. Using four-terminal sensing, we measured T C and R N of AlMn 2000 ppm films and devices of varying thicknesses fabricated at Argonne National Laboratory to motivate device geometries and fabrication processes to tune T C to 150–200 mK and R N to ~10 mΩ. Furthermore, measurements of IV curves and time constants for the resulting devices of varying leg length were made using time-division SQUID multiplexing and determined TC, G, k, f 3db , and RN. We present the results of these tests along with the geometries and fabrication steps used to tune the device parameters to the desired limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Stability Analysis of Interface Conditions for Ocean-Atmosphere Coupling

In this paper, we analyze the stability of different coupling strategies for multidomain PDEs that arise in general circulation models used in climate simulations. We focus on fully coupled ocean–atmosphere models that are needed to represent and understand the complicated interactions of these two systems, becoming increasingly important in climate change assessment in recent years. Numerical stability issues typically arise because of different time-stepping strategies applied to the coupled PDE system. In particular, the contributing factors include using large time steps, lack of accurate interface flux, and single-iteration coupling. We investigate the stability of the coupled ocean–atmosphere models for various interface conditions such as the Dirichlet–Neumann condition and the bulk interface condition, which is unique to climate modeling. By analyzing a simplified model, we demonstrate here how the parameterization of the bulk condition and other numerical and physical parameters affect the coupling stability and establish stability conditions for different coupling strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Unveiling the microstructural evolution of carbon fibers derived from polyamide-6

Polyacrylonitrile-based carbon fibers have dominated the industry for decades, but the high cost of polyacrylonitrile has prevented the widespread adoption of carbon fiber in high-volume structural applications. As such, a significant amount of research has been dedicated to finding an alternative, low-cost carbon fiber precursor. In this work, carbon fibers were produced from polyamide-6 using metal salt impregnation and a thermo-oxidative stabilization step. To gain further insight into the carbonization process and microstructural transformation, the morphologies, crystallinities, elemental compositions, and thermal stabilities of the fibers were characterized at various stages of processing. The stabilization step resulted in a significant increase in carbon yield, indicating a dramatic increase in thermal stability. This is due to the crosslinking of polyamide-6 chains, which was confirmed by functional group analysis. The crystallinity of the fibers was also significantly altered during processing, as the produced carbon fibers consisted of pseudo-amorphous carbon with two distinct regions of metal salt impregnation. Finally, the findings and microstructural evolution mechanisms provide guidelines for further research into carbon fiber produced from polyamide-6.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molar Absorptivities of U(VI), U(IV), and Pu(III) in Nitric Acid Solutions of Various Concentrations Relevant to Developing Nuclear Fuel Recycling Flowsheets

Testing of a co-decontamination (CoDCon) tributyl phosphate (TBP) based solvent extraction flowsheet is being performed to optimize conditions for achieving a target U/Pu ratio in the Pu-containing product. The flowsheet is designed to directly produce a mixture of U and Pu with a U/Pu mass ratio of 7/3. The system is monitored in real time using optical spectroscopy, which allows for control of the U/Pu ratio in the aqueous U/Pu nitrate product. On-line spectrophotometric monitoring of the aqueous stream at the reductive stripping step involves determination of U(VI), U(IV) and Pu(III) concentrations and relies on molar absorptivities of these species in nitric acid at a number of wavelengths corresponding to their peaks’ maxima. The magnitudes of these molar absorptivities are a sensitive function of nitric acid concentration and have to be determined as precisely as possible using an off-line spectrophotometry under well controlled conditions. This step (called training set acquisition) is typically performed not in a high contamination area of radiological glovebox where flow-through cells and other on-line equipment are installed but in a less aggressive environment of a radiological fume hood with much better control of stock solutions quality, dilution factors, optical absorbance scale calibration, etc. This is followed by the calibration transfer between the off-line instrument and the instrument deployed during extraction and stripping runs. This paper reports values of molar absorptivities of U(VI), U(IV), and Pu(III) in nitric acid solution of 0.5 M to 4 M concentration and compares them with available technical literature data. Results of off-line spectrophotometric analysis of selected aqueous grab samples from one of CoDCon runs show satisfactory agreement with total uranium and plutonium concentrations determined in the same samples by ICP-MS.

Sinkov, Sergey I.↗
DOE OSTI2022

Column separation of tetravalent cerium fission products from trivalent rare earth radio-isotopes

Rapid and efficient isolation of individual rare earth element (REE) radioisotopes from complex mixtures is necessary to support the fields of nuclear forensics, medical isotope production and nuclear physics measurements. The separations must be robust and generate sufficiently high purity samples for subsequent radiological analysis. Current methodologies utilize a laborious two step Ce(IV)-iodate precipitation followed by Ce(IV)-nitrate di-(2-ethylhexyl)phosphoric acid/n-heptane extraction. This work reports an alternative method to isolate Ce fission product isotopes equilibrated with milligrams of cold carrier, from the remaining REE’s using a Ce(IV)/LN-resin (Eichrom Technologies) based separation. Optimization of Ce(IV) loading, REE(III) elution volume and subsequent Ce(III) stripping step utilized 166m Ho tracer. Mixed fission product experiments provided a further test of this LN-resin method, with analysis of purified 141/144 Ce, 147 Nd and 91 Y. Overall, while decontamination factors were lower for the LN-resin method, (vs. the combined precipitation/extraction) the simplicity of the technique and the satisfactory fission product analysis results point to the efficacy of this method.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Use of a Pressurized Water Treatment to Prevent Cracking of Internal Gelation Sol-Gel Microspheres

Pressurized water treatments were investigated as a means to prevent cracking of internal gelation sol-gel microspheres upon drying and heating. A batch of gelled cerium oxide microspheres was prepared, split into two samples that were processed with and without a pressurized water treatment, and dried in air at room temperature. Thermal gravimetric analysis was used to compare the mass loss profiles between the two samples and to understand volatile emissions from cerium oxide microspheres during heat treatments up to 600 °C. Microspheres that were processed without a pressurized water treatment lost approximately 4.6 times more mass upon heating than microspheres that were subjected to a pressurized water treatment. Comparison of the cracking behavior of microspheres prepared with and without a pressurized water treatment indicated that microspheres subjected to a pressurized water treatment were virtually crack-free after heat treatments. Conversely, microspheres processed without a pressurized water treatment, which is typical of traditional washing processes, exhibited extensive cracking. Analysis of residues dried from pressurized water treatment effluents indicated that hexamethylenetetramine and urea, which are excess reactants from the internal gelation process, were removed from gelled microspheres during the pressurized water treatment step. Traditional processing methods for preparing internal gelation microspheres use specific feed solution and gelation parameters to reduce subsequent cracking upon drying and heating. Results from this study indicate that more aggressive gel-forming conditions may be used to produce crack-free microspheres, if a pressurized water treatment step is included after washing and before air-drying.

cerium oxide, internal gelation, microsphere, TRIS↗
DOE OSTI2021

Effects of a Pressurized Water Treatment on Internal Gelation Sol-Gel Microspheres

Based on previous observations that a pressurized water treatment (PWT) prevented cracking of sol-gel microspheres, we investigated the effects of a PWT on microsphere crystallinity, density, and specific surface area (SSA). Results were used to determine how a PWT alters the properties of microspheres upon drying and heating. Microspheres with diameters near 100-200 µm were prepared with and without a PWT and measured using x-ray diffraction (XRD), transmission electron microscopy (TEM), nitrogen adsorption (BET), and pycnometry. Properties of air-dried microspheres processed with and without a PWT are compared. Additionally, the properties of microspheres processed using a PWT are reported after heating to 150 °C, 450 °C, 750 °C, 1050 °C, and 1350 °C. X-ray diffraction measurements indicate that the PWT step improves the crystallinity of air-dried microspheres. XRD data was also used to calculate crystallite size, which increases with higher heat treatment temperatures. TEM images support crystallite size calculations from XRD data and provide an indication of the range of crystallite sizes, particularly for samples processed at higher temperatures where crystallite sizes are too large for estimation using the Scherrer formula. Density and SSA measurements performed as a function of heat treatment temperature indicate that a PWT increases the density of air-dried microspheres, creates a pore network, and that significant densification occurs between 450 - 750 °C. These results may be used to inform decisions on internal gelation flowsheet parameters to optimize the microsphere formation and gelation step, prevent microsphere cracking, and produce microspheres suitable for subsequent coating operations or pressing into pellets.

cerium oxide, TRISO, internal gelation, nuclear fu↗
DOE OSTI2024

Assessing compounding climate-related stresses and development pathways on the power sector in the central U.S.

Future configurations of the power system in the central region of the USA are dependent on relative costs of alternative power generation technologies, energy and environmental policies, and multiple climate-induced stresses. Higher demand in the summer months combined with compounding supply shocks in several power generation technologies can potentially cause a “perfect storm” leading to failure of the power system. Potential future climate stress must be incorporated in investment decisions and energy system planning and operation. We assess how projected future climate impacts on the power system would affect alternative pathways for the electricity sector considering a broad range of generation technologies and changes in demand. We calculate a “potential supply gap” metric for each pathway, system component, and sub-region of the US Heartland due to climate-induced effects on electricity demand and power generation. Potential supply gaps range from 5% in the North Central region under mild changes in climate to 21% in the Lakes-Mid Atlantic region under more severe climate change. We find increases in electricity demand to be more important in determining the size of the potential supply gap than stresses on power generation, while larger shares of renewables in the power system contribute to lower supply gaps. Our results provide a first step toward considering systemic climate impacts that may require changes in managing the grid or on potential additional capacity/reserves that may be needed.

Gurgel, Angelo Costa↗
DOE OSTI2020

Capabilities of multivariate Bayesian inference toward seismic hazard assessment

Multivariate Bayesian analysis can bring significant benefits to seismic hazard analysis: Its multivariate feature enables computing scalar and vector hazard without making any approximations; Correlations between intensity measures are implicitly modeled, permitting direct simulation of ground motion selection tools such as the conditional mean spectrum and the generalized conditioning intensity measure; and Its updating feature enables a seamless integration of new ground motion data into the hazard results. Here, we first develop a multivariate Bayesian ground motion model through the NGA-West2 database. The model functional form considers fault-type, magnitude, and distance dependencies, and also the linear and the rock intensity dependent site response. We use a hybrid Markov Chain Monte Carlo sampling to perform Bayesian inference consisting of Gibbs step and a multilevel Metropolis-Hastings step. We then perform several checks on the model and note that its performance is satisfactory. Finally, we illustrate the merits of this multivariate Bayesian analysis, which include: ground motion model updating with ground motion data recorded in the last four years not part of the NGA-West2 database; computation of scalar and vector seismic hazard using the un-updated and updated ground motion models for Los Angeles, CA; and simulation of the conditional mean spectrum under scalar and vector IM conditioning while accounting for different sources of aleatoric and epistemic uncertainties.

58 GEOSCIENCES↗
DOE OSTI2024

Large-scale harmonic balance simulations with Krylov subspace and preconditioner recycling

The multi-harmonic balance method combined with numerical continuation provides an efficient framework to compute a family of time-periodic solutions, or response curves, for large-scale, nonlinear mechanical systems. The predictor and corrector steps repeatedly solve a sequence of linear systems that scale by the model size and number of harmonics in the assumed Fourier series approximation. In this paper, a novel Newton–Krylov iterative method is embedded within the multi-harmonic balance and continuation algorithm to efficiently compute the approximate solutions from the sequence of linear systems that arise during the prediction and correction steps. Further, the method recycles, or reuses, both the preconditioner and the Krylov subspace generated by previous linear systems in the solution sequence. A delayed frequency preconditioner refactorizes the preconditioner only when the performance of the iterative solver deteriorates. The GCRO-DR iterative solver recycles a subset of harmonic Ritz vectors to initialize the solution subspace for the next linear system in the sequence. The performance of the iterative solver is demonstrated on two exemplars with contact-type nonlinearities and benchmarked against a direct solver with traditional Newton–Raphson iterations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter 0 < ε &NestedLessLess; 1 inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at O(ε), while in space, it is at O(1) with the wave speed at O(ε -1 ). We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step τ, mesh size h, and the small parameter ε. The error bounds indicate that, to obtain the “correct” numerical solution in the massless and nonrelativistic regime, i.e., 0 < ε &NestedLessLess; 1, all these finite difference methods share the same ε-scalability as time step τ = O(ε 3/2 ) and mesh size h = O(ε 1/2 ). A large number of numerical results are reported to verify the error estimates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Enhancing scalability and accuracy of quantum poisson solver

The Poisson equation has many applications across the broad areas of science and engineering. Most quantum algorithms for the Poisson solver presented so far either suffer from lack of accuracy and/or are limited to very small sizes of the problem and thus have no practical usage. In this regard, our previous work showed a proof-of-concept demonstration in advancing quantum Poisson solver algorithm and validated preliminary results for a simple case of 3 x 3 problem. In this work, we delve into comprehensive research details, presenting the results on up to 15 x 15 problems that include step-by-step improvements in Poisson equation solutions, scaling performance, and experimental exploration. In particular, we demonstrate the implementation of eigenvalue amplification by a factor of up to 2 8 , achieving a significant improvement in the accuracy of our quantum Poisson solver and comparing that to the exact solution. Additionally, we present success probability results, highlighting the reliability of our quantum Poisson solver. Moreover, we explore the scaling performance of our algorithm against the circuit depth and width, demonstrating how our approach scales with larger problem sizes and thus further solidifies the practicality of easy adaptation of this algorithm in real-world applications. We also discuss a multilevel strategy for how this algorithm might be further improved to explore much larger problems with greater performance. Finally, through our experiments on the IBM quantum hardware, we conclude that though overall results on the existing NISQ hardware are dominated by the error in the CNOT gates, this work opens a path to realizing a multidimensional Poisson solver on near-term quantum hardware.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Scalable multiscale modeling of platelets with 100 million particles

Here, we developed the core components of the AI-aided multiple time stepping algorithm for multiscale modeling of cell dynamics. This algorithm was implemented and analyzed on two supercomputer architectures with an application of simulating the aggregation of 250 platelets, or 102 million particles. To scale on these computers with complex memory and network architectures with GPUs, we devised a biomechanics-informed task mapping scheme to optimize load imbalance, communications, and memory utilization. Our simulations, scaling well up to 192 nodes on a Summit-like supercomputer with a peak speed of 11 petaflops, achieved a rate of 423 μs/day which is 500 times faster than the conventional algorithm using static time step and this has enabled studies of record size blood clots at record spatial–temporal resolutions. Additionally, we discovered the sensitive dependence of the scalability and execution time on the methods of decomposition, CPU–GPU coupling, and task mapping.

97 MATHEMATICS AND COMPUTING↗