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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Parallel Cut-Cell Algorithm for the Free-Boundary Grad--Shafranov Problem

A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad--Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. Here the algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad--Shafranov solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Observations and Modeling of Fiber Optic Strain on Hydraulic Fracture Height Growth in Hydraulic Fracturing Test Site 2 (HFTS-2)

Summary Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced data set to allow us to observe and understand fracture geometries rigorously. Low-frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1,000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. The simulation results improved the understating of the HFTS-2 LF-DAS data, and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

Engineering↗
DOE OSTI2021

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗
NASA NTRS1998

Integration of a CAD System Into an MDO Framework

NASA Langley has developed a heterogeneous distributed computing environment, called the Framework for Inter-disciplinary Design Optimization, or FIDO. Its purpose has been to demonstrate framework technical feasibility and usefulness for optimizing the preliminary design of complex systems and to provide a working environment for testing optimization schemes. Its initial implementation has been for a simplified model of preliminary design of a high-speed civil transport. Upgrades being considered for the FIDO system include a more complete geometry description, required by high-fidelity aerodynamics and structures codes and based on a commercial Computer Aided Design (CAD) system. This report presents the philosophy behind some of the decisions that have shaped the FIDO system and gives a brief case study of the problems and successes encountered in integrating a CAD system into the FEDO framework.

Townsend, J. C.↗
NASA NTRS2017

Structural Health Management of Damaged Aircraft Structures Using the Digital Twin Concept

The development of multidisciplinary integrated Structural Health Management (SHM) tools will enable accurate detection, and prognosis of damaged aircraft under normal and adverse conditions during flight. As part of the digital twin concept, methodologies are developed by using integrated multiphysics models, sensor information and input data from an in-service vehicle to mirror and predict the life of its corresponding physical twin. SHM tools are necessary for both damage diagnostics and prognostics for continued safe operation of damaged aircraft structures. The adverse conditions include loss of control caused by environmental factors, actuator and sensor faults or failures, and structural damage conditions. A major concern in these structures is the growth of undetected damage/cracks due to fatigue and low velocity foreign object impact that can reach a critical size during flight, resulting in loss of control of the aircraft. To avoid unstable, catastrophic propagation of damage during a flight, load levels must be maintained that are below a reduced load-carrying capacity for continued safe operation of an aircraft. Hence, a capability is needed for accurate real-time predictions of damage size and safe load carrying capacity for structures with complex damage configurations. In the present work, a procedure is developed that uses guided wave responses to interrogate damage. As the guided wave interacts with damage, the signal attenuates in some directions and reflects in others. This results in a difference in signal magnitude as well as phase shifts between signal responses for damaged and undamaged structures. Accurate estimation of damage size, location, and orientation is made by evaluating the cumulative signal responses at various pre-selected sensor locations using a genetic algorithm (GA) based optimization procedure. The damage size, location, and orientation is obtained by minimizing the difference between the reference responses and the responses obtained by wave propagation finite element analysis of different representative cracks, geometries, and sizes.

Seshadri, Banavara R.↗
NASA NTRSDate not supplied

Langley Automated Sensor Inter-calibration System (LASICS): Open Access Tools for Satellite Imager Inter- Calibration

Satellite imager calibration teams are tasked with maintaining stable measurement records to facilitate reliable monitoring of geophysical parameters and ensure dependable input for forecast models. Satellite imagers are neither uniformly calibrated nor radiometrically scaled to a common reference standard. Consistent inter-calibration between various earth-orbiting satellite imager pairs is a critical step in the creation of seamless earth-scene reflectance data records over time for input to higher level algorithms that retrieve earth system climate-sensitive properties. Each imager inherently by virtue of its optics (and associated properties like spectral response etc.) and orbit will have a unique measurement of the same earth-scene reflected signal. A key part of this is the computationally efficient and optimal identification and prediction of science opportunities where the imager pairs from the irrespective earth-orbits near-simultaneously view the same stable terrestrial targets with nearly identical viewing and solar geometry.

Arun Gopalan↗
DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2022

Understanding resistance oscillation in the CsV 3 Sb 5 superconductor

