Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometry optimization”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Suppressing trapped-electron-mode-driven turbulence via optimization of three-dimensional shaping

Abstract Turbulent transport driven by trapped electron modes (TEMs) is believed to drive significant heat and particle transport in quasihelically symmetric stellarators. Two three-dimensionally-shaped magnetic configurations with suppressed TEM-driven turbulence were generated through optimization that targeted quasihelical symmetry and the available energy of trapped electrons. Initial equilibria have flux surface shapes with a helically rotating negative triangularity (NT) and positive triangularity (PT). In gyrokinetic simulations, TEMs are suppressed in the reduced-TEM NT and PT configurations, showing that NT does not have the same beneficial turbulence properties over PT as seen in tokamaks. Heat fluxes from TEMs are also suppressed. Without temperature gradients and with a strong density gradient, the most unstable modes at low k y were consistent with toroidal universal instabilities (UIs) in the NT case and slab UIs in the PT case. Nonlinear simulations show that UIs drive substantial heat flux in both the NT and PT configurations. A moderate increase in β halves the heat flux in the NT configuration, while suppressing the heat flux in the PT geometry. Based on the present work, future optimizations aimed at reducing electrostatic drift wave-driven turbulent transport will need to consider UIs if β is sufficiently small.

Duff, J. M. (ORCID:0000000153086694)↗
DOE OSTI2023

Unified, Geometric Framework for Nonequilibrium Protocol Optimization

Controlling thermodynamic cycles to minimize the dissipated heat is a long-standing goal in thermodynamics, and more recently, a central challenge in stochastic thermodynamics for nanoscale systems. Here, we introduce a theoretical and computational framework for optimizing nonequilibrium control protocols that can transform a system between two distributions in a minimally dissipative fashion. These protocols optimally transport a system along paths through the space of probability distributions that minimize the dissipative cost of a transformation. Furthermore, we show that the thermodynamic metric—determined via a linear response approach—can be directly derived from the same objective function that is optimized in the optimal transport problem, thus providing a unified perspective on thermodynamic geometries. As a result, we investigate this unified geometric framework in two model systems and observe that our procedure for optimizing control protocols is robust beyond linear response.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Blade and Rim Seal Design of a First Stage High Pressure Turbine for a 300 MWe Supercritical CO2 Power Cycle

A first stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design 3D with objectives of increasing efficiency and decreasing heat load to the blade. The National Institute of Standards and Technology Reference Fluid Thermodynamic and Transport Properties Database (NIST REFPROP) [1] data for supercritical CO2 is formatted into tables of bicubic polynomial coefficients for use in condensable gas simulations in FINE/Turbo. Nearly 3000 unique shapes are evaluated via three-dimensional Reynolds Averaged Navier Stokes simulations, yielding increases in efficiency of up to 0.85 percentage points and decreases in heat load of 14%. A final blade, deemed the advanced blade, is chosen for future experimental analysis. Following this, a squealer tip optimization is performed on both the baseline and advanced blade designs. This optimization resulted in a performance gain of 1.25 points in efficiency and 15% reduction in tip heat load compared to the baseline flat tip design at the same clearance. In tandem, an optimization of the rotor-stator platform rim seal is performed using a parametrized geometry allowing for straight, meandering, and knife seal cavities. This multi-objective optimization focuses on decreasing the cooling mass flow and increasing the heat flux from the rotor and stator disk. The optimization resulted in cooling mass flow decreases of up to 26% while maintaining the average heat flux on the rim seal.

Tuite, Logan↗
DOE OSTI2024

Adaptive Mesh Refinement Simulations for Turbulent Reacting Flow

With the increased availability of exascale computing hardware, detailed simulations of realistic devices can be performed at practically relevant time and length scales. Insights into the multiscale driving mechanisms in compressible reacting flow systems with complex geometry, such as combustors, can be used for design optimization and technology improvements. However, to effectively perform these simulations, advanced numerical algorithms must be used to maintain solution accuracy without incurring undue computational costs. PeleC, part of the Pele suite of codes, leverages block-structured adaptive mesh refinement (AMR) through the AMReX library to capture fine-scale flow features in compressible reacting flows. In this talk, we discuss recent improvements to the numerical algorithms, particularly in regard to describing flows at complex boundary structures, and PeleC's performance on exascale computing hardware. We will demonstrate that PeleC is well-suited for modern, extreme-scale, heterogenous compute platforms.

