Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time 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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct reconstruction of isolated XUV or soft x-ray attosecond pulses from high-harmonic generation streaking spectra

Characterization of an isolated attosecond pulse (IAP) in the extreme ultraviolet (XUV) or soft x-ray (SXR) region is essential for its applications. Here we propose to retrieve an IAP in the time domain directly through the modulation of high-harmonic generation (HHG) spectra in the presence of a time-delayed intense few-cycle infrared or mid-infrared laser. The retrieval algorithm is derived based on the strong-field approximation and an extended quantitative rescattering model. We show that both isolated XUV pulses with a narrow spectral bandwidth and isolated SXR pulses with a broad bandwidth can be well characterized through the HHG streaking spectra. Such an all-optical method for characterizing the IAP differs from the commonly used approach based on the streaked photoelectron spectra that would require electron spectrometers. We check the robustness of the retrieval method by changing the dressing laser or by adjusting the steps of time delay. We also show that the XUV pulse can be accurately retrieved by treating the HHG streaking spectra calculated from solving the time-dependent Schrödinger equation for single atoms as the ‘experimental’ data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Optimal Electric Grid Black Start Restoration Subject to Intentional Threats

Efficient restoration of the electric grid from significant disruptions – both natural and manmade – that lead to the grid entering a failed state is essential to maintaining resilience under a wide range of threats. Restoration follows a set of black start plans, allowing operators to select among these plans to meet the constraints imposed on the system by the disruption. Restoration objectives aim to restore power to a maximum number of customers in the shortest time. Current state-of-the-art for restoration modeling breaks the problem into multiple parts, assuming a known network state and full observability and control by grid operators. These assumptions are not guaranteed under some threats. This paper focuses on a novel integration of modeling and analysis capabilities to aid operators during restoration activities. A power flow-informed restoration framework, comprised of a restoration mixed-integer program informed by power flow models to identify restoration alternatives, interacts with a dynamic representation of the grid through a cognitive model of operator decision-making, to identify and prove an optimal restoration path. Application of this integrated approach is illustrated on exemplar systems. Validation of the restoration is performed for one of these exemplars using commercial solvers, and comparison is made between the steps and time involved in the commercial solver, and that required by the restoration optimization in and of itself, and by the operator model in acting on the restoration optimization output. Publications and proposals developed under this work, along with a path forward for additional expansion of the work, and summary of what was achieved, are also documented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Constraints on Future Analysis Metadata Systems in High Energy Physics

In high energy physics (HEP), analysis metadata comes in many forms—from theoretical cross-sections, to calibration corrections, to details about file processing. Correctly applying metadata is a crucial and often time-consuming step in an analysis, but designing analysis metadata systems has historically received little direct attention. Among other considerations, an ideal metadata tool should be easy to use by new analysers, should scale to large data volumes and diverse processing paradigms, and should enable future analysis reinterpretation. This document, which is the product of community discussions organised by the HEP Software Foundation, categorises types of metadata by scope and format and gives examples of current metadata solutions. Important design considerations for metadata systems, including sociological factors, analysis preservation efforts, and technical factors, are discussed. A list of best practices and technical requirements for future analysis metadata systems is presented. These best practices could guide the development of a future cross-experimental effort for analysis metadata tools.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗
DOE OSTI2023

Synthesis of Na + -gated nanochannel membranes for the ammonia (NH 3 ) separation

Na + -gated nanochannel membranes were synthesized on α-Al 2 O 3 supports by secondary growth for ammonia (NH 3 ) separation from a mixture gas containing NH 3 , hydrogen (H 2 ) and nitrogen (N 2 ) at high temperature and pressure. Here, the effects of synthesis parameters (number of seeding steps, aging time, synthesis time, and Si/Al ratio) on membrane properties were investigated by characterization and gas permeation testing. Under the optimized condition (two-time seeding, 8 h aging time, 5 h synthesis time, and gel composition of 1.05 Al 2 O 3 : 5 SiO 2 : 50 Na 2 O: 100 H 2 O), the best membranes achieved NH 3 /H 2 and NH 3 /N 2 selectivities of 328 ± 60 and 1106 ± 207, respectively at 200 °C and 21 bar. Separation performance of the membrane was shown to be highly dependent on the permeances of N 2 and H 2 , which were greatly influenced by the membrane quality.

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

Photocurrent in carbon nanotube bundle: Graded Seebeck coefficient phenomenon