A recent demonstration of the periodic oscillation of resistance in the thin film of the CsV 3 ⁢Sb 5 superconductor with a hole in the film suggests that charge-4⁢ e and charge-6⁢⁢ e Cooper pairs may have condensed in this compound. While exciting, such interpretation calls for a precise determination of the effective area for the passage of Cooper pairs from one end of the lead to the other. Unlike the traditional Little-Parks effect where the rim around the hole is thin, the effective hole area is not obviously defined for the “thick-rim geometry” adopted in the experiment. Here, we note that the experiment was conducted in a regime where the superconductivity is strongly fluctuating, which motivates an analysis based on the spacetime formulation of the time-dependent Ginzburg-Landau theory. Further, we argue that under appropriate conditions, the optimal semiclassical path is not the geometrically shortest one, but the one that moves along the edge of the hole and takes advantage of the reduced fluctuations at the boundary of the hole. The condition for the crossover from the geometrically shortest path to the path that sticks to the wall is clarified. In such a scenario, the geometric area of the hole indeed emerges as the effective area for the flux, providing a theoretical justification to the interpretation given by J. Ge et al. (arXiv:2201.10352). The conclusion of our analysis may have implications for similar experiments in other superconductors where the geometry of the device is not obviously that of the Little-Parks experiment employing the thin wall, but of the thick-rim type such as used in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

Aeroacoustic sensitivity analysis and optimal aeroacoustic design of turbomachinery blades

During the first year of the project, we have developed a theoretical analysis - and wrote a computer code based on this analysis - to compute the sensitivity of unsteady aerodynamic loads acting on airfoils in cascades due to small changes in airfoil geometry. The steady and unsteady flow though a cascade of airfoils is computed using the full potential equation. Once the nominal solutions have been computed, one computes the sensitivity. The analysis takes advantage of the fact that LU decomposition is used to compute the nominal steady and unsteady flow fields. If the LU factors are saved, then the computer time required to compute the sensitivity of both the steady and unsteady flows to changes in airfoil geometry is quite small. The results to date are quite encouraging, and may be summarized as follows: (1) The sensitivity procedure has been validated by comparing the results obtained by 'finite difference' techniques, that is, computing the flow using the nominal flow solver for two slightly different airfoils and differencing the results. The 'analytic' solution computed using the method developed under this grant and the finite difference results are found to be in almost perfect agreement. (2) The present sensitivity analysis is computationally much more efficient than finite difference techniques. We found that using a 129 by 33 node computational grid, the present sensitivity analysis can compute the steady flow sensitivity about ten times more efficiently that the finite difference approach. For the unsteady flow problem, the present sensitivity analysis is about two and one-half times as fast as the finite difference approach. We expect that the relative efficiencies will be even larger for the finer grids which will be used to compute high frequency aeroacoustic solutions. Computational results show that the sensitivity analysis is valid for small to moderate sized design perturbations. (3) We found that the sensitivity analysis provided important insight into how airfoils should be modified to improve their aeroelastic stability. Using this insight, we redesigned an aeroelastically unstable cascade to produce a stable cascade.

Hall, Kenneth C.↗
DOE OSTI2024

Predictive binder jet additive manufacturing enabled by clean burn-off binder design

Binder jet additive manufacturing (BJAM) presents an avenue for advanced manufacturing of various high-value materials due to high deposition rates, scalability, and geometric flexibility. However, conventional organic binders in BJAM introduce residual carbon upon pyrolysis, often leading to imprecise alloy composition in the final sintered part. The undesirable residual carbon from binder burn-off limits the application of BJAM for high-performance alloys due to their high sensitivity to carbon addition. In this study, we have designed poly(vinylpyrrolidone-co-vinyl acetate) (PVP-VAc) as a clean burn-off binder for BJAM, where excess oxygen groups in VAc enable cleaner burn-off and reduce residual carbon retention. Compared to a widely used commercial binder, the optimized PVP-VAc binder reduced residual carbon retention by over 90% in H13 tool steel. The significant reduction in residual carbon enables predictable printing and subsequent sintering of complex H13 tool steel geometries, an alloy known to have substantial challenges around distortion due to the sintering window shifting from carbon addition. The design of a clean burnout binder provides a major path forward for BJAM by enabling new AM designs and applications for composition-sensitive high-performance alloys, such as nickel-based superalloys, titanium alloys, and high alloyed steels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nature-inspired lotus-shaped fins combined with hybrid nanoparticles and metal foam for high-performance latent heat thermal energy storage