combustion↗
DOE OSTI2024

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Piston Bowl Geometry Effects on Gasoline Compression Ignition in a Heavy-Duty Diesel Engine

A design optimization campaign was conducted to search for improved combustion profiles that enhance gasoline compression ignition in a heavy-duty diesel engine with a geometric compression ratio of 17.3. A large-scale design of experiments approach was used for the optimization, employing three-dimensional computational fluid dynamics simulations. The main parameters explored include geometric features, injector specifications, and swirl motion. Both stepped-lip and re-entrant bowls were included in order to assess their respective performance implications. A total of 256 design candidates were prepared using the software package CAESES for automated and simultaneous geometry generation and combustion recipe perturbation. The design optimization was conducted for three engine loads representing light to medium load conditions. The design candidates were evaluated for fuel efficiency, emissions, fuel-air mixing, and global combustion behavior. Simulation results showed that the optimum designs were all stepped-lip bowls, due to improvements in fuel-air mixing, as well as reduced heat loss and emissions formation. Improvements in indicated specific fuel consumption of up to 3.2% were achieved while meeting engine-out NOx emission targets of 1-1.5 g/kW center dot h. Re-entrant bowls performed worse compared to the baseline design, and significant performance variations occurred across the load points. Specifically, the re-entrant bowls were on par with the stepped-lip bowls under light load conditions, but significant deteriorations occurred under higher load conditions. As a final task, selected optimized designs were then evaluated under full-load conditions.

CFD↗
DOE OSTI2022

Krylov subspace recycling for evolving structures

Krylov subspace recycling is a powerful tool when solving a long series of large, sparse linear systems that change only slowly over time. In PDE constrained shape optimization, these series appear naturally, as typically hundreds or thousands of optimization steps are needed with only small changes in the geometry. In this setting, however, applying Krylov subspace recycling can be a difficult task. As the geometry evolves, in general, so does the finite element mesh defined on or representing this geometry, including the numbers of nodes and elements and element connectivity. This is especially the case if re-meshing techniques are used. As a result, the number of algebraic degrees of freedom in the system changes, and in general the linear system matrices resulting from the finite element discretization change size from one optimization step to the next. Changes in the mesh connectivity also lead to structural changes in the matrices. In the case of re-meshing, even if the geometry changes only a little, the corresponding mesh might differ substantially from the previous one. Obviously, this prevents any straightforward mapping of the approximate invariant subspace of the linear system matrix (the focus of recycling in this work) from one optimization step to the next; similar problems arise for other selected subspaces. In this paper, we present an algorithm to map an approximate invariant subspace of the linear system matrix for the previous optimization step to an approximate invariant subspace of the linear system matrix for the current optimization step, for general meshes. This is achieved by exploiting the map from coefficient vectors to finite element functions on the mesh, combined with interpolation or approximation of functions on the finite element mesh. We demonstrate the effectiveness of our approach numerically with several proof of concept studies for a specific meshing technique.

42 ENGINEERING↗
DOE OSTI2021

Joining AA7099 to Ni-Cr-Mo Steel Using Single Pass Friction Stir Dovetailing and AA6061 Butter Layer

A novel single pass friction stir dovetailing (FSD) was developed to join thick plate AA7099 to Ni-Cr-Mo steel in a lap configuration through the use of an AA6061 butter layer. The butter layer approach can be shown to significantly reduce interfacial failures by inhibiting the formation of brittle Fe–Al intermetallic compounds (IMCs) which are typically formed during conventional friction stir welding (FSW) of AA7099 directly to steel. This study evaluates the strength and microstructural evolution of the aluminum/steel joint interface between single pass FSD of AA7099/AA6061(butter)/steel fabricated using three different FSD tool configurations. The butter layer was achieved by placement of a free-floating strip of AA6061 within the dovetail groove that was pre-machined in the steel substrate. Friction stir welding was then used to simultaneously join the thick plate AA7099, the AA6061 strip and the AA7099 in a single weld pass. To optimize the joining process, three different pin geometries (full thread, half smooth/half thread, no thread) were evaluated; it was demonstrated that the half smooth/half thread tool produced the strongest weld joint due to the absence of a Zn rich layer at the Al/steel interface. The efficacy of this unique processing technique was tested using lap shear tensile testing as well as high resolution transmission electron microscopy (TEM). A discussion of the new FSD technique, joint configurations, and process parameters are provided along with joint microstructural analyses and mechanical performance.

