Engineering Papers⌕ Search

Engineering topics

Liu, Yang

Publications and source records attributed to Liu, Yang.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Linear complexity

We present factorization and solution phases for a new linear complexity direct solver designed for concurrent batch operations on fine-grained parallel architectures, for matrices amenable to hierarchical representation. We focus on the strong-admissibility-based $\mathscr{H}^{2}$ format, where strong recursive skeletonization factorization compresses remote interactions. We build upon previous implementations of $\mathscr{H}^{2}$ matrix construction for efficient factorization and solution algorithm design, which are illustrated graphically in stepwise detail. The algorithms are ‘blackbox’ in the sense that the only inputs are the matrix and right-hand side, without analytical or geometrical information about the origin of the system. We demonstrate linear complexity scaling in both time and memory on four representative families of dense matrices up to one million in size. Parallel scaling up to 16 threads is enabled by a multi-level matrix graph coloring and avoidance of dynamic memory allocations thanks to prefix-sum memory management. An experimental backward error analysis is included. We break down the timings of different phases, identify phases that are memory-bandwidth limited, and discuss alternatives for phases that may be sensitive to the trend to employ lower precisions for performance.

Boukaram, Wajih↗
DOE OSTI2026

Complete Phase Transformation of Ir Nanowire Network into Defect-Rich Oxide Catalyst for High-Performance PEM Water Electrolysis

Iridium-based catalysts remain the most reliable option for the oxygen evolution reaction (OER) in proton exchange membrane water electrolyzers (PEMWEs). However, their high cost and limited performance represent critical barriers to the commercialization of this green hydrogen production technology. Herein, we report the creation of a metallic Ir nanowire network (IrNWN), which exhibits superior OER performance through its in-situ transition into an oxide structure with high intrinsic activity. At a low loading of 0.25 mgIr/cm2 in PEMWEs, IrNWN achieved a current density of 3.13 A/cm2 at a cell voltage of 1.8 V, outperforming the commercial Ir-based catalyst and surpassing the Department of Energy (DOE) 2026 technical target. Moreover, the high activity of IrNWN was maintained for 900 hours in a durability test at 2 A/cm2, showing a low degradation rate of 0.042 mV/hour. Structural analysis of the electrochemically oxidized IrNWN revealed the presence of mixed Ir oxidation states and a high density of surface terminal oxygen groups (μ1-O), which contributed to a reduced energy barrier for the rate-determining O-O coupling step.

oxygen evolution reaction↗
DOE OSTI2025

Parallelizing autotuning for HPC applications: Unveiling the potential of the speculation strategy in Bayesian optimization

In the exascale computing era, tuning High-Performance Computing (HPC) applications has become a significant computational challenge. Although Bayesian optimization (BO) has emerged as a promising tool for HPC performance tuning, the BO workflow is inherently sequential (i.e., one function evaluation at a time) and cannot leverage the huge amount of parallel resources present in modern supercomputers, resulting in a considerable underutilization of their computational capabilities. This paper explores the trade-off between search quality and parallelism in BO, investigating a diverse set of methods. Building upon both previous approaches from the literature and novel methodologies introduced in this work, our study provides a deep analysis to accelerate BO performance tuning. By examining a set of synthetic functions and practical HPC applications, our exploration analyzes the interaction among various BO methods for parallelization, the quantity of parallel resources, the runtime distribution of target HPC applications, and the costs associated with different search orchestration mechanisms that have been overlooked in previous studies. Compared to sequential BO, our novel methodology achieves comparable quality while demonstrating robust scalability in search time as the amount of parallel resources increases; it also outperforms a state-of-the-art tuner, which supports parallelization, achieving up to 3.67x faster search time. We provide high-value insights for practitioners seeking to leverage the power of parallel computing for efficient HPC application tuning. Additionally, to further assist researchers in accelerating the performance tuning of their HPC applications, we provide an extension of an existing open-source tuning framework that incorporates our methods.

Bayesian optimization↗
DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗
DOE OSTI2025