Here, polarized photovoltage of a suspended aligned carbon nanotube (CNT) bundle under uniform optical irradiation is discovered without additional structural modification or bias voltage. Such a phenomenon is very surprising considering the metallic behavior of the overall bundle and zero temperature difference between ends. The photovoltage characteristic time is found similar to the thermal response time under step Joule heating and implies a relation to the thermal behavior of the CNT bundle. A similar thermoelectric voltage is also observed during step Joule heating. Localized laser heating and scanning along the axial direction of the bundle uncovers a linear spatial variation of the local Seebeck coefficient. The Seebeck coefficient linearly decreases from root to tip of the CNT bundle with a rate of a few - µV∙K –1 ·mm –1 . Deep investigation in both the microscopic and macroscopic structures of the CNT bundle reveals that the local alignment of CNT assemblies rather than the minor defects in individual CNTs brings about this linear distribution of Seebeck coefficient in space. The finding presents a new way for direct photon-to-electric energy conversion via Seebeck coefficient grading in CNT structures.

42 ENGINEERING↗
DOE OSTI2021

Elucidating Electron Storage and Distribution within the Pentaheme Scaffold of Cytochrome c Nitrite Reductase (NrfA)

Cytochrome c nitrite reductases (CcNIR or NrfA) play important roles in the global nitrogen cycle by conserving the usable nitrogen in the soil. Here, the electron storage and distribution properties within the pentaheme scaffold of Geobacter lovleyi NrfA were investigated via EPR spectroscopy coupled with chemical titration experiments. Initially, a chemical reduction method was established to sequentially add electrons to the fully oxidized protein, one equivalent at a time. The step–by–step reduction of the hemes was then followed using UV–Vis and EPR spectroscopy. EPR–spectral simulations elucidate the sequence of heme reduction within the pentaheme scaffold, and identify the signals of all five hemes in the EPR spectra. Electrochemical experiments ascertain the reduction potentials for each heme, observed in a narrow +11 (Heme 5) to –212 mV (Heme 3) range (vs SHE). Based on quantitative analysis/simulation of the EPR data, we demonstrate that Hemes 4 and 5 are reduced first (before the active site Heme 1) and serve the purpose of an electron storage unit within the protein. To probe the role of the central Heme 3, a H108M NrfA variant was generated where the reduction potential of Heme 3 is shifted positively (from –212 mV to +29 mV). Here, the H108M mutation significantly impacts the distribution of electrons within the pentaheme scaffold and the reduction potentials of the hemes, reducing the catalytic activity of the enzyme to 3% compared to wild–type. We propose that this is due to Heme 3’s important role as an electron gateway in the wild–type enzyme.

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

Photochemical and Photophysical Dynamics of the Aqueous Ferrate(VI) Ion

Ferrate(VI) has the potential to play a key role in future water supplies. Furthermore, its salts have been suggested as “green” alternatives to current advanced oxidation and disinfection methods in water treatment, especially when combined with ultraviolet light to stimulate generation of highly oxidizing Fe(V) and Fe(IV) species. However, the nature of these intermediates, the mechanisms by which they form, and their roles in downstream oxidation reactions remain unclear. Here, we use a combination of optical and X-ray transient absorption spectroscopies to study the formation, interconversion, and relaxation of several excited-state and metastable high-valent iron species following excitation of aqueous potassium ferrate(VI) by ultraviolet and visible light. Branching from the initially populated ligand-to-metal charge transfer state into independent photophysical and photochemical pathways occurs within tens of picoseconds, with the quantum yield for the generation of reactive Fe(V) species determined by relative rates of the competing intersystem crossing and reverse electron transfer processes. Relaxation of the metal-centered states then occurs within 4 ns, while the formation of metastable Fe(V) species occurs in several steps with time constants of 250 ps and 300 ns. Results here improve the mechanistic understanding of the formation and fate of Fe(V) and Fe(IV), which will accelerate the development of novel advanced oxidation processes for water treatment applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Scanning Transmission Electron Microscopy–Atom Probe Tomography Correlative Analysis for the Characterization of Solute-defect Interactions

Atom probe tomography (APT) and (scanning) transmission electron microscopy ((S)TEM) are complementary techniques that provide spatially resolved chemical and structural information at the atomic scale. Here, in this study, we employ two different STEM/APT correlative analysis methods to investigate Cr segregation at dislocation loops in ultra-high purity Fe–Cr alloys. APT needles for the correlative analysis were extracted either from bulk material or from thinned TEM lamellae. STEM analysis was used to determine the Burgers vectors of ion-irradiation-induced dislocation loops, while APT reconstruction of the same region revealed the Cr segregation to these loops. We extended the g•b = 0 invisibility criterion of dislocation loops from TEM mode in a lamella to STEM mode in a needle-shaped specimen. STEM and APT analysis on the same needle provide straightforward correlative analysis, although it is limited by a small observation volume. In contrast, iterative STEM analysis of TEM lamellae, followed by the selective extraction of specific regions of interest for APT analysis, expands the observation area by up to 100 times but requires additional time-consuming steps for APT needle extraction from the lamellae.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Quantum algorithms for disordered physics

