Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “topological gradient”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluating Material Design Principles for Calcium-Ion Mobility in Intercalation Cathodes

Multivalent-ion batteries offer an alternative to Li-based technologies, with the potential for greater sustainability, improved safety, and higher energy density, primarily due to their rechargeable system featuring a passivating metal anode. Although a system based on the Ca 2+ /Ca couple is particularly attractive given the low electrochemical plating potential of Ca 2+ , the remaining challenge for a viable rechargeable Ca battery is to identify Ca cathodes with fast ion transport. In this work, a high-throughput computational pipeline is adapted to (1) discover novel Ca cathodes in a largely unexplored space of empty intercalation hosts and (2) develop material design rules for Ca-ion mobility. One candidate from the screening, W 2 O 3 (PO 4 ) 2 , is confirmed to have a low Nudged Elastic Band (NEB) barrier of 168 meV within a one-dimensional (1D) ion percolation topology. This candidate is subsequently synthesized and electrochemically tested, achieving reversible Ca cycling with a capacity of 25 mA h/g. To further accelerate the screening for promising Ca intercalation electrodes, machine learning (ML) Random Forest (RF) and Extreme Gradient Boosting (XGB) classification models are created with local environment descriptors based on a large, structurally and chemically diverse dataset of minimum energy pathways, spanning over 5,000 density functional theory (DFT) site energy calculations. Accuracies of 92% are achieved, material design metrics are quantified, ML force-fields are leveraged in an accelerated iteration of the screening, and a total of 27 novel Ca cathode materials are highlighted for further investigation.

25 ENERGY STORAGE↗
DOE OSTI2025

Double-Crucible Vertical Bridgman Technique for Stoichiometry-Controlled Chalcogenide Crystal Growth

Precise stoichiometry control in single-crystal growth is essential for both technological applications and fundamental research. However, conventional growth methods often face challenges such as non-stoichiometry, compositional gradients, and phase impurities, particularly in non-congruent melting systems. Even in congruent melting systems like Bi₂Se₃, deviations from the ideal stoichiometric composition can lead to significant property degradation, such as excessive bulk conductivity, which limits its topological applications. In this study, we introduce the double-crucible vertical Bridgman (DCVB) method, a novel approach that enhances stoichiometry control through the combined use of continuous source material feeding, traveling-solvent growth, and liquid encapsulation, which suppresses volatile element loss under high pressure. Using Bi₂Se₃ as a model system, we demonstrate that crystals grown via DCVB exhibit enhanced stoichiometric control, significantly reducing defect density and achieving much lower carrier concentrations compared to those produced by conventional Bridgman techniques. Moreover, the continuous feeding of source material enables the growth of large crystals. As a result, this approach presents a promising strategy for synthesizing high-quality, large-scale crystals, particularly for metal chalcogenides and pnictides that exhibit challenging non-congruent melting behaviors.

Crystal structure↗
DOE OSTI2023

Multi-objective, robust constraints enforced global topology optimizer for optical devices

A method for optimization of photonic devices is disclosed. The method includes receiving a set of unconstrained latent variables; mapping the set of unconstrained latent variables to a constrained space to generate a constrained device; calculating the permittivity across each element of the constrained device; determining a permittivity-constrained width gradient based at least partially on the permittivity across each element; and optimizing the set of unconstrained latent variables by at least partially using the permittivity-constrained width gradient.

Fan, Jonathan↗
DOE OSTI2021

Magnetic field-induced non-trivial electronic topology in Fe3− x GeTe2