Olszta, Matthew J.↗
DOE OSTI2024

Closed-loop electron-beam-induced spectroscopy and nanofabrication around individual quantum emitters

Color centers in diamond play a central role in the development of quantum photonic technologies, and their importance is only expected to grow in the near future. For many quantum applications, high collection efficiency from individual emitters is required, but the refractive index mismatch between diamond and air limits the optimal collection efficiency with conventional diamond device geometries. While different out-coupling methods with near-unity efficiency exist, many have yet to be realized due to current limitations in nanofabrication methods, especially for mechanically hard materials like diamond. Here, we leverage electron-beam-induced etching to modify Sn-implanted diamond quantum microchiplets containing integrated waveguides with a width and thickness of 280 nm and 200 nm, respectively. This approach allows for simultaneous high-resolution imaging and modification of the host matrix with an open geometry and direct writing. When coupled with the cathodoluminescence signal generated from the electron–emitter interactions, we can monitor the enhancement of the quantum emitters in real-time with nanoscale spatial resolution. The operando cathodoluminescence measurement and fabrication around single photon emitters demonstrated here provide a new foundation for the potential control of emitter–cavity interactions in integrated quantum photonics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toward Qualifications of HB and LB 650 MHz Cavities for the Prototype Cryomodules for the PIP-II Project

High-beta (HB) and low-beta (LB) 650 MHz cryomodules are key components of the Proton Improvement Plan II (PIP-II) project. In this contribution we present the results of several 5-cell HB650 cavities that have been processed and tested with the purpose of qualifying them for the prototype cryomodule assembly, which will take place later this year. We also present the first results obtained in LB650 single-cell cavities process optimization. Taking advantage of their very similar geometry, we are also analyzing the effect of different surface treatments in FRIB’s 5-cell medium-beta 644MHz cavities. Cavities processed with N-doping and mid-T baking showed very promising results in term of both Q-factors and accelerating gradient for these low-beta structures.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Additive Manufacturing of Heat Pipes for Microreactor Applications

Heat pipes are highly effective devices used to passively transport heat by two-phase capillary action. Their small footprint, light weight and lack of moving parts make them ideal for transportable nuclear microreactors, where size and weight are limited, minimal to no operation and maintenance are required, and high reliability is desired. The use of additive manufacturing allows for the fabrication of heat pipes with performance enhancements, such as microchannels, grooves, arteries, or tailored porosity, which are difficult, expensive or impossible to achieve with traditional fabrication techniques. We are evaluating the capabilities of the Admaflex 130 Digital Light Printer to produce these geometries by characterizing feature shrinkage, developing an optimal heat treatment recipe, assessing the effects of heat treatment atmosphere and varying print orientation. Computational simulations are being used to estimate performance of the various 3D printed configurations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Analysis of the ratio of ℓ 1 and ℓ 2 norms in compressed sensing

We study the ratio of ℓ 1 and ℓ 2 norms ( ℓ 1 / ℓ 2 ) as a sparsity-promoting objective in compressed sensing. We first propose a novel criterion that guarantees that an s-sparse signal is the local minimizer of the ℓ 1 / ℓ 2 objective; our criterion is interpretable and useful in practice. We also give the first uniform recovery condition using a geometric characterization of the null space of the measurement matrix, and show that this condition is satisfied for a class of random matrices. We also present analysis on the robustness of the procedure when noise pollutes data. Numerical experiments are provided that compare ℓ 1 / ℓ 2 with some other popular non-convex methods in compressed sensing. Finally, we propose a novel initialization approach to accelerate the numerical optimization procedure. We call this initialization approach support selection, and we demonstrate that it empirically improves the performance of existing ℓ 1 / ℓ 2 algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Enabling the broader adoption of fusion simulation on complex geometry