Natural nanoparticle complexes at water-water interfaces

Aqueous two-phase systems (ATPSs) with multifunctional attributes have significant promise as biomimetic materials, but current approaches do not harness their full potential. Here, we show that ATPSs can be stabilized to form microcapsules by the interfacial assembly of chitin nanofibers (ChNF) and rod-like cellulose nanocrystals (CNC). The high structural stability of the ChNF/CNC complexes integrates permeability and transport across the membrane at the liquid/liquid interface. Driven by density and osmotic stress gradients, the microcapsules show switchable motility, including cyclic meniscus-climbing and subsurface transport. These observations demonstrate a self-regulating system with potential for cargo transfer, cell biomimicry, as well as ATPSs-based microreactors and microrobots.

Bioinspired materials↗
DOE OSTI2024

Batched sparse direct solver design and evaluation in SuperLU_DIST

Over the course of interactions with various application teams, the need for batched sparse linear algebra functions has emerged in order to make more efficient use of the GPUs for many small and sparse linear algebra problems. In this paper, we present our recent work on a batched sparse direct solver for GPUs. The sparse LU factorization is computed by the levels of the elimination tree, leveraging the batched dense operations at each level and a new batched Scatter GPU kernel. The sparse triangular solve is computed by the level sets of the directed acyclic graph (DAG) of the triangular matrix. Batched operations overcome the large overhead associated with launching many small kernels. For medium sized matrix batches with not-so-small bandwidth, using an NVIDIA A100 GPU, our new batched sparse direct solver is orders of magnitude faster than a batched banded solver and uses less than one-tenth of the memory.

Boukaram, Wajih↗
DOE OSTI2024

Development of Whole System Digital Twins for Advanced Reactors: Leveraging Graph Neural Networks and SAM Simulations

Here, in this work, we introduce a novel method to develop whole system digital twins (DTs) for advanced nuclear reactors. This method treats a complex reactor system as a heterogeneous graph: with the system components as different types of graph nodes and their physical interconnections as edges. Based on the heterogeneous graph, a graph neural network combining graph convolution and temporal node attention is developed as the DT, facilitating a comprehensive understanding of the system's dynamic behavior. By utilizing the System Analysis Module (SAM) code for simulating various operational transients, we develop a graph-based database that trains the DT. This DT is characterized by two primary functions: It can infer the entire system's status using sparse node information, and it can predict the progress of transients based on current and historical system information. Our approach is validated through case studies on the Experimental Breeder Reactor II (EBR-II) system and a generic Fluoride-salt-cooled High-temperature Reactor (gFHR), demonstrating the DT's accuracy in forecasting operational transients. The DT's rapid computation capabilities enhance its potential for supporting advanced reactor operations, offering benefits in intelligent simulation, autonomous control, and anomaly detection, paving the way for improved safety analysis and intelligent component health management for advanced reactor systems and reducing their operations and maintenance cost.

EBR-II↗
DOE OSTI2024

A model for oblique impacts on material surfaces

Many practical situations of material damage, wear, and erosion involve collisions between small particles and surfaces at inclined angles. While there are many well-validated models of normal incidence impact situations, elastic-plastic models for oblique incidence impact events are lacking. Here the finite element method is used to predict the normal and tangential coefficient of restitution in oblique impacts for hard, elastic spheres impacting an elastic-perfectly plastic material surface. The proposed model covers various impact angles ranging from 0° to 45°, within a limiting impact velocity below which the effects of heating are negligible. The normal coefficient of restitution follows power-laws with respect to normalized values of the impact velocity. Interestingly, the tangential coefficient of restitution follows a linear relationship with impact velocity. Together, these results provide a semi-empirical set of equations predicting oblique impact rebounds (both velocity and trajectory) for a wide range of conditions and material properties, with which experimental results can be rapidly interpreted. As a result, Laser-Induced Particle Impact Test (LIPIT) data are also presented for aluminum particles impacting aluminum substrates, at impact angles of 25° and 40°; the results compare favorably with the model and validate the general use of such models for the analysis of experimental data.