We show how a quantum computer may efficiently simulate a disordered Hamiltonian, by incorporating a pseudo-random number generator directly into the time evolution circuit. This technique is applied to quantum simulation of few-body disordered systems in the large volume limit; in particular, Anderson localization. The method requires a number of (error corrected) qubits proportional to the logarithm of the volume of the system, and each time evolution step requires a number of gates polylogarithmic in the volume. We simulate the method to observe the metal-insulator transition on a three-dimensional lattice. Additionally, we demonstrate the algorithm on a one-dimensional lattice, using physical quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Multilevel Parareal Algorithm with Averaging for Oscillatory Problems

The present study is an extension of the work done by Peddle, Haut, and Wingate and Haut and Wingate, where a two-level Parareal method with mapping and averaging is examined. The method proposed in this paper is a multilevel Parareal method with arbitrarily many levels, which is not restricted to the two-level case. We give an asymptotic error estimate which reduces to the two-level estimate for the case when only two levels are considered. Introducing more than two levels has important consequences for the averaging procedure, as we choose separate averaging windows for each of the different levels, which is an additional new feature of the present study. The different averaging windows make the proposed method especially appropriate for nonlinear multiscale problems, because we can introduce a level for each intrinsic scale of the problem and adapt the averaging procedure such that we reproduce the behavior of the model on the particular scale resolved by the level. The method is applied to nonlinear differential equations. The nonlinearities can generate a range of frequencies in the problem. The computational cost of the new method is investigated and studied on several examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Transmission Innovation Symposium: Modernizing the U.S. Electrical Grid

The foundation of the United States Department of Energy (DOE) Transmission Reliability research program was established 20 years ago by a series of commissioned white papers. Those white papers described the dramatic institutional and regulatory changes that the U.S. electricity transmission grid was undergoing at the time and articulated the technical challenges that these changes created. The challenges outlined in the white papers were the basis for the initial research goals of the DOE Transmission Reliability program. To a large extent, the reliability research needs outlined in the original white papers have now been met. As a result, now is an appropriate time to step back and review the technical challenges that the industry currently faces and to use those challenges as the basis for identifying the next set of targets for DOE’s transmission-related research and development (R&D) programs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Generalizing mkFit and its Application to HL-LHC

mkFit is an implementation of the Kalman filter-based track reconstruction algorithm that exploits both thread- and data-level parallelism. In the past few years the project transitioned from the R&D phase to deployment in the Run-3 offline workflow of the CMS experiment. The CMS tracking performs a series of iterations, targeting reconstruction of tracks of increasing difficulty after removing hits associated to tracks found in previous iterations. mkFit has been adopted for several of the tracking iterations, which contribute to the majority of reconstructed tracks. When tested in the standard conditions for production jobs, speedups in track pattern recognition are on average of the order of 3.5x for the iterations where it is used (3-7x depending on the iteration). Multiple factors contribute to the observed speedups, including vectorization and a lightweight geometry description, as well as improved memory management and single precision. Efficient vectorization is achieved with both the icc and the gcc (default in CMSSW) compilers and relies on a dedicated library for small matrix operations, Matriplex, which has recently been released in a public repository. While the mkFit geometry description already featured levels of abstraction from the actual Phase-1 CMS tracker, several components of the implementations were still tied to that specific geometry. We have further generalized the geometry description and the configuration of the run-time parameters, in order to enable support for the Phase-2 upgraded tracker geometry for the HL-LHC and potentially other detector configurations. The implementation strategy and high-level code changes required for the HL-LHC geometry are presented. Speedups in track building from mkFit imply that track fitting becomes a comparably time consuming step of the tracking chain.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Three-step colloidal gelation revealed by time-resolved x-ray photon correlation spectroscopy

The gelation of PEGylated gold nanoparticles dispersed in a glycerol–water mixture is probed in situ by x-ray photon correlation spectroscopy. Following the evolution of structure and dynamics over 10 4 s, a three-step gelation process is found. First, a simultaneous increase of the Ornstein–Zernike length ξ and slowdown of dynamics is characterized by an anomalous q-dependence of the relaxation times of τ ∝q -6 and strongly stretched intermediate scattering functions. After the structure of the gel network has been established, evidenced by a constant ξ, the dynamics show aging during the second gelation step accompanied by a change toward ballistic dynamics with τ ∝ q -1 and compressed correlation functions. In the third step, aging continues after the arrest of particle motion. Our observations further suggest that gelation is characterized by stress release as evidenced by anisotropic dynamics once gelation sets in.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