This project addressed a key barrier to advanced fusion and nuclear simulation: the difficulty of performing high-fidelity Monte Carlo neutronics directly on complex, real-world CAD geometry. Traditional workflows require engineers to rebuild CAD models as simplified constructive solid geometry, a time-consuming and error-prone process that limits design iteration and broader adoption of simulation tools. The goal of this Phase I SBIR was to make CAD-based neutronics practical, accessible, and robust for industrial and research users. During the project, Coreform significantly enhanced the Direct Accelerated Geometry Monte Carlo (DAGMC) workflow and fully integrated it into Coreform Cubit as a first-class capability. Major achievements include optimized material assignment and surface meshing workflows, substantial performance improvements to geometry imprinting and preparation, native export of DAGMC models, and new visualization tools to support OpenMC source definition and lost-particle debugging. Coreform also expanded Cubit’s capabilities as a full OpenMC preprocessor, including the ability to convert OpenMC constructive solid geometry models back into CAD for visualization, multiphysics coupling, and debugging. In collaboration with Argonne National Laboratory, the project delivered comprehensive new DAGMC documentation and training materials, transforming DAGMC from a research-oriented tool into a production-ready workflow. Results were disseminated through tutorials, conference training, and multiple well-attended webinars demonstrating integrated CAD-based neutronics and multiphysics workflows. Overall, this project demonstrated that high-fidelity Monte Carlo simulations can be performed directly on complex CAD geometry, reducing setup time, improving usability, and enabling faster, more informed design decisions for fusion and nuclear energy systems.

42 ENGINEERING↗
DOE OSTI2021

A low-frequency buncher field-variation study on a 750 keV H – beam to increase Drift Tube Linac capture

We report low-frequency (such as 16.77 MHz) RF bunchers are widely used in RF accelerator systems for longitudinal compression of pulses into a single RF bucket, which increases instantaneous beam intensity for time-dependent studies. In this study, the dependency of capture into a 201.25 MHz Drift Tube Linac (DTL) was measured as a function of gap voltage for a 16.77 MHz buncher on chopped H - beam (approximately 25 ns at 750 keV, 10 mA peak current). The multiparticle code PARMILA was used to simulate the phase-space distribution of the 10 mA, 750 keV, H - beam at the entrance to DTL with a wide range of the Low-Frequency Buncher (LFB) field (10 kV to 35 kV). The measurement and simulation indicated that the DTL capture could be dilute (reduced) for a non-optimized buncher field to a pre-configured beamline geometry. The data shows that changing the bunch field while keeping the incoming beam current and energy constant does not significantly alter the beam’s emittance. However, downstream beam capture into the DTL is changed for a non-optimized phase-space bunching distance with the buncher field.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Annual Performance of the Second-Generation Variable-Geometry Oscillating Surge Wave Energy Converter

Recent studies in wave energy have highlighted the need for a structured innovation approach in wave energy converter (WEC) design because cost-of-energy estimates have remained high. One specific innovation being investigated by the National Renewable Energy Laboratory is WEC geometry control, which uses control surfaces in combination with an oscillating surge WEC (OSWEC) to increase device availability and power generation while limiting structural costs. This study performs the first analysis of annual performance for a novel OSWEC with geometry control to understand how the geometry control affects availability, annual power generation, and structural loadings like the power-take-off (PTO) torque or surge foundation force. Device hydrodynamic coefficients are calculated using linear potential theory for six geometry configurations. A nonideal PTO system is assumed and quadratic viscous damping losses are considered. Annual performance is evaluated and compared for three U.S. wave energy sites. The WEC geometry and PTO system are controlled on a sea-state basis to optimize for power capture while remaining under limits set on motion amplitude and structural loads. Results indicate that geometry control can increase availability up to 25 days in an average year depending on design limits, increase average power generation, and significantly reduce peak structural loads.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