Coefficient of restitution↗
DOE OSTI2024

Effects of irradiation damage on the hardness and elastic properties of quaternary and high entropy transition metal diborides

Multi-principal component transition metal (TM) diborides represent a class of high-entropy ceramics (HECs) that have received considerable interest in recent years owing to their promising properties for extreme environment applications that include thermal/ environmental barriers, hypersonic vehicles, turbine engines, and next-generation nuclear reactors. While the addition of chemical disorder through the random distribution of TM elements on the cation sublattice has offered opportunities to tailor elastic stiffness and hardness, the effects of irradiation-induced structural damage on the physical properties of these complex materials have remained largely unexplored. To this end, changes in the hardness and elastic moduli of a high-entropy TM diboride (Hf 0.2 Nb 0.2 Ta 0.2 Ti 0.2 Zr 0.2 )B 2 and three of its quaternary subsets following irradiation with 10 MeV gold (Au) ions to fluences of up to 6 × 10 15 Au cm -2 are investigated at the micrometer and sub-micrometer length-scales via the dispersion of laser-generated surface acoustic waves (SAW) and nanoindentation, respectively. The nanoindentation measurements show that the TM diborides exhibit an initial increase in hardness following irradiation with energetic Au ions, with a subsequent decrease in hardness following further irradiation. One quaternary composition, (Hf 1/3 Ta 1/3 Ti 1/3 )B 2 , exhibits a notable exception to the trend and continues to exhibit an increase in hardness with ion irradiation fluence. Although differences in the absolute values of the effective elastic moduli obtained from the measured SAW dispersion and nanoindentation are observed (and attributed to microstructural variations at the measurement length-scale), both techniques yield similar trends in the form of an initial reduction and subsequent saturation in the elastic modulus with increasing ion irradiation fluence. The quaternary TM diboride (Hf 1/3 Ta 1/3 Ti 1/3 )B 2 again exhibits a departure from this trend. The high-entropy TM diboride (Hf 0.2 Nb 0.2 Ta 0.2 Ti 0.2 Zr 0.2 )B 2 exhibits the greatest recovery in hardness and modulus when irradiated to high ion fluences following initial changes at low fluence, indicating superior resistance to radiation-induced damage over its quaternary counterparts. Opportunities for designing HECs with superior hardness and modulus for enhanced radiation resistance (compared to their single constituent counterparts) by tailoring chemical disorder and bond character in the lattice are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Non-smooth Bayesian optimization in tuning scientific applications

Tuning algorithmic parameters to optimize the performance of large, complicated computational codes is an important problem involving finding the optima and identifying regimes defined by non-smooth boundaries in black-box functions. Within the Bayesian optimization framework, the Gaussian process surrogate model produces smooth mean functions, but functions in the tuning problem are often non-smooth, which is exacerbated by the fact that we usually have limited sequential samples from the black-box function. Here, motivated by these issues encountered in tuning, we propose a novel Gaussian process model called a clustered Gaussian process (cGP), where the components are dynamically updated by clustering. In our studies, the performance of cGP can be better than stationary GPs in nearly 90% of the experiments and better than non-stationary GPs in nearly 70% of the repeated experiments while requiring less computational cost. cGP provides a novel approach for dynamic GP, computes more efficiently than recursive partitioning, and discovers non-smoothness regimes. We provide extensive experiments including high-performance computing (HPC) and industrial simulation functions to show the effectiveness of our methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Semi-Analytical Approach for State-Space Electromagnetic Transient Simulation

Here, this paper proposes a semi-analytical approach for efficient and accurate electromagnetic transient (EMT) simulation of a power grid. The approach first derives a high-order semi-analytical solution (SAS) of the grid’s state-space EMT model using the differential transformation (DT), and then evaluates the solution over enlarged, variable time steps to significantly accelerate the simulations while maintaining its high accuracy on detailed fast EMT dynamics. The approach also addresses switches during large time steps by using a limit violation detection algorithm with a binary search-enhanced quadratic interpolation. Case studies are conducted on EMT models of the IEEE 39-bus system and large-scale systems to demonstrate the merits of the new simulation approach against traditional numerical methods.