The anomalous Hall, Nernst, and thermal Hall coefficients of the itinerant ferromagnet Fe3−xGeTe2 display anomalies upon cooling that are consistent with a topological transition that could induce deviations with respect to the Wiedemann–Franz (WF) law. This law has not yet been validated for the anomalous transport variables, with recent experimental studies yielding material-dependent results. Nevertheless, the anomalous Hall and thermal Hall coefficients of Fe3−xGeTe2 are found, within our experimental accuracy, to satisfy the WF law for magnetic fields μ0H applied along its c axis. Remarkably, large anomalous transport is also observed for μ0H||a axis with the field aligned along the gradient of the chemical potential generated by thermal gradients or electrical currents, a configuration that should not lead to their observation. These anomalous planar quantities are found to not scale with the component of the planar magnetization (M||), showing instead a sharp decrease beyond μ0H||= 4 T or the field required to align the magnetic moments along μ0H||. We argue that chiral spin structures associated with Bloch domain walls lead to a field-dependent spin chirality that produces a novel type of topological transport in the absence of interaction between the magnetic field and electrical or thermal currents. Locally chiral spin structures are captured by our Monte Carlo simulations incorporating small Dzyaloshinskii–Moriya and biquadratic exchange interactions. These observations reveal not only a new way to detect and expose topological excitations, but also a new configuration for heat conversion that expands the current technological horizon for thermoelectric energy applications.

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

Mapping structures and dynamics with frequency-correlated diffusion exchange

Understanding molecular motion in diffusion-driven complex environments is critical for designing sustainable materials and improving chemical processes. Here, we introduce a multidimensional nuclear magnetic resonance (NMR) method that captures how molecular populations exchange across different dynamic regimes. By extending the modulated gradient spin-echo technique to include frequency-frequency correlations, our approach reveals diffusion pathways that are otherwise obscured in heterogeneous systems. Implemented on a unilateral NMR magnet, the method eliminates gradient pulsing constraints and accesses dynamics in the kilohertz regime. We apply this technique to swelling and acid-catalyzed deconstruction of cross-linked and linear polymers to observe how structural heterogeneity evolves over time. By linking molecular motion to topology and chemical state, we extract physical metrics such as fractal surface dimensionality and reaction wavefront velocity, properties inaccessible with standard diffusion measurements. This work expands the capabilities of NMR for probing soft matter, with implications for polymer recycling and materials design.

