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.

159 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Laser‐Powder Melt Pool Solidification Dynamics and Microstructural Engineering of Ti‐5553 Microlattices

The fine geometric and topological control afforded by additive manufacturing technologies has enabled the manufacture of architected materials across length scales, and enabling tunable mechanical performance as a function of local and global design. Progress has been made to tune the mechanical response of architected materials through geometry, but understanding how the geometry and processing conditions will inform the microstructure remains a challenge due to the rapid solidification in laser powder bed fusion. This study uses in situ X-ray imaging and electron backscatter diffraction microscopy to demonstrate that the melt pool size, microstructure morphology, and elastic strain distribution is influenced by a combination of lattice geometry and laser processing conditions. These results indicate that within larger melt pools the local thermal gradients are sufficient to enable a columnar-to-equiaxed transition across the melt pool. Furthermore, the solidification mechanisms producing these microstructures are examined across the first 5 ms of melting and solidification, described via in situ high-speed X-ray imaging and mirrored via multiphysics simulation.

additive manufacturing↗
NASA NTRS1983

Force balance in the magnetospheres of Jupiter and Saturn

Spacecraft measurements of the plasma populations and magnetic fields near Jupiter and Saturn have revealed that large magnetospheres surround both planets. Magnetic field measurements have indicated closed field line topologies in the dayside magnetospheres of both planets while plasma instruments have shown these regions to be populated by both hot and cold plasma components convected azimuthally in the sense of planetary rotation. By using published data from the Voyager Plasma Science (PLS), Low Energy Charged Particle (LECP), and Magnetometer (MAG) instruments, it is possible to investigate the validity of the time stationary MHD momentum equation in the middle magnetospheres of Jupiter and Saturn. At Saturn, the hot plasma population is negligible in the dynamic sense and the centrifugal force of the cold rotating plasma appears to balance the Lorentz force. At Jupiter, the centrifugal force balances about 25 percent of the Lorentz force. The remaining inward Lorentz force is balanced by pessure gradients in the hot, high-beta plasma of the Jovian magnetodisk.

Mcnutt, R. L., Jr.↗
DOE OSTI2024

Real space iterative reconstruction for vector tomography (RESIRE-V)

Tomography has had an important impact on the physical, biological, and medical sciences. To date, most tomographic applications have been focused on 3D scalar reconstructions. However, in some crucial applications, vector tomography is required to reconstruct 3D vector fields such as the electric and magnetic fields. Over the years, several vector tomography methods have been developed. Here, we present the mathematical foundation and algorithmic implementation of REal Space Iterative REconstruction for Vector tomography, termed RESIRE-V. RESIRE-V uses multiple tilt series of projections and iterates between the projections and a 3D reconstruction. Each iteration consists of a forward step using the Radon transform and a backward step using its transpose, then updates the object via gradient descent. Incorporating with a 3D support constraint, the algorithm iteratively minimizes an error metric, defined as the difference between the measured and calculated projections. The algorithm can also be used to refine the tilt angles and further improve the 3D reconstruction. To validate RESIRE-V, we first apply it to a simulated data set of the 3D magnetization vector field, consisting of two orthogonal tilt series, each with a missing wedge. Our quantitative analysis shows that the three components of the reconstructed magnetization vector field agree well with the ground-truth counterparts. We then use RESIRE-V to reconstruct the 3D magnetization vector field of a ferromagnetic meta-lattice consisting of three tilt series. Our 3D vector reconstruction reveals the existence of topological magnetic defects with positive and negative charges. We expect that RESIRE-V can be incorporated into different imaging modalities as a general vector tomography method. To make the algorithm accessible to a broad user community, we have made our RESIRE-V MATLAB source codes and the data freely available at https://github.com/minhpham0309/RESIRE-V.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Divertor shaping with neutral baffling as a solution to the tokamak power exhaust challenge