electromagnetic transient↗
DOE OSTI2024

A small secreted protein serves as a plant-derived effector mediating symbiosis between Populus and Laccaria bicolor

Beneficial symbiotic fungi colonize plant tissues, delivering crucial ecosystem services such as carbon sequestration and plant fertilization. Specifically, trees that form a nutrient-acquiring symbiosis with mutualistic ectomycorrhizal (ECM) fungi gain advantages from these associations by experiencing enhanced growth rates and increased resilience to both biotic and abiotic stresses. Despite the vital role ECM fungi play in the nutrition and well-being of trees, identifying key regulators participating in the molecular communication between plant and fungal cells is still in its early stages. The mutualistic relationship between Laccaria bicolor and Populus spp. has been utilized as a model system for investigating ECM symbiosis at the molecular level. It has been demonstrated that the fungus L. bicolor secretes Mycorrhiza-induced Small Secreted Proteins (MiSSPs) required for ECM development. Meanwhile, we have previously shown that P. trichocarpa small, secreted proteins (PtSSPs) are highly induced during mutualistic symbiosis and some of them can enter, via in-vitro feeding, L. bicolor hyphae affecting their growth and morphology. However, the exact role and mode of action of PtSSPs in mutualistic symbiosis remain unknown. Because previous study showed that PtSSP1 is taking up by fungal cell and then localize in fungal cells, we decide to dig further on its putative role in fungal cells and also by overexpressing it in poplar because it is not technically possible yet to overexpress poplar protein in Laccaria bicolor hyphae. Here, we further characterized the function of PtSSP1(Potri.009G063200) in ectomycorrhization, which accumulates in the nucleus of L. bicolor in an in-vitro feeding experiment. Our results provide new knowledge for the genetic engineering of plants to control associated microbes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

On the Emergence of Ferromagnetism in LaCoO 3 Ultrathin Films

It is well known that the properties of a crystal evolve as it increases in size from a single atomic plane to that of the bulk. Such size-dependent transitions can stem from many different origins and depend on minute changes to crystal bonding and composition. Here, a model example is that of LaCoO 3 , which is non-magnetic in the bulk but can display ferromagnetism at the nanoscale. Here, the evolution of structure-property relationships is studied in the LaCoO 3-δ /SrTiO 3 (001) system as the thickness of LaCoO 3-δ is increased from a single plane to 10 unit cells. In situ synchrotron X-ray studies are performed during and post-deposition to probe changes in the interactions between structure, stoichiometry, and magnetic behavior. Structural quantification indicates that the oxygen octahedral rotation pattern evolves with thickness, due to inherent differences in crystal symmetry between the film and substrate. The change in rotation modifies the required energy barrier for the spin state transition via the Co–O bond length and Co–O–Co bond angle, affecting the appearance of ferromagnetism. Our results highlight the contributions of high spin Co 2+ and/or high spin Co 3+ to respective weak and robust ferromagnetism and the evolution of properties with size in ultrathin LaCoO 3-δ heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomic-scale identification of active sites of oxygen reduction nanocatalysts

Heterogeneous nanocatalysts play a crucial role in both the chemical and energy industries. Despite substantial advancements in theoretical, computational and experimental studies, identifying their active sites remains a major challenge. Here we utilize atomic electron tomography to determine the three-dimensional atomic structure of PtNi and Mo-doped PtNi nanocatalysts for the electrochemical oxygen reduction reaction. We then employ the experimental atomic structures as input to first-principles-trained machine learning to identify the active sites of the nanocatalysts. Through the analysis of the structure–activity relationships, we formulate an equation termed the local environment descriptor, which balances the strain and ligand effects to provide physical and chemical insights into active sites in the oxygen reduction reaction. The ability to determine the three-dimensional atomic structure and chemical composition of realistic nanoparticles, combined with machine learning, could transform our fundamental understanding of the active sites of catalysts and guide the rational design of optimal nanocatalysts.

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

