Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “block structure”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗
DOE OSTI2025

A conserved viral RNA fold enables nuclease resistance across kingdoms of life

Abstract Viral exoribonuclease-resistant RNA (xrRNA) structures block cellular nucleases to produce subgenomic viral RNAs during infection. High sequence variability among xrRNAs from distantly related viruses raises questions about the shared molecular features that enable these RNAs to withstand the strong unwinding forces of exoribonucleases. Here, we present the first structure of a plant-virus xrRNA in its active conformation and uncover universal principles of xrRNA folding. Comparison with the structure of a human-pathogenic flavivirus xrRNA reveals that both share a core structural motif—a protective ring encircling the RNA’s 5′ end—despite lacking sequence similarity. Disrupting this core motif through targeted mutagenesis eliminates exoribonuclease-resistance and attenuates viral infection. We identify hundreds of related structures across multiple virus families, supporting the conservation of this mechanism. Our study demonstrates how distantly related RNA viruses have converged on a common structural strategy to inhibit cellular nucleases, with a universal ring topology as the defining feature of viral xrRNAs.

Biochemistry & Molecular Biology↗
DOE OSTI2022

Renewable low viscosity polyester-polyols for biodegradable thermoplastic polyurethanes

In the transition to renewably sourced, biodegradable polymers, the preparation of low viscosity polyester-polyols has posed a challenge for renewable polyurethane (PU) development. Low viscosity polyols not only reduce the requirement for high process temperatures but also decrease manufacturing time. In our efforts to incorporate increasing ratios of bio-based monomers into renewable PUs, we mixed diacids such as even carbon sebacic acid and odd carbon azelaic acid along with a renewable diol. This provided library of 2000 g/mol molecular weight polyester-polyols, and structures were established by 1 H and 13 C NMR analysis. The prepared polyester-polyols offered lower viscosity and enable lower fabrication temperatures to make TPUs, and their structure and material metrics were evaluated. The formation of TPUs is ascertained from FTIR and NMR analysis. The final TPUs displayed good physical and mechanical properties. Further, these TPUs exhibited T g in the range of -56.5 to -39.7°C, corresponding to TPU soft block structure, and T m between 98.3 and 105.1°C originating from the hard segment. Prepared TPUs exhibit excellent biodegradation under compost environmental conditions. These TPUs showed up to 57% decrease in molecular weight by GPC analysis after 9 weeks of biodegradation, and respirometer analysis displayed up to 97% biodegradation over 120 days.

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

Metallicity, Atomic Disorder, and Li-Ion Storage in Fast-Charging Anodes

Oxides of Nb with Wadsley-Roth shear structures comprise a family of stable, high-rate anode materials for Li-ion batteries. A particular pair of them offers the unusual opportunity to test how important metallic conduction of the starting electrode is for electrode performance. The selected pair of compounds with similar 4 × 3 Wadsley-Roth block structures are insulating Ti 2 Nb 10 O 29 and metallic Nb 12 O 29 . A combination of diffraction, electrochemistry, magnetic measurements, and entropic potential measurements is employed to establish key findings for these two anode materials. We find that starting with a metallic oxide is not especially advantageous over a comparable material that readily transitions into a metallic state upon lithiation. Second, the rate performance appears to be dictated by ion mobility, and atomic Ti/Nb disorder in Ti 2 Nb 10 O 29 contributes to improved capacity retention at high rates by suppressing Li-ion ordering. However, subtle details in the nature of redox processes make Nb 12 O 29 a slightly better electrode material for long-term cycling at slower rates.

crystal structure↗
DOE OSTI2021

Island scanning pattern optimization for residual deformation mitigation in laser powder bed fusion via sequential inherent strain method and sensitivity analysis