Exhausting power from the hot fusion core to the plasma-facing components is one fusion energy’s biggest challenges. The MAST Upgrade tokamak uniquely integrates strong containment of neutrals within the exhaust area (divertor) with extreme divertor shaping capability. By systematically altering the divertor shape, this study shows the strongest evidence to date to our knowledge that long-legged divertors with a high magnetic field gradient (total flux expansion) deliver key power exhaust benefits without adversely impacting the hot fusion core. These benefits are already achieved with relatively modest geometry adjustments that are more feasible to integrate in reactor designs. Benefits include reduced target heat loads and improved access to, and stability of, a neutral gas buffer that ‘shields’ the target and enhances power exhaust (detachment). Analysis and model comparisons shows these benefits are obtained by combining multiple shaping aspects: long-legged divertors have expanded plasma-neutral interaction volume that drive reductions in particle and power loads, while total flux expansion enhances detachment access and stability. Containing the neutrals in the exhaust area with physical structures further augments these shaping benefits. These results demonstrate strategic variation in the divertor geometry and magnetic topology is a potential solution to one of fusion’s power exhaust challenge.

magnetically confined plasmas↗
DOE OSTI2022

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Scalable Deep-Learning-Accelerated Topology Optimization for Additively Manufactured Materials

Topology optimization (TO) is a popular and powerful computational approach for designing novel structures, materials, and devices. Two computational challenges have limited the applicability of TO to a variety of industrial applications. First, a TO problem often involves a large number of design variables to guarantee sufficient expressive power. Second, many TO problems require a large number of expensive physical model simulations, and those simulations cannot be parallelized. To address these issues, we propose a general scalable deep-learning (DL) based TO framework, referred to as SDL-TO, which utilizes parallel schemes in high performance computing (HPC) to accelerate the TO process for designing additively manufactured (AM) materials. Unlike the existing studies of DL for TO, our framework accelerates TO by learning the iterative history data and simultaneously training on the mapping between the given design and its gradient. The surrogate gradient is learned by utilizing parallel computing on multiple CPUs incorporated with a distributed DL training on multiple GPUs. The learned TO gradient enables a fast online update scheme instead of an expensive update based on the physical simulator or solver. Using a local sampling strategy, we achieve to reduce the intrinsic high dimensionality of the design space and improve the training accuracy and the scalability of the SDL-TO framework. The method is demonstrated by benchmark examples and AM materials design for heat conduction. The proposed SDL-TO framework shows competitive performance compared to the baseline methods but significantly reduces the computational cost by a speed up of around 8.6x over the standard TO implementation.

Bi, Sirui↗
NASA NTRS1997

Pioneer 10 and 11 Data Analysis

This report finishes the work of NASA Grant NAS2-153, which supported data analysis for the UCSD instruments on Pioneers 10 and 11. The data analyzed under this grant span 22 years of interplanetary measurements in the inner and outer heliosphere. The UCSD instruments made their mark in cosmic ray research based upon their high energy thresholds, directional responses, and reliable data streams. one of their primary scientific objectives concerns the size, configuration, and time behavior of the heliosphere. The size scale is inferred from the radial intensity gradient, which is measured between the two spacecraft and extrapolated to interstellar intensity levels at the cosmic ray modulation boundary. This boundary still eludes us, and its position, motion, and the best method of extrapolation are ongoing problems. Current projections place the boundary beyond 100 AU, which may be beyond the termination shock, and raises the question of possible modulation in the heliosheath. Probably our only hope of seeing this region in the immediate future rides on the possibility that the boundary will move inward. Our instruments have recorded many Forbush, or transient, decreases in the outer heliosphere. These observations led us to a model that attributes many of the decreases to solar wind stream-stream interactions, and relates the cosmic ray variations to the locally observed magnetic field magnitude. As the cosmic ray variations in this model result only from topological changes in the modulation integral, the model is a tool for studying the possibility that the 11 year cosmic ray modulation cycle can be accounted for by a superposition of Forbush decreases. The cosmic ray angular distribution function is measurable, given a good telemetry rate, by the UCSD Cerenkov detector which counts particles of energy greater than 500 MeV/n. We obtained statistically significant samples from 1 to 9 AU, at 13 AU, and at 34 AU. The anisotropy tends to be a few tenths of a per cent at all radial distances. A quasiperiodic variation in the east-west anisotropy with period of about 50 days remains unexplained.

Fillius, Walker↗
DOE OSTI2020

Generalized Grad–Shafranov equation for non-axisymmetric MHD equilibria