Fricke, Sophia N. [University of California, Berke↗
DOE OSTI2022

Vapor transport growth of MnBi 2 Te 4 and related compounds

Motivated by fine tuning of the magnetic and topological properties of MnBi 2 Te 4 via defect engineering, in this work, we report the crystal growth of MnBi 2 Te 4 and related compounds using vapor transport method and crystal characterization by measuring elemental ratio, magnetic and transport properties, and scanning tunneling microscopy. For the growth of MnBi 2 Te 4 single crystals, I 2 , MnI 2 , MnCl 2 , TeCl 4 , and MoCl 5 are all effective transport agents; chemical transportation occurs faster in the presence of iodides than chlorides. We further successfully grow MnSb 2 Te 4 , MnBi 2-x Sb x Te 4 , and Sb-doped MnBi 4 Te 7 crystals. A small temperature gradient< 20°C between the hot and cold ends of the growth ampoule is critical for the successful crystal growth of MnBi 2 Te 4 and related compounds. Compared to flux grown crystals, vapor transported crystals tend to be Mn stoichiometric, and Sb-bearing compositions have more Mn/Sb site mixing. The vapor transport growth provides a new materials synthesis approach to fine tune the magnetic and topological properties of these intrinsic magnetic topological insulators where controlling defects is vital.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electric polarization in inhomogeneous crystals

In this work, we derive the charge density up to second order in spatial gradient in inhomogeneous crystals using the semiclassical coarse graining procedure based on the wave packet method. It can be recast as divergence of polarization, whose first-order contribution consists of three parts, a perturbative correction to the original Berry connection expression, a topological part that can be written as an integral of the Chern-Simons 3-form, and a previously-unknown, quadrupolelike contribution. The topological part can be related to the quantized fractional charge carried by a vortex in two-dimensional systems. We then generalize our results to the multiband case and show that the quadrupolelike contribution plays an important role, as it makes the total polarization gauge-independent. Finally, we verify our theory in several model systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Localized Evaluation for Constructing Discrete Vector Fields

Topological abstractions offer a method to summarize the behavior of vector fields, but computing them robustly can be challenging due to numerical precision issues. One alternative is to represent the vector field using a discrete approach, which constructs a collection of pairs of simplices in the input mesh that satisfies criteria introduced by Forman's discrete Morse theory. While numerous approaches exist to compute pairs in the restricted case of the gradient of a scalar field, state-of-the-art algorithms for the general case of vector fields require expensive optimization procedures. This paper introduces a fast, novel approach for pairing simplices of two-dimensional, triangulated vector fields that do not vary in time. The key insight of our approach is that we can employ a local evaluation, inspired by the approach used to construct a discrete gradient field, where every simplex in a mesh is considered by no more than one of its vertices. Specifically, we observe that for any edge in the input mesh, we can uniquely assign an outward direction of flow. We can further expand this consistent notion of outward flow at each vertex, which corresponds to the concept of a downhill flow in the case of scalar fields. Working with outward flow enables a linear-time algorithm that processes the (outward) neighborhoods of each vertex one-by-one, similar to the approach used for scalar fields. Here, we couple our approach to constructing discrete vector fields with a method to extract, simplify, and visualize topological features. Empirical results on analytic and simulation data demonstrate drastic improvements in running time, produce features similar to the current state-of-the-art, and show the application of simplification to large, complex flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fracture-based shape optimization built upon the topological derivative

In Silva et al. (2011) and Alidoost et al. (2020), the authors developed an approximation of the energy release rate field associated with a small edge or surface crack at any boundary location and with any orientation using the topological derivative. The approximation is computationally attractive because it requires only a single analysis on the non-cracked domain in contrast with conventional boundary-element and finite-element-based methods, which require a separate and costlier analysis for each crack length-location-orientation combination. Here, a shape optimization scheme for fracture-resistant structures is developed using the energy release rate approximation. In the gradient-based optimization scheme, the domain and its boundary are defined implicitly using level-set functions. The level-set functions of arbitrary geometries are constructed using Boolean operations from the level-set functions of simple primitives. This geometrical representation has the dual advantage of (i) allowing shapes to intersect and/or separate during the optimization and (ii) simplifying the computation of the shape sensitivities.

42 ENGINEERING↗
DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗
DOE OSTI2020

Deep Generative Models that Solve PDEs: Distributed Computing for Training Large Data-Free Models

Recent progress in scientific machine learning (SciML) has opened up the possibility of training novel neural network architectures that solve complex partial differential equations (PDEs). Several (nearly data free) approaches have been recently reported that successfully solve PDEs, with examples including deep feed forward networks, generative networks, and deep encoder-decoder networks. However, practical adoption of these approaches is limited by the difficulty in training these models, especially to make predictions at large output resolutions (≥1024×1024). Here we report on a software framework for data parallel distributed deep learning that resolves the twin challenges of training these large SciML models - training in reasonable time as well as distributing the storage requirements. Our framework provides several out of the box functionality including (a) loss integrity independent of number of processes, (b) synchronized batch normalization, and (c) distributed higher-order optimization methods. We show excellent scalability of this framework on both cloud as well as HPC clusters, and report on the interplay between bandwidth, network topology and bare metal vs cloud. We deploy this approach to train generative models of sizes hitherto not possible, showing that neural PDE solvers can be viably trained for practical applications. We also demonstrate that distributed higher-order optimization methods are 2-3× faster than stochastic gradient-based methods and provide minimal convergence drift with higher batch-size.

PDEs↗
DOE OSTI2023

Protection Against Graph-Based False Data Injection Attacks on Power Systems

Graph signal processing (GSP) has emerged as a powerful tool for practical network applications, including power system monitoring. By representing power system voltages as smooth graph signals, recent research has focused on developing GSP-based methods for state estimation, attack detection, and topology identification. Included, efficient methods have been developed for detecting false data injection (FDI) attacks, which until now were perceived as non-smooth with respect to the graph Laplacian matrix. Consequently, these methods may not be effective against smooth FDI attacks. In this paper, we propose a graph FDI (GFDI) attack that minimizes the Laplacian-based graph total variation (TV) under practical constraints. In addition, we develop a low-complexity algorithm that solves the non-convex GDFI attack optimization problem using ell_1-norm relaxation, the projected gradient descent (PGD) algorithm, and the alternating direction method of multipliers (ADMM). We then propose a protection scheme that identifies the minimal set of measurements necessary to constrain the GFDI output to high graph TV, thereby enabling its detection by existing GSP-based detectors. Our numerical simulations on the IEEE-57 bus test case reveal the potential threat posed by well-designed GSP-based FDI attacks. Moreover, we demonstrate that integrating the proposed protection design with GSP-based detection can lead to significant hardware cost savings compared to previous designs of protection methods against FDI attacks.

Morgenstern, Gal↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2024

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗
DOE OSTI2024

Electron‐Assisted Generation and Straight Movement of Skyrmion Bubble in Kagome TbMn 6 Sn 6

Abstract Topological magnetic textures are promising candidates as binary data units for the next‐generation memory device. The precise generation and convenient control of nontrivial spin topology at zero field near room temperature endows the critical advantages in skyrmionic devices but is not simultaneously integrated into one material. Here, in the Kagome plane of quantum TbMn 6 Sn 6 , the expedient generation of the skyrmion bubbles in versatile forms of lattice, chain, and isolated one by converging the electron beam, where the electron intensity gradient contributes to the dynamic generation from local anisotropy variation near spin reorientation transition (SRT) is reported. Encouragingly, by utilizing the dynamic shift of the SRT domain interface, the straight movement is actualized with the skyrmion bubble slave to the SRT domain interface forming an elastic composite object, avoiding the usual deflection from the skyrmion Hall effect. The critical contribution of the SRT domain interface via conveniently electron‐assisted heating is further theoretically validated in micromagnetic simulation, highlighting the compatible application possibility in advanced devices.

Li, Zhuolin↗
DOE OSTI2021

Multiscale Chirping Modes Driven by Thermal Ions in a Plasma with Reactor-Relevant Ion Temperature

A thermal ion driven bursting instability with rapid frequency chirping, considered as an Alfv´enic ion temperature gradient mode, has been observed in plasmas having reactor-relevant temperature in the DIII-D tokamak. The modes are excited over a wide spatial range from macroscopic device size to micro-turbulence size and the perturbation energy propagates across multiple spatial scales. The radial mode structure is able to expand from local to global in ~ 0.1 ms and it causes magnetic topology changes in the plasma edge, which can lead to a minor disruption event. Since the mode is typically observed in the high ion temperature ≳ 10 keV and high-! plasma regime, the manifestation of the mode in future reactors should be studied with development of mitigation strategies, if needed. In conclusion, this is the first observation of destabilization of the Alfv´en continuum caused by the compressibility of ions with reactor-relevant ion temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Measurement of phonon angular momentum

In condensed matter, angular momentum is intimately related to the emergence of topological quantum states, including chiral superconductivity, quantum spin liquids and various chiral quasiparticles. Recently, it has been predicted that microscopic lattice excitations, known as phonons, can carry finite angular momentum, leading to specific physical properties of materials. However, phonon angular momentum has not yet been observed directly. Here, in this study, we demonstrate that angular momentum conservation results in a macroscopic mechanical torque when applying a time-reversal symmetry-breaking thermal gradient along the chiral axis of single-crystal tellurium. We probe this torque using a cantilever-based device and establish that it changes sign by flipping the chirality or thermal gradient. This behavior disappears in polycrystalline samples that lack a preferred chirality. Our experimental results align well with theoretical calculations. We provide compelling evidence for phonon angular momentum, which might enable quantum states with potential applications in microelectronics.

36 MATERIALS SCIENCE↗