Laser powder bed fusion (L-PBF) has emerged as one of the mainstream additive manufacturing approaches for fabricating metal parts with complex geometries and intricate internal structures. However, large deformation associated with rapid heating and cooling can lead to build failure and requires post-processing which may increase manufacturing cost and prolong the production period. Here in this work, an island scanning pattern design method is proposed to optimize the scanning direction of each island in order to reduce part deformation after cutting off the build platform. The objective of this optimization is to minimize the upward bending of the part after sectioning, which allows the part deformation to satisfy the tolerance requirement or reduce the post heat treatment time. Inherent strain method is employed in the sequential finite element analysis consisting of layer-by-layer activations and sectioning for fast residual distortion prediction. Full sequential sensitivity analysis for the formulated optimization is provided to update the island scanning directions. To show the feasibility and effectiveness of the proposed method, the scanning patterns of a block structure and a connecting rod were designed and parts were fabricated using an open architecture L-PBF machine. The fabrication experiments demonstrated that the residual deformation of both parts fabricated by optimized scanning pattern can be reduced by over 50% compared to the initial scanning patterns, which demonstrate the effectiveness of the proposed method.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Early experiences on the OLCF Frontier system with AthenaPK and Parthenon–Hydro

The Oak Ridge Leadership Computing Facility (OLCF) has been preparing the nation's first exascale system, Frontier, for production and end users. Frontier is based on HPE Cray's new EX architecture and Slingshot interconnect and features 74 cabinets of optimized 3rd Gen AMD EPYC CPUs for HPC and AI and AMD Instinct 250X accelerators. As a part of this preparation, “real-world” user codes have been selected to help assess the functionality, performance, and usability of the system. This article describes early experiences using the system in collaboration with the Hamburg Observatory for two selected codes, which have since been adopted in the OLCF test harness. Experiences discussed include efforts to resolve performance variability and per-cycle slowdowns. Results are shown for a performance portable astrophysical magnetohydronamics code, AthenaPK, and a mini-application stressing the core functionality of a performance portable block-structured adaptive mesh refinement framework, Parthenon-Hydro. Here, these results show good scaling characteristics to the full system. At the largest scale, the Parthenon-Hydro miniapp reaches a total of $1.7$ $\times$ $10^{13}$ zone-cycles/s on 9216 nodes (73,728 logical GPUs) at ≈92% weak scaling parallel efficiency (starting from a single node using a second-order, finite-volume method).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗
DOE OSTI2025

AMR-Wind: A Performance-Portable, High-Fidelity Flow Solver for Wind Farm Simulations

We present AMR-Wind, a verified and validated high-fidelity computational-fluid-dynamics code for wind farm flows. AMR-Wind is a block-structured, adaptive-mesh, incompressible-flow solver that enables predictive simulations of the atmospheric boundary layer and wind plants. It is a highly scalable code designed for parallel high-performance computing with a specific focus on performance portability for current and future computing architectures, including graphical processing units (GPUs). In this paper, we detail the governing equations, the numerical methods, and the turbine models. Establishing a foundation for the correctness of the code, we present the results of formal verification and validation. The verification studies, which include a novel actuator line test case, indicate that AMR-Wind is spatially and temporally second-order accurate. The validation studies demonstrate that the key physics capabilities implemented in the code, including actuator disk models, actuator line models, turbulence models, and large eddy simulation (LES) models for atmospheric boundary layers, perform well in comparison to reference data from established computational tools and theory. We conclude with a demonstration simulation of a 12-turbine wind farm operating in a turbulent atmospheric boundary layer, detailing computational performance and realistic wake interactions.

17 WIND ENERGY↗
DOE OSTI2023

An implicit-in-time DPG formulation of the 1D1V Vlasov-Poisson equations

Efficient solution of the Vlasov equation, which can be up to six-dimensional, is key to the simulation of many difficult problems in plasma physics. The discontinuous Petrov-Galerkin (DPG) finite element methodology provides a framework for the development of stable (in the sense of Ladyzhenskaya–Babuška–Brezzi conditions) finite element formulations, with built-in mechanisms for adaptivity. While DPG has been studied extensively in the context of steady-state problems and to a lesser extent with space-time discretizations of transient problems, relatively little attention has been paid to time-marching approaches. In the present work, we study a first application of time-marching DPG to the Vlasov equation, using backward Euler for a Vlasov-Poisson discretization. We demonstrate adaptive mesh refinement for two problems: the two-stream instability problem, and a cold diode problem. Furthermore, we believe the present work is novel both in its application of unstructured adaptive mesh refinement (as opposed to block-structured adaptivity, which has been studied previously) in the context of Vlasov-Poisson, as well as in its application of DPG to the Vlasov-Poisson system. We also discuss extensive additions to the Camellia library in support of both the present formulation as well as extensions to higher dimensions, Maxwell equations, and space-time formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗
DOE OSTI2020

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A hybrid adaptive multiresolution approach for the efficient simulation of reactive flows