The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad–Shafranov equation, which is a single elliptic partial differential equation in two space dimensions. By revealing a hidden symmetry, here we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad–Shafranov equation. In contrast to solutions of the classical Grad–Shafranov equation, solutions of the generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad–Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A cryogenically cooled 200 kV DC photoemission electron gun for ultralow emittance photocathodes

Novel photocathode materials like ordered surfaces of single crystal metals, epitaxially grown high quantum efficiency thin films, and topologically non-trivial materials with dirac cones show great promise for generating brighter electron beams for various accelerator and ultrafast electron scattering applications. Despite several materials being identified as brighter photocathodes, none of them have been tested in electron guns to extract electron beams due to technical and logistical challenges. In this paper, we present the design and commissioning of a cryocooled 200 kV DC electron gun that is capable of testing a wide variety of novel photocathode materials over a broad range of temperatures from 298 to 35 K for bright electron beam generation. This gun is designed to enable easy transfer of the photocathode to various standard ultra-high-vacuum surface diagnostics and preparation techniques, allowing a full characterization of the dependence of beam brightness on the photocathode material and surface properties. We demonstrate the development of such a high-voltage, high-gradient gun using materials and equipment that are easily available in any standard university lab, making the development of such 200 kV electron guns more accessible.

Instruments & Instrumentation↗
DOE OSTI2023

An adaptive sampling augmented Lagrangian method for stochastic optimization with deterministic constraints

The primary goal of this paper is to provide an efficient solution algorithm based on the augmented Lagrangian framework for optimization problems with a stochastic objective function and deterministic constraints. Our main contribution is combining the augmented Lagrangian framework with adaptive sampling, resulting in an efficient optimization methodology validated with practical examples. To achieve the presented efficiency, here we consider inexact solutions for the augmented Lagrangian subproblems, and through an adaptive sampling mechanism, we control the variance in the gradient estimates. Furthermore, we analyze the theoretical performance of the proposed scheme by showing equivalence to a gradient descent algorithm on a Moreau envelope function, and we prove sublinear convergence for convex objectives and linear convergence for strongly convex objectives with affine equality constraints. The worst-case sample complexity of the resulting algorithm, for an arbitrary choice of penalty parameter in the augmented Lagrangian function, is $\mathscr{O}$(ϵ -3-δ ) , where ϵ > 0 is the expected error of the solution and δ > 0 is a user-defined parameter. If the penalty parameter is chosen to be $\mathscr{O}$(ϵ -1 ), we demonstrate that the result can be improved to $\mathscr{O}$(ϵ -2 ) , which is competitive with the other methods employed in the literature. Moreover, if the objective function is strongly convex with affine equality constraints, we obtain $\mathscr{O}$(ϵ -1 log(1/ϵ)) complexity. Finally, we empirically verify the performance of our adaptive sampling augmented Lagrangian framework in machine learning optimization and engineering design problems, including topology optimization of a heat sink with environmental uncertainty.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Global flow structure and exact formal transseries of the Gubser flow in kinetic theory

In this work we introduce the generic conditions for the existence of a non-equilibrium attractor that is an invariant manifold determined by the long-wavelength modes of the physical system. We investigate the topological properties of the global flow structure of the Gubser flow for the Israel-Stewart theory and a kinetic model for the Boltzmann equation by employing Morse-Smale theory. We present a complete classification of the invariant submanifolds of the flow and determine all the possible flow lines connecting any pair of UV/IR fixed points. The formal transseries solutions to the Gubser dynamical system around the early-time (UV) and late-time (IR) fixed points are constructed and analyzed. It is proven that these solutions are purely perturbative (or power-law asymptotic) series with a finite radius of convergence. Based on these analyses, we find that Gubser-like expanding kinetic systems do not hydrodynamize owing to the failure of the hydrodynamization process which heavily relies on the classification of (non)hydrodynamic modes in the IR regime. This is in contrast to longitudinal boost-invariant plasmas where the asymptotic dynamics is described by a few terms of the hydrodynamic gradient expansion. We finally compare our results for both Bjorken and Gubser conformal kinetic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A shift from inorganic to organic nitrogen-dominance shapes soil microbiome composition and co-occurrence networks