DNA parts and gene constructs for plant biodesign

Plant biodesign requires the knowledge of DNA parts (e.g., genes, promoters, terminators), along with their combinations (as gene constructs) linked to engineered traits. DNA parts with validated or predicted functions in plants have been deposited in various online databases. However, these existing databases focus on basic biological functions of individual DNA parts, leaving a gap between basic knowledge and bioengineering applications. To fill this knowledge gap, we have created a user-friendly, open-ended database as a knowledge graph linking DNA parts to gene constructs to traits. This database contains experimentally validated DNA parts and gene constructs documented in peer-reviewed publications. The DNA parts include 1) molecular components with biological functions, such as genes involved in various biological processes (e.g., metabolic and signal transduction pathways) and 2) molecular components with technical functions, such as gene expression, genome engineering and sequence splicing. The gene constructs deposited in this database include both single-gene and multi-gene constructs. This database allows users to submit DNA parts and gene construct compositions linked to engineered traits described in peer-reviewed publications, providing a public digital repository for sharing the biodesign information among the researchers in the fields of plant biotechnology and plant synthetic biology.

plant biodesign synthetic biology gene constructs ↗
DOE OSTI2024

Atomic dynamics of electrified solid–liquid interfaces in liquid-cell TEM

Electrified solid–liquid interfaces (ESLIs) play a key role in various electrochemical processes relevant to energy, biology and geochemistry. The electron and mass transport at the electrified interfaces may result in structural modifications that markedly influence the reaction pathways. For example, electrocatalyst surface restructuring during reactions can substantially affect the catalysis mechanisms and reaction products. Despite its importance, direct probing the atomic dynamics of solid–liquid interfaces under electric biasing is challenging owing to the nature of being buried in liquid electrolytes and the limited spatial resolution of current techniques for in situ imaging through liquids. Here, with our development of advanced polymer electrochemical liquid cells for transmission electron microscopy (TEM), we are able to directly monitor the atomic dynamics of ESLIs during copper (Cu)-catalysed CO 2 electroreduction reactions (CO 2 ERs). Our observation reveals a fluctuating liquid-like amorphous interphase. It undergoes reversible crystalline–amorphous structural transformations and flows along the electrified Cu surface, thus mediating the crystalline Cu surface restructuring and mass loss through the interphase layer. Furthermore, the combination of real-time observation and theoretical calculations unveils an amorphization-mediated restructuring mechanism resulting from charge-activated surface reactions with the electrolyte. Our results open many opportunities to explore the atomic dynamics and its impact in broad systems involving ESLIs by taking advantage of the in situ imaging capability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring Cr and molten salt interfacial interactions for molten salt applications

Molten salts play an important role in various energy-related applications such as high-temperature heat transfer fluids and reaction media. However, the extreme molten salt environment causes the degradation of materials, raising safety and sustainability challenges. A fundamental understanding of material–molten salt interfacial evolution is needed. This work studies the transformation of metallic Cr in molten 50/50 mol% KCl–MgCl 2 via multi-modal in situ synchrotron X-ray nano-tomography, diffraction and spectroscopy combined with density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. Notably, in addition to the dissolution of Cr in the molten salt to form porous structures, a δ-A15 Cr phase was found to gradually form as a result of the metal–salt interaction. This phase change of Cr is associated with a change in the coordination environment of Cr at the interface. DFT and AIMD simulations provide a basis for understanding the enhanced stability of δ-A15 Cr vs. bcc Cr, by revealing their competitive phase thermodynamics at elevated temperatures and probing the interfacial behavior of the molten salt at relevant facets. This study provides critical insights into the morphological and chemical evolution of metal–molten salt interfaces. Finally, the combination of multimodal synchrotron analysis and atomic simulation also offers an opportunity to explore a broader range of systems critical to energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