Computational studies that use block-structured adaptive mesh refinement (AMR) approaches suffer from unnecessarily high mesh resolution in regions adjacent to important solution features. This deficiency limits the performance of AMR codes. In this work a novel hybrid adaptive multiresolution (HAMR) approach to AMR-based calculations is introduced to address this issue. The multiresolution (MR) smoothness indicators are used to identify regions of smoothness on the mesh where the computational cost of individual physics solvers may be decreased by replacing direct calculations with interpolation. We suggest an approach to balance the errors due to the adaptive discretization and the interpolation of physics quantities such that the overall accuracy of the HAMR solution is consistent with that of the MR-driven AMR solution. The performance of the HAMR scheme is evaluated for a range of test problems, from pure hydrodynamics to turbulent combustion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗
DOE OSTI2022

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗
DOE OSTI2021

Neutron diffraction discriminates between models for the nanoarchitecture of graphene sheets in glassy carbon

Glassy carbon is a chemically inert, biocompatible, disordered material with the graphene sheet as its basic building block. Structural characterisation techniques have so far been unable to provide definitive distinction between the many proposed models for its structure. Computer based simulation methods have made a step forward by predicting structure from interaction potentials, but the results are sensitive to the choice of the potential. Here we use the white neutron beam of the Spallation Neutron Source at Oak Ridge National Laboratory coupled with the SNAP time-of-flight diffractometer to sample the reciprocal space of glassy carbon and calculate accurate radial distribution functions. From the radial distribution information, we determine graphene sheet dimensions, registration, curvature, and assess defect content to rank the proposed structures in terms of their agreement with experimental data. In this work, we find that the recent models derived from empirical potentials give the best agreement.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Theoretical Treatment of the Coriolis Effect Using Hyperspherical Coordinates, with Application to the Ro-Vibrational Spectrum of Ozone

Here, several alternative methods for the description of the interaction between rotation and vibration are compared and contrasted using hyperspherical coordinates for a triatomic molecule. These methods differ by the choice of the z -axis and by the assumption of a prolate or oblate rotor shape of the molecule. For each case, a block-structure of the rotational–vibrational Hamiltonian matrix is derived and analyzed, and the advantages and disadvantages of each method are made explicit. This theory is then employed to compute ro-vibrational spectra of singly substituted ozone; roughly, 600 vibrational states of 16 O 18 O 16 O and 16 O 16 O 18 O isomers combined, with rotational excitations up to J = 5 and both inversion parities (21600 coupled ro-vibrational states in total). Splittings between the states of different parities, so-called K -doublings, are calculated and analyzed. The roles of the asymmetric-top rotor term and the Coriolis coupling term are determined individually, and it is found that they both affect these splittings, but in the opposite directions. Thus, the two effects partially cancel out, and the residual splittings are relatively small. Energies of the ro-vibrational states reported in this work for 16 O 18 O 16 O and 16 O 16 O 18 O are in excellent agreement with literature (available for low-vibrational excitation). New data obtained here for the highly excited vibrational states enable the first systematic study of the Coriolis effect in symmetric and asymmetric isotopomers of ozone.

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

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Data-driven Distributed Learning of Multi-agent Systems: A Koopman Operator Approach

Koopman operator theory provides a model-free technique for studying nonlinear dynamical systems purely from data. Since the Koopman operator is infinite-dimensional, researchers have developed several methods that provide a finite-dimensional approximation of the Koopman operator so that it can be applied for practical use cases. One common thing with most of the methods is that their solutions are obtained by solving a centralized minimization problem. In this work, we treat the dynamical system to be a multi-agent system and propose an algorithm to compute the finite-dimensional approximation of the Koopman operator in a distributed manner using the knowledge of the topology of the underlying multi-agent system. The proposed distributed approach is shown to be equivalent to the centralized learning problem and results in a sparse Koopman whose block structure mimics the Laplacian of the multi-agent system. Extensive simulation studies illustrate the proposed framework on the network of oscillators and the IEEE 68 bus system.

Nandanoori, Sai Pushpak↗