Soil microbiomes are characterized by their composition and networks, which are linked to soil nitrogen (N) availability. In nature, inorganic N dominates at one end and organic N dominates at the other end along soil N gradients; however, little is known about how this shift influences soil microbiome composition and co-occurrence networks, as well as their controls. To this end, we conducted an experiment with the host plant Solidago canadensis , which was subject to three N regimes: inorganic N-dominated, co-dominated by inorganic and organic N (CIO), and organic N-dominated. Organic N dominance exhibited stronger effects on the composition and co-occurrence networks of soil microbiomes than inorganic N dominance. The predominant control was plant traits for bacterial and fungal richness, and soil pH for keystone species. Relative to the CIO regime, inorganic N dominance did not affect fungal richness and increased keystone species; organic N dominance decreased fungal richness and keystone species. Pathogenic fungi and arbuscular mycorrhizal fungi were suppressed by organic N dominance but not by inorganic N dominance. These findings suggest that the shift from soil inorganic N-dominance to soil organic N-dominance could strongly shape soil microbiome composition and co-occurrence networks by altering species diversity and topological properties.

Xin, Yue↗
DOE OSTI2022

Exploring the Basic Physical Mechanisms of Cathode- and Anode-Initiated High-Voltage Surface Flashover

Surface flashover in vacuum imposes a substantial physical limit on modern, large-scale pulsed power. One of the ramifications is a minimum size requirement for new machines, which in itself becomes a hard barrier to the modernization and improvement of existing infrastructure. Pulsed power topologies require the physical mechanisms of both anode- and cathode-initiated flashover to be considered. Originally, the geometrical implications of field emission at the cathode triple junction (CTJ) motivated the usage of configurations that avoid electrons impinging on the insulating material. This will largely suppress the cathode-initiated flashover, which is best described by the secondary electron avalanche mechanism, gas desorption, and final breakdown in the desorbed gas. It depends on the cascade growth of a conducting plasma along the length of the insulator from the cathode. Mitigating the cathode-initiated flashover typically comes at the cost of a significant field enhancement at the anode triple junction (ATJ). In a typical implementation, the anode field may be three times higher than the cathode field for a given voltage, making the corresponding anode-initiated flashover much more common than otherwise. In the case of pulsed, anode-initiated flashover, experimental evidence suggests that charge is directly extracted from the insulator resulting in the insulator taking on a net positive charge advancing the anode potential. Furthermore, along with accompanying gas desorption from the surface, the potential will then propagate from the anode toward the cathode until the effective length of the gap is sufficiently reduced to support flashover. The underlying physical mechanisms of cathode- and anode-directed flashover are discussed in light of previously gathered experimental data and recent experiments with pulsed, high-gradient, anode-initiated flashover.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗
DOE OSTI2025

Differentiable lagrangian shock hydrodynamics with application to stable shock acceleration of density interfaces

We develop a gradient based optimization approach for the equations of compressible, Lagrangian hydrodynamics and demonstrate how it can be employed to automatically uncover strategies to control hydrodynamic instabilities arising from shock acceleration of density interfaces. Strategies for controlling the Richtmyer-Meshkov instability (RMI) are of great benefit for inertial confinement fusion (ICF) where shock interactions with many small imperfections in the density interface lead to instabilities which rapidly grow over time. These instabilities lead to mixing which, in the case of laser driven ICF, quenches the runaway fusion process ruining the potential for positive energy return. Here, we demonstrate that control of these instabilities can be achieved by optimization of initial conditions with ( > 100) parameters. Optimizing over a large parameter space like this is not possible with gradient-free optimization strategies. This requires computation of the gradient of the outputs of a numerical solution to the equations of Lagrangian hydrodynamics with respect to the inputs. We show that the efficient computation of these gradients is made possible via a judicious application of (i) adjoint methods, the exact formal representation of sensitivities involving partial differential equations, and (ii) automatic differentiation (AD), the algorithmic calculation of derivatives of functions. Careful regularization of multiple operators including artificial viscosity and timestep control is required. We perform design optimization of > 100 parameter energy field driving the Richtmyer Meshkov instability showing significant suppression while simultaneously enhancing the acceleration of the interface relative to a nominal baseline case.

Hydrophysics↗