Latent heat thermal energy storage (LHTES) systems play a critical role in renewable energy integration by providing high energy density and nearly isothermal operation during phase transitions. However, their performance is often limited by slow melting/charging rates, which motivates the search for enhanced heat transfer designs. This study investigates the melting behavior of RT-82 phase change material (PCM) using novel lotus-shaped fins combined with copper metal foam and conductive graphene nanoparticles and carbon nanotubes. A two-dimensional enthalpy-porosity model in ANSYS Fluent was developed to simulate the charging/melting process, capturing non-thermal equilibrium between the foam and PCM/nano-PCM. In this study, effects of fin geometry, nanoparticle concentration, and foam porosity on melting dynamics and cost-performance trade-offs were investigated. Results showed that natural convection accelerated melting by ~12% compared to conduction-only scenarios. Optimized lotus-shaped fins with higher fin density (T3F4 and T3F10) achieved up to 63% faster melting relative to sparse configurations. Graphene nanoparticles improved thermal conductivity, with a 6% volume fraction, by reducing melting time by ~6.9%, while their combination with 75% porosity foam achieved a maximum reduction in the melting time of ~51% compared to pure PCM. Cost-performance analysis identified T3F4 as the most balanced design, offering rapid thermal response without excessive material costs, while moderate-density designs like T3S6 provided economical alternatives with acceptable performance. These results highlight the performance enhancement that can be achieved by integrating bio-inspired fins, nanoparticles, and foams, into compact and efficient LHTES for solar heating, building thermal management, and industrial waste-heat recovery applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗
DOE OSTI2022

Design of inertial fusion implosions reaching the burning plasma regime

Abstract In a burning plasma state 1–7 , alpha particles from deuterium–tritium fusion reactions redeposit their energy and are the dominant source of heating. This state has recently been achieved at the US National Ignition Facility 8 using indirect-drive inertial-confinement fusion. Our experiments use a laser-generated radiation-filled cavity (a hohlraum) to spherically implode capsules containing deuterium and tritium fuel in a central hot spot where the fusion reactions occur. We have developed more efficient hohlraums to implode larger fusion targets compared with previous experiments 9,10 . This delivered more energy to the hot spot, whereas other parameters were optimized to maintain the high pressures required for inertial-confinement fusion. We also report improvements in implosion symmetry control by moving energy between the laser beams 11–16 and designing advanced hohlraum geometry 17 that allows for these larger implosions to be driven at the present laser energy and power capability of the National Ignition Facility. These design changes resulted in fusion powers of 1.5 petawatts, greater than the input power of the laser, and 170 kJ of fusion energy 18,19 . Radiation hydrodynamics simulations 20,21 show energy deposition by alpha particles as the dominant term in the hot-spot energy balance, indicative of a burning plasma state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

CHESS: The future direct geometry spectrometer at the second target station

CHESS, chopper spectrometer examining small samples, is a planned direct geometry neutron chopper spectrometer designed to detect and analyze weak signals intrinsic to small cross sections (e.g., small mass, small magnetic moments, or neutron absorbing materials) in powders, liquids, and crystals. CHESS is optimized to enable transformative investigations of quantum materials, spin liquids, thermoelectrics, battery materials, and liquids. The broad dynamic range of the instrument is also well suited to study relaxation processes and excitations in soft and biological matter. The 15 Hz repetition rate of the Second Target Station at the Spallation Neutron Source enables the use of multiple incident energies within a single source pulse, greatly expanding the information gained in a single measurement. Furthermore, the high flux grants an enhanced capability for polarization analysis. This enables the separation of nuclear from magnetic scattering or coherent from incoherent scattering in hydrogenous materials over a large range of energy and momentum transfer. This paper presents optimizations and technical solutions to address the key requirements envisioned in the science case and the anticipated uses of this instrument.

47 OTHER INSTRUMENTATION↗
NASA NTRS1972

Design of microstrip components by computer

A number of computer programs are presented for use in the synthesis of microwave components in microstrip geometries. The programs compute the electrical and dimensional parameters required to synthesize couplers, filters, circulators, transformers, power splitters, diode switches, multipliers, diode attenuators and phase shifters. Additional programs are included to analyze and optimize cascaded transmission lines and lumped element networks, to analyze and synthesize Chebyshev and Butterworth filter prototypes, and to compute mixer intermodulation products. The programs are written in FORTRAN and the emphasis of the study is placed on the use of these programs and not on the theoretical aspects of the structures.

Cisco, T. C